Enhancing Important Fluctuations: Rare Events and Metadynamics from a Conceptual Viewpoint.

Atomistic simulations play a central role in many fields of science. However, their usefulness is often limited by the fact that many systems are characterized by several metastable states separated by high barriers, leading to kinetic bottlenecks. Transitions between metastable states are thus rare events that occur on significantly longer timescales than one can simulate in practice. Numerous enhanced sampling methods have been introduced to alleviate this timescale problem, including methods based on identifying a few crucial order parameters or collective variables and enhancing the sampling of these variables. Metadynamics is one such method that has proven successful in a great variety of fields. Here we review the conceptual and theoretical foundations of metadynamics. As demonstrated, metadynamics is not just a practical tool but can also be considered an important development in the theory of statistical mechanics.

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

[2]  Eric F Darve,et al.  Calculating free energies using average force , 2001 .

[3]  M. Buehler,et al.  A robust nanoscale experimental quantification of fracture energy in a bilayer material system , 2014, Proceedings of the National Academy of Sciences.

[4]  Eric T. Kim,et al.  How does a drug molecule find its target binding site? , 2011, Journal of the American Chemical Society.

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

[6]  B. Berne,et al.  Replica exchange with solute tempering: a method for sampling biological systems in explicit water. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[7]  A. Laio,et al.  Assessing the accuracy of metadynamics. , 2005, The journal of physical chemistry. B.

[8]  W. E,et al.  Finite temperature string method for the study of rare events. , 2002, Journal of Physical Chemistry B.

[9]  T. Bredow,et al.  Molecular dynamics investigation of oxygen vacancy diffusion in rutile. , 2005, Physical chemistry chemical physics : PCCP.

[10]  W. Andreoni,et al.  Capture and Release of CO₂ in Monoethanolamine Aqueous Solutions: New Insights from First-Principles Reaction Dynamics. , 2015, Journal of chemical theory and computation.

[11]  Michele Parrinello,et al.  Well-tempered metadynamics converges asymptotically. , 2014, Physical review letters.

[12]  Giovanni Bussi,et al.  Enhanced Conformational Sampling Using Replica Exchange with Collective-Variable Tempering , 2015, Journal of chemical theory and computation.

[13]  Giovanni Bussi,et al.  Enhanced Sampling in Molecular Dynamics Using Metadynamics, Replica-Exchange, and Temperature-Acceleration , 2013, Entropy.

[14]  Fabio Pietrucci,et al.  Graph theory meets ab initio molecular dynamics: atomic structures and transformations at the nanoscale. , 2011, Physical review letters.

[15]  Marino Arroyo,et al.  Modeling and enhanced sampling of molecular systems with smooth and nonlinear data-driven collective variables. , 2013, The Journal of chemical physics.

[16]  Rosalind J Allen,et al.  Forward flux sampling for rare event simulations , 2009, Journal of physics. Condensed matter : an Institute of Physics journal.

[17]  D. Haar,et al.  Statistical Physics , 1971, Nature.

[18]  R. Elber,et al.  Computing time scales from reaction coordinates by milestoning. , 2004, The Journal of chemical physics.

[19]  Andrew E. Torda,et al.  Local elevation: A method for improving the searching properties of molecular dynamics simulation , 1994, J. Comput. Aided Mol. Des..

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

[21]  C. Dellago,et al.  Reaction coordinates of biomolecular isomerization. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[22]  Axel van de Walle,et al.  Accelerated molecular dynamics through stochastic iterations and collective variable based basin identification , 2013 .

[23]  B. Berne,et al.  Replica exchange with solute scaling: a more efficient version of replica exchange with solute tempering (REST2). , 2011, The journal of physical chemistry. B.

[24]  Michele Parrinello,et al.  Variational approach to enhanced sampling and free energy calculations. , 2014, Physical review letters.

[25]  Juan J. de Pablo,et al.  Flux Tempered Metadynamics , 2011 .

[26]  Michele Parrinello,et al.  Path Integral Metadynamics. , 2015, Journal of chemical theory and computation.

[27]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[28]  M. Parrinello,et al.  A perturbative solution to metadynamics ordinary differential equation. , 2015, The Journal of chemical physics.

