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.

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.

Cooperative mechanics of PR65 scaffold underlies the allosteric regulation of the phosphatase PP2A

PR65, a horseshoe-shaped scaffold composed of 15 HEAT (observed in Huntingtin, elongation factor 3, protein phosphatase 2A, and the yeast kinase TOR1) repeats, forms, together with catalytic and regulatory subunits, the heterotrimeric protein phosphatase PP2A. We examined the role of PR65 in enabling PP2A enzymatic activity with computations at various levels of complexity, including hybrid approaches that combine full-atomic and elastic network models. Our study points to the high flexibility of this scaffold allowing for end-to-end distance fluctuations of 40-50 Å between compact and extended conformations. Notably, the intrinsic dynamics of PR65 facilitates complexation with the catalytic subunit and is retained in the PP2A complex enabling PR65 to engage the two domains of the catalytic subunit and provide the mechanical framework for enzymatic activity, with support from the regulatory subunit. In particular, the intra-repeat coils at the C-terminal arm play an important role in allosterically mediating the collective dynamics of PP2A, pointing to target sites for modulating PR65 function.

Parameterizing elastic network models to capture the dynamics of proteins

Coarse-grained normal mode analyses of protein dynamics rely on the idea that the geometry of a protein structure contains enough information for computing its fluctuations around its equilibrium conformation. This geometry is captured in the form of an elastic network (EN), namely a network of edges between its residues. The normal modes of a protein are then identified with the normal modes of its EN. Different approaches have been proposed to construct ENs, focusing on the choice of the edges that they are comprised of, and on their parameterizations by the force constants associated with those edges. Here we propose new tools to guide choices on these two facets of EN. We study first different geometric models for ENs. We compare cutoff-based ENs, whose edges have lengths that are smaller than a cutoff distance, with Delaunay-based ENs and find that the latter provide better representations of the geometry of protein structures. We then derive an analytical method for the parameterization of the EN such that its dynamics leads to atomic fluctuations that agree with experimental B-factors. To limit overfitting, we attach a parameter referred to as flexibility constant to each atom instead of to each edge in the EN. The parameterization is expressed as a non-linear optimization problem whose parameters describe both rigid-body and internal motions. We show that this parameterization leads to improved ENs, whose dynamics mimic MD simulations better than ENs with uniform force constants, and reduces the number of normal modes needed to reproduce functional conformational changes.

The amounts of thermal vibrations and static disorder in protein X-ray crystallographic B-factors

Crystallographic B-factors provide direct dynamical information on the internal mobility of proteins that is closely linked to function, and are also widely used as a benchmark in assessing elastic network models. A significant question in the field is: what is the exact amount of thermal vibrations in protein crystallographic B-factors? This work sets out to answer this question. First, we carry out a thorough, statistically sound analysis of crystallographic B-factors of over 10 000 structures. Second, by employing a highly accurate all-atom model based on the well-known CHARMM force field, we obtain computationally the magnitudes of thermal vibrations of nearly 1000 structures. Our key findings are: (i) the magnitude of thermal vibrations, surprisingly, is nearly protein-independent, as a corollary to the universality for the vibrational spectra of globular proteins established earlier; (ii) the magnitude of thermal vibrations is small, less than 0.1 Å² at 100 K; (iii) the percentage of thermal vibrations in B-factors is the lowest at low resolution and low temperature (<10%) but increases to as high as 60% for structures determined at high resolution and at room temperature. The significance of this work is that it provides for the first time, using an extremely large dataset, a thorough analysis of B-factors and their thermal and static disorder components. The results clearly demonstrate that structures determined at high resolution and at room temperature have the richest dynamics information. Since such structures are relatively rare in the PDB database, the work naturally calls for more such structures to be determined experimentally.

Coupling dynamics and evolutionary information with structure to identify protein regulatory and functional binding sites

