Identification of slow molecular order parameters for Markov model construction

Abstract

A goal in the kinetic characterization of a macromolecular system is the description of its slow relaxation processes, involving (i) identification of the structural changes involved in these processes, and (ii) estimation of the rates or timescales at which these slow processes occur. Most of the approaches to this task, including Markov models, Master-equation models, and kinetic network models, start by discretizing the high-dimensional state space and then characterize relaxation processes in terms of the eigenvectors and eigenvalues of a discrete transition matrix. The practical success of such an approach depends very much on the ability to finely discretize the slow order parameters. How can this task be achieved in a high-dimensional configuration space without relying on subjective guesses of the slow order parameters? In this paper, we use the variational principle of conformation dynamics to derive an optimal way of identifying the "slow subspace" of a large set of prior order parameters - either generic internal coordinates (distances and dihedral angles), or a user-defined set of parameters. It is shown that a method to identify this slow subspace exists in statistics: the time-lagged independent component analysis (TICA). Furthermore, optimal indicators-order parameters indicating the progress of the slow transitions and thus may serve as reaction coordinates-are readily identified. We demonstrate that the slow subspace is well suited to construct accurate kinetic models of two sets of molecular dynamics simulations, the 6-residue fluorescent peptide MR121-GSGSW and the 30-residue natively disordered peptide KID. The identified optimal indicators reveal the structural changes associated with the slow processes of the molecular system under analysis.

Cited by in corpus (117)

