Thermodynamics and mechanism of alpha helix initiation in alanine and valine peptides

Abstract

We used molecular dynamics simulations to study the folding/unfolding of one of turn of an alpha helix in Ac-(Ala)3-NHMe and Ac-(Val)3-NHMe. Using specialized sampling techniques, we computed free energy surfaces as functions of a conformational coordinate that corresponds to alpha helices at small values and to extended conformations at large values. Analysis of the peptide conformations populated during the simulations showed that alpha helices, reverse turns, and extended conformations correspond to minima on the free energy surfaces of both peptides. The free energy difference between alpha helix and extended conformations, determined from the equilibrium constants for helix unfolding, is approximately -1 kcal/mol for Ac-(Ala)3-NHMe and -5 kcal/mol for Ac-(Val)3-NHMe. The mechanism observed in our simulations, which includes reverse turns as important intermediates along the helix folding/unfolding pathway, is consistent with a mechanism proposed previously. Our results predict that both peptides (but especially the Ala peptide) have a much larger equilibrium constant for helix initiation than is predicted by the helix-coil transition theory with the host-guest parameters. We also predict a much greater difference in the equilibrium constants than the theory predicts. Insofar as helix initiation is concerned, our results suggest that the large difference between the helical propensities of Ala and Val cannot be explained by simple concepts such as side-chain rotamer restriction or unfavorable steric interactions. Rather, the origin of the difference appears to be quite complicated because it involves subtle differences in the solvation of the two peptides. The two peptides have similar turn-extended equilibria but very different helix-turn equilibria, and the difference in helical propensities reflects the fact that the helix-turn equilibrium strongly favors the turns in Ac-(Val)3-NHMe, while it favors the helices in Ac-(Ala)3-NHMe. We also computed thermodynamic decompositions of the free energy surfaces, and these revealed that the helix-turn equilibria are vastly different primarily because the changes in peptide-water interactions that accompany helix-to-turn conformational changes are qualitatively different for the two peptides.