Integrating Path Sampling with Enhanced Sampling for Rare-event Kinetics

Dhiman Ray

DOI: https://doi.org/10.26434/chemrxiv-2024-qr6kv

2024-09-18

Abstract:Studying the kinetics of long-timescale rare events is a fundamental challenge in molecular simulation. To address this problem, we propose an integration of two different rare-event sampling philosophies: biased enhanced sampling and unbiased path sampling. Enhanced sampling methods e.g. metadynamics can facilitate enthalpic barrier crossing by applying an external bias potential. On the contrary, path sampling methods like weighted ensemble (WE) lack explicit mechanisms to overcome energetic barriers. However, they can accelerate the exploration of rugged free energy surfaces through trajectory resampling. We show that a judicious combination of the weighted ensemble with a metadynamics-like algorithm, can synergize the strengths and mitigate the deficiencies of path sampling and enhanced sampling approaches. The resulting integrated sampling (IS) algorithm improves the computational efficiency of calculating the kinetics of peptide conformational transitions, protein unfolding, and the dissociation of a ligand-receptor complex. Furthermore, the IS approach can direct sampling along the minimum free energy pathway even when the collective variable used for biasing is suboptimal. These advantages make the integrated sampling algorithm suitable for studying the kinetics of complex molecular systems of biological and pharmaceutical relevance.

Chemistry

What problem does this paper attempt to address?

The paper attempts to address the problem of efficiently studying the dynamics of rare events over long time scales in molecular dynamics simulations. Specifically, the authors propose a new algorithm that combines two different sampling methods (enhanced sampling and path sampling) to overcome the limitations of traditional methods in dealing with high energy barriers and complex free energy landscapes.

The main issues discussed in the paper include:
1. **High energy barrier crossing**: In molecular dynamics simulations, high energy barriers prevent the system from transitioning from one metastable state to another within computationally accessible time scales.
2. **Acceleration of diffusion processes**: For complex biomolecular systems (such as protein folding and ligand-receptor binding), slow diffusion processes also increase transition times.

To address these issues, the authors combine two different methods:
- **Enhanced sampling methods** (such as metadynamics): Accelerate barrier crossing by applying external bias potentials.
- **Path sampling methods** (such as weighted ensemble): Accelerate exploration of the free energy landscape through trajectory resampling but lack a clear mechanism to overcome energy barriers.

By integrating the advantages of these two methods, the newly proposed Integrated Sampling (IS) algorithm can significantly improve computational efficiency without sacrificing dynamical accuracy. This method is suitable for studying the dynamic properties of complex molecular systems with biological and pharmacological relevance.