Contour forward flux sampling: Sampling rare events along multiple collective variables.

Many rare event transitions involve multiple collective variables (CVs), and the most appropriate combination of CVs is generally unknown a priori. We thus introduce a new method, contour forward flux sampling (cFFS), to study rare events with multiple CVs simultaneously. cFFS places nonlinear interfaces on-the-fly from the collective progress of the simulations, without any prior knowledge of the energy landscape or appropriate combination of CVs. We demonstrate cFFS on analytical potential energy surfaces and a conformational change in alanine dipeptide.

[1]  C. Dellago,et al.  A string reaction coordinate for the folding of a polymer chain , 2015, Journal of physics. Condensed matter : an Institute of Physics journal.

[2]  Carmeline J. Dsilva,et al.  Systematic characterization of protein folding pathways using diffusion maps: application to Trp-cage miniprotein. , 2015, The Journal of chemical physics.

[3]  Lydia E Kavraki,et al.  Low-dimensional, free-energy landscapes of protein-folding reactions by nonlinear dimensionality reduction , 2006, Proc. Natl. Acad. Sci. USA.

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

[5]  Ernesto E. Borrero,et al.  Kinetics and mechanism of the unfolding native-to-loop transition of Trp-cage in explicit solvent via optimized forward flux sampling simulations. , 2010, The Journal of chemical physics.

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

[7]  Frank Noé,et al.  Time-lagged autoencoders: Deep learning of slow collective variables for molecular kinetics , 2017, The Journal of chemical physics.

[8]  A. Michaelides,et al.  Microscopic Mechanism and Kinetics of Ice Formation at Complex Interfaces: Zooming in on Kaolinite , 2016, The journal of physical chemistry letters.

[9]  Aaron R Dinner,et al.  Umbrella sampling for nonequilibrium processes. , 2007, The Journal of chemical physics.

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

[11]  Ryan S. DeFever,et al.  Nucleation mechanism of clathrate hydrates of water-soluble guest molecules. , 2017, The Journal of chemical physics.

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

[13]  P. Bolhuis,et al.  Rate constant and reaction coordinate of Trp-cage folding in explicit water. , 2008, Biophysical journal.

[14]  D. Frenkel,et al.  Simulating rare events in equilibrium or nonequilibrium stochastic systems. , 2005, The Journal of chemical physics.

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

[16]  T. S. V. Erp,et al.  Reaction rate calculation by parallel path swapping. , 2007, Physical review letters.

[17]  Hajime Tanaka,et al.  Crystal nucleation as the ordering of multiple order parameters. , 2016, The Journal of chemical physics.

[18]  Elaborating transition interface sampling methods , 2005, cond-mat/0405116.

[19]  Peter G. Bolhuis,et al.  A novel path sampling method for the calculation of rate constants , 2003 .

[20]  Pablo G Debenedetti,et al.  Direct calculation of ice homogeneous nucleation rate for a molecular model of water , 2015, Proceedings of the National Academy of Sciences.

[21]  Gerhard Hummer,et al.  Transition path sampling of rare events by shooting from the top. , 2017, The Journal of chemical physics.

[22]  Efficiency analysis of reaction rate calculation methods using analytical models I: The two-dimensional sharp barrier , 2006, cond-mat/0605260.

[23]  Ao Ma,et al.  Recent developments in methods for identifying reaction coordinates , 2014, Molecular simulation.

[24]  Christoph Dellago,et al.  Optimizing transition interface sampling simulations. , 2011, The Journal of chemical physics.

[25]  Wei Chen,et al.  Molecular enhanced sampling with autoencoders: On‐the‐fly collective variable discovery and accelerated free energy landscape exploration , 2017, J. Comput. Chem..

[26]  Amit Singer,et al.  Detecting intrinsic slow variables in stochastic dynamical systems by anisotropic diffusion maps , 2009, Proceedings of the National Academy of Sciences.

[27]  Jianhua Lin,et al.  Divergence measures based on the Shannon entropy , 1991, IEEE Trans. Inf. Theory.

[28]  E. Vanden-Eijnden,et al.  String method for the study of rare events , 2002, cond-mat/0205527.