Binding sites in proteins can be either specifically functional binding sites (active sites) that bind specific substrates with high affinity or regulatory binding sites (allosteric sites), that modulate the activity of functional binding sites through effector molecules. Owing to their significance in determining protein function, the identification of protein functional and regulatory binding sites is widely acknowledged as an important biological problem. In this work, we present a novel binding site prediction method, Active and Regulatory site Prediction (AR-Pred), which supplements protein geometry, evolutionary, and physicochemical features with information about protein dynamics to predict putative active and allosteric site residues. As the intrinsic dynamics of globular proteins plays an essential role in controlling binding events, we find it to be an important feature for the identification of protein binding sites. We train and validate our predictive models on multiple balanced training and validation sets with random forest machine learning and obtain an ensemble of discrete models for each prediction type. Our models for active site prediction yield a median area under the curve (AUC) of 91% and Matthews correlation coefficient (MCC) of 0.68, whereas the less well-defined allosteric sites are predicted at a lower level with a median AUC of 80% and MCC of 0.48. When tested on an independent set of proteins, our models for active site prediction show comparable performance to two existing methods and gains compared to two others, while the allosteric site models show gains when tested against three existing prediction methods. AR-Pred is available as a free downloadable package at https://github.com/sambitmishra0628/AR-PRED_source.

A deep learning approach to the structural analysis of proteins

Deep learning (DL) algorithms hold great promise for applications in the field of computational biophysics. In fact, the vast amount of available molecular structures, as well as their notable complexity, constitutes an ideal context in which DL-based approaches can be profitably employed. To express the full potential of these techniques, though, it is a prerequisite to express the information contained in a molecule’s atomic positions and distances in a set of input quantities that the network can process. Many of the molecular descriptors devised so far are effective and manageable for relatively small structures, but become complex and cumbersome for larger ones. Furthermore, most of them are defined locally, a feature that could represent a limit for those applications where global properties are of interest. Here, we build a DL architecture capable of predicting non-trivial and intrinsically global quantities, that is, the eigenvalues of a protein’s lowest-energy fluctuation modes. This application represents a first, relatively simple test bed for the development of a neural network approach to the quantitative analysis of protein structures, and demonstrates unexpected use in the identification of mechanically relevant regions of the molecule.

Folding cooperativity and allosteric function in the tandem-repeat protein class

The term allostery was originally developed to describe structural changes in one binding site induced by the interaction of a partner molecule with a distant binding site, and it has been studied in depth in the field of enzymology. Here, we discuss the concept of action at a distance in relation to the folding and function of the solenoid class of tandem-repeat proteins such as tetratricopeptide repeats (TPRs) and ankyrin repeats. Distantly located repeats fold cooperatively, even though only nearest-neighbour interactions exist in these proteins. A number of repeat-protein scaffolds have been reported to display allosteric effects, transferred through the repeat array, that enable them to direct the activity of the multi-subunit enzymes within which they reside. We also highlight a recently identified group of tandem-repeat proteins, the RRPNN subclass of TPRs, recent crystal structures of which indicate that they function as allosteric switches to modulate multiple bacterial quorum-sensing mechanisms. We believe that the folding cooperativity of tandem-repeat proteins and the biophysical mechanisms that transform them into allosteric switches are intimately intertwined. This opinion piece aims to combine our understanding of the two areas and develop ideas on their common underlying principles.This article is part of a discussion meeting issue 'Allostery and molecular machines'.

Shared dynamics of LeuT superfamily members and allosteric differentiation by structural irregularities and multimerization

The LeuT-fold superfamily includes secondary active transporters from different functional families, which share a common tertiary structure, despite having a remarkably low sequence similarity. By identifying the common structural and dynamical features upon principal component analysis of a comprehensive ensemble of 90 experimentally resolved structures and anisotropic network model evaluation of collective motions, we provide a unified point of view for understanding the reasons why this particular fold has been selected by evolution to accomplish such a broad spectrum of functions. The parallel identification of conserved sequence features, localized at specific sites of transmembrane helices, sheds light on the role of broken helices (TM1 and TM6 in LeuT) in promoting ion/substrate binding and allosteric interconversion between the outward- and inward-facing conformations of transporters. Finally, the determination of the dynamics landscape for the structural ensemble provides a promising framework for the classification of transporters based on their dynamics, and the characterization of the collective movements that favour multimerization.This article is part of a discussion meeting issue 'Allostery and molecular machines'.

Allostery in Its Many Disguises: From Theory to Applications

Allosteric regulation plays an important role in many biological processes, such as signal transduction, transcriptional regulation, and metabolism. Allostery is rooted in the fundamental physical properties of macromolecular systems, but its underlying mechanisms are still poorly understood. A collection of contributions to a recent interdisciplinary CECAM (Center Européen de Calcul Atomique et Moléculaire) workshop is used here to provide an overview of the progress and remaining limitations in the understanding of the mechanistic foundations of allostery gained from computational and experimental analyses of real protein systems and model systems. The main conceptual frameworks instrumental in driving the field are discussed. We illustrate the role of these frameworks in illuminating molecular mechanisms and explaining cellular processes, and describe some of their promising practical applications in engineering molecular sensors and informing drug design efforts.

