Reaction Coordinates for Conformational Transitions using Linear Discriminant Analysis on Positions
arXiv:2301.03479 · doi:10.1021/acs.jctc.3c00051
Abstract
In this work, we demonstrate that Linear Discriminant Analysis (LDA) applied to atomic positions in two different states of a biomolecule produces a good reaction coordinate between those two states. Atomic coordinates of a macromolecule are a direct representation of a macromolecular configuration, and yet they have not been used in enhanced sampling studies due to a lack of rotational and translational invariance. We resolve this issue using the technique of our prior work, where a molecular configuration is considered a member of an equivalence class in size-and-shape space, which is the set of all configurations that can be translated and rotated to a single point within a reference multivariate Gaussian distribution characterizing a single molecular state. The reaction coordinates produced by LDA applied to positions are shown to be good reaction coordinates both in terms of characterizing the transition between two states of a system within a long molecular dynamics (MD) simulation, and also ones that allow us to readily produce free energy estimates along that reaction coordinate using enhanced sampling MD techniques.
arXiv comment11 pages, 5 figures (plus 2 SI pages and 3 SI figures)
References in corpus (23)
- PLUMED 2: New feathers for an old bird
- Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method
- Machine learning for molecular simulation
- Enhanced Sampling in Molecular Dynamics Using Metadynamics, Replica-Exchange, and Temperature-Acceleration
- Time-lagged autoencoders: Deep learning of slow collective variables for molecular kinetics
- Rethinking Metadynamics: from bias potentials to probability distributions
- Enhanced sampling methods for molecular dynamics simulations
- Caliber based spectral gap optimization of order parameters (SGOOP) for sampling complex molecular systems
- Molecular enhanced sampling with autoencoders: On-the-fly collective variable discovery and accelerated free energy landscape exploration
- Data-Driven Collective Variables for Enhanced Sampling
- Nonlinear Discovery of Slow Molecular Modes using State-Free Reversible VAMPnets
- A self-learning algorithm for biased molecular dynamics
- A Unified Approach to Enhanced Sampling
- Enhanced Conformational Sampling using Replica Exchange with Collective-Variable Tempering
- Exploring High Dimensional Free Energy Landscapes: Temperature Accelerated Sliced Sampling
- Folding a Small Protein Using Harmonic Linear Discriminant Analysis
- Size-and-Shape Space Gaussian Mixture Models for Structural Clustering of Molecular Dynamics Trajectories
- Finding Efficient Collective Variables: The Case of Crystallization
- Collective Variables for the Study of Crystallization
- Accurate multiple time step in biased molecular simulations
- Assessing models of force-dependent unbinding rates via infrequent metadynamics
- Reweighted Manifold Learning of Collective Variables from Enhanced Sampling Simulations
- Infinite Switch Simulated Tempering in Force (FISST)