Shared Signature Dynamics Tempered by Local Fluctuations Enables Fold Adaptability and Specificity

WatFinder: a ProDy tool for protein-water interactions

SUMMARY

We introduce WatFinder, a tool designed to identify and visualize protein-water interactions (water bridges, water-mediated associations, or water channels, fluxes, and clusters) relevant to protein stability, dynamics, and function. WatFinder is integrated into ProDy, a Python API broadly used for structure-based prediction of protein dynamics. WatFinder provides a suite of functions for generating raw data as well as outputs from statistical analyses. The ProDy framework facilitates comprehensive automation and efficient analysis of the ensembles of structures resolved for a given protein or the time-evolved conformations from simulations in explicit water, as illustrated in five case studies presented in the Supplementary Material.

AVAILABILITY AND IMPLEMENTATION

ProDy is open-source and freely available under MIT License from https://github.com/ProDy/ProDy.

Genome structural dynamics: insights from Gaussian network analysis of Hi-C data

Characterization of the spatiotemporal properties of the chromatin is essential to gaining insights into the physical bases of gene co-expression, transcriptional regulation and epigenetic modifications. The Gaussian network model (GNM) has proven in recent work to serve as a useful tool for modeling chromatin structural dynamics, using as input high-throughput chromosome conformation capture data. We focus here on the exploration of the collective dynamics of chromosomal structures at hierarchical levels of resolution, from single gene loci to topologically associating domains or entire chromosomes. The GNM permits us to identify long-range interactions between gene loci, shedding light on the role of cross-correlations between distal regions of the chromosomes in regulating gene expression. Notably, GNM analysis performed across diverse cell lines highlights the conservation of the global/cooperative movements of the chromatin across different types of cells. Variations driven by localized couplings between genomic loci, on the other hand, underlie cell differentiation, underscoring the significance of the four-dimensional properties of the genome in defining cellular identity. Finally, we demonstrate the close relation between the cell type-dependent mobility profiles of gene loci and their gene expression patterns, providing a clear demonstration of the role of chromosomal 4D features in defining cell-specific differential expression of genes.

Repurposing of SARS-CoV-2 compounds against Marburg Virus using MD simulation, mm/GBSA, PCA analysis, and free energy landscape

The significant mortality rate associated with Marburg virus infection made it the greatest hazard among infectious diseases. Drug repurposing using in silico methods has been crucial in identifying potential compounds that could prevent viral replication by targeting the virus's primary proteins. This study aimed at repurposing the drugs of SARS-CoV-2 for identifying potential candidates against the matrix protein VP40 of the Marburg virus. Virtual screening was performed where the control compound, Nilotinib, showed a binding score of -9.99 kcal/mol. Based on binding scores, hit compounds 9549298, 11960895, 44545852, 51039094, and 89670174 were selected that had a lower binding score than the control. Subsequent molecular dynamics (MD) simulation revealed that compound 9549298 consistently formed a hydrogen bond with the residue Gln290. This was observed both in molecular docking and MD simulation poses, indicating a strong and significant interaction with the protein. 11960895 had the most stable and consistent RMSD pattern exhibited in 100 ns simulation, while 9549298 had the most identical RMSD plot compared to the control molecule. MM/PBSA analysis showed that the binding free energy (ΔG) of 9549298 and 11960895 was lower than the control, with -30.84 and -38.86 kcal/mol, respectively. It was observed by the PCA (principal component analysis) and FEL (free energy landscape) analysis that compounds 9549298 and 11960895 had lesser conformational variation. Overall, this study proposed 9549298 and 11960895 as potential binders of VP40 MARV that can cause its inhibition, however it inherently lacks experimental validation. Furthermore, the study proposes in-vitro experiments as the next step to validate these computational findings, offering a practical approach to further explore these compounds' potential as antiviral agents.Communicated by Ramaswamy H. Sarma.

Tandem-repeat proteins conformational mechanics are optimized to facilitate functional interactions and complexations

The architectures of tandem-repeat proteins are distinct from those of globular proteins. Individual modules, each comprising small structural motifs of 20-40 residues, are arrayed in a quasi one-dimensional fashion to form striking, elongated, horseshoe-like, and superhelical architectures, stabilized solely by short-range interaction. The spring-like shapes of repeat arrays point to elastic modes of action, and these proteins function as adapter molecules or 'hubs,' propagating signals within multi-subunit assemblies in diverse biological contexts. This flexibility is apparent in the dramatic variability observed in the structures of tandem-repeat proteins in different complexes. Here, using computational analysis, we demonstrate the striking ability of just one or a few global motions to recapitulate these structures. These findings show how the mechanics of repeat arrays are robustly enabled by their unique architecture. Thus, the repeating architecture has been optimized by evolution to favor functional modes of motions. The global motions enabling functional transitions can be fully visualized at http://bahargroup.org/tr_web.

Scipion-EM-ProDy: A Graphical Interface for the ProDy Python Package within the Scipion Workflow Engine Enabling Integration of Databases, Simulations and Cryo-Electron Microscopy Image Processing

