Computing ensembles of transitions with molecular dynamics simulations.

A molecular understanding of conformational change is important for connecting structure and function. Without the ability to sample on the meaningful large-scale conformational changes, the ability to infer biological function and to understand the effect of mutations and changes in environment is not possible. Our Dynamic Importance Sampling method (DIMS), part of the CHARMM simulation package, is a method that enables sampling over ensembles of transition intermediates. This chapter outlines the context for the method and the usage within the program.

[1]  A. Seth,et al.  Granger causality and transfer entropy are equivalent for Gaussian variables. , 2009, Physical review letters.

[2]  V. Pande,et al.  Can conformational change be described by only a few normal modes? , 2006, Biophysical journal.

[3]  H. Kamberaj,et al.  Extracting the causality of correlated motions from molecular dynamics simulations. , 2009, Biophysical journal.

[4]  D. Leahy,et al.  Molecular dynamics simulations of transitions for ECD epidermal growth factor receptors show key differences between human and drosophila forms of the receptors , 2013, Proteins.

[5]  Lars Onsager,et al.  Fluctuations and Irreversible Processes , 1953 .

[6]  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.

[7]  R. Elber,et al.  Reaction path study of conformational transitions in flexible systems: Applications to peptides , 1990 .

[8]  Michael R. Shirts,et al.  Statistically optimal analysis of samples from multiple equilibrium states. , 2008, The Journal of chemical physics.

[9]  Oliver F. Lange,et al.  Can principal components yield a dimension reduced description of protein dynamics on long time scales? , 2006, The journal of physical chemistry. B.

[10]  Mark Gerstein,et al.  MolMovDB: analysis and visualization of conformational change and structural flexibility , 2003, Nucleic Acids Res..

[11]  T. Schreiber,et al.  Information transfer in continuous processes , 2002 .

[12]  T. Woolf,et al.  Rapid Determination of Multiple Reaction Pathways in Molecular Systems: The Soft-Ratcheting Algorithm , 2002, physics/0209098.

[13]  García,et al.  Large-amplitude nonlinear motions in proteins. , 1992, Physical review letters.

[14]  G. Chirikjian,et al.  Elastic models of conformational transitions in macromolecules. , 2002, Journal of molecular graphics & modelling.

[15]  Thomas B. Woolf,et al.  Multiple pathways in conformational transitions of the alanine dipeptide: An application of dynamic importance sampling , 2006, J. Comput. Chem..

[16]  Eric Vanden-Eijnden,et al.  Transition pathways in complex systems: Application of the finite-temperature string method to the alanine dipeptide. , 2005, The Journal of chemical physics.

[17]  A. Laio,et al.  Escaping free-energy minima , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[18]  Thomas B. Woolf,et al.  Path corrected functionals of stochastic trajectories: towards relative free energy and reaction coordinate calculations , 1998 .

[19]  Bradley M Dickson,et al.  Pitfalls of choosing an order parameter for rare event calculations. , 2009, The Journal of chemical physics.

[20]  A. Seth,et al.  Multivariate Granger causality and generalized variance. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[21]  Juan R Perilla,et al.  Computing ensembles of transitions from stable states: Dynamic importance sampling , 2011, J. Comput. Chem..

[22]  Gennady V Miloshevsky,et al.  Open-state conformation of the KcsA K+ channel: Monte Carlo normal mode following simulations. , 2007, Structure.

[23]  Christophe Chipot,et al.  Good practices in free-energy calculations. , 2010, The journal of physical chemistry. B.

[24]  Elizabeth J. Denning,et al.  Cooperative nature of gating transitions in K+ channels as seen from dynamic importance sampling calculations , 2010, Proteins.

[25]  Martin Karplus,et al.  Simulation of conformational transitions by the restricted perturbation-targeted molecular dynamics method. , 2005, The Journal of chemical physics.

[26]  D M Zuckerman,et al.  Efficient dynamic importance sampling of rare events in one dimension. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[27]  He Huang,et al.  A Comparison of Three Perturbation Molecular Dynamics Methods for Modeling Conformational Transitions. , 2009, Journal of chemical theory and computation.

