Abstract:Molecular dynamics simulations have become essential in many areas of atomistic modelling from drug discovery to materials science. They provide critical atomic-level insights into key dynamical events experiments cannot easily capture. However, their impact often falls short as the timescales of the important processes are inaccessible using standard molecular dynamics. Enhanced sampling methods provided avenues to access such crucial rare events, for example key slow conformational changes of biomolecules. However, the bias in enhanced sampling simulations is rarely optimized, and even if they are, the optimization criteria is based on the thermodynamics or Hamiltonian of the system, but do not directly consider molecular kinetics. Here, we introduce a novel enhanced sampling algorithm that adaptively optimizes the bias based on the kinetics of the system for the first time. We identify the optimal bias that minimizes a key physical observable, the mean first passage time (MFPT) from a starting state to a target state. Our algorithm makes use of the relation between biased and unbiased kinetics obtained from discretized Markov state models (MSMs), as established in the dynamic histogram analysis method (DHAM). We demonstrate the applicability of the method for different 1D and 2D analytical potential-based model examples, NaCl dissociation in explicit water, and phosphate unbinding in Ras GTPase. Our algorithm has excellent performance compared with state-of-art enhanced sampling methods in terms of the timescales required to reach the final state in the benchmarking systems. Our findings provide a novel, kinetics-driven enhanced sampling strategy, signatured by a targeted approach to facilitate mapping rare events, with the potential for breakthrough applications in drug discovery and materials science.

Short-Time Infrequent Metadynamics for Improved Kinetics Inference

Short-Time Infrequent Metadynamics for Improved Kinetics Inference

Frequency Adaptive Metadynamics for the Calculation of Rare-Event Kinetics.

Good rates from bad coordinates: the exponential average time-dependent rate approach

Transition rates, survival probabilities, and quality of bias from time-dependent biased simulations

Inference of non-exponential kinetics through stochastic resetting

Combining stochastic resetting with Metadynamics to speed-up molecular dynamics simulations

Metadynamics to Enhance Sampling in Biomolecular Simulations

Kinetics-Optimized Enhanced Sampling Using Mean First Passage Times

Resetting Metadynamics

Accelerating inference for stochastic kinetic models

Inference for Stochastic Chemical Kinetics Using Moment Equations and System Size Expansion

Metadynamics: a method to simulate rare events and reconstruct the free energy in biophysics, chemistry and material science

From Thermodynamics to Kinetics: Enhanced Sampling of Rare Events

Acceleration of Molecular Simulations by Parametric Time-Lagged tSNE Metadynamics

Enhanced sampling of chemical and biochemical reactions with metadynamics

Enhanced Sampling in Molecular Dynamics Using Metadynamics, Replica-Exchange, and Temperature-Acceleration

Bayesian Inference for Kinetic Models of Biotransformation Using A Generalized Rate Equation

A Vine-Copula Based Adaptive Mcmc Sampler for Efficient Inference of Dynamical Systems

Acceleration and sensitivity analysis of lattice kinetic Monte Carlo simulations using parallel processing and rate constant rescaling.

Accelerating Molecular Dynamics through Informed Resetting