Macromolecular assemblies, such as protein complexes, undergo continuous structural dynamics, including global reconfigurations critical for their function. Two fast analytical methods are widely used to study these global dynamics, namely elastic network model normal mode analysis and principal component analysis of ensembles of structures. These approaches have found wide use in various computational studies, driving the development of complex pipelines in several software packages. One common theme has been conformational sampling through hybrid simulations incorporating all-atom molecular dynamics and global modes of motion. However, wide functionality is only available for experienced programmers with limited capabilities for other users. We have, therefore, integrated one popular and extensively developed software for such analyses, the ProDy Python application programming interface, into the Scipion workflow engine. This enables a wider range of users to access a complete range of macromolecular dynamics pipelines beyond the core functionalities available in its command-line applications and the normal mode wizard in VMD. The new protocols and pipelines can be further expanded and integrated into larger workflows, together with other software packages for cryo-electron microscopy image analysis and molecular simulations. We present the resulting plugin, Scipion-EM-ProDy, in detail, highlighting the rich functionality made available by its development.

Global Survey of Protein Dynamic Properties

Using tools developed to study the dynamic bioinformatics of proteins, we are able to study the dynamic characteristics of very large numbers of protein sequences simultaneously. We study herein the distribution of protein sequences in a space determined by sequence mobility. It is shown that there are statistically significant differences in mobility distribution between folded sequences of different structural classes and between those and sequences of intrinsically disordered proteins. It is also shown that the several regions of mobility space differ significantly with respect to structural makeup. Helical proteins are shown to have distinctive dynamic characteristics at both extremes of the mobility spectrum.

Functional insight into Cordyceps militaris sugar transporters by structure modeling, network analysis and allosteric regulation

Insights into the structures, functions and dynamics of Cordyceps militaris (C. militaris) sugar transporters are necessary for understanding their versatile metabolic capability for fungal growth. The sequence-function relationship study of 85 C. militaris sugar transporters showed that there is a gap between phylogenetic-based subfamily classification and their functions. Beyond protein sequences, structural modeling and principal component analysis of the structural ensemble revealed the different folds of the Car and Org subfamilies. Performing channel detection and network analysis found that the Alp and Hex subfamilies can be specifically distinguished from others by the betweenness of channel residues. Signature dynamics analysis further suggested that the Hex subfamily demonstrates different dynamics, with high flexibility at the H1 region in TM11. Furthermore, the H1 region as an allosteric site was examined by network parameter calculations that guided allosteric pathways between this region and the channel cavity. Together with gene expression data of C. militaris, e.g., Hex06741 in the Hex subfamily, it was promisingly expressed when sugar utilization was altered. This work demonstrates an in silico framework for investigating C. militaris sugar transporters as an example case study of the allosteric activity of the Hex subfamily and can facilitate sugar transporter engineering design that can further optimize the preferable sugar utilization and fermentation process of C. militaris.

Structural Dynamics Predominantly Determine the Adaptability of Proteins to Amino Acid Deletions

The insertion or deletion (indel) of amino acids has a variety of effects on protein function, ranging from disease-forming changes to gaining new functions. Despite their importance, indels have not been systematically characterized towards protein engineering or modification goals. In the present work, we focus on deletions composed of multiple contiguous amino acids (mAA-dels) and their effects on the protein (mutant) folding ability. Our analysis reveals that the mutant retains the native fold when the mAA-del obeys well-defined structural dynamics properties: localization in intrinsically flexible regions, showing low resistance to mechanical stress, and separation from allosteric signaling paths. Motivated by the possibility of distinguishing the features that underlie the adaptability of proteins to mAA-dels, and by the rapid evaluation of these features using elastic network models, we developed a positive-unlabeled learning-based classifier that can be adopted for protein design purposes. Trained on a consolidated set of features, including those reflecting the intrinsic dynamics of the regions where the mAA-dels occur, the new classifier yields a high recall of 84.3% for identifying mAA-dels that are stably tolerated by the protein. The comparative examination of the relative contribution of different features to the prediction reveals the dominant role of structural dynamics in enabling the adaptation of the mutant to mAA-del without disrupting the native fold.

Correlated Evolution of Low-Frequency Vibrations and Function in Enzymes

Previous studies of the flexibility of ancestral proteins suggest that proteins evolve their function by altering their native state ensemble. Here, we propose a more direct method to analyze such changes during protein evolution by comparing thermally activated vibrations at frequencies below 6 THz, which report on the dynamics of collective protein modes. We analyzed the backbone vibrational density of states of ancestral and extant β-lactamases and thioredoxins and observed marked changes in the vibrational spectrum in response to evolution. Coupled with previously observed changes in protein flexibility, the observed shifts of vibrational mode densities suggest that protein dynamics and dynamical allostery are critical factors for the evolution of enzymes with specialized catalytic and biophysical properties.

Activation and Speciation Mechanisms in Class A GPCRs

