α/310‐Helix transitions in α‐methylalanine homopeptides: Conformational transition pathway and potential of mean force

Abstract

The conformational transition between the α‐ and 3₁₀‐helical states of α‐methylalanine homopeptides is studied with molecular mechanics. Conformational transition pathways for Ace‐(MeA)_n‐NMe with n = 7, 9, and 11 are obtained with the algorithms of Elber and co‐workers [R. Czerminski & R. Elber (1990) International Journal of Quantum Chemistry, Vol. 24, pp. 167–186; A. Ulitsky & R. Elber (1990) Journal of Chemical Physics, Vol. 92, pp. 1510–1511]. The free energy surface, or potential of mean force, for the conformational transition of Ace‐(MeA)₉‐NMe is calculated from molecular dynamics simulations, and a method is presented for the decomposition of the free energy surface into the constituent energetic and entropic terms, via the calculation of the required temperature derivatives in situ. For the AMBER/OPLS model employed here, the conformational transition pathways each contain a single 3₁₀‐helical‐like transition state, and the transition state potential energy relative to the 3₁₀‐conformation is 3 kcal/mol, independent of peptide length. Entropic stabilization in the barrier region significantly lowers the activation free energies for the forward and reverse transitions from the estimates of the barrier heights based simply on potential energy alone. © 1994 John Wiley & Sons, Inc.