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 Data-Driven Accelerated Sampling Method for Searching Functional States of Proteins

A Two-ended DAta-Driven Accelerated Sampling Method for Exploring the Transition Pathways Between Two Known States of Protein.

Efficient Conformational Search Based on Structural Dissimilarity Sampling: Applications for Reproducing Structural Transitions of Proteins.

A Discard-and-Restart MD algorithm for the sampling of realistic protein transition states and enhance structure-based drug discovery

Digital Net Conformational Sampler (DNCS) – An Enhanced Sampling Technique to Explore Conformational Space of Intrinsically Disordered Peptides

Enhanced molecular dynamics sampling of drug target conformations

Exploration of Multi-State Conformational Dynamics and Underlying Global Functional Landscape of Maltose Binding Protein

Scaled Alternating Steepest Descent Algorithm Applied for Protein Structure Determination from Nuclear Magnetic Resonance Data.

Monte Carlo Sampling of Protein Folding by Combining an All-Atom Physics-Based Model with a Native State Bias

LAST: Latent Space Assisted Adaptive Sampling for Protein Trajectories

Automated Collective Variable Discovery for MFSD2A transporter from molecular dynamics simulations

Advances in Molecular Dynamics Simulations and Enhanced Sampling Methods for the Study of Protein Systems

Structural Refinement of Proteins by Restrained Molecular Dynamics Simulations with Non-Interacting Molecular Fragments

Fast Sampling of Protein Conformational Dynamics

Accurate Conformation Sampling via Protein Structural Diffusion

Coarse-Grained Conformational Sampling of Protein Structure Improves the Fit to Experimental Hydrogen-Exchange Data

Conformation search across multiple-level potential-energy- surfaces (CSAMP): A strategy for accurate prediction of protein-ligand binding structures.

Deep-Learning-Assisted Enhanced Sampling for Exploring Molecular Conformational Changes

Exploring Protein Conformational Changes Using a Large‐Scale Biophysical Sampling Augmented Deep Learning Strategy

Real-time structural signatures of protein allostery by electrically- and mechanically-coupled sensing