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.

Prediction of protein–protein complexes using replica exchange with repulsive scaling

Refinement of Docked Protein-Protein Complexes Using Repulsive Scaling Replica Exchange Simulations

Induced fit with replica exchange improves protein complex structure prediction

Replica Exchange Simulations of Transient Encounter Complexes in Protein–protein Association

Replica exchange molecular dynamics method for protein folding simulation.

Predicting Protein Interactions by Brownian Dynamics Simulations

Replica Exchange Molecular Dynamics: A Practical Application Protocol with Solutions to Common Problems and a Peptide Aggregation and Self-Assembly Example.

Exploring the Binding Mechanism Between Human Profilin (PFN1) and Polyproline-10 Through Binding Mode Screening

Computational prediction of protein–protein binding affinities

Computational Prediction of Heteromeric Protein Complex Disassembly Order with Hybrid Monte Carlo/Molecular Dynamics Simulation

Replica exchange with guided annealing for accelerated sampling of disordered protein conformations

Replica Exchange with Solute Tempering for Protein Conformational Sampling

Accurate Prediction for Protein-Peptide Binding Based on High-Temperature Molecular Dynamics Simulations

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

Reliable protein-protein docking with AlphaFold, Rosetta and replica-exchange

Protein-Protein Docking Using A Brownian Dynamics Simulations Approach

Replica exchange of expanded ensembles: A generalized ensemble approach with enhanced flexibility and parallelizability

Dynamic Protein-Protein Complexes: How Alternative Interactions Create Ensembles and How Solution Nmr and Md Simulations Can Characterize Them

Learning effective amino acid interactions through iterative stochastic techniques

Prediction of protein complex structure: methods and progress

Predicting protein dynamics from structural ensembles