[29]  M. Parrinello,et al.  A time-independent free energy estimator for metadynamics. , 2015, The journal of physical chemistry. B.

[30]  Pratyush Tiwary,et al.  Role of water and steric constraints in the kinetics of cavity–ligand unbinding , 2015, Proceedings of the National Academy of Sciences.

[31]  John L. Klepeis,et al.  Anton, a special-purpose machine for molecular dynamics simulation , 2007, ISCA '07.

[32]  Juan de Pablo,et al.  A boundary correction algorithm for metadynamics in multiple dimensions. , 2013, The Journal of chemical physics.

[33]  Michele Parrinello,et al.  A self-learning algorithm for biased molecular dynamics , 2010, Proceedings of the National Academy of Sciences.

[34]  Eric Vanden-Eijnden,et al.  Transition-path theory and path-finding algorithms for the study of rare events. , 2010, Annual review of physical chemistry.

[35]  John E. Straub,et al.  Classical and modern methods in reaction rate theory , 1988 .

[36]  Preston Moore,et al.  Metadynamics as a tool for exploring free energy landscapes of chemical reactions. , 2006, Accounts of chemical research.

[37]  R. Zwanzig Nonequilibrium statistical mechanics , 2001, Physics Subject Headings (PhySH).

[38]  Oliver F. Lange,et al.  Application of Enhanced Sampling Monte Carlo Methods for High-Resolution Protein-Protein Docking in Rosetta , 2015, PloS one.

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

[40]  J. Pfaendtner,et al.  Structural Effects of Methionine Oxidation on Isolated Subdomains of Human Fibrin D and αC Regions , 2014, PloS one.

[41]  A. Voter A method for accelerating the molecular dynamics simulation of infrequent events , 1997 .

[42]  Axel van de Walle,et al.  Hybrid deterministic and stochastic approach for efficient atomistic simulations at long time scales , 2011 .

[43]  Juan J de Pablo,et al.  Density of states-based molecular simulations. , 2012, Annual review of chemical and biomolecular engineering.

[44]  A. Laio,et al.  Characterization of the free-energy landscapes of proteins by NMR-guided metadynamics , 2013, Proceedings of the National Academy of Sciences.

[45]  H. Kushner,et al.  Stochastic Approximation and Recursive Algorithms and Applications , 2003 .

[46]  A. Voter,et al.  Chapter 4 Accelerated Molecular Dynamics Methods: Introduction and Recent Developments , 2009 .

[47]  P. Hänggi,et al.  Reaction-rate theory: fifty years after Kramers , 1990 .

[48]  B. Berne,et al.  A rapid method for determining rate constants by molecular dynamics , 1985 .

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

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

[51]  C. Landes,et al.  Measuring the Hydrodynamic Size of Nanoparticles Using Fluctuation Correlation Spectroscopy. , 2016, Annual review of physical chemistry.

[52]  Michele Parrinello,et al.  Using sketch-map coordinates to analyze and bias molecular dynamics simulations , 2012, Proceedings of the National Academy of Sciences.

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

[54]  Michele Parrinello,et al.  Well-Tempered Variational Approach to Enhanced Sampling. , 2015, Journal of chemical theory and computation.

[55]  D. Landau,et al.  Efficient, multiple-range random walk algorithm to calculate the density of states. , 2000, Physical review letters.

[56]  Shaohui Zheng,et al.  Enhanced sampling of chemical and biochemical reactions with metadynamics , 2015 .

[57]  M. Maggioni,et al.  Determination of reaction coordinates via locally scaled diffusion map. , 2011, The Journal of chemical physics.

[58]  M. Parrinello,et al.  From metadynamics to dynamics. , 2013, Physical review letters.

[59]  A. Cavalli,et al.  Molecular basis of cyclooxygenase enzymes (COXs) selective inhibition , 2010, Proceedings of the National Academy of Sciences.

[60]  R. Needs,et al.  Continuum variational and diffusion quantum Monte Carlo calculations , 2010, Journal of physics. Condensed matter : an Institute of Physics journal.

[61]  M. Klein,et al.  Free-energy landscape of ion-channel voltage-sensor–domain activation , 2014, Proceedings of the National Academy of Sciences of the United States of America.

