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.

Approximating First Hitting Point Distribution in Milestoning for Rare Event Kinetics

Comparison of Accuracy and Efficiency of Milestoning Variants: Introducing Buffer Milestoning

A Mathematical Theory of Optimal Milestoning (with a Detour Via Exact Milestoning)

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

Advancements in Milestoning I: Accelerated Milestoning via "Wind" Assisted Re-weighted Milestoning (WARM)

Weighted ensemble milestoning (WEM): A combined approach for rare event simulations

Milestoning network refinement by incorporating experimental thermodynamic and kinetic data

Kinetics-Optimized Enhanced Sampling Using Mean First Passage Times

First-Passage Kinetic Monte Carlo method

Rapid state-recrossing kinetics in non-Markovian systems

A Random-Batch Monte Carlo Method for Many-Body Systems with Singular Kernels.

Mean First-Passage Time in a Bistable Kinetic Model with Stochastic Potentials

A mathematical framework for exact milestoning

A Particle Fokker-Planck Algorithm with Multiscale Temporal Discretization for Rarefied and Continuum Gas Flows

On Committor Functions in Milestoning

Parallelization, processor communication and error analysis in lattice kinetic Monte Carlo

A Multilevel Monte Carlo Asymptotic-Preserving Particle Method for Kinetic Equations in the Diffusion Limit

Unified perspective on exponential tilt and bridge algorithms for rare trajectories of discrete Markov processes

Adaptive tau-leaping methods for microscopic-lattice kinetic Monte Carlo simulations

An Efficient Quasi-3D Particle Tracking-Based Approach for Transport Through Fractures with Application to Dynamic Dispersion Calculation