• Machine learning for molecular simulation
• VAMPnets: Deep learning of molecular kinetics
• Time-lagged autoencoders: Deep learning of slow collective variables for molecular kinetics
• Caliber based spectral gap optimization of order parameters (SGOOP) for sampling complex molecular systems
• Data-driven model reduction and transfer operator approximation
• TorchMD: A deep learning framework for molecular simulations
• Variational cross-validation of slow dynamical modes in molecular kinetics
• Machine learning force fields and coarse-grained variables in molecular dynamics: application to materials and biological systems
• Coarse Graining Molecular Dynamics with Graph Neural Networks
• Variational Encoding of Complex Dynamics
• Unsupervised machine learning in atomistic simulations, between predictions and understanding
• Projected and Hidden Markov Models for calculating kinetics and metastable states of complex molecules
• 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
• A variational conformational dynamics approach to the selection of collective variables in metadynamics
• Machine Learning Coarse-Grained Potentials of Protein Thermodynamics
• Eigendecompositions of Transfer Operators in Reproducing Kernel Hilbert Spaces
• Identification of simple reaction coordinates from complex dynamics
• Learning Molecular Dynamics with Simple Language Model built upon Long Short-Term Memory Neural Network
• Galerkin Approximation of Dynamical Quantities using Trajectory Data
• Perspective: Markov Models for Long-Timescale Biomolecular Dynamics
• Coarse-Grained Modelling Out of Equilibrium
• A unified framework for machine learning collective variables for enhanced sampling simulations: $\text{mlcolvar}$
• Water regulates the residence time of Benzamidine in Trypsin
• Special Topic: Markov Models of Molecular Kinetics
• Machine Learning Implicit Solvation for Molecular Dynamics
• Flow-matching -- efficient coarse-graining of molecular dynamics without forces
• Variational Selection of Features for Molecular Kinetics
• Ensemble Learning of Coarse-Grained Molecular Dynamics Force Fields with a Kernel Approach
• Transition manifolds of complex metastable systems: Theory and data-driven computation of effective dynamics
• Deep Learning Collective Variables from Transition Path Ensemble
• GraphVAMPNet, using graph neural networks and variational approach to markov processes for dynamical modeling of biomolecules
• Time-Lagged t-Distributed Stochastic Neighbor Embedding (t-SNE) of Molecular Simulation Trajectories
• Using Markov State Models to Study Self-Assembly
• Adversarial-Residual-Coarse-Graining: Applying machine learning theory to systematic molecular coarse-graining
• Adaptive Markov State Model estimation using short reseeding trajectories
• Simulations meet Machine Learning in Structural Biology
• Unfolding Hidden Barriers by Active Enhanced Sampling
• Correlation-based feature selection to identify functional dynamics in proteins
• Efficient implementation of atom-density representations
• Finding Efficient Collective Variables: The Case of Crystallization
• Enhanced Sampling in the Age of Machine Learning: Algorithms and Applications
• Computationally driven discovery of SARS-CoV-2 Mpro inhibitors: from design to experimental validation
• Efficient Bayesian estimation of Markov model transition matrices with given stationary distribution
• Human learning for molecular simulations: the Collective Variables Dashboard in VMD
• Weighted Ensemble Milestoning (WEM): A Combined Approach for Rare Event Simulations
• Relaxation Mode Analysis and Markov State Relaxation Mode Analysis for Chignolin in Aqueous Solution near a Transition Temperature
• Ensuring thermodynamic consistency with invertible coarse-graining
• Identification and Analysis of Transition and Metastable Markov States
• Integrated VAC: A robust strategy for identifying eigenfunctions of dynamical operators
• Identification of kinetic order parameters for non-equilibrium dynamics
• Dynamical reweighting methods for Markov models
• How wet should be the reaction coordinate for ligand unbinding?
• Capabilities and Limitations of Time-lagged Autoencoders for Slow Mode Discovery in Dynamical Systems
• Predicting the kinetics of RNA oligonucleotides using Markov state models
• The mechanism of RNA base fraying: molecular dynamics simulations analyzed with core-set Markov state models
• MSM/RD: Coupling Markov state models of molecular kinetics with reaction-diffusion simulations
• Optimal data-driven estimation of generalized Markov state models for non-equilibrium dynamics
• A Deep Autoencoder Framework for Discovery of Metastable Ensembles in Biomacromolecules
• Temporally coherent backmapping of molecular trajectories from coarse-grained to atomistic resolution
• xTRAM: Estimating equilibrium expectations from time-correlated simulation data at multiple thermodynamic states
• Neural Canonical Transformation with Symplectic Flows
• Recent advances in describing and driving crystal nucleation using machine learning and artificial intelligence
• Kernel methods for detecting coherent structures in dynamical data
• TimeSOAP: Tracking high-dimensional fluctuations in complex molecular systems via time-variations of SOAP spectra
• Characterizing metastable states with the help of machine learning
• Dimensional Reduction of Markov State Models from Renormalization Group Theory
• Selecting Features for Markov Modeling: A Case Study on HP35
• Maximally predictive states: from partial observations to long timescales
• Dimensionality reduction to maximize prediction generalization capability
• Data-driven Computation of Molecular Reaction Coordinates
• Learning Clustered Representation for Complex Free Energy Landscapes
• Multidimensional minimum-work control of a 2D Ising model
• Path probability ratios for Langevin dynamics -- exact and approximate
• Towards a Benchmark for Markov State Models: The Folding of HP35
• An implementation of the maximum-caliber principle by replica-averaged time-resolved restrained simulations
• Data-driven construction of stochastic reduced dynamics encoded with non-Markovian features
• Manifold Learning in Atomistic Simulations: A Conceptual Review
• Optimal reaction coordinates and kinetic rates from the projected dynamics of transition paths
• Solving eigenvalue PDEs of metastable diffusion processes using artificial neural networks
• Deeptime: a Python library for machine learning dynamical models from time series data
• Efficient Approximation of Molecular Kinetics using Random Fourier Features
• MSM lag time cannot be used for variational model selection
• Deep learning the slow modes for rare events sampling
• Inexact iterative numerical linear algebra for neural network-based spectral estimation and rare-event prediction
• Deflation reveals dynamical structure in nondominant reaction coordinates
• Non-equilibrium Markov state modeling of periodically driven biomolecules
• Learning Collective Variables with Synthetic Data Augmentation through Physics-Inspired Geodesic Interpolation
• Non-Equilibrium Markov State Modeling of the Globule-Stretch Transition
• On the Role of Solvent in Hydrophobic Cavity-ligand Recognition Kinetics
• Flow Matching for Optimal Reaction Coordinates of Biomolecular System
• Error bounds for dynamical spectral estimation
• Thermodynamically Optimized Machine-learned Reaction Coordinates for Hydrophobic Ligand Dissociation
• Learning Markovian Dynamics with Spectral Maps
• AMARO: All Heavy-Atom Transferable Neural Network Potentials of Protein Thermodynamics
• Coarse-graining conformational dynamics with multi-dimensional generalized Langevin equation: how, when, and why
• Markov modeling of peptide folding in the presence of protein crowders
• WeTICA: A directed search weighted ensemble based enhanced sampling method to estimate rare event kinetics in a reduced dimensional space
• Identifying the leading dynamics of ubiquitin: a comparison between the tICA and the LE4PD slow fluctuations in amino acids' position
• Computing Long Timescale Biomolecular Dynamics using Quasi-Stationary Distribution Kinetic Monte Carlo (QSD-KMC)
• Markov State Modeling of Sliding Friction
• Relevant, hidden, and frustrated information in high-dimensional analyses of complex dynamical systems with internal noise
• Collective variables between large-scale states in turbulent convection
• Large Scale Training of Graph Neural Networks for Optimal Markov-Chain Partitioning Using the Kemeny Constant
• Progress in deep Markov State Modeling: Coarse graining and experimental data restraints
• A Sinking Approach to Explore Arbitrary Areas in Free Energy Landscapes
• The Effect of Hydration and Dynamics on the Mass Density of Single Proteins
• Recovering Hidden Degrees of Freedom Using Gaussian Processes
• Exploring the free energy gain of phase separation via Markov State Modeling
• Volume-Scaled Common Nearest Neighbor Clustering Algorithm with Free-Energy Hierarchy
• A Markov State Modeling analysis of sliding dynamics of a 2D model
• Dimensional Reduction of Dynamical Systems by Machine Learning: Automatic Generation of the Optimum Extensive Variables and Their Time-Evolution Map
• Improving Estimation of the Koopman Operator with Kolmogorov-Smirnov Indicator Functions
• MDIntrinsicDimension: Dimensionality-Based Analysis of Collective Motions in Macromolecules from Molecular Dynamics Trajectories
• Markov State Models for Tracking Reaction Dynamics on Catalytic Nanoparticles
• Adaptive tensor train metadynamics for high-dimensional free energy exploration
• Variational embedding of protein folding simulations using gaussian mixture variational autoencoders