Discrepancies between conformational distributions of a polyalanine peptide in solution obtained from molecular dynamics force fields and amide I' band profiles.

Structural preferences in the unfolded state of peptides determined by molecular dynamics still contradict experimental data. A remedy in this regard has been suggested by MD simulations with an optimized Amber force field ff03* ( Best, R. Hummer, G. J. Phys. Chem. B 2009 , 113 , 9004 - 9015 ). The simulations yielded a statistical coil distribution for alanine which is at variance with recent experimental results. To check the validity of this distribution, we investigated the peptide H-A(5)W-OH, which with the exception of the additional terminal tryptophan is analogous to the peptide used to optimize the force fields ff03*. Electronic circular dichroism, vibrational circular dichroism, and infrared spectroscopy as well as J-coupling constants obtained from NMR experiments were used to derive the peptide's conformational ensemble. Additionally, Förster resonance energy transfer between the terminal chromophores of the fluorescently labeled peptide analogue H-Dbo-A(5)W-OH was used to determine its average length, from which the end-to-end distance of the unlabeled peptide was estimated. Qualitatively, the experimental (3)J(H(N),C(α)), VCD, and ECD indicated a preference of alanine for polyproline II-like conformations. The experimental (3)J(H(N),C(α)) for A(5)W closely resembles the constants obtained for A(5). In order to quantitatively relate the conformational distribution of A(5) obtained with the optimized AMBER ff03* force field to experimental data, the former was used to derive a distribution function which expressed the conformational ensemble as a mixture of polyproline II, β-strand, helical, and turn conformations. This model was found to satisfactorily reproduce all experimental J-coupling constants. We employed the model to calculate the amide I' profiles of the IR and vibrational circular dichroism spectrum of A(5)W, as well as the distance between the two terminal peptide carbonyls. This led to an underestimated negative VCD couplet and an overestimated distance between terminal carbonyl groups. In order to more accurately account for the experimental data, we changed the distribution parameters based on results recently obtained for the alanine-based tripeptides. The final model, which satisfactorily reproduced amide I' profiles, J-coupling constant, and the end-to-end distance of A(5)W, reinforces alanine's high structural preference for polyproline II. Our results suggest that distributions obtained from MD simulations suggesting a statistical coil-like distribution for alanine are still based on insufficiently accurate force fields.