Molecular simulation of multistate peptide dynamics: A comparison between microsecond timescale sampling and multiple shorter trajectories

Molecular dynamics simulations of the RN24 peptide, which includes a diverse set of structurally heterogeneous states, are carried out in explicit solvent. Two approaches are employed and compared directly under identical simulation conditions. Specifically, we examine sampling by two individual long trajectories (microsecond timescale) and many shorter (MS) uncoupled trajectories. Statistical analysis of the structural properties indicates a qualitative agreement between these approaches. Microsecond timescale sampling gives large uncertainties on most structural metrics, while the shorter timescale of MS simulations results in slight structural memory for beta‐structure starting states. Additionally, MS sampling detects numerous transitions on a relatively short timescale that are not observed in microsecond sampling, while long simulations allow for detection of a few transitions on significantly longer timescales. A correlation between the complex free energy landscape and the kinetics of the equilibrium is highlighted by principal component analysis on both simulation sets. This report highlights the increased precision of the MS approach when studying the kinetics of complex conformational change, while revealing the complementary insight and qualitative agreement offered by far fewer individual simulations on significantly longer timescales. © 2008 Wiley Periodicals, Inc. J Comput Chem, 2008

[1]  Eric J. Sorin,et al.  β-hairpin folding simulations in atomistic detail using an implicit solvent model1 , 2001 .

[2]  A. Fersht.  On the simulation of protein folding by short time scale molecular dynamics and distributed computing , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[3]  V. Pande,et al.  Calculation of the distribution of eigenvalues and eigenvectors in Markovian state models for molecular dynamics. , 2007, The Journal of chemical physics.

[4]  J Tirado-Rives,et al.  Molecular dynamics simulations of the unfolding of an alpha-helical analogue of ribonuclease A S-peptide in water. , 1991, Biochemistry.

[5]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[6]  P E Wright,et al.  1H NMR studies of the solution conformations of an analogue of the C-peptide of ribonuclease A. , 1989, Biochemistry.

[7]  J. Onuchic,et al.  Folding a protein in a computer: An atomic description of the folding/unfolding of protein A , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[8]  Eric J. Sorin,et al.  Insights into nucleic acid conformational dynamics from massively parallel stochastic simulations. , 2003, Biophysical journal.

[9]  J. E. Brown,et al.  Helix-coil transition of the isolated amino terminus of ribonuclease. , 1971, Biochemistry.

[10]  Vijay S. Pande,et al.  Empirical force‐field assessment: The interplay between backbone torsions and noncovalent term scaling , 2005, J. Comput. Chem..

[11]  C. Brooks,et al.  Exploring Flory's isolated-pair hypothesis: Statistical mechanics of helix–coil transitions in polyalanine and the C-peptide from RNase A , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[12]  Vijay S. Pande,et al.  Screen Savers of the World Unite! , 2000, Science.

[13]  W. Kabsch,et al.  Dictionary of protein secondary structure: Pattern recognition of hydrogen‐bonded and geometrical features , 1983, Biopolymers.

[14]  B. Hess.  Determining the shear viscosity of model liquids from molecular dynamics simulations , 2002 .

[15]  F. Young.  Biochemistry , 1955, The Indian Medical Gazette.

[16]  Jeremy C. Smith,et al.  Hierarchical analysis of conformational dynamics in biomolecules: transition networks of metastable states. , 2007, The Journal of chemical physics.

[17]  R. Rosenfeld.  Nature , 2009, Otolaryngology--head and neck surgery : official journal of American Academy of Otolaryngology-Head and Neck Surgery.

[18]  Berk Hess,et al.  LINCS: A linear constraint solver for molecular simulations , 1997 .

[19]  Eric J. Sorin,et al.  Exploring the helix-coil transition via all-atom equilibrium ensemble simulations. , 2005, Biophysical journal.

[20]  Alexander D. MacKerell,et al.  Force field influence on the observation of π-helical protein structures in molecular dynamics simulations , 2003 .

[21]  V. Pande,et al.  Does water play a structural role in the folding of small nucleic acids? , 2005, Biophysical journal.

[22]  R Reichelt,et al.  Scanning force microscopy at the air-water interface of an air bubble coated with pulmonary surfactant. , 2002, Biophysical journal.

[23]  Vijay S Pande,et al.  Validation of Markov state models using Shannon's entropy. , 2006, The Journal of chemical physics.

[24]  P. Kollman,et al.  Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models , 1992 .

[25]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[26]  M. Karplus,et al.  Solution conformations and thermodynamics of structured peptides: molecular dynamics simulation with an implicit solvation model. , 1998, Journal of molecular biology.

[27]  X. Daura,et al.  Parametrization of aliphatic CHn united atoms of GROMOS96 force field , 1998 .

[28]  Berk Hess,et al.  Improving efficiency of large time‐scale molecular dynamics simulations of hydrogen‐rich systems , 1999, Journal of computational chemistry.

[29]  Vijay S Pande,et al.  Using path sampling to build better Markovian state models: predicting the folding rate and mechanism of a tryptophan zipper beta hairpin. , 2004, The Journal of chemical physics.

[30]  William Swope,et al.  Describing Protein Folding Kinetics by Molecular Dynamics Simulations. 1. Theory , 2004 .

[31]  Eric J. Sorin,et al.  Simulations of the role of water in the protein-folding mechanism. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[32]  P. Kollman,et al.  Pathways to a protein folding intermediate observed in a 1-microsecond simulation in aqueous solution. , 1998, Science.

[33]  A. Cavalli,et al.  Analysis of the distributed computing approach applied to the folding of a small β peptide , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[34]  Lisa J. Lapidus,et al.  Fast kinetics and mechanisms in protein folding. , 2000, Annual review of biophysics and biomolecular structure.

[35]  P. S. Kim,et al.  A salt bridge stabilizes the helix formed by isolated C-peptide of RNase A. , 1982, Proceedings of the National Academy of Sciences of the United States of America.

[36]  D. van der Spoel,et al.  GROMACS: A message-passing parallel molecular dynamics implementation , 1995 .

[37]  V. Muñoz,et al.  Kinetics and Dynamics of Loops, α-Helices, β-Hairpins, and Fast-Folding Proteins , 1998 .

[38]  Michael R. Shirts,et al.  Atomistic protein folding simulations on the submillisecond time scale using worldwide distributed computing. , 2003, Biopolymers.

[39]  V. Pande,et al.  Absolute comparison of simulated and experimental protein-folding dynamics , 2002, Nature.

[40]  Luca Monticelli,et al.  Mechanism of helix nucleation and propagation: microscopic view from microsecond time scale MD simulations. , 2005, The journal of physical chemistry. B.

[41]  K. Dill,et al.  Automatic discovery of metastable states for the construction of Markov models of macromolecular conformational dynamics. , 2007, The Journal of chemical physics.

[42]  Robert L. Baldwin,et al.  Tests of the helix dipole model for stabilization of α-helices , 1987, Nature.

[43]  M. Shirts,et al.  Corrigendum to “Simulation of Folding of a Small Alpha-helical Protein in Atomistic Detail using Worldwide-distributed Computing” , 2002 .

[44]  Christian Bartels,et al.  Solution conformations of structured peptides: continuum electrostatics versus distance-dependent dielectric functions , 1999 .

[45]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

[46]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .