Tackling sampling challenges in biomolecular simulations.

Molecular dynamics (MD) simulations are a powerful tool to give an atomistic insight into the structure and dynamics of proteins. However, the time scales accessible in standard simulations, which often do not match those in which interesting biological processes occur, limit their predictive capabilities. Many advanced sampling techniques have been proposed over the years to overcome this limitation. This chapter focuses on metadynamics, a method based on the introduction of a time-dependent bias potential to accelerate sampling and recover equilibrium properties of a few descriptors that are able to capture the complexity of a process at a coarse-grained level. The theory of metadynamics and its combination with other popular sampling techniques such as the replica exchange method is briefly presented. Practical applications of these techniques to the study of the Trp-Cage miniprotein folding are also illustrated. The examples contain a guide for performing these calculations with PLUMED, a plugin to perform enhanced sampling simulations in combination with many popular MD codes.

[1]  M. Parrinello,et al.  Well-tempered metadynamics: a smoothly converging and tunable free-energy method. , 2008, Physical review letters.

[2]  Massimiliano Bonomi,et al.  Efficient Simulation of Explicitly Solvated Proteins in the Well-Tempered Ensemble. , 2012, Journal of chemical theory and computation.

[3]  G. Hummer,et al.  Thermostat artifacts in replica exchange molecular dynamics simulations. , 2009, Journal of chemical theory and computation.

[4]  Conrad C. Huang,et al.  UCSF Chimera—A visualization system for exploratory research and analysis , 2004, J. Comput. Chem..

[5]  M. Parrinello,et al.  Canonical sampling through velocity rescaling. , 2007, The Journal of chemical physics.

[6]  Vijay S. Pande,et al.  Folding@home: Lessons from eight years of volunteer distributed computing , 2009, 2009 IEEE International Symposium on Parallel & Distributed Processing.

[7]  Michele Parrinello,et al.  Simplifying the representation of complex free-energy landscapes using sketch-map , 2011, Proceedings of the National Academy of Sciences.

[8]  Sören Doose,et al.  A microscopic view of miniprotein folding: enhanced folding efficiency through formation of an intermediate. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[9]  George I Makhatadze,et al.  Unfolding thermodynamics of Trp-cage, a 20 residue miniprotein, studied by differential scanning calorimetry and circular dichroism spectroscopy. , 2007, Biochemistry.

[10]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[11]  John D. Hunter,et al.  Matplotlib: A 2D Graphics Environment , 2007, Computing in Science & Engineering.

[12]  C. Dellago,et al.  Transition Path Sampling and Other Advanced Simulation Techniques for Rare Events , 2009 .

[13]  U. Hansmann Parallel tempering algorithm for conformational studies of biological molecules , 1997, physics/9710041.

[14]  Alessandro Laio,et al.  A Kinetic Model of Trp-Cage Folding from Multiple Biased Molecular Dynamics Simulations , 2009, PLoS Comput. Biol..

[15]  Oliver Riordan,et al.  The inefficiency of re-weighted sampling and the curse of system size in high-order path integration , 2011, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

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

[17]  Michele Ceriotti,et al.  Solid-liquid interface free energy through metadynamics simulations , 2009, 0911.2615.

[18]  Vojtěch Spiwok,et al.  Metadynamics in the conformational space nonlinearly dimensionally reduced by Isomap. , 2011, The Journal of chemical physics.

[19]  A. Laio,et al.  Free-energy landscape for beta hairpin folding from combined parallel tempering and metadynamics. , 2006, Journal of the American Chemical Society.

[20]  J. W. Neidigh,et al.  Designing a 20-residue protein , 2002, Nature Structural Biology.

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

[22]  Massimiliano Bonomi,et al.  Linking well-tempered metadynamics simulations with experiments. , 2010, Biophysical journal.

[23]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[24]  R. A. Leibler,et al.  On Information and Sufficiency , 1951 .

[25]  P. Bolhuis,et al.  Sampling the multiple folding mechanisms of Trp-cage in explicit solvent , 2006, Proceedings of the National Academy of Sciences.

[26]  R. Zhou Trp-cage: Folding free energy landscape in explicit water , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[27]  Michele Parrinello,et al.  Replica Temperatures for Uniform Exchange and Efficient Roundtrip Times in Explicit Solvent Parallel Tempering Simulations. , 2011, Journal of chemical theory and computation.

[28]  H. Nymeyer,et al.  Replica exchange simulation of reversible folding/unfolding of the Trp-cage miniprotein in explicit solvent: on the structure and possible role of internal water. , 2007, Journal of structural biology.

[29]  A. Roitberg,et al.  Smaller and faster: the 20-residue Trp-cage protein folds in 4 micros. , 2002, Journal of the American Chemical Society.

[30]  S. Nosé A unified formulation of the constant temperature molecular dynamics methods , 1984 .

[31]  Daniel J. Sindhikara,et al.  Exchange Often and Properly in Replica Exchange Molecular Dynamics. , 2010, Journal of chemical theory and computation.

[32]  A. Pohorille,et al.  Free energy calculations : theory and applications in chemistry and biology , 2007 .

[33]  Motonori Ota,et al.  Phylogeny of protein-folding trajectories reveals a unique pathway to native structure. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[34]  Massimiliano Bonomi,et al.  PLUMED: A portable plugin for free-energy calculations with molecular dynamics , 2009, Comput. Phys. Commun..

[35]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[36]  M. Parrinello,et al.  Enhanced sampling in the well-tempered ensemble. , 2009, Physical review letters.

[37]  Massimiliano Bonomi,et al.  PLUMED 2: New feathers for an old bird , 2013, Comput. Phys. Commun..

[38]  Massimiliano Bonomi,et al.  Reconstructing the equilibrium Boltzmann distribution from well‐tempered metadynamics , 2009, J. Comput. Chem..

[39]  R. Dror,et al.  How Fast-Folding Proteins Fold , 2011, Science.

[40]  M. Parrinello,et al.  Metadynamics with Adaptive Gaussians. , 2012, Journal of chemical theory and computation.

[41]  Michael W Deem,et al.  Parallel tempering: theory, applications, and new perspectives. , 2005, Physical chemistry chemical physics : PCCP.

[42]  Francesco Luigi Gervasio,et al.  From A to B in free energy space. , 2007, The Journal of chemical physics.

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

[44]  S. Asher,et al.  UV-resonance raman thermal unfolding study of Trp-cage shows that it is not a simple two-state miniprotein. , 2005, Journal of the American Chemical Society.

[45]  Francesco Luigi Gervasio,et al.  Comparing the Efficiency of Biased and Unbiased Molecular Dynamics in Reconstructing the Free Energy Landscape of Met-Enkephalin , 2010 .

[46]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[47]  A. Laio,et al.  A bias-exchange approach to protein folding. , 2007, The journal of physical chemistry. B.

[48]  Michele Parrinello,et al.  Effect of urea on the β-hairpin conformational ensemble and protein denaturation mechanism. , 2011, Journal of the American Chemical Society.