Locally accessible conformations of proteins: Multiple molecular dynamics simulations of crambin

Citation for published article

Constitutive activation of MEK5 promotes a mesenchymal and migratory cell phenotype in triple negative breast cancer



Document Type

Journal Article

Publication Date


Publication Title

Protein Science


Conformational space, Dimensionality reduction, Ergodicity, Molecular dynamics simulation, Multiple minima, Potential energy surface, Principal components analysis, Sensitivity to initial conditions


Bayer School of Natural and Environmental Sciences


Multiple molecular dynamics (MD) simulations of crambin with different initial atomic velocities are used in sample conformations in the vicinity of the native structure. Individual trajectories of length up to 5 ns sample only a fraction of the conformational distribution generated by ten independent 120 ps trajectories at 300 K. The backbone atom conformational space distribution is analyzed using principal components analysis (PCA). Four different major conformational regions are found. In general, a trajectory samples only one region and few transitions between the regions are observed. Consequently, the averages of structural and dynamic properties over the ten trajectories differ significantly from those obtained from individual trajectories. The nature of the conformational sampling has important consequences for the utilization of MD simulations for a wide range of problems, such as comparisons with X-ray or NMR data. The overall average structure is significantly closer to the X-ray structure than any of the individual trajectory average structures. The high frequency (less than 10 ps) atomic fluctuations from the ten trajectories tend to be similar, but the lower frequency (100 ps) motions are different. To improve conformational sampling in molecular dynamics simulations of proteins, as in nucleic acids, multiple trajectories with different initial conditions should be used rather than a single long trajectory.