Abstract:Molecular dynamics (MD) simulations are ideally suited to describe conformational ensembles of biomolecules such as proteins and nucleic acids. Microsecond-long simulations are now routine, facilitated by the emergence of graphical processing units. Processing such ensembles on the basis of statistical mechanics can bring insights about different biologically relevant states, their representative structures, states, and even dynamics between states. Clustering, which groups objects based on structural similarity, is typically used to process ensembles, leading to different states, their populations, and the identification of representative structures. For some purposes, such as in protein structure prediction, we are interested in identifying the representative structure that is more similar to the native state of the protein. The traditional pipeline combines hierarchical clustering for clustering and selecting the cluster centroid as representative of the cluster. However, even when the first cluster represents the native basin, the centroid can be several angstroms away in RMSD from the native state – and many other structures inside this cluster could be better choices of representative structures, reducing the need for protein structure refinement. In this study, we developed a module—Protein Retrieval via Integrative Molecular Ensemble (PRIME), that consists of tools to determine the most prevalent states in an ensemble using extended continuous similarity. PRIME is integrated with our Molecular Dynamics Analysis with -ary Clustering Ensembles (MDANCE) package and can be used as a post-processing tool for arbitrary clustering algorithms, compatible with several MD suites. PRIME was validated with ensembles of different protein and protein complex systems for their ability to reliably identify the most native-like state, which we compare to their experimental structure, and to the traditional approach. Systems were chosen to represent different degrees of difficulty such as folding processes and binding which require large conformational changes. PRIME predictions produced structures that when aligned to the experimental structure were better superposed (lower RMSD). A further benefit of PRIME is its linear scaling – rather than the traditional O( ) traditionally associated to comparisons of elements in a set.

Reverse mapping of coarse grained polyglutamine conformations from PRIME20 sampling

Learning Effective Molecular Models from Experimental Observables.

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

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

Enhanced Sampling Method with Coarse Graining of Conformational Space.

Generative Coarse-Graining of Molecular Conformations

Adaptive Resolution Force Probe Simulations: Coarse Graining in the Ideal Gas Approximation

Chemically Transferable Generative Backmapping of Coarse-Grained Proteins

Adversarial reverse mapping of equilibrated condensed-phase molecular structures

Protocol for the development of coarse-grained structures for macromolecular simulation using GROMACS

Multiscale Simulation of Polyglutamine and the Effect of Neighboring Amino Acids on Oligomerization

Formalizing Coarse-Grained Representations of Anisotropic Interactions at Multimeric Protein Interfaces Using Virtual Sites

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

Using a Coarse-Grained Modeling Framework to Identify Oligomeric Motifs with Tunable Secondary Structure

Evaluation of Essential Dynamics and Fixed-Length Coarse Graining for Multidomain Proteins

Dynamic coarse-graining fills the gap between atomistic simulations and experimental investigations of mechanical unfolding

Protein Retrieval via Integrative Molecular Ensembles (PRIME) through extended similarity indices

MD-assisted Refinement of X-Ray Coordinates

A coarse‐grained model for disordered and multi‐domain proteins

Iterative derivation of effective potentials to sample the conformational space of proteins at atomistic scale

Efficient and Reproducible Folding Simulations of the Trp-cage Protein with Multiscale Molecular Dynamics