Structural dynamics is a determinant of the functional significance of missense variants

Discrimination of clinically relevant mutations from neutral mutations is of paramount importance in precision medicine and pharmacogenomics. Our study shows that current computational predictions of pathogenicity, mostly based on analysis of sequence conservation, may be improved by considering the changes in the structural dynamics of the protein due to point mutations. We introduce and demonstrate the utility of a classifier that takes advantage of efficient evaluation of structural dynamics by elastic network models.

Accurate evaluation of the effect of point mutations on protein function is essential to assessing the genesis and prognosis of many inherited diseases and cancer types. Currently, a wealth of computational tools has been developed for pathogenicity prediction. Two major types of data are used to this aim: sequence conservation/evolution and structural properties. Here, we demonstrate in a systematic way that another determinant of the functional impact of missense variants is the protein’s structural dynamics. Measurable improvement is shown in pathogenicity prediction by taking into consideration the dynamical context and implications of the mutation. Our study suggests that the class of dynamics descriptors introduced here may be used in conjunction with existing features to not only increase the prediction accuracy of the impact of variants on biological function, but also gain insight into the physical basis of the effect of missense variants.

Dynamic Modulation of Binding Affinity as a Mechanism for Regulating Interferon Signaling

How structural dynamics affects cytokine signaling is under debate. Here, we investigated the dynamics of the type I interferon (IFN) receptor, IFNAR1, and its effect on signaling upon binding IFN and IFNAR2 using a combination of structure-based mechanistic studies, in situ binding, and gene induction assays. Our study reveals that IFNAR1 flexibility modulates ligand-binding affinity, which, in turn, regulates biological signaling. We identified the hinge sites and key interactions implicated in IFNAR1 inter-subdomain (SD1-SD4) movements. We showed that the predicted cooperative movements are essential to accommodate intermolecular interactions. Engineered disulfide bridges, computationally predicted to interfere with IFNAR1 dynamics, were experimentally confirmed. Notably, introducing disulfide bonds between subdomains SD2 and SD3 modulated IFN binding and activity in accordance with the relative attenuation of cooperative movements with varying distance from the hinge center, whereas locking the SD3-SD4 interface flexibility in favor of an extended conformer increased activity.

Altered dynamics upon oligomerization corresponds to key functional sites

It is known that over half of the proteins encoded by most organisms function as oligomeric complexes. Oligomerization confers structural stability and dynamics changes in proteins. We investigate the effects of oligomerization on protein dynamics and its functional significance for a set of 145 multimeric proteins. Using coarse-grained elastic network models, we inspect the changes in residue fluctuations upon oligomerization and then compare with residue conservation scores to identify the functional significance of these changes. Our study reveals conservation of about ½ of the fluctuations, with ¼ of the residues increasing in their mobilities and ¼ having reduced fluctuations. The residues with dampened fluctuations are evolutionarily more conserved and can serve as orthosteric binding sites, indicating their importance. We also use triosephosphate isomerase as a test case to understand why certain enzymes function only in their oligomeric forms despite the monomer including all required catalytic residues. To this end, we compare the residue communities (groups of residues which are highly correlated in their fluctuations) in the monomeric and dimeric forms of the enzyme. We observe significant changes to the dynamical community architecture of the catalytic core of this enzyme. This relates to its functional mechanism and is seen only in the oligomeric form of the protein, answering why proteins are oligomeric structures. Proteins 2017; 85:1422-1434. © 2017 Wiley Periodicals, Inc.

Comparative Normal Mode Analysis of the Dynamics of DENV and ZIKV Capsids

Key steps in the life cycle of a virus, such as the fusion event as the virus infects a host cell and its maturation process, relate to an intricate interplay between the structure and the dynamics of its constituent proteins, especially those that define its capsid, much akin to an envelope that protects its genomic material. We present a comprehensive, comparative analysis of such interplay for the capsids of two viruses from the flaviviridae family, Dengue (DENV) and Zika (ZIKV). We use for that purpose our own software suite, DD-NMA, which is based on normal mode analysis. We describe the elements of DD-NMA that are relevant to the analysis of large systems, such as virus capsids. In particular, we introduce our implementation of simplified elastic networks and justify their parametrization. Using DD-NMA, we illustrate the importance of packing interactions within the virus capsids on the dynamics of the E proteins of DENV and ZIKV. We identify differences between the computed atomic fluctuations of the E proteins in DENV and ZIKV and relate those differences to changes observed in their high resolution structures. We conclude with a discussion on additional analyses that are needed to fully characterize the dynamics of the two viruses.

