A molecular simulation protocol to avoid sampling redundancy and discover new states

作者: Marco Bacci , Andreas Vitalis , Amedeo Caflisch

DOI: 10.1016/J.BBAGEN.2014.08.013

关键词:

摘要: Abstract Background For biomacromolecules or their assemblies, experimental knowledge is often restricted to specific states. Ambiguity pervades simulations of these complex systems because there no prior relevant phase space domains, and sampling recurrence difficult achieve. In molecular dynamics methods, ruggedness the free energy surface exacerbates this problem by slowing down unbiased exploration space. Sampling inefficient if dwell times in metastable states are large. Methods We suggest a heuristic algorithm terminate reseed trajectories run multiple copies parallel. It uses recent method order snapshots, which provides notions “interesting” “unique” for individual simulations. define criteria guide reseeding runs from more points they sample overlapping regions Results Using pedagogical example an α-helical peptide, approach demonstrated amplify rate discover not found conventional schemes. Evidence provided that accurate kinetics pathways can be extracted Conclusions The method, termed PIGS Progress Index Guided Sampling, proceeds unsupervised fashion, scalable, benefits synergistically larger numbers replicas. confirm underlying ideas appropriate sufficient enhance sampling. General Significance simulations, errors caused exploring domains always unquantifiable arbitrarily Our protocol adds toolkit available researchers reducing types errors. This article part Special Issue entitled “Recent Developments Molecular Dynamics.”

参考文章(71)
Andreas Vitalis, Rohit V. Pappu, Chapter 3 Methods for Monte Carlo Simulations of Biomacromolecules Annual Reports in Computational Chemistry. ,vol. 5, pp. 49- 76 ,(2009) , 10.1016/S1574-1400(09)00503-9
Alan Grossfield, Daniel M. Zuckerman, Chapter 2 Quantifying Uncertainty and Sampling Quality in Biomolecular Simulations Annual Reports in Computational Chemistry. ,vol. 5, pp. 23- 48 ,(2009) , 10.1016/S1574-1400(09)00502-7
S. Muff, A. Caflisch, ETNA: Equilibrium Transitions Network and Arrhenius Equation for Extracting Folding Kinetics from REMD Simulations Journal of Physical Chemistry B. ,vol. 113, pp. 3218- 3226 ,(2009) , 10.1021/JP807261H
Andreas Vitalis, Rohit V. Pappu, ABSINTH: A new continuum solvation model for simulations of polypeptides in aqueous solutions Journal of Computational Chemistry. ,vol. 30, pp. 673- 699 ,(2009) , 10.1002/JCC.21005
Eric Vanden-Eijnden, Maddalena Venturoli, Markovian milestoning with Voronoi tessellations. Journal of Chemical Physics. ,vol. 130, pp. 194101- 194101 ,(2009) , 10.1063/1.3129843
Alan Grossfield, Scott E. Feller, Michael C. Pitman, Convergence of molecular dynamics simulations of membrane proteins Proteins: Structure, Function, and Bioinformatics. ,vol. 67, pp. 31- 40 ,(2007) , 10.1002/PROT.21308
Nina Singhal, Christopher D. Snow, Vijay S. Pande, Using path sampling to build better Markovian state models: Predicting the folding rate and mechanism of a tryptophan zipper beta hairpin Journal of Chemical Physics. ,vol. 121, pp. 415- 425 ,(2004) , 10.1063/1.1738647
Régis Pomès, J.Andrew McCammon, Mass and step length optimization for the calculation of equilibrium properties by molecular dynamics simulation Chemical Physics Letters. ,vol. 166, pp. 425- 428 ,(1990) , 10.1016/0009-2614(90)85055-H
Frank Noé, Christof Schütte, Eric Vanden-Eijnden, Lothar Reich, Thomas R. Weikl, Constructing the equilibrium ensemble of folding pathways from short off-equilibrium simulations Proceedings of the National Academy of Sciences of the United States of America. ,vol. 106, pp. 19011- 19016 ,(2009) , 10.1073/PNAS.0905466106
C. Komalapriya, M. Thiel, M. C. Romano, N. Marwan, U. Schwarz, J. Kurths, Reconstruction of a system's dynamics from short trajectories Physical Review E. ,vol. 78, pp. 066217- ,(2008) , 10.1103/PHYSREVE.78.066217