Accurate development of allosteric modulators of GPCRs require a thorough assessment of their sequence, structure, and dynamics, toward gaining insights into their mechanisms of actions shared by family members, as well as dynamic features that distinguish subfamilies. Building on recent progress in the characterization of the signature dynamics of proteins, we analyzed here a dataset of 160 Class A GPCRs to determine their sequence similarities, structural landscape, and dynamic features across different species (human, bovine, mouse, squid, and rat), different activation states (active/inactive), and different subfamilies. The two dominant directions of variability across experimentally resolved structures, identified by principal component analysis of the dataset, shed light to cooperative mechanisms of activation, subfamily differentiation, and speciation of Class A GPCRs. The analysis reveals the functional significance of the conformational flexibilities of specific structural elements, including: the dominant role of the intracellular loop 3 (ICL3) together with the cytoplasmic ends of the adjoining helices TM5 and TM6 in enabling allosteric activation; the role of particular structural motifs at the extracellular loop 2 (ECL2) connecting TM4 and TM5 in binding ligands specific to different subfamilies; or even the differentiation of the N-terminal conformation across different species. Detailed analyses of the modes of motions accessible to the members of the dataset and their variations across members demonstrate how the active and inactive states of GPCRs obey distinct conformational dynamics. The collective fluctuations of the GPCRs are robustly defined in the active state, while the inactive conformers exhibit broad variance among members.

In Search of a Dynamical Vocabulary: A Pipeline to Construct a Basis of Shared Traits in Large-Scale Motions of Proteins

The paradigmatic sequence–structure–dynamics–function relation in proteins is currently well established in the scientific community; in particular, a large effort has been made to probe the first connection, indeed providing convincing evidence of its strength and rationalizing it in a quantitative and general framework. In contrast, however, the role of dynamics as a link between structure and function has eluded a similarly clear-cut verification and description. In this work, we propose a pipeline aimed at building a basis for the quantitative characterization of the large-scale dynamics of a set of proteins, starting from the sole knowledge of their native structures. The method hinges on a dynamics-based clusterization, which allows a straightforward comparison with structural and functional protein classifications. The resulting basis set, obtained through the application to a group of related proteins, is shown to reproduce the salient large-scale dynamical features of the dataset. Most interestingly, the basis set is shown to encode the fluctuation patterns of homologous proteins not belonging to the initial dataset, thus highlighting the general applicability of the pipeline used to build it.

Protein dynamics developments for the large scale and cryoEM: case study of ProDy 2.0

Cryo-electron microscopy (cryoEM) has become a well established technique with the potential to produce structures of large and dynamic supramolecular complexes that are not amenable to traditional approaches for studying structure and dynamics. The size and low resolution of such molecular systems often make structural modelling and molecular dynamics simulations challenging and computationally expensive. This, together with the growing wealth of structural data arising from cryoEM and other structural biology methods, has driven a trend in the computational biophysics community towards the development of new pipelines for analysing global dynamics using coarse-grained models and methods. At the centre of this trend has been a return to elastic network models, normal mode analysis (NMA) and ensemble analyses such as principal component analysis, and the growth of hybrid simulation methodologies that make use of them. Here, this field is reviewed with a focus on ProDy, the Python application programming interface for protein dynamics, which has been developed over the last decade. Two key developments in this area are highlighted: (i) ensemble NMA towards extracting and comparing the signature dynamics of homologous structures, aided by the recent SignDy pipeline, and (ii) pseudoatom fitting for more efficient global dynamics analyses of large and low-resolution supramolecular assemblies from cryoEM, revisited in the CryoDy pipeline. It is believed that such a renewal and extension of old models and methods in new pipelines will be critical for driving the field forward into the next cryoEM revolution.

A Comparative Evaluation of the Structural and Dynamic Properties of Insect Odorant Binding Proteins

Insects devote a major part of their metabolic resources to the production of odorant binding proteins (OBPs). Although initially, these proteins were implicated in the solubilisation, binding and transport of semiochemicals to olfactory receptors, it is now recognised that they may play diverse, as yet uncharacterised, roles in insect physiology. The structures of these OBPs, the majority of which are known as “classical” OBPs, have shed some light on their potential functional roles. However, the dynamic properties of these proteins have received little attention despite their functional importance. Structural dynamics are encoded in the native protein fold and enable the adaptation of proteins to substrate binding. This paper provides a comparative review of the structural and dynamic properties of OBPs, making use of sequence/structure analysis, statistical and theoretical physics-based methods. It provides a new layer of information and additional methodological tools useful in unravelling the relationship between structure, dynamics and function of insect OBPs. The dynamic properties of OBPs, studied by means of elastic network models, reflect the similarities/dissimilarities observed in their respective structures and provides insights regarding protein motions that may have important implications for ligand recognition and binding. Furthermore, it was shown that the OBPs studied in this paper share conserved structural ‘core’ that may be of evolutionary and functional importance.

