Abstract:Biomolecules often exhibit complex free energy landscapes in which long-lived metastable states are separated by large energy barriers. Overcoming these barriers to robustly sample transitions between the metastable states with classical molecular dynamics (MD) simulations presents a challenge. To circumvent this issue, collective variable (CV)-based enhanced sampling MD approaches are often employed. Traditional CV selection relies on intuition and prior knowledge of the system. This approach introduces bias, which can lead to incomplete mechanistic insights. Thus, automated CV detection is desired to gain a deeper understanding of the system/process. Analysis of MD data with various machine learning algorithms, such as Principal Component Analysis (PCA), Support Vector Machine (SVM), and Linear Discriminant Analysis (LDA)-based approaches have been implemented for automated CV detection. However, their performance has not been systematically evaluated on structurally and mechanistically complex biological systems. Here, we applied these methods to MD simulations of the MFSD2A (Major Facilitator Superfamily Domain 2A) lysolipid transporter in multiple functionally relevant metastable states with the goal of identifying optimal CVs that would structurally discriminate these states. Specific emphasis was on the automated detection and interpretive power of LDA-based CVs. We found that LDA methods, which included a novel gradient descent-based multiclass harmonic variant, termed GDHLDA, we developed here, outperform PCA in class separation, exhibiting remarkable consistency in extracting CVs critical for distinguishing metastable states. Furthermore, the identified CVs included features previously associated with conformational transitions in MFSD2A. Specifically, conformational shifts in transmembrane helix 7 and in residue Y294 on this helix emerged as critical features discriminating the metastable states in MFSD2A. This highlights the effectiveness of LDA-based approaches in automatically extracting from MD trajectories CVs of functional relevance that can be used to drive biased MD simulations to efficiently sample conformational transitions in the molecular system.

A Mode Evolution Metric to Extract Reaction Coordinates for Biomolecular Conformational Transitions

A Mode Evolution Metric to Extract Converged Reaction Coordinates for Biomolecular Conformational Transitions

Learning Effective Molecular Models from Experimental Observables.

MLCV: Bridging Machine-Learning-Based Dimensionality Reduction and Free-Energy Calculation

Reaction Coordinates for Conformational Transitions using Linear Discriminant Analysis on Positions

Effective Data-Driven Collective Variables for Free Energy Calculations from Metadynamics of Paths

Data-driven Computation of Molecular Reaction Coordinates

Molecular dynamics-based approaches for enhanced sampling of long-time, large-scale conformational changes in biomolecules

Mapping transiently formed and sparsely populated conformations on a complex energy landscape

Enhanced Sampling of Protein Conformational Changes via True Reaction Coordinates from Energy Relaxation

Identifying a Feasible Transition Pathway between Two Conformational States for a Protein

Machine Learning Deciphered Molecular Mechanistics with Accurate Kinetic and Thermodynamic Prediction

Learning reaction coordinates via cross-entropy minimization: Application to alanine dipeptide

Automated Collective Variable Discovery for MFSD2A transporter from molecular dynamics simulations

One Descriptor to Fold Them All: Harnessing Intuition and Machine Learning to Identify Transferable Lasso Peptide Reaction Coordinates

Biomolecular transitions: efficient computation of pathways, free energies, and rates

Discriminant Analysis Optimizes Progress Coordinate in Weighted Ensemble Simulations of Rare Event Kinetics

Acceleration with Interpretability: Surrogate Model Based Collective Variable for Enhanced Sampling

Reaction Coordinates by Nonlinear Dimensionality Reduction

Screening Fast-Mode Motion in Collective Variable Discovery for Biochemical Processes

Spectral Map for Slow Collective Variables, Markovian Dynamics, and Transition State Ensembles