A multi-resolution model to capture both global fluctuations of an enzyme and molecular recognition in the ligand-binding site

In multi-resolution simulations, different system components are simultaneously modeled at different levels of resolution, these being smoothly coupled together. In the case of enzyme systems, computationally expensive atomistic detail is needed in the active site to capture the chemistry of ligand binding. Global properties of the rest of the protein also play an essential role, determining the structure and fluctuations of the binding site; however, these can be modeled on a coarser level. Similarly, in the most computationally efficient scheme only the solvent hydrating the active site requires atomistic detail. We present a methodology to couple atomistic and coarse-grained protein models, while solvating the atomistic part of the protein in atomistic water. This allows a free choice of which protein and solvent degrees of freedom to include atomistically. This multi-resolution methodology can successfully model stable ligand binding, and we further confirm its validity by exploring the reproduction of system properties relevant to enzymatic function. In addition to a computational speedup, such an approach can allow the identification of the essential degrees of freedom playing a role in a given process, potentially yielding new insights into biomolecular function. Proteins 2016; 84:1902-1913. © 2016 Wiley Periodicals, Inc.

Cooperative Dynamics of Intact AMPA and NMDA Glutamate Receptors: Similarities and Subfamily-Specific Differences

Ionotropic glutamate receptors (iGluRs) are tetrameric ion channels that mediate excitatory neurotransmission. Recent structures of α-amino-3-hydroxy-5-methyl-4-isoxazole propionate (AMPA) and N-methyl-D-aspartate (NMDA) receptors permit a comparative analysis of whole-receptor dynamics for the first time. Despite substantial differences in the packing of their two-domain extracellular region, the two iGluRs share similar dynamics, elucidated by elastic network models. Motions accessible to either structure enable conformational interconversion, such as compression of the AMPA receptor toward the more tightly packed NMDA receptor conformation, which has been linked to allosteric regulation. Pivoting motions coupled to concerted rotations of the transmembrane ion channel are prominent between dimers of distal N-terminal domains in the loosely packed AMPA receptor. The occurrence and functional relevance of these motions is verified by cross-linking experiments designed to probe the computationally predicted distance changes. Together with the identification of hotspot residues acting as mediators of allosteric communication, our data provide a glimpse into the dynamic spectrum of iGluRs.

WEBnm@ v2.0: Web server and services for comparing protein flexibility

BACKGROUND

Normal mode analysis (NMA) using elastic network models is a reliable and cost-effective computational method to characterise protein flexibility and by extension, their dynamics. Further insight into the dynamics-function relationship can be gained by comparing protein motions between protein homologs and functional classifications. This can be achieved by comparing normal modes obtained from sets of evolutionary related proteins.

RESULTS

We have developed an automated tool for comparative NMA of a set of pre-aligned protein structures. The user can submit a sequence alignment in the FASTA format and the corresponding coordinate files in the Protein Data Bank (PDB) format. The computed normalised squared atomic fluctuations and atomic deformation energies of the submitted structures can be easily compared on graphs provided by the web user interface. The web server provides pairwise comparison of the dynamics of all proteins included in the submitted set using two measures: the Root Mean Squared Inner Product and the Bhattacharyya Coefficient. The Comparative Analysis has been implemented on our web server for NMA, WEBnm@, which also provides recently upgraded functionality for NMA of single protein structures. This includes new visualisations of protein motion, visualisation of inter-residue correlations and the analysis of conformational change using the overlap analysis. In addition, programmatic access to WEBnm@ is now available through a SOAP-based web service. Webnm@ is available at http://apps.cbu.uib.no/webnma .

CONCLUSION

WEBnm@ v2.0 is an online tool offering unique capability for comparative NMA on multiple protein structures. Along with a convenient web interface, powerful computing resources, and several methods for mode analyses, WEBnm@ facilitates the assessment of protein flexibility within protein families and superfamilies. These analyses can give a good view of how the structures move and how the flexibility is conserved over the different structures.