Sampling of Protein Conformational Space Using Hybrid Simulations: A Critical Assessment of Recent Methods

Recent years have seen several hybrid simulation methods for exploring the conformational space of proteins and their complexes or assemblies. These methods often combine fast analytical approaches with computationally expensive full atomic molecular dynamics (MD) simulations with the goal of rapidly sampling large and cooperative conformational changes at full atomic resolution. We present here a systematic comparison of the utility and limits of four such hybrid methods that have been introduced in recent years: MD with excited normal modes (MDeNM), collective modes-driven MD (CoMD), and elastic network model (ENM)-based generation, clustering, and relaxation of conformations (ClustENM) as well as its updated version integrated with MD simulations (ClustENMD). We analyzed the predicted conformational spaces using each of these four hybrid methods, applied to four well-studied proteins, triosephosphate isomerase (TIM), 3-phosphoglycerate kinase (PGK), HIV-1 protease (PR) and HIV-1 reverse transcriptase (RT), which provide extensive ensembles of experimental structures for benchmarking and comparing the methods. We show that a rigorous multi-faceted comparison and multiple metrics are necessary to properly assess the differences between conformational ensembles and provide an optimal protocol for achieving good agreement with experimental data. While all four hybrid methods perform well in general, being especially useful as computationally efficient methods that retain atomic resolution, the systematic analysis of the same systems by these four hybrid methods highlights the strengths and limitations of the methods and provides guidance for parameters and protocols to be adopted in future studies.

Uncovering the Dominant Motion Modes of Allosteric Regulation Improves Allosteric Site Prediction

Allostery is an important mechanism that biological systems use to regulate function at a distant site. Allosteric drugs have attracted much attention in recent years due to their high specificity and the possibility of overcoming existing drug-resistant mutations. However, the discovery of allosteric drugs remains challenging as allosteric regulation mechanisms are not clearly understood and allosteric sites cannot be accurately predicted. In this study, we analyzed the dominant modes that determine motion correlations between allosteric and orthosteric sites using the Gaussian network model and found that motion correlations between allosteric and orthosteric sites are dominated by either fast or slow vibrational modes. This dependence of modes results from the relative locations of the two sites and local secondary structures. Based on these analyses, we developed CorrSite2.0 to predict allosteric sites by taking the maximum of the Z-scores calculated from using either slow or fast modes. The prediction accuracy of CorrSite2.0 outperformed other commonly used allosteric site prediction methods with prediction accuracy over 90.0%. Our study uncovers the relationship of protein structure, dynamics, and allosteric regulation and demonstrates that using the dominant motion modes greatly improves allosteric site prediction accuracy. CorrSite2.0 has been integrated into the CavityPlus web server, which can be accessed at http://www.pkumdl.cn/cavityplus. CorrSite2.0 provides a powerful and user-friendly tool for allosteric drug and protein design.

Structural dynamics in the evolution of a bilobed protein scaffold

Proteins conduct numerous complex biological functions by use of tailored structural dynamics. The molecular details of how these emerged from ancestral peptides remains mysterious. How does nature utilize the same repertoire of folds to diversify function? To shed light on this, we analyzed bilobed proteins with a common structural core, which is spread throughout the tree of life and is involved in diverse biological functions such as transcription, enzymatic catalysis, membrane transport, and signaling. We show here that the structural dynamics of the structural core differentiate predominantly via terminal additions during a long-period evolution. This diversifies substrate specificity and, ultimately, biological function.

Novel biophysical tools allow the structural dynamics of proteins and the regulation of such dynamics by binding partners to be explored in unprecedented detail. Although this has provided critical insights into protein function, the means by which structural dynamics direct protein evolution remain poorly understood. Here, we investigated how proteins with a bilobed structure, composed of two related domains from the periplasmic-binding protein–like II domain family, have undergone divergent evolution, leading to adaptation of their structural dynamics. We performed a structural analysis on ∼600 bilobed proteins with a common primordial structural core, which we complemented with biophysical studies to explore the structural dynamics of selected examples by single-molecule Förster resonance energy transfer and Hydrogen–Deuterium exchange mass spectrometry. We show that evolutionary modifications of the structural core, largely at its termini, enable distinct structural dynamics, allowing the diversification of these proteins into transcription factors, enzymes, and extracytoplasmic transport-related proteins. Structural embellishments of the core created interdomain interactions that stabilized structural states, reshaping the active site geometry, and ultimately altered substrate specificity. Our findings reveal an as-yet-unrecognized mechanism for the emergence of functional promiscuity during long periods of evolution and are applicable to a large number of domain architectures.

Evolution of dynamical networks enhances catalysis in a designer enzyme

