Reliable Approximation of Long Relaxation Timescales in Molecular Dynamics

Many interesting rare events in molecular systems, like ligand association, protein folding or conformational changes, occur on timescales that often are not accessible by direct numerical simulation. Therefore, rare event approximation approaches like interface sampling, Markov state model building, or advanced reaction coordinate-based free energy estimation have attracted huge attention recently. In this article we analyze the reliability of such approaches. How precise is an estimate of long relaxation timescales of molecular systems resulting from various forms of rare event approximation methods? Our results give a theoretical answer to this question by relating it with the transfer operator approach to molecular dynamics. By doing so we also allow for understanding deep connections between the different approaches.

[1]  Peter Deuflhard,et al.  Transfer Operator Approach to Conformational Dynamics in Biomolecular Systems , 2001 .

[2]  Frank Noé,et al.  Markov state models based on milestoning. , 2011, The Journal of chemical physics.

[3]  Vijay S. Pande,et al.  Everything you wanted to know about Markov State Models but were afraid to ask. , 2010, Methods.

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

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

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

[7]  C. Schütte,et al.  Effective dynamics along given reaction coordinates, and reaction rate theory. , 2016, Faraday discussions.

[8]  T. Lelièvre,et al.  Effective dynamics using conditional expectations , 2009, 0906.4865.

[9]  Christof Schütte,et al.  Estimating the Eigenvalue Error of Markov State Models , 2012, Multiscale Model. Simul..

[10]  Gabriel Stoltz,et al.  Langevin dynamics with constraints and computation of free energy differences , 2010, Math. Comput..

[11]  C. Schütte,et al.  Supplementary Information for “ Constructing the Equilibrium Ensemble of Folding Pathways from Short Off-Equilibrium Simulations ” , 2009 .

[12]  Frank Noé,et al.  Variational Approach to Molecular Kinetics. , 2014, Journal of chemical theory and computation.

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

[14]  C. Schütte,et al.  Model reduction algorithms for optimal control and importance sampling of diffusions , 2016 .

[15]  P. Deuflhard,et al.  A Direct Approach to Conformational Dynamics Based on Hybrid Monte Carlo , 1999 .

[16]  S. Thompson Advances in experimental medicine and biology , 1996 .

[17]  G. Ciccotti,et al.  Blue moon sampling, vectorial reaction coordinates, and unbiased constrained dynamics. , 2005, Chemphyschem : a European journal of chemical physics and physical chemistry.

[18]  手塚明則 Blue‐moon samplingによる変形経路の解析 , 2005 .

[19]  Christof Schütte,et al.  Metastability and Markov State Models in Molecular Dynamics Modeling, Analysis , 2016 .

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

[21]  I. Gyöngy Mimicking the one-dimensional marginal distributions of processes having an ito differential , 1986 .

[22]  Ioannis G. Kevrekidis,et al.  Equation-free: The computer-aided analysis of complex multiscale systems , 2004 .

[23]  Vicente Hernández,et al.  SLEPc: A scalable and flexible toolkit for the solution of eigenvalue problems , 2005, TOMS.

[24]  Frank Noé,et al.  An Introduction to Markov State Models and Their Application to Long Timescale Molecular Simulation , 2014, Advances in Experimental Medicine and Biology.

[25]  Toni Giorgino,et al.  Identification of slow molecular order parameters for Markov model construction. , 2013, The Journal of chemical physics.

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

[27]  Giovanni Samaey,et al.  Equation-free multiscale computation: algorithms and applications. , 2009, Annual review of physical chemistry.

[28]  Frank Noé,et al.  Variational tensor approach for approximating the rare-event kinetics of macromolecular systems. , 2016, The Journal of chemical physics.

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

[30]  G. Ciccotti,et al.  Projection of diffusions on submanifolds: Application to mean force computation , 2008 .

[31]  Eric Darve,et al.  Adaptive biasing force method for scalar and vector free energy calculations. , 2008, The Journal of chemical physics.

[32]  Jonathan C. Mattingly,et al.  Ergodicity for SDEs and approximations: locally Lipschitz vector fields and degenerate noise , 2002 .

[33]  Gary Froyland,et al.  A Computational Method to Extract Macroscopic Variables and Their Dynamics in Multiscale Systems , 2013, SIAM J. Appl. Dyn. Syst..

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

[35]  Peter G. Bolhuis,et al.  Investigating Rare Events by Transition Interface Sampling , 2004 .

[36]  Frank Noé,et al.  A Variational Approach to Modeling Slow Processes in Stochastic Dynamical Systems , 2012, Multiscale Model. Simul..

[37]  P. R. ten Wolde,et al.  Non-stationary forward flux sampling. , 2012, The Journal of chemical physics.

[38]  Frank Noé,et al.  On the Approximation Quality of Markov State Models , 2010, Multiscale Model. Simul..

[39]  C. W. Gear,et al.  Equation-Free, Coarse-Grained Multiscale Computation: Enabling Mocroscopic Simulators to Perform System-Level Analysis , 2003 .

[40]  Christof Schütte,et al.  Approximating Selected Non-dominant Timescales by Markov State Models , 2012 .