[29]  J. Tenenbaum,et al.  A global geometric framework for nonlinear dimensionality reduction. , 2000, Science.

[30]  Charles H. Bennett,et al.  Molecular Dynamics and Transition State Theory: The Simulation of Infrequent Events , 1977 .

[31]  J. Araque,et al.  Transition path sampling and forward flux sampling. Applications to biological systems , 2009, Journal of physics. Condensed matter : an Institute of Physics journal.

[32]  Ericka Stricklin-Parker,et al.  Ann , 2005 .

[33]  Ioannis G. Kevrekidis,et al.  Nonlinear dimensionality reduction in molecular simulation: The diffusion map approach , 2011 .

[34]  A. Dinner,et al.  Nonequilibrium umbrella sampling in spaces of many order parameters. , 2009, The Journal of chemical physics.

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

[36]  Mohammad M. Sultan,et al.  Automated design of collective variables using supervised machine learning. , 2018, The Journal of chemical physics.

[37]  B. Peters.  Reaction Coordinates and Mechanistic Hypothesis Tests. , 2016, Annual review of physical chemistry.

[38]  T. S. V. Erp.  Dynamical Rare Event Simulation Techniques for Equilibrium and Nonequilibrium Systems , 2011, 1101.0927.

[39]  Ryan S. DeFever,et al.  Suppression of sub-surface freezing in free-standing thin films of a coarse-grained model of water. , 2014, Physical chemistry chemical physics : PCCP.

[40]  P. R. ten Wolde,et al.  Sampling rare switching events in biochemical networks. , 2004, Physical review letters.

[41]  Wolfgang Lechner,et al.  Nonlinear reaction coordinate analysis in the reweighted path ensemble. , 2010, The Journal of chemical physics.

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

[43]  Exact rate calculations by trajectory parallelization and tilting. , 2009, The Journal of chemical physics.

[44]  Dependence of the Rate of LiF Ion-Pairing on the Description of Molecular Interaction. , 2016, The journal of physical chemistry. B.

[45]  P. R. ten Wolde,et al.  Rate constants for proteins binding to substrates with multiple binding sites using a generalized forward flux sampling expression. , 2018, The Journal of chemical physics.

[46]  Christoph Dellago,et al.  Kinetic pathways of ion pair dissociation in water , 1999 .

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

[48]  Michele Parrinello,et al.  Collective Variables from Local Fluctuations. , 2018, The journal of physical chemistry letters.

[49]  B. Trout,et al.  Obtaining reaction coordinates by likelihood maximization. , 2006, The Journal of chemical physics.

[50]  P ? ? ? ? ? ? ? % ? ? ? ? , 1991 .

[51]  Ernesto E. Borrero,et al.  Kinetics and reaction coordinate for the isomerization of alanine dipeptide by a forward flux sampling protocol. , 2009, The Journal of chemical physics.

[52]  Tianshu Li,et al.  Probing methane hydrate nucleation through the forward flux sampling method. , 2014, The journal of physical chemistry. B.

[53]  Pratyush Tiwary,et al.  Reweighted autoencoded variational Bayes for enhanced sampling (RAVE). , 2018, The Journal of chemical physics.

[54]  Aaron R Dinner,et al.  Automatic method for identifying reaction coordinates in complex systems. , 2005, The journal of physical chemistry. B.

[55]  Ernesto E. Borrero,et al.  Reaction coordinates and transition pathways of rare events via forward flux sampling. , 2007, The Journal of chemical physics.

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

[57]  Ann B. Lee,et al.  Geometric diffusions as a tool for harmonic analysis and structure definition of data: diffusion maps. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[58]  B. Alsberg,et al.  Local initiation conditions for water autoionization , 2018, Proceedings of the National Academy of Sciences.

[59]  Ioannis G Kevrekidis,et al.  Intrinsic map dynamics exploration for uncharted effective free-energy landscapes , 2016, Proceedings of the National Academy of Sciences.

[60]  Christoph Dellago,et al.  Trajectory‐Based Rare Event Simulations , 2010 .