[62]  Francesco Luigi Gervasio,et al.  New advances in metadynamics , 2012 .

[63]  Ronald R. Coifman,et al.  Diffusion Maps, Reduction Coordinates, and Low Dimensional Representation of Stochastic Systems , 2008, Multiscale Model. Simul..

[64]  G. Ciccotti,et al.  Constrained reaction coordinate dynamics for the simulation of rare events , 1989 .

[65]  Peter M. Kasson,et al.  GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit , 2013, Bioinform..

[66]  Eric Moulines,et al.  Non-strongly-convex smooth stochastic approximation with convergence rate O(1/n) , 2013, NIPS.

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

[68]  Michele Parrinello,et al.  Demonstrating the Transferability and the Descriptive Power of Sketch-Map. , 2013, Journal of chemical theory and computation.

[69]  Michele Parrinello,et al.  Probing the Unfolded Configurations of a β-Hairpin Using Sketch-Map. , 2015, Journal of chemical theory and computation.

[70]  Alessandro Laio,et al.  Reconstructing the density of states by history-dependent metadynamics. , 2004, Physical review letters.

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

[72]  E. Vanden-Eijnden,et al.  A temperature accelerated method for sampling free energy and determining reaction pathways in rare events simulations , 2006 .

[73]  Joan-Emma Shea,et al.  Transmission Coefficients, Committors, and Solvent Coordinates in Ion-Pair Dissociation. , 2014, Journal of chemical theory and computation.

[74]  M. Parrinello,et al.  The allosteric communication pathways in KIX domain of CBP , 2013, Proceedings of the National Academy of Sciences.

[75]  Bernd Ensing,et al.  Formamide hydrolysis in alkaline aqueous solution: insight from Ab initio metadynamics calculations. , 2006, Angewandte Chemie.

[76]  F L Gervasio,et al.  Efficient numerical reconstruction of protein folding kinetics with partial path sampling and pathlike variables. , 2013, Physical review letters.

[77]  G. Henkelman,et al.  Long time scale kinetic Monte Carlo simulations without lattice approximation and predefined event table , 2001 .

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

[79]  B. Ensing,et al.  Complex Reaction Environments and Competing Reaction Mechanisms in Zeolite Catalysis: Insights from Advanced Molecular Dynamics. , 2015, Chemistry.

[80]  Marissa G. Saunders,et al.  Unraveling the Mystery of ATP Hydrolysis in Actin Filaments , 2014, Journal of the American Chemical Society.

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

[82]  J. Iglesias-Fernández,et al.  The complete conformational free energy landscape of β-xylose reveals a two-fold catalytic itinerary for β-xylanases† †Electronic supplementary information (ESI) available. See DOI: 10.1039/c4sc02240h Click here for additional data file. , 2014, Chemical science.

[83]  Berend Smit,et al.  Understanding Molecular Simulation , 2001 .

[84]  Michele Parrinello,et al.  Variationally Optimized Free-Energy Flooding for Rate Calculation. , 2015, Physical review letters.

[85]  Giovanni Bussi,et al.  Hamiltonian replica exchange in GROMACS: a flexible implementation , 2013, 1307.5144.

[86]  Carsten Kutzner,et al.  Tackling Exascale Software Challenges in Molecular Dynamics Simulations with GROMACS , 2015, EASC.

[87]  James T. Hynes,et al.  The stable states picture of chemical reactions. II. Rate constants for condensed and gas phase reaction models , 1980 .

[88]  Alessandro Laio,et al.  Protein Folding and Ligand-Enzyme Binding from Bias-Exchange Metadynamics Simulations , 2012 .

[89]  Ilias Bilionis,et al.  Free energy computations by minimization of Kullback-Leibler divergence: An efficient adaptive biasing potential method for sparse representations , 2010, J. Comput. Phys..

[90]  Marco Mazzotti,et al.  Uncovering molecular details of urea crystal growth in the presence of additives. , 2012, Journal of the American Chemical Society.

[91]  Juan J de Pablo,et al.  Efficient Free Energy Calculation of Biomolecules from Diffusion-Biased Molecular Dynamics. , 2012, Journal of chemical theory and computation.

