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.

Advantages of multi-dimensional biasing in accelerated dynamics: application to the calculation of the acid $pK_a$ for acetic acid

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

Enhanced Sampling of Multidimensional Free-Energy Landscapes Using Adaptive Biasing Forces

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

Ridge-based Bias Potentials to Accelerate Molecular Dynamics.

Improving the Exploration of High Dimensional Free Energy Landscape by a Combination of Temperature Accelerated Sliced Sampling and Parallel Biasing

Free energy calculations: An efficient adaptive biasing potential method

Improved Sampling of Adaptive Path Collective Variables by Stabilized Extended-System Dynamics

Reinforced Dynamics for Enhanced Sampling in Large Atomic and Molecular Systems

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

Data‐Efficient Active Learning for Thermodynamic Integration: Acidity Constants of BiVO4 in Water

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

Biasing the center of charge in molecular dynamics simulations with empirical valence bond models: free energetics of an excess proton in a water droplet

Enhanced sampling of chemical and biochemical reactions with metadynamics

A microscopic description of acid-base equilibrium

In silico calculation of acidity constants of carbonic acid conformers.

Reweighted Jarzynski sampling: Acceleration of rare events and free energy calculation with a bias potential learned from nonequilibrium work

Towards the "Eldorado" of pKa Determination: A Reliable and Rapid DFT Model

Achieving an Order of Magnitude Speed-up in Hybrid Functional and Plane Wave based Ab Initio Molecular Dynamics: Applications to Proton Transfer Reactions in Enzymes and in Solution

Kinetics-Optimized Enhanced Sampling Using Mean First Passage Times

Short-Time Infrequent Metadynamics for Improved Kinetics Inference