Path sampling challenges in large biomolecular systems: RETIS and REPPTIS for ABL-imatinib kinetics
arXiv:2504.14722 · doi:10.1016/j.bpj.2025.04.020
Abstract
Predicting the kinetics of drug-protein interactions is crucial for understanding drug efficacy, particularly in personalized medicine, where protein mutations can significantly alter drug residence times. This study applies Replica Exchange Transition Interface Sampling (RETIS) and its Partial Path variant (REPPTIS) to investigate the dissociation kinetics of imatinib from Abelson nonreceptor tyrosine kinase (ABL) and mutants relevant to chronic myeloid leukemia therapy. These path-sampling methods offer a bias-free alternative to conventional approaches requiring qualitative predefined reaction coordinates. Nevertheless, the complex free-energy landscape of ABL-imatinib dissociation presents significant challenges. Multiple metastable states and orthogonal barriers lead to parallel unbinding pathways, complicating convergence in TIS-based methods. Despite employing computational efficiency strategies such as asynchronous replica exchange, full convergence remained elusive. This work provides a critical assessment of path sampling in high-dimensional biological systems, discussing the need for enhanced initialization strategies, advanced Monte Carlo path generation moves, and machine learning-derived reaction coordinates to improve kinetic predictions of drug dissociation with minimal prior knowledge.
Preprint
References in corpus (13)
- Canonical sampling through velocity-rescaling
- PLUMED 2: New feathers for an old bird
- String Method for the Study of Rare Events
- A Novel Path Sampling Method for the Calculation of Rate Constants
- Forward Flux Sampling for rare event simulations
- Rate constants for diffusive processes by partial path sampling
- Reaction rate calculation by parallel path swapping
- A multiple replica approach to simulate reactive trajectories
- Precision shooting: Sampling long transition pathways
- Efficiency analysis of reaction rate calculation methods using analytical models I: The 2D sharp barrier
- Preserving Correlations Between Trajectories for Efficient Path Sampling
- Exchanging replicas with unequal cost, infinitely and permanently
- Enhanced path sampling using subtrajectory Monte Carlo moves