Activation heat capacity is emerging as a crucial factor in enzyme thermoadaptation, as shown by the non-Arrhenius behaviour of many natural enzymes. However, its physical origin and relationship to the evolution of catalytic activity remain uncertain. Here we show that directed evolution of a computationally designed Kemp eliminase reshapes protein dynamics, which gives rise to an activation heat capacity absent in the original design. These changes buttress transition-state stabilization. Extensive molecular dynamics simulations show that evolution results in the closure of solvent-exposed loops and a better packing of the active site. Remarkably, this gives rise to a correlated dynamical network that involves the transition state and large parts of the protein. This network tightens the transition-state ensemble, which induces a negative activation heat capacity and non-linearity in the activity-temperature dependence. Our results have implications for understanding enzyme evolution and suggest that selectively targeting the conformational dynamics of the transition-state ensemble by design and evolution will expedite the creation of novel enzymes.

Probing Universal Protein Dynamics Using Hydrogen-Deuterium Exchange Mass Spectrometry-Derived Residue-Level Gibbs Free Energy

Hydrogen-deuterium exchange mass spectrometry (HDX-MS) is a powerful technique to monitor protein intrinsic dynamics. The technique provides high-resolution information on how protein intrinsic dynamics are altered in response to biological signals, such as ligand binding, oligomerization, or allosteric networks. However, identification, interpretation, and visualization of such events from HDX-MS data sets is challenging as these data sets consist of many individual data points collected across peptides, time points, and experimental conditions. Here, we present PyHDX, an open-source Python package and webserver, that allows the user to batch extract the universal quantity Gibbs free energy at residue levels over multiple protein conditions and homologues. The output is directly visualized on a linear map or 3D structures or is exported as .csv files or PyMOL scripts.

Dynamics-Evolution Correspondence in Protein Structures

The genotype-phenotype mapping of proteins is a fundamental question in structural biology. In this Letter, with the analysis of a large dataset of proteins from hundreds of protein families, we quantitatively demonstrate the correlations between the noise-induced protein dynamics and mutation-induced variations of native structures, indicating the dynamics-evolution correspondence of proteins. Based on the investigations of the linear responses of native proteins, the origin of such a correspondence is elucidated. It is essential that the noise- and mutation-induced deformations of the proteins are restricted on a common low-dimensional subspace, as confirmed from the data. These results suggest an evolutionary mechanism of the proteins gaining both dynamical flexibility and evolutionary structural variability.

Computational approaches to predict protein functional families and functional sites

Understanding the mechanisms of protein function is indispensable for many biological applications, such as protein engineering and drug design. However, experimental annotations are sparse, and therefore, theoretical strategies are needed to fill the gap. Here, we present the latest developments in building functional subclassifications of protein superfamilies and using evolutionary conservation to detect functional determinants, for example, catalytic-, binding- and specificity-determining residues important for delineating the functional families. We also briefly review other features exploited for functional site detection and new machine learning strategies for combining multiple features.

Gating and modulation of a hetero-octameric AMPA glutamate receptor

AMPA receptors (AMPARs) mediate the majority of excitatory transmission in the brain and enable the synaptic plasticity that underlies learning¹. A diverse array of AMPAR signalling complexes are established by receptor auxiliary subunits, which associate with the AMPAR in various combinations to modulate trafficking, gating and synaptic strength². However, their mechanisms of action are poorly understood. Here we determine cryo-electron microscopy structures of the heteromeric GluA1-GluA2 receptor assembled with both TARP-γ8 and CNIH2, the predominant AMPAR complex in the forebrain, in both resting and active states. Two TARP-γ8 and two CNIH2 subunits insert at distinct sites beneath the ligand-binding domains of the receptor, with site-specific lipids shaping each interaction and affecting the gating regulation of the AMPARs. Activation of the receptor leads to asymmetry between GluA1 and GluA2 along the ion conduction path and an outward expansion of the channel triggers counter-rotations of both auxiliary subunit pairs, promoting the active-state conformation. In addition, both TARP-γ8 and CNIH2 pivot towards the pore exit upon activation, extending their reach for cytoplasmic receptor elements. CNIH2 achieves this through its uniquely extended M2 helix, which has transformed this endoplasmic reticulum-export factor into a powerful AMPAR modulator that is capable of providing hippocampal pyramidal neurons with their integrative synaptic properties.

Fast computational mutation-response scanning of proteins

Studying the effect of perturbations on protein structure is a basic approach in protein research. Important problems, such as predicting pathological mutations and understanding patterns of structural evolution, have been addressed by computational simulations that model mutations using forces and predict the resulting deformations. In single mutation-response scanning simulations, a sensitivity matrix is obtained by averaging deformations over point mutations. In double mutation-response scanning simulations, a compensation matrix is obtained by minimizing deformations over pairs of mutations. These very useful simulation-based methods may be too slow to deal with large proteins, protein complexes, or large protein databases. To address this issue, I derived analytical closed formulas to calculate the sensitivity and compensation matrices directly, without simulations. Here, I present these derivations and show that the resulting analytical methods are much faster than their simulation counterparts.

ProDy 2.0: Increased Scale and Scope after 10 Years of Protein Dynamics Modelling with Python

Summary