[28]  Sebastian Doniach,et al.  Simulation of protein folding by reaction path annealing , 2001 .

[29]  G. Chirikjian,et al.  Efficient generation of feasible pathways for protein conformational transitions. , 2002, Biophysical journal.

[30]  C. Chennubhotla,et al.  Intrinsic dynamics of enzymes in the unbound state and relation to allosteric regulation. , 2007, Current opinion in structural biology.

[31]  Martin J. Field,et al.  A temperature-dependent nudged-elastic-band algorithm , 2003 .

[32]  J. Andrew McCammon,et al.  Substrate Induced Population Shifts and Stochastic Gating in the PBCV-1 mRNA Capping Enzyme , 2009, Journal of the American Chemical Society.

[33]  Giovanni Ciccotti,et al.  Book Review: Classical and Quantum Dynamics in Condensed Phase Simulations , 1998 .

[34]  Hyun-soo Cho,et al.  EGF activates its receptor by removing interactions that autoinhibit ectodomain dimerization. , 2003, Molecular cell.

[35]  Kevin L. Shaw,et al.  Linear extrapolation method of analyzing solvent denaturation curves , 2000, Proteins.

[36]  David Chandler,et al.  Statistical mechanics of chemical equilibria and intramolecular structures of nonrigid molecules in condensed phases , 1976 .

[37]  M. Karplus,et al.  Dynamics of folded proteins , 1977, Nature.

[38]  Ron Elber,et al.  Temperature dependent reaction coordinates , 2000 .

[39]  Milan Palus,et al.  Detecting couplings between interacting oscillators with time-varying basic frequencies: instantaneous wavelet bispectrum and information theoretic approach. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[40]  Ron Elber,et al.  Yet another look at the steepest descent path , 1997 .

[41]  Joseph A. Bank,et al.  Supporting Online Material Materials and Methods Figs. S1 to S10 Table S1 References Movies S1 to S3 Atomic-level Characterization of the Structural Dynamics of Proteins , 2022 .

[42]  G. Ciccotti,et al.  String method in collective variables: minimum free energy paths and isocommittor surfaces. , 2006, The Journal of chemical physics.

[43]  Schreiber,et al.  Nonlinear noise reduction: A case study on experimental data. , 1993, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[44]  Ron Elber,et al.  A new technique to calculate steepest descent paths in flexible polyatomic systems , 1990 .

[45]  D. Wales,et al.  A doubly nudged elastic band method for finding transition states. , 2004, The Journal of chemical physics.

[46]  Massimo Marchi,et al.  Adiabatic bias molecular dynamics: A method to navigate the conformational space of complex molecular systems , 1999 .

[47]  Ron Elber,et al.  A method for determining reaction paths in large molecules: application to myoglobin , 1987 .

[48]  Ron Elber,et al.  Reaction path study of helix formation in tetrapeptides: Effect of side chains , 1991 .

[49]  N. Go,et al.  Investigating protein dynamics in collective coordinate space. , 1999, Current opinion in structural biology.

[50]  Kazuyuki Aihara,et al.  Identifying hidden common causes from bivariate time series: a method using recurrence plots. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[51]  J. Andrew McCammon,et al.  Generalized Langevin dynamics simulations with arbitrary time‐dependent memory kernels , 1983 .

[52]  Ron Elber,et al.  Kinetics of helix unfolding: molecular dynamics simulations with milestoning. , 2009, The journal of physical chemistry. A.

[53]  Y. Sanejouand,et al.  Building‐block approach for determining low‐frequency normal modes of macromolecules , 2000, Proteins.

[54]  Bin W. Zhang,et al.  Efficient and verified simulation of a path ensemble for conformational change in a united-residue model of calmodulin , 2007, Proceedings of the National Academy of Sciences.

[55]  Lawrence R. Pratt,et al.  A statistical method for identifying transition states in high dimensional problems , 1986 .

