Virtual atom representation of hydrogen bonds in minimal off-lattice models of α helices: effect on stability, cooperativity and kinetics

Background: The most conspicuous feature of a right-handed α helix is the presence of hydrogen bonds between the backbone carbonyl oxygen and NH groups along the chain. A simple off-lattice model that includes hydrogen bond interactions using virtual atoms is used to examine the stability, cooperativity and kinetics of the helix–coil transition.

Results: We have studied the thermodynamics (using multiple histogram method) and kinetics (by Brownian dynamics simulations) of 16-mer minimal off-lattice models of four-turn α-helix sequences. The carbonyl and NH groups are represented as virtual moieties located between two α-carbon atoms along the polypeptide chain. The characteristics of the native conformations of the model helices, such as the helical pitch and angular correlations, coincide with those found in real proteins. The transition from coil to helix is quite broad, which is typical of these finite-sized systems. The cooperativity, as measured by a dimensionless parameter, Ω _c, that takes into account the width and the slope of the transition curves, is enhanced when hydrogen bonds are taken into account. The value of Ω _c for our model is consistent with that inferred from experiment for an alanine-based helix-forming peptide. The folding time τ_F ranges from 6 to 1000 ns in the temperature range 0.7–1.9 T_F, where T_F is the helix–coil transition temperature. These values are in excellent agreement with the results from recent fast folding experiments. The temperature dependence of τ_F exhibits a nearly Arrhenius behavior. Thermally induced unfolding occurs on a time scale that is less than 40–170 ps depending on the final temperature. Our calculations also predict that, although τ_F can be altered by changes in the sequence, the dynamic range over which such changes take place is not as large as that predicted for β-turn formation.

Conclusions:Hydrogen bonds not only affect the stability of α-helix formation but also have profound influence on the kinetics. The excellent agreement between our calculations and experiments suggests that these models can be used to investigate the effects of sequence, temperature and viscosity on the helix–coil transition.