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.

Phyletic Relationships of Protein Structures Based on Spatial Preference of Residues

A Multiple Criteria Framework for 3D Protein Structure Similarity Retrieval

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

Dynamics of the tuning process between singers

Muscle-3D: scalable multiple protein structure alignment

A Method for Measuring Protein Structure Similarity Based on the Molecular Inner Spatial Density Distribution

Revealing divergent evolution, identifying circular permutations and detecting active-sites by protein structure comparison

Phylogenetic Profiles as a Unified Framework for Measuring Protein Structure, Function and Evolution

Grouping of Amino Acids and Recognition of Protein Structurally Conserved Regions by Reduced Alphabets of Amino Acids

Local sequence-structure relationships in proteins

A new method on reconstructing protein structure from NOESY distances

Predicting Protein Interaction Sites from Residue Spatial Sequence Profile and Evolution Rate

Protein Structure Determination Using a Riemannian Approach

A general substitution matrix for structural phylogenetics.

Mining Protein Sequence Motifs Representing Common 3D Structures.

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

Comparing Graph Representations of Protein Structure for Mining Family-Specific Residue-Based Packing Motifs

MRFalign: Protein Homology Detection through Alignment of Markov Random Fields

MM-align: a quick algorithm for aligning multiple-chain protein complex structures using iterative dynamic programming

Protein Structure Prediction by Protein Alignments

SeqMapPDB: A Standalone Pipeline to Identify Representative Structures of Protein Sequences and Mapping Residue Indices in Real-Time at Proteome Scale