ProDy, an integrated application programming interface developed for modelling and analysing protein dynamics, has significantly evolved in recent years in response to the growing data and needs of the computational biology community. We present major developments that led to ProDy 2.0: (i) improved interfacing with databases and parsing new file formats, (ii) SignDy for signature dynamics of protein families, (iii) CryoDy for collective dynamics of supramolecular systems using cryo-EM density maps and (iv) essential site scanning analysis for identifying sites essential to modulating global dynamics.

Availability and implementation

ProDy is open-source and freely available under MIT License from https://github.com/prody/ProDy.

Supplementary information

Supplementary data are available at Bioinformatics online.

State-dependent sequential allostery exhibited by chaperonin TRiC/CCT revealed by network analysis of Cryo-EM maps

The eukaryotic chaperonin TRiC/CCT plays a major role in assisting the folding of many proteins through an ATP-driven allosteric cycle. Recent structures elucidated by cryo-electron microscopy provide a broad view of the conformations visited at various stages of the chaperonin cycle, including a sequential activation of its subunits in response to nucleotide binding. But we lack a thorough mechanistic understanding of the structure-based dynamics and communication properties that underlie the TRiC/CCT machinery. In this study, we present a computational methodology based on elastic network models adapted to cryo-EM density maps to gain a deeper understanding of the structure-encoded allosteric dynamics of this hexadecameric machine. We have analysed several structures of the chaperonin resolved in different states toward mapping its conformational landscape. Our study indicates that the overall architecture intrinsically favours cooperative movements that comply with the structural variabilities observed in experiments. Furthermore, the individual subunits CCT1-CCT8 exhibit state-dependent sequential events at different states of the allosteric cycle. For example, in the ATP-bound state, subunits CCT5 and CCT4 selectively initiate the lid closure motions favoured by the overall architecture; whereas in the apo form of the heteromer, the subunit CCT7 exhibits the highest predisposition to structural change. The changes then propagate through parallel fluxes of allosteric signals to neighbours on both rings. The predicted state-dependent mechanisms of sequential activation provide new insights into TRiC/CCT intra- and inter-ring signal transduction events.

Equally Weighted Multiscale Elastic Network Model and Its Comparison with Traditional and Parameter-Free Models

Dynamical properties of proteins play an essential role in their function exertion. The elastic network model (ENM) is an effective and efficient tool in characterizing the intrinsic dynamical properties encoded in biomacromolecule structures. The Gaussian network model (GNM) and anisotropic network model (ANM) are the two often-used ENM models. Here, we introduce an equally weighted multiscale ENM (equally weighted mENM) based on the original mENM (denoted as mENM), in which fitting weights of Kirchhoff/Hessian matrixes in mENM are removed since they neglect the details of pairwise interactions. Then, we perform its comparison with the mENM, traditional ENM, and parameter-free ENM (pfENM) in reproducing dynamical properties for the six representative proteins whose molecular dynamics (MD) trajectories are available in http://mmb.pcb.ub.es/MoDEL/. In the results, for B-factor prediction, mENM performs best, while the equally weighted mENM performs also well, better than the traditional ENM and pfENM models. As to the dynamical cross-correlation map calculation, mENM performs worst, while the results produced from the equally weighted mENM and pfENM models are close to those from MD trajectories with the latter a little better than the former. Furthermore, encouragingly, the equally weighted mANM displays the best performance in capturing the functional motional modes, followed by pfANM and traditional ANM models, while the mANM fails in all the cases. This work is helpful for strengthening the understanding of the elastic network model and provides a valuable guide for researchers to utilize the model to explore protein dynamics.

Rhapsody: predicting the pathogenicity of human missense variants

MOTIVATION

The biological effects of human missense variants have been studied experimentally for decades but predicting their effects in clinical molecular diagnostics remains challenging. Available computational tools are usually based on the analysis of sequence conservation and structural properties of the mutant protein. We recently introduced a new machine learning method that demonstrated for the first time the significance of protein dynamics in determining the pathogenicity of missense variants.

RESULTS

Here, we present a new interface (Rhapsody) that enables fully automated assessment of pathogenicity, incorporating both sequence coevolution data and structure- and dynamics-based features. Benchmarked against a dataset of about 20 000 annotated variants, the methodology is shown to outperform well-established and/or advanced prediction tools. We illustrate the utility of Rhapsody by in silico saturation mutagenesis studies of human H-Ras, phosphatase and tensin homolog and thiopurine S-methyltransferase.

AVAILABILITY AND IMPLEMENTATION

The new tool is available both as an online webserver at http://rhapsody.csb.pitt.edu and as an open-source Python package (GitHub repository: https://github.com/prody/rhapsody; PyPI package installation: pip install prody-rhapsody). Links to additional resources, tutorials and package documentation are provided in the 'Python package' section of the website.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Adaptability and specificity: how do proteins balance opposing needs to achieve function?

