Rare Event Kinetics from Adaptive Bias Enhanced Sampling

Abstract

We introduce a novel enhanced sampling approach named OPES flooding for calculating the kinetics of rare events from atomistic molecular dynamics simulation. This method is derived from the On-the-fly-Probability-Enhanced-Sampling (OPES) approach [Invernizzi and Parrinello, JPC Lett. 2020], which has been recently developed for calculating converged free energy surfaces for complex systems. In this paper, we describe the theoretical details of the OPES flooding technique and demonstrate the application on three systems of increasing complexity: barrier crossing in a two-dimensional double well potential, conformational transition in the alanine dipeptide in gas phase, and the folding and unfolding of the chignolin polypeptide in aqueous environment. From extensive tests, we show that the calculation of accurate kinetics not only requires the transition state to be bias-free, but the amount of bias deposited should also not exceed the effective barrier height measured along the chosen collective variables. In this vein, the possibility of computing rates from biasing suboptimal order parameters has also been explored. Furthermore, we describe the choice of optimum parameter combinations for obtaining accurate results from limited computational effort.

References in corpus (14)

• Escaping free-energy minima
• PLUMED 2: New feathers for an old bird
• Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method
• From Metadynamics to Dynamics
• Simulating rare events in equilibrium or non-equilibrium stochastic systems
• Rethinking Metadynamics: from bias potentials to probability distributions
• Data-Driven Collective Variables for Enhanced Sampling
• A Unified Approach to Enhanced Sampling
• Folding a Small Protein Using Harmonic Linear Discriminant Analysis
• Gaussian Mixture Based Enhanced Sampling
• Variationally Optimized Free Energy Flooding for Rate Calculation
• Making the best of a bad situation: a multiscale approach to free energy calculation
• Weighted Ensemble Milestoning (WEM): A Combined Approach for Rare Event Simulations
• Conformational Entropy as Collective Variable for Proteins

Cited by in corpus (9)

• Water regulates the residence time of Benzamidine in Trypsin
• Deep Learning Collective Variables from Transition Path Ensemble
• Enhanced Sampling in the Age of Machine Learning: Algorithms and Applications
• Stochastic Resetting for Enhanced Sampling
• Intramolecular and water mediated tautomerism of solvated glycine
• Learning Collective Variables with Synthetic Data Augmentation through Physics-Inspired Geodesic Interpolation
• Weighted Active Space Protocol for Multireference Machine-Learned Potentials
• WeTICA: A directed search weighted ensemble based enhanced sampling method to estimate rare event kinetics in a reduced dimensional space
• Guiding Peptide Kinetics via Collective-Variable Tuning of Free-Energy Barriers