Kinetics-Optimized Enhanced Sampling Using Mean First Passage Times

Tiejun Wei,Balint Dudas,Edina Rosta

2024-06-13

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.

Biological Physics

What problem does this paper attempt to address?

The paper proposes a novel enhanced sampling algorithm called KOMBI (Kinetics-Optimized Markovian Bias Integration) to address the sampling problem of important long time-scale processes in molecular dynamics simulations. In traditional molecular dynamics simulations, capturing important dynamic events can be challenging since they occur very slowly. To overcome this problem, various enhanced sampling methods have been developed, but these methods often focus on thermodynamic optimization rather than molecular dynamics.

The core of the KOMBI algorithm is to optimize biases in the system by minimizing the average mean first passage time (MFPT) from the starting state to the target state. MFPT measures the average time required for the system to transition from one state to another. In the paper, the authors establish the relationship between biases and unbiased dynamics by using discretized Markov state models (MSMs) and dynamic histogram analysis method (DHAM) to estimate unbiased Markov matrices.

By iteratively optimizing the parameters of Gaussian functions, the algorithm can adaptively adjust biases to make transitions from the starting state to the target state more efficient. In 1D and 2D models, as well as molecular dynamics simulations of actual NaCl ion pairs and Ras GTPase systems, the KOMBI algorithm demonstrates faster time scales to reach the target state compared to existing enhanced sampling methods, proving its potential breakthrough application in rare event studies in fields such as drug discovery and materials science.

In conclusion, this paper addresses how to more effectively simulate key slow transitions in biological molecular systems by optimizing dynamics-based biases, providing new strategies for understanding and simulating complex biological processes.