Bernhard Jaun

Folding and Unfolding of Two Mixed α/β Peptides

We present a molecular dynamics simulation study of two peptides containing α‐ and β‐amino acid residues. According to experiment, the two peptides differ in the dominant fold when solvated in methanol: one shows a helical fold, the other a β hairpin. The simulations at 300 and 340 K were done by starting from a NMR spectroscopic model structure and from an extended (denatured) structure. The typical structural features of the two peptides are reproduced and a folding/unfolding equilibrium is observed on the nanosecond timescale at 300 K. Analysis of proton–proton NOE distance bounds and backbone 3J coupling constants gives results consistent with the experimental data. We conclude that our simulations are complementary to the experiments by providing detailed information on the conformational distributions.

Interpreting Experimental Data by Using Molecular Simulation Instead of Model Building

A proper description of the conformational equilibrium of polypeptides or proteins is essential for a correct description of their function. The conformational ensembles from 16 molecular dynamic simulations of two beta- heptapeptides were used to interpret the primary NMR data, which were also compared to a set of NMR model structures (see graphic).One of the most used spectroscopic techniques for resolving the structure of a biomolecule, such as a protein or peptide, is NMR spectroscopy. Because only NMR signal intensities and frequencies are measured in the experiment, a conformational interpretation of the primary data, that is, measured data, is not straightforward, especially for flexible molecules. It is hampered by the occurrence of conformational and/or time-averaging, by insufficient number of experimental data and by insufficient accuracy of experimental data. All three problematic aspects of structure refinement based on NMR nuclear Overhauser effect (NOE) intensities and (3)J coupling data are illustrated by using two beta-heptapeptides in methanol as an example. We have performed 16 molecular dynamics (MD) simulations between 20 to 100 ns in length of unrestrained and NOE distance-restrained cases (instantaneous and time-averaged) of two beta-heptapeptides with a central beta-HAla(alpha-OH) amino acid in methanol at two different temperatures using two different GROMOS force-field parameter sets, 45 A3 and 53 A6. The created conformational ensembles were used to interpret the primary NMR data on these molecules. They also were compared to a set of NMR model structures derived by single-structure refinement in vacuo by using standard techniques. It is shown that the conformational interpretation of measured experimental data can be significantly improved by using unrestrained, instantaneous and time-averaged restrained MD simulations of the peptides by using a thermodynamically calibrated force field and by explicitly including solvent degrees of freedom.