Many proteins select from a small repertoire of 3-dimensional folds retained over evolutional timescales and recruited for different functions, with changes in local structure and sequence to enable specificity. Recent studies have revealed the evolutionary constraints on protein dynamics to achieve function. The significance of protein dynamics in simultaneously satisfying conformational flexibility/malleability and stability/precision requirements becomes clear upon dissecting the spectrum of equilibrium motions accessible to fold families. Accessibility to highly conserved global modes of motions shared by family members, to low-to-intermediate-frequency modes that distinguish subfamilies and confer specificity, and to conserved high-frequency modes ensuring chemical precision and core stability underlies functional specialization while exploiting highly versatile folds. These design principles are illustrated for the family of PDZ domains.

Towards gaining sight of multiscale events: utilizing network models and normal modes in hybrid methods

With the explosion of normal mode analyses (NMAs) based on elastic network models (ENMs) in the last decade, and the proven precision of MD simulations for visualizing interactions at atomic scale, many hybrid methods have been proposed in recent years. These aim at exploiting the best of both worlds: the atomic precision of MD that often fall short of exploring time and length scales of biological interest, and the capability of ENM-NMA to predict the cooperative and often functional rearrangements of large structures and assemblies, albeit at low resolution. We present an overview of recent progress in the field with examples of successful applications highlighting the utility of such hybrid methods and pointing to emerging future directions guided by advances in experimental characterization of biomolecular systems structure and dynamics.

Recent advances suggest increased influence of selective pressure in allostery

Allosteric regulation of protein functions is ubiquitous in organismal biology, but the principles governing its evolution are not well understood. Here we discuss recent studies supporting the large-scale existence of latent allostery in ancestor proteins of superfamilies. As suggested, the evolution of allostery could be driven by the need for specificity in paralogs of slow evolving protein complexes with conserved active sites. The same slow evolution is displayed by purifying selection exhibited in allosteric proteins with somatic mutations involved in cancer, where disease-associated mutations are enriched in both orthosteric and allosteric sites. Consequently, disease-associated variants can be used to identify druggable allosteric sites that are specific for paralogs in protein superfamilies with otherwise similar functions.

Differences in the intrinsic spatial dynamics of the chromatin contribute to cell differentiation

Advances in chromosome conformation capture techniques as well as computational characterization of genomic loci structural dynamics open new opportunities for exploring the mechanistic aspects of genome-scale differences across different cell types. We examined here the dynamic basis of variabilities between different cell types by investigating their chromatin mobility profiles inferred from Hi-C data using an elastic network model representation of the chromatin. Our comparative analysis of sixteen cell lines reveals close similarities between chromosomal dynamics across different cell lines on a global scale, but notable cell-specific variations emerge in the detailed spatial mobilities of genomic loci. Closer examination reveals that the differences in spatial dynamics mainly originate from the difference in the frequencies of their intrinsically accessible modes of motion. Thus, even though the chromosomes of different types of cells have access to similar modes of collective movements, not all modes are deployed by all cells, such that the effective mobilities and cross-correlations of genomic loci are cell-type-specific. Comparison with RNA-seq expression data reveals a strong overlap between highly expressed genes and those distinguished by high mobilities in the present study, in support of the role of the intrinsic spatial dynamics of chromatin as a determinant of cell differentiation.

Dynamics-function relationship in the catalytic domains of N-terminal acetyltransferases

N-terminal acetyltransferases (NATs) belong to the superfamily of acetyltransferases. They are enzymes catalysing the transfer of an acetyl group from acetyl coenzyme A to the N-terminus of polypeptide chains. N-terminal acetylation is one of the most common protein modifications. To date, not much is known on the molecular basis for the exclusive substrate specificity of NATs. All NATs share a common fold called GNAT. A characteristic of NATs is the β6β7 hairpin loop covering the active site and forming with the α1α2 loop a narrow tunnel surrounding the catalytic site in which cofactor and polypeptide meet and exchange an acetyl group. We investigated the dynamics-function relationships of all available structures of NATs covering the three domains of Life. Using an elastic network model and normal mode analysis, we found a common dynamics pattern conserved through the GNAT fold; a rigid V-shaped groove formed by the β4 and β5 strands and splitting the fold in two dynamical subdomains. Loops α1α2, β3β4 and β6β7 all show clear displacements in the low frequency normal modes. We characterized the mobility of the loops and show that even limited conformational changes of the loops along the low-frequency modes are able to significantly change the size and shape of the ligand binding sites. Based on the fact that these movements are present in most low-frequency modes, and common to all NATs, we suggest that the α1α2 and β6β7 loops may regulate ligand uptake and the release of the acetylated polypeptide.

The Role of Conformational Dynamics and Allostery in Modulating Protein Evolution