[92]  G Saladino,et al.  Assessing the Performance of Metadynamics and Path Variables in Predicting the Binding Free Energies of p38 Inhibitors. , 2012, Journal of chemical theory and computation.

[93]  David Chandler,et al.  Transition path sampling: throwing ropes over rough mountain passes, in the dark. , 2002, Annual review of physical chemistry.

[94]  D. Marx,et al.  The Janus-faced role of external forces in mechanochemical disulfide bond cleavage. , 2013, Nature chemistry.

[95]  Francesco Luigi Gervasio,et al.  The role of the peripheral anionic site and cation-pi interactions in the ligand penetration of the human AChE gorge. , 2005, Journal of the American Chemical Society.

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

[97]  Alessandro Laio,et al.  Efficient exploration of reactive potential energy surfaces using Car-Parrinello molecular dynamics. , 2003, Physical review letters.

[98]  M. Parrinello,et al.  Thermodynamical Description of a Quasi-First-Order Phase Transition from the Well-Tempered Ensemble. , 2013, Journal of chemical theory and computation.

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

[100]  N. Destainville,et al.  DNA denaturation bubbles: free-energy landscape and nucleation/closure rates. , 2014, The Journal of chemical physics.

[101]  Martin Karplus,et al.  Gaussian-mixture umbrella sampling. , 2009, The journal of physical chemistry. B.

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

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

[104]  A. Laio,et al.  Efficient reconstruction of complex free energy landscapes by multiple walkers metadynamics. , 2006, The journal of physical chemistry. B.

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

[106]  Jacques G. Amar,et al.  Accelerated molecular dynamics methods: introduction and recent developments , 2009 .

[107]  Julius Jellinek,et al.  Energy Landscapes: With Applications to Clusters, Biomolecules and Glasses , 2005 .

[108]  Jesse G. McDaniel,et al.  Next-Generation Force Fields from Symmetry-Adapted Perturbation Theory. , 2016, Annual review of physical chemistry.

[109]  Donald Hamelberg,et al.  Towards fast, rigorous and efficient conformational sampling of biomolecules: Advances in accelerated molecular dynamics. , 2015, Biochimica et biophysica acta.

[110]  A. Laio,et al.  Metadynamics: a method to simulate rare events and reconstruct the free energy in biophysics, chemistry and material science , 2008 .

[111]  W. Thiel,et al.  Solvent influence on cellulose 1,4-β-glycosidic bond cleavage: a molecular dynamics and metadynamics study. , 2015, Chemistry.

[112]  A. Laio,et al.  Equilibrium free energies from nonequilibrium metadynamics. , 2006, Physical Review Letters.

[113]  Xevi Biarnés,et al.  Molecular motions in drug design: the coming age of the metadynamics method , 2011, J. Comput. Aided Mol. Des..

[114]  Albert C. Pan,et al.  Molecular determinants of drug-receptor binding kinetics. , 2013, Drug discovery today.

[115]  Jacek Komasa,et al.  Explicitly Correlated Functions in Variational Calculations , 2003 .

[116]  R. Broglia,et al.  Exploring the protein G helix free‐energy surface by solute tempering metadynamics , 2007, Proteins.

[117]  Giovanni Bussi,et al.  Free‐Energy Calculations with Metadynamics: Theory and Practice , 2015 .

[118]  S. Takada,et al.  On the Hamiltonian replica exchange method for efficient sampling of biomolecular systems: Application to protein structure prediction , 2002 .

[119]  Michele Parrinello,et al.  Assessing the Reliability of the Dynamics Reconstructed from Metadynamics. , 2014, Journal of chemical theory and computation.

[120]  A. Cavalli,et al.  Investigating drug-target association and dissociation mechanisms using metadynamics-based algorithms. , 2015, Accounts of chemical research.

[121]  C. Clementi,et al.  Discovering mountain passes via torchlight: methods for the definition of reaction coordinates and pathways in complex macromolecular reactions. , 2013, Annual review of physical chemistry.

[122]  Artem R. Oganov,et al.  Systematic search for low-enthalpy sp3 carbon using evolutionary metadynamics , 2012, 1203.0706.

[123]  M. Parrinello,et al.  Targeting biomolecular flexibility with metadynamics. , 2010, Current opinion in structural biology.

