Positional preference of proline in α‐helices

Abstract
Conformational free energy calculations have been carried out for proline‐containing alanine‐based pentadecapeptides with the sequence Ac‐(Ala)n‐Pro‐(Ala)m‐NHMe, where n + m = 14, to figure out the positional preference of proline in α‐helices. The relative free energy of each peptide was calculated by subtracting the free energy of the extended conformation from that of the α‐helical one, which is used here as a measure of preference. The highest propensity is found for the peptide with proline at the N‐terminus (i.e., Ncap + 1 position), and the next propensities are found at Ncap, N' (Ncap — 1), and C' (Ccap + 1) positions. These computed results are reasonably consistent with the positional propensities estimated from X‐ray structures of proteins. The breaking in hydrogen bonds around proline is found to play a role in destabilizing α‐helical conformations, which, however, provides the favored hydration of the corresponding N‐H and C=O groups. The highest preference of proline at the beginning of α‐helix appears to be due to the favored electrostatic and nonbonded energies between two residues preceding proline and the intrinsic stability of α‐helical conformation of the proline residue itself as well as no disturbance in hydrogen bonds of α‐helix by proline. The average free energy change for the substitution of Ala by Pro in a α‐helix is computed to be 4.6 kcal/mol, which is in good agreement with the experimental value of ∼4 kcal/mol estimated for an oligopeptide dimer and proteins of barnase and T4 lysozyme.