Analysis of the adaptive multilevel splitting method on the isomerization of alanine dipeptide

We apply the adaptive multilevel splitting (AMS) method to the C eq → C ax transition of alanine dipeptide in vacuum. Some properties of the algorithm are numerically illustrated, such as the unbiasedness of the probability estimator and the robustness of the method with respect to the reaction coordinate. We also calculate the transition time obtained via the probability estimator, using an appropriate ensemble of initial conditions. Finally, we show how the AMS method can be used to compute an approximation of the committor function. © 2019 Wiley Periodicals, Inc.

[1]  Jianfeng Lu,et al.  PR ] 7 M ar 2 01 3 REACTIVE TRAJECTORIES AND THE TRANSITION PATH , 2014 .

[2]  T. Lelièvre,et al.  ANALYSIS OF ADAPTIVE MULTILEVEL SPLITTING ALGORITHMS IN AN IDEALIZED CASE , 2014, 1405.1352.

[3]  R. Copeland,et al.  Drug–target residence time and its implications for lead optimization , 2007, Nature Reviews Drug Discovery.

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

[5]  Klaus Schulten,et al.  Adaptive Multilevel Splitting Method for Molecular Dynamics Calculation of Benzamidine-Trypsin Dissociation Time. , 2016, Journal of chemical theory and computation.

[6]  Rémi Bardenet,et al.  Monte Carlo Methods , 2013, Encyclopedia of Social Network Analysis and Mining. 2nd Ed..

[7]  Christophe Chipot,et al.  Exploring Multidimensional Free Energy Landscapes Using Time-Dependent Biases on Collective Variables. , 2010, Journal of chemical theory and computation.

[8]  A Rojnuckarin,et al.  Brownian dynamics simulations of protein folding: access to milliseconds time scale and beyond. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

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

[10]  Arnaud Guyader,et al.  A multiple replica approach to simulate reactive trajectories. , 2010, The Journal of chemical physics.

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

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

[13]  Ludovic Goudenège,et al.  Unbiasedness of some generalized Adaptive Multilevel Splitting algorithms , 2016 .

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

[15]  F. Cérou,et al.  Adaptive Multilevel Splitting for Rare Event Analysis , 2007 .