Flow Matching for Optimal Reaction Coordinates of Biomolecular System
arXiv:2408.17139 · doi:10.1021/acs.jctc.4c01139
Abstract
We present flow matching for reaction coordinates (FMRC), a novel deep learning algorithm designed to identify optimal reaction coordinates (RC) in biomolecular reversible dynamics. FMRC is based on the mathematical principles of lumpability and decomposability, which we reformulate into a conditional probability framework for efficient data-driven optimization using deep generative models. While FMRC does not explicitly learn the well-established transfer operator or its eigenfunctions, it can effectively encode the dynamics of leading eigenfunctions of the system transfer operator into its low-dimensional RC space. We further quantitatively compare its performance with several state-of-the-art algorithms by evaluating the quality of Markov state models (MSM) constructed in their respective RC spaces, demonstrating the superiority of FMRC in three increasingly complex biomolecular systems. In addition, we successfully demonstrated the efficacy of FMRC for bias deposition in the enhanced sampling of a simple model system. Finally, we discuss its potential applications in downstream applications such as enhanced sampling methods and MSM construction.
For Supporting Information, please see https://pubs.acs.org/doi/full/10.1021/acs.jctc.4c01139
References in corpus (13)
- Escaping free-energy minima
- The information bottleneck method
- Identification of slow molecular order parameters for Markov model construction
- Rethinking Metadynamics: from bias potentials to probability distributions
- Nonlinear Discovery of Slow Molecular Modes using State-Free Reversible VAMPnets
- Variational Koopman models: slow collective variables and molecular kinetics from short off-equilibrium simulations
- Autonomous artificial intelligence discovers mechanisms of molecular self-organization in virtual experiments
- State Predictive Information Bottleneck
- Identification of simple reaction coordinates from complex dynamics
- Computing the Committor with the Committor: an Anatomy of the Transition State Ensemble
- Transition manifolds of complex metastable systems: Theory and data-driven computation of effective dynamics
- Data-driven Computation of Molecular Reaction Coordinates
- Optimal Low-dimensional Approximation of Transfer Operators via Flow Matching: Computation and Error Analysis