Advances in sequencing techniques and statistical methods have made it possible not only to predict sequences of ancestral proteins but also to identify thousands of mutations in the human exome, some of which are disease associated. These developments have motivated numerous theories and raised many questions regarding the fundamental principles behind protein evolution, which have been traditionally investigated horizontally using the tip of the phylogenetic tree through comparative studies of extant proteins within a family. In this article, we review a vertical comparison of the modern and resurrected ancestral proteins. We focus mainly on the dynamical properties responsible for a protein's ability to adapt new functions in response to environmental changes. Using the Dynamic Flexibility Index and the Dynamic Coupling Index to quantify the relative flexibility and dynamic coupling at a site-specific, single-amino-acid level, we provide evidence that the migration of hinges, which are often functionally critical rigid sites, is a mechanism through which proteins can rapidly evolve. Additionally, we show that disease-associated mutations in proteins often result in flexibility changes even at positions distal from mutational sites, particularly in the modulation of active site dynamics.

Modulation of Toroidal Proteins Dynamics in Favor of Functional Mechanisms upon Ligand Binding

Toroidal proteins serve as molecular machines and play crucial roles in biological processes such as DNA replication and RNA transcription. Despite progress in the structural characterization of several toroidal proteins, we still lack a mechanistic understanding of the significance of their architecture, oligomerization states, and intermolecular interactions in defining their biological function. In this work, we analyze the collective dynamics of toroidal proteins with different oligomerization states, namely, dimeric and trimeric DNA sliding clamps, nucleocapsid proteins (4-, 5-, and 6-mers) and Trp RNA-binding attenuation proteins (11- and 12-mers). We observe common global modes, among which cooperative rolling stands out as a mechanism enabling DNA processivity, and clamshell motions as those underlying the opening/closure of the sliding clamps. Alterations in global dynamics due to complexation with DNA or the clamp loader are shown to assist in enhancing motions to enable robust function. The analysis provides new insights into the differentiation and enhancement of functional motions upon intersubunit and intermolecular interactions.

Domain-mediated interactions for protein subfamily identification

Within a protein family, proteins with the same domain often exhibit different cellular functions, despite the shared evolutionary history and molecular function of the domain. We hypothesized that domain-mediated interactions (DMIs) may categorize a protein family into subfamilies because the diversified functions of a single domain often depend on interacting partners of domains. Here we systematically identified DMI subfamilies, in which proteins share domains with DMI partners, as well as with various functional and physical interaction networks in individual species. In humans, DMI subfamily members are associated with similar diseases, including cancers, and are frequently co-associated with the same diseases. DMI information relates to the functional and evolutionary subdivisions of human kinases. In yeast, DMI subfamilies contain proteins with similar phenotypic outcomes from specific chemical treatments. Therefore, the systematic investigation here provides insights into the diverse functions of subfamilies derived from a protein family with a link-centric approach and suggests a useful resource for annotating the functions and phenotypic outcomes of proteins.

Intrinsic dynamics is evolutionarily optimized to enable allosteric behavior

Allosteric behavior is central to the function of many proteins, enabling molecular machinery, metabolism, signaling and regulation. Recent years have shown that the intrinsic dynamics of allosteric proteins defined by their 3-dimensional architecture or by the topology of inter-residue contacts favors cooperative motions that bear close similarity to structural changes they undergo during their allosteric actions. These conformational motions are usually driven by energetically favorable or soft modes at the low frequency end of the mode spectrum, and they are evolutionarily conserved among orthologs. These observations brought into light evolutionary adaptation mechanisms that help maintain, optimize or regulate allosteric behavior as the evolution from bacterial to higher organisms introduces sequential heterogeneities and structural complexities.

Trimerization of dopamine transporter triggered by AIM-100 binding: Molecular mechanism and effect of mutations

Recent work demonstrated the propensity of dopamine transporters (DATs) to form trimers or higher oligomers, enhanced upon binding a furopyrimidine, AIM-100. AIM-100 binding promotes DAT endocytosis and thereby moderates dopaminergic transmission. Despite the neurobiological significance of these events, the molecular mechanisms that underlie the stabilization of DAT trimer and the key interactions that modulate the trimerization of DAT, and not serotonin transporter SERT, remain unclear. In the present study, we determined three structural models, termed trimer-W238, -C306 and -Y303, for possible trimerization of DATs . To this aim, we used structural data resolved for DAT and its structural homologs that share the LeuT fold, advanced computational modeling and simulations, site-directed mutagenesis experiments and live-cell imaging assays. The models are in accord with the versatility of LeuT fold to stabilize dimeric or higher order constructs. Selected residues show a high propensity to occupy interfacial regions. Among them, D231-W238 in the extracellular loop EL2, including the intersubunit salt-bridge forming pair D231/D232-R237 (not present in SERT) (in trimer-W238), the loop EL3 (trimers-C306 and -Y303), and W497 on the intracellularly exposed IL5 loop (trimer-C306) and its spatial neighbors (e.g. K525) near the C-terminus are computationally predicted and experimentally confirmed to play important roles in enabling the correct folding and/or oligomerization of DATs in the presence of AIM-100. The study suggests the possibility of controlling the effective transport of dopamine by altering the oligomerization state of DAT upon small molecule binding, as a possible intervention strategy to modulate dopaminergic signaling. This article is part of the issue entitled 'Special Issue on Neurotransmitter Transporters'.