[124]  Grant M. Rotskoff,et al.  Transition-Tempered Metadynamics: Robust, Convergent Metadynamics via On-the-Fly Transition Barrier Estimation. , 2014, Journal of chemical theory and computation.

[125]  Davide Provasi,et al.  Exploring molecular mechanisms of ligand recognition by opioid receptors with metadynamics. , 2009, Biochemistry.

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

[127]  Michele Ceriotti,et al.  Recognizing molecular patterns by machine learning: an agnostic structural definition of the hydrogen bond. , 2014, The Journal of chemical physics.

[128]  Rudolph A. Marcus,et al.  Electron transfer reactions in chemistry. Theory and experiment , 1993 .

[129]  Grubmüller,et al.  Predicting slow structural transitions in macromolecular systems: Conformational flooding. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[130]  D. Nesbitt,et al.  Biophysical Insights from Temperature-Dependent Single-Molecule Förster Resonance Energy Transfer. , 2016, Annual review of physical chemistry.

[131]  Michele Parrinello,et al.  Enhanced, targeted sampling of high-dimensional free-energy landscapes using variationally enhanced sampling, with an application to chignolin , 2016, Proceedings of the National Academy of Sciences.

[132]  Michele Parrinello,et al.  Metadynamics studies of crystal nucleation , 2015, IUCrJ.

[133]  Mark E Tuckerman,et al.  Locating landmarks on high-dimensional free energy surfaces , 2015, Proceedings of the National Academy of Sciences.

[134]  Shie Mannor,et al.  A Tutorial on the Cross-Entropy Method , 2005, Ann. Oper. Res..

[135]  Alessandro Laio,et al.  Metadynamics convergence law in a multidimensional system. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[136]  M. Parrinello,et al.  Kinetics of protein–ligand unbinding: Predicting pathways, rates, and rate-limiting steps , 2015, Proceedings of the National Academy of Sciences.

[137]  A. Voter,et al.  Extending the Time Scale in Atomistic Simulation of Materials Annual Re-views in Materials Research , 2002 .

[138]  D. Huse,et al.  Optimizing the ensemble for equilibration in broad-histogram Monte Carlo simulations. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[139]  Matthias U. Bohner,et al.  Nudged-elastic band used to find reaction coordinates based on the free energy. , 2014, The Journal of chemical physics.

[140]  M. Tuckerman Statistical Mechanics: Theory and Molecular Simulation , 2010 .

[141]  A. Voter Parallel replica method for dynamics of infrequent events , 1998 .

[142]  Berg,et al.  Multicanonical ensemble: A new approach to simulate first-order phase transitions. , 1992, Physical review letters.

[143]  J. VandeVondele,et al.  Canonical Adiabatic Free Energy Sampling (CAFES): A Novel Method for the Exploration of Free Energy Surfaces , 2002 .

[144]  M. Parrinello,et al.  Funnel metadynamics as accurate binding free-energy method , 2013, Proceedings of the National Academy of Sciences.

[145]  M. Parrinello,et al.  Combustion chemistry via metadynamics: benzyl decomposition revisited. , 2015, The journal of physical chemistry. A.

[146]  Frank Noé,et al.  Markov models of molecular kinetics: generation and validation. , 2011, The Journal of chemical physics.

[147]  F. Gervasio,et al.  New Insights into the Molecular Mechanism of E-Cadherin-Mediated Cell Adhesion by Free Energy Calculations. , 2015, Journal of chemical theory and computation.

[148]  Ulrich H E Hansmann,et al.  Global optimization by energy landscape paving. , 2002, Physical review letters.

[149]  M. Tuckerman,et al.  On the use of the adiabatic molecular dynamics technique in the calculation of free energy profiles , 2002 .

[150]  G. Hummer,et al.  Reaction coordinates and rates from transition paths. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[151]  A. Voter Hyperdynamics: Accelerated Molecular Dynamics of Infrequent Events , 1997 .

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

[153]  M Scott Shell,et al.  Coarse-graining errors and numerical optimization using a relative entropy framework. , 2011, The Journal of chemical physics.

[154]  G. Henkelman,et al.  A climbing image nudged elastic band method for finding saddle points and minimum energy paths , 2000 .

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

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