Entropic Stabilization of Proteins by TMAO

To understand the mechanism of trimethylamine N-oxide (TMAO) induced stabilization of folded protein states, we systematically investigated the action of TMAO on several model dipeptides (Leucine, L2, Serine, S2, Glutamine, Q2, Lysine, K2, and Glycine, G2) in order to elucidate the effect of residue-specific TMAO interactions on small fragments of solvent-exposed conformations of the denatured states of proteins. We find that TMAO preferentially hydrogen bonds with the exposed dipeptide backbone, but generally not with nonpolar or polar side chains. However, interactions with the positively charged Lys are substantially greater than with the backbone. The dipeptide G2, is a useful model of pure amide backbone, interacts with TMAO by forming a hydrogen bond between the amide nitrogen and the oxygen in TMAO. In contrast, TMAO is depleted from the protein backbone in the hexapeptide G6, which shows that the length of the polypeptide chain is relevant in aqueous TMAO solutions. These simulations lead to the hypothesis that TMAO-induced stabilization of proteins and peptides is a consequence of depletion of the solute from the protein surface provided intramolecular interactions are more favorable than those between TMAO and the backbone. To test our hypothesis we performed additional simulations of the action of TMAO on an intrinsically disordered A{\beta}16-22 (KLVFFAE) monomer. In the absence of TMAO A{\beta}16-22 is a disordered random coil. However, in aqueous TMAO solution A{\beta}16-22 monomer samples compact conformations. A transition from random coil to {\alpha}-helical secondary structure is observed at high TMAO concentrations. Our work highlights the potential similarities between the action of TMAO on long polypeptide chains and entropic stabilization of proteins in a crowded environment due to excluded volume interactions. In this sense TMAO is a nano-crowding particle.

💡 Research Summary

The paper investigates the molecular mechanism by which the osmolyte trimethylamine N‑oxide (TMAO) stabilizes folded protein states. Using all‑atom molecular dynamics simulations with the CHARMM22/CMAP force field and NAMD, the authors examined a series of model peptides in explicit water at 1 M TMAO concentration: five dipeptides representing non‑polar (Leu), polar (Ser, Gln), basic (Lys) residues, as well as diglycine (G₂) and hexaglycine (G₆). Additional simulations were performed on the intrinsically disordered Aβ₁₆‑₂₂ (KLVFFAE) peptide at 0, 1, 2.5, and 5 M TMAO.

Key findings include: (1) TMAO preferentially forms hydrogen bonds with the peptide backbone nitrogen (g(r≈3 Å)≈1.5) across all dipeptides, while interactions with side‑chain atoms are generally weak (g(r≈4 Å)≈1.0). The basic Lys side‑chain nitrogen shows a markedly stronger interaction (g(r≈3 Å)≈2.5), indicating that charged side‑chains can attract TMAO more effectively than neutral ones. (2) Diglycine, which lacks side‑chains, readily hydrogen‑bonds to TMAO via its amide nitrogen, confirming that the backbone is a primary binding site for short peptides. (3) In contrast, hexaglycine exhibits a pronounced depletion of TMAO from its vicinity; the radial distribution function between TMAO oxygen and backbone nitrogen is much lower than for G₂. This depletion coincides with an increased population of α‑helical conformations, as shown by Ramachandran free‑energy surfaces and a shift in the radius of gyration toward values characteristic of an ideal α‑helix. (4) The Aβ₁₆‑₂₂ peptide, initially a random coil, becomes progressively more compact with increasing TMAO concentration (average R_g decreases from 6.9 Å at 0 M to 5.7 Å at 5 M, a 17 % reduction). At ≥2.5 M TMAO, β‑sheet basins disappear from the Ramachandran landscape, and right‑handed α‑helical basins dominate, indicating a sharp coil‑to‑helix transition driven by TMAO.

The authors interpret these results through the lens of “nano‑crowding.” TMAO, being excluded from the peptide surface as the chain length increases, generates an osmotic pressure that favors intramolecular hydrogen bonding and compaction, analogous to the entropic stabilization observed in macromolecular crowding. For short peptides (G₂) direct hydrogen bonding to TMAO dominates, but for longer chains (G₆, Aβ₁₆‑₂₂) depletion becomes the governing effect, leading to structural collapse and secondary‑structure formation. This depletion‑induced mechanism contrasts with the direct interaction model proposed for denaturants such as urea, highlighting that protective osmolytes act by reducing solvent‑accessible volume rather than by forming stabilizing contacts.

In conclusion, the study provides compelling computational evidence that TMAO stabilizes proteins primarily via entropic exclusion from the protein surface, acting as a nano‑crowding particle that promotes compact, α‑helical conformations when intramolecular interactions outweigh TMAO‑backbone contacts. This insight advances our understanding of osmolyte‑mediated protein stability and may inform the design of novel stabilizing agents for biotechnological applications.