[56]  Ming Lei,et al.  Sampling protein conformations and pathways , 2004, J. Comput. Chem..

[57]  W. Wagner Unbiased Monte Carlo evaluation of certain functional integrals , 1987 .

[58]  H. Kantz,et al.  Nonlinear time series analysis , 1997 .

[59]  Stefan Fischer,et al.  Mechanism of a molecular valve in the halorhodopsin chloride pump. , 2005, Structure.

[60]  Peter G. Bolhuis,et al.  Transition path sampling of cavitation between molecular scale solvophobic surfaces , 2000 .

[61]  J. Straub,et al.  The MaxFlux algorithm for calculating variationally optimized reaction paths for conformational transitions in many body systems at finite temperature , 1997 .

[62]  G. Huber,et al.  Weighted-ensemble Brownian dynamics simulations for protein association reactions. , 1996, Biophysical journal.

[63]  M. Karplus,et al.  Conjugate peak refinement: an algorithm for finding reaction paths and accurate transition states in systems with many degrees of freedom , 1992 .

[64]  Examining the limits of time reweighting and Kramers' rate theory to obtain correct kinetics from accelerated molecular dynamics. , 2010, The Journal of chemical physics.

[65]  R. Elber,et al.  Self‐avoiding walk between two fixed points as a tool to calculate reaction paths in large molecular systems , 1990 .

[66]  T. Schlick,et al.  Biomolecular modeling and simulation: a field coming of age , 2011, Quarterly Reviews of Biophysics.

[67]  Martin Karplus,et al.  Large amplitude conformational change in proteins explored with a plastic network model: adenylate kinase. , 2005, Journal of molecular biology.

[68]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[69]  C. Dellago,et al.  Transition path sampling and the calculation of rate constants , 1998 .

[70]  Thomas Schreiber,et al.  Detecting and Analyzing Nonstationarity in a Time Series Using Nonlinear Cross Predictions , 1997, chao-dyn/9909044.

[71]  H. Jónsson,et al.  Nudged elastic band method for finding minimum energy paths of transitions , 1998 .

[72]  Sebastian Doniach,et al.  A comparative study of motor-protein motions by using a simple elastic-network model , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[73]  Jianpeng Ma,et al.  Usefulness and limitations of normal mode analysis in modeling dynamics of biomolecular complexes. , 2005, Structure.

[74]  A. Chakraborty,et al.  A growing string method for determining transition states: comparison to the nudged elastic band and string methods. , 2004, The Journal of chemical physics.

[75]  Schreiber,et al.  Measuring information transfer , 2000, Physical review letters.

[76]  M. Magnasco,et al.  Forced thermal ratchets. , 1993, Physical review letters.

[77]  M. Karplus,et al.  Determination of a transition state at atomic resolution from protein engineering data. , 2002, Journal of molecular biology.

[78]  T. Woolf,et al.  Towards the prediction of order parameters from molecular dynamics simulations in proteins. , 2012, The Journal of chemical physics.

[79]  R Elber,et al.  Reaction path study of conformational transitions and helix formation in a tetrapeptide. , 1989, Proceedings of the National Academy of Sciences of the United States of America.

[80]  Hernan F. Stamati,et al.  Application of nonlinear dimensionality reduction to characterize the conformational landscape of small peptides , 2010, Proteins.

[81]  David H Mathews,et al.  Nudged elastic band calculation of minimal energy paths for the conformational change of a GG non-canonical pair. , 2006, Journal of molecular biology.

[82]  T. Creighton Proteins: Structures and Molecular Properties , 1986 .

[83]  A. Lesk,et al.  Structural mechanisms for domain movements in proteins. , 1994, Biochemistry.

[84]  T. Woolf,et al.  DYNAMIC REACTION PATHS AND RATES THROUGH IMPORTANCE-SAMPLED STOCHASTIC DYNAMICS , 1999 .

[85]  Bernard R Brooks,et al.  Normal-modes-based prediction of protein conformational changes guided by distance constraints. , 2005, Biophysical journal.