Structural dynamics flexibility informs function and evolution at a proteome scale

Protein structures are dynamic entities with a myriad of atomic fluctuations, side-chain rotations, and collective domain movements. Although the importance of these dynamics to proper functioning of proteins is emerging in the studies of many protein families, there is a lack of broad evidence for the critical role of protein dynamics in shaping the biological functions of a substantial fraction of residues for a large number of proteins in the human proteome. Here, we propose a novel dynamic flexibility index (dfi) to quantify the dynamic properties of individual residues in any protein and use it to assess the importance of protein dynamics in 100 human proteins. Our analyses involving functionally critical positions, disease-associated and putatively neutral population variations, and the rate of interspecific substitutions per residue produce concordant patterns at a proteome scale. They establish that the preservation of dynamic properties of residues in a protein structure is critical for maintaining the protein/biological function. Therefore, structural dynamics needs to become a major component of the analysis of protein function and evolution. Such analyses will be facilitated by the dfi, which will also enable the integrative use of structural dynamics with evolutionary conservation in genomic medicine as well as functional genomics investigations.

Conformational Sampling of Maltose-transporter Components in Cartesian Collective Variables is Governed by the Low-frequency Normal Modes

We have studied large-scale conformational transitions in the maltose-binding protein, and the nucleotide binding domains of a maltose-transporter using enhanced conformational sampling in Cartesian collective variables (CVs) with temperature-accelerated molecular dynamics (TAMD), and C(α)-based elastic network normal mode analysis. Significantly, every functional displacement in the TAMD-generated pathways of each protein could be rationalized via a single low-frequency soft mode, while a combination of 2 to 3 low-frequency modes were found to describe the entire conformational change suggesting that collective functional movement in TAMD trajectories is facilitated by the intrinsically accessible low-frequency normal modes. By applying a harmonic potential to facilitate functional motion in TAMD simulations, we also provide a recipe to reproducibly generate structural transitions in both proteins, which can be used to characterize large-scale conformational changes in other biomolecules.

Pre-existing soft modes of motion uniquely defined by native contact topology facilitate ligand binding to proteins

Modeling protein flexibility constitutes a major challenge in accurate prediction of protein-ligand and protein-protein interactions in docking simulations. The lack of a reliable method for predicting the conformational changes relevant to substrate binding prevents the productive application of computational docking to proteins that undergo large structural rearrangements. Here, we examine how coarse-grained normal mode analysis has been advantageously applied to modeling protein flexibility associated with ligand binding. First, we highlight recent studies that have shown that there is a close agreement between the large-scale collective motions of proteins predicted by elastic network models and the structural changes experimentally observed upon ligand binding. Then, we discuss studies that have exploited the predicted soft modes in docking simulations. Two general strategies are noted: pregeneration of conformational ensembles that are then utilized as input for standard fixed-backbone docking and protein structure deformation along normal modes concurrent to docking. These studies show that the structural changes apparently "induced" upon ligand binding occur selectively along the soft modes accessible to the protein prior to ligand binding. They further suggest that proteins offer suitable means of accommodating/facilitating the recognition and binding of their ligand, presumably acquired by evolutionary selection of the suitable three-dimensional structure.

Structural changes involved in protein binding correlate with intrinsic motions of proteins in the unbound state

Protein-protein binding usually involves structural changes that may extend beyond the rearrangements on a local scale, and cannot be explained by a classical lock-and-key mechanism. Several models have been advanced to explain the flexible binding of proteins such as the induced fit mechanism where the ligand is postulated to induce a conformational change at the interaction site upon binding, or the preexisting equilibrium hypothesis that assumes that protein samples an ensemble of conformations at equilibrium conditions and that the ligand binds selectively to an active conformation. We explored the equilibrium motions of proteins that exhibit relatively large (nonlocal) conformational changes upon protein binding using the Gaussian network model and the anisotropic network model of protein dynamics. For four complexes, LIR-1/HLA-A2, Actin/DNase I, CDK2/cyclin, and CDK6/p16(INK4a), the motions calculated for the monomer exhibiting the largest conformational change, in its unbound (free) form, correlate with the experimentally observed structural changes upon binding. This study emphasizes the preexisting equilibrium/conformational selection as a mechanism for protein-protein interaction and lends support the concept that proteins, in their native conformation, are predisposed to undergo conformational fluctuations that are relevant to, or even required for, their biological functions.