Interaction energy based protein structure networks

Computational Modeling of the Hsp90 Interactions with Cochaperones and Small-Molecule Inhibitors

Allosteric interactions of the molecular chaperone Hsp90 with a diverse array of cochaperones and client proteins, such as protein kinases and transcription factors, allow for efficient molecular communication in signal transduction networks. Deregulation of pathways involving these proteins is commonly associated with cancer pathologies and allosteric inhibition of oncogenic clients by targeting Hsp90 provides a powerful therapeutic strategy in cancer research. We review several validated computational approaches and tools used in the studies of the Hsp90 interactions with proteins and small molecules. These methods include experimentally guided docking to predict Hs90-protein interactions, molecular and binding free energy simulations to analyze Hsp90 binding with small molecules, and structure-based network modeling to evaluate allosteric interactions and communications in the Hsp90 regulatory complexes. Through the lens of allosteric-centric view on Hsp90 function and regulation, we discuss newly emerging computational tools that link protein structure modeling with biophysical simulations and network-based systems biology approaches.

Uncovering GPCR and G Protein Function by Protein Structure Network Analysis

Protein structure network (PSN) analysis is one of the graph theory-based approaches currently used for investigating structural communication in biomolecular systems. Information on the system's dynamics can be provided by atomistic molecular dynamics (MD) simulations or coarse grained elastic network models paired with normal mode analysis (ENM-NMA). This chapter reports on selected applications of PSN analysis to uncover the structural communication in G protein coupled receptors (GPCRs) and G proteins. Strategies to highlight changes in structural communication caused by mutations, ligand and protein binding are described. Conserved amino acids, sites of misfolding mutations, or ligands acting as functional switches tend to behave as hubs in the native structure networks. Densely linked regions in the protein structure graphs could be identified as playing central roles in protein stability and function. Changes in the communication pathway fingerprints depending on the bound ligand or following amino acid mutation could be highlighted as well. A bridge between misfolding and misrouting could be established in rhodopsin mutants linked to inherited blindness. The analysis of native network perturbations by misfolding mutations served to infer key structural elements of protein responsiveness to small chaperones with implications for drug discovery.

Ensemble-based modeling and rigidity decomposition of allosteric interaction networks and communication pathways in cyclin-dependent kinases: Differentiating kinase clients of the Hsp90-Cdc37 chaperone

The overarching goal of delineating molecular principles underlying differentiation of protein kinase clients and chaperone-based modulation of kinase activity is fundamental to understanding activity of many oncogenic kinases that require chaperoning of Hsp70 and Hsp90 systems to attain a functionally competent active form. Despite structural similarities and common activation mechanisms shared by cyclin-dependent kinase (CDK) proteins, members of this family can exhibit vastly different chaperone preferences. The molecular determinants underlying chaperone dependencies of protein kinases are not fully understood as structurally similar kinases may often elicit distinct regulatory responses to the chaperone. The regulatory divergences observed for members of CDK family are of particular interest as functional diversification among these kinases may be related to variations in chaperone dependencies and can be exploited in drug discovery of personalized therapeutic agents. In this work, we report the results of a computational investigation of several members of CDK family (CDK5, CDK6, CDK9) that represented a broad repertoire of chaperone dependencies—from nonclient CDK5, to weak client CDK6, and strong client CDK9. By using molecular simulations of multiple crystal structures we characterized conformational ensembles and collective dynamics of CDK proteins. We found that the elevated dynamics of CDK9 can trigger imbalances in cooperative collective motions and reduce stability of the active fold, thus creating a cascade of favorable conditions for chaperone intervention. The ensemble-based modeling of residue interaction networks and community analysis determined how differences in modularity of allosteric networks and topography of communication pathways can be linked with the client status of CDK proteins. This analysis unveiled depleted modularity of the allosteric network in CDK9 that alters distribution of communication pathways and leads to impaired signaling in the client kinase. According to our results, these network features may uniquely define chaperone dependencies of CDK clients. The perturbation response scanning and rigidity decomposition approaches identified regulatory hotspots that mediate differences in stability and cooperativity of allosteric interaction networks in the CDK structures. By combining these synergistic approaches, our study revealed dynamic and network signatures that can differentiate kinase clients and rationalize subtle divergences in the activation mechanisms of CDK family members. The therapeutic implications of these results are illustrated by identifying structural hotspots of pathogenic mutations that preferentially target regions of the increased flexibility to enable modulation of activation changes. Our study offers a network-based perspective on dynamic kinase mechanisms and drug design by unravelling relationships between protein kinase dynamics, allosteric communications and chaperone dependencies.

Hidden electrostatic basis of dynamic allostery in a PDZ domain

Allosteric effect implies ligand binding at one site leading to structural and/or dynamical changes at a distant site. PDZ domains are classic examples of dynamic allostery without conformational changes, where distal side-chain dynamics is modulated on ligand binding and the origin has been attributed to entropic effects. In this work, we unearth the energetic basis of the observed dynamic allostery in a PDZ3 domain protein using molecular dynamics simulations. We demonstrate that electrostatic interaction provides a highly sensitive yardstick to probe the allosteric modulation in contrast to the traditionally used structure-based parameters. There is a significant population shift in the hydrogen-bonded network and salt bridges involving side chains on ligand binding. The ligand creates a local energetic perturbation that propagates in the form of dominolike changes in interresidue interaction pattern. There are significant changes in the nature of specific interactions (nonpolar/polar) between interresidue contacts and accompanied side-chain reorientations that drive the major redistribution of energy. Interestingly, this internal redistribution and rewiring of side-chain interactions led to large cancellations resulting in small change in the overall enthalpy of the protein, thus making it difficult to detect experimentally. In contrast to the prevailing focus on the entropic or dynamic effects, we show that the internal redistribution and population shift in specific electrostatic interactions drive the allosteric modulation in the PDZ3 domain protein.

The adaptive nature of protein residue networks

Protein residue networks PRNs are used to describe proteins. These networks are usually based on an average structure for the protein. However, proteins are dynamic entities that are affected by their surroundings. In this work, we study the effect of temperatures above and below the protein dynamical transition temperature(≈200 K), on three important network parameters gleaned from weighted PRNs for the solvated β-lactamase inhibitory protein BLIP: the betweenness centrality B, the closeness centrality C, and the clustering coefficient CC. The B and C values will be extracted for each node from PRNs at six different temperatures: 150 K, 180 K, 200 K, 220 K, 250 K, and 310 K respectively. The average value for the CC for each PRN will also be calculated at each temperature, respectively. We find that at temperatures ≤200 K, the network nodes with the most significant B and C values tend to have lower relative solvent accessibility RSA values, and to fall within the protein secondary structure elements (α helices and β sheets). At temperatures >200 K, the significant nodes in terms of B and C tend to have larger RSA values, and to fall on the connecting loops in the protein. The average CC decreases in value for the PRNs up to 200 K, and then remains basically constant above 200 K. This clearly shows that any conclusions based on static PRNs should be handled with care. The dynamic nature of proteins and its coupling to the surrounding environment should be taken into consideration when using the PRN paradigm. Proteins 2017; 85:917-923. © 2016 Wiley Periodicals, Inc.

Computational Analysis of Residue Interaction Networks and Coevolutionary Relationships in the Hsp70 Chaperones: A Community-Hopping Model of Allosteric Regulation and Communication

Allosteric interactions in the Hsp70 proteins are linked with their regulatory mechanisms and cellular functions. Despite significant progress in structural and functional characterization of the Hsp70 proteins fundamental questions concerning modularity of the allosteric interaction networks and hierarchy of signaling pathways in the Hsp70 chaperones remained largely unexplored and poorly understood. In this work, we proposed an integrated computational strategy that combined atomistic and coarse-grained simulations with coevolutionary analysis and network modeling of the residue interactions. A novel aspect of this work is the incorporation of dynamic residue correlations and coevolutionary residue dependencies in the construction of allosteric interaction networks and signaling pathways. We found that functional sites involved in allosteric regulation of Hsp70 may be characterized by structural stability, proximity to global hinge centers and local structural environment that is enriched by highly coevolving flexible residues. These specific characteristics may be necessary for regulation of allosteric structural transitions and could distinguish regulatory sites from nonfunctional conserved residues. The observed confluence of dynamics correlations and coevolutionary residue couplings with global networking features may determine modular organization of allosteric interactions and dictate localization of key mediating sites. Community analysis of the residue interaction networks revealed that concerted rearrangements of local interacting modules at the inter-domain interface may be responsible for global structural changes and a population shift in the DnaK chaperone. The inter-domain communities in the Hsp70 structures harbor the majority of regulatory residues involved in allosteric signaling, suggesting that these sites could be integral to the network organization and coordination of structural changes. Using a network-based formalism of allostery, we introduced a community-hopping model of allosteric communication. Atomistic reconstruction of signaling pathways in the DnaK structures captured a direction-specific mechanism and molecular details of signal transmission that are fully consistent with the mutagenesis experiments. The results of our study reconciled structural and functional experiments from a network-centric perspective by showing that global properties of the residue interaction networks and coevolutionary signatures may be linked with specificity and diversity of allosteric regulation mechanisms.

The diversity of allosteric mechanisms in the Hsp70 proteins could range from modulation of the inter-domain interactions and conformational dynamics to fine-tuning of the Hsp70 interactions with co-chaperones. The goal of this study is to present a systematic computational analysis of the dynamic and evolutionary factors underlying allosteric structural transformations of the Hsp70 proteins. We investigated the relationship between functional dynamics, residue coevolution, and network organization of residue interactions in the Hsp70 proteins. The results of this study revealed that conformational dynamics of the Hsp70 proteins may be linked with coevolutionary propensities and mutual information dependencies of the protein residues. Modularity and connectivity of allosteric interactions in the Hsp70 chaperones are coordinated by stable functional sites that feature unique coevolutionary signatures and high network centrality. The emergence of the inter-domain communities that are coordinated by functional centers and include highly coevolving residues could facilitate structural transitions through cooperative reorganization of the local interacting modules. We determined that the differences in the modularity of the residue interactions and organization of coevolutionary networks in DnaK may be associated with variations in their allosteric mechanisms. The network signatures of the DnaK structures are characteristic of a population-shift allostery that allows for coordinated structural rearrangements of local communities. A dislocation of mediating centers and insufficient coevolutionary coupling between functional regions may render a reduced cooperativity and promote a limited entropy-driven allostery in the Sse1 chaperone that occurs without structural changes. The results of this study showed that a network-centric framework and a community-hopping model of allosteric communication pathways may provide novel insights into molecular and evolutionary principles of allosteric regulation in the Hsp70 proteins.

Network-based modelling and percolation analysis of conformational dynamics and activation in the CDK2 and CDK4 proteins: dynamic and energetic polarization of the kinase lobes may determine divergence of the regulatory mechanisms

Network modeling and percolation analysis of conformational dynamics and energetics of regulatory mechanisms in cyclin-dependent kinases.

Exploring Molecular Mechanisms of Paradoxical Activation in the BRAF Kinase Dimers: Atomistic Simulations of Conformational Dynamics and Modeling of Allosteric Communication Networks and Signaling Pathways

The recent studies have revealed that most BRAF inhibitors can paradoxically induce kinase activation by promoting dimerization and enzyme transactivation. Despite rapidly growing number of structural and functional studies about the BRAF dimer complexes, the molecular basis of paradoxical activation phenomenon is poorly understood and remains largely hypothetical. In this work, we have explored the relationships between inhibitor binding, protein dynamics and allosteric signaling in the BRAF dimers using a network-centric approach. Using this theoretical framework, we have combined molecular dynamics simulations with coevolutionary analysis and modeling of the residue interaction networks to determine molecular determinants of paradoxical activation. We have investigated functional effects produced by paradox inducer inhibitors PLX4720, Dabrafenib, Vemurafenib and a paradox breaker inhibitor PLX7904. Functional dynamics and binding free energy analyses of the BRAF dimer complexes have suggested that negative cooperativity effect and dimer-promoting potential of the inhibitors could be important drivers of paradoxical activation. We have introduced a protein structure network model in which coevolutionary residue dependencies and dynamic maps of residue correlations are integrated in the construction and analysis of the residue interaction networks. The results have shown that coevolutionary residues in the BRAF structures could assemble into independent structural modules and form a global interaction network that may promote dimerization. We have also found that BRAF inhibitors could modulate centrality and communication propensities of global mediating centers in the residue interaction networks. By simulating allosteric communication pathways in the BRAF structures, we have determined that paradox inducer and breaker inhibitors may activate specific signaling routes that correlate with the extent of paradoxical activation. While paradox inducer inhibitors may facilitate a rapid and efficient communication via an optimal single pathway, the paradox breaker may induce a broader ensemble of suboptimal and less efficient communication routes. The central finding of our study is that paradox breaker PLX7904 could mimic structural, dynamic and network features of the inactive BRAF-WT monomer that may be required for evading paradoxical activation. The results of this study rationalize the existing structure-functional experiments by offering a network-centric rationale of the paradoxical activation phenomenon. We argue that BRAF inhibitors that amplify dynamic features of the inactive BRAF-WT monomer and intervene with the allosteric interaction networks may serve as effective paradox breakers in cellular environment.

Mechanisms for allosteric activation of protease DegS by ligand binding and oligomerization as revealed from molecular dynamics simulations

A local perturbation of a protein may lead to functional changes at some distal site, a phenomenon denoted as allostery. Here, we study the allosteric control of a protease using molecular dynamics simulations. The system considered is the bacterial protein DegS which includes a protease domain activated on ligand binding to an adjacent PDZ domain. Starting from crystallographic structures of DegS homo-trimers, we perform simulations of the ligand-free and -bound state of DegS at equilibrium. Considering a single protomer only, the trimeric state was mimicked by applying restraints on the residues in contact with other protomers in the DegS trimer. In addition, the bound state was also simulated without any restraints to mimic the monomer. Our results suggest that not only ligand release but also disassembly of a DegS trimer inhibits proteolytic activity. Considering various observables for structural changes, we infer allosteric pathways from the interface with other protomers to the active site. Moreover, we study how ligand release leads to (i) catalytically relevant changes involving residues 199-201 and (ii) a transition from a stretched to a bent conformation for residues 217-219 (which prohibits proper substrate binding). Finally, based on ligand-induced C_α shifts we identify residues in contact with other protomers in the DegS trimer that likely transduce the perturbation from ligand release from a given protomer to adjacent protomers. These residues likely play a key role in the experimentally known effect of ligand release from a protomer on the proteolytic activity of the other protomers. Proteins 2016; 84:1690-1705. © 2016 Wiley Periodicals, Inc.

Rigid Residue Scan Simulations Systematically Reveal Residue Entropic Roles in Protein Allostery

Intra-protein information is transmitted over distances via allosteric processes. This ubiquitous protein process allows for protein function changes due to ligand binding events. Understanding protein allostery is essential to understanding protein functions. In this study, allostery in the second PDZ domain (PDZ2) in the human PTP1E protein is examined as model system to advance a recently developed rigid residue scan method combining with configurational entropy calculation and principal component analysis. The contributions from individual residues to whole-protein dynamics and allostery were systematically assessed via rigid body simulations of both unbound and ligand-bound states of the protein. The entropic contributions of individual residues to whole-protein dynamics were evaluated based on covariance-based correlation analysis of all simulations. The changes of overall protein entropy when individual residues being held rigid support that the rigidity/flexibility equilibrium in protein structure is governed by the La Châtelier's principle of chemical equilibrium. Key residues of PDZ2 allostery were identified with good agreement with NMR studies of the same protein bound to the same peptide. On the other hand, the change of entropic contribution from each residue upon perturbation revealed intrinsic differences among all the residues. The quasi-harmonic and principal component analyses of simulations without rigid residue perturbation showed a coherent allosteric mode from unbound and bound states, respectively. The projection of simulations with rigid residue perturbation onto coherent allosteric modes demonstrated the intrinsic shifting of ensemble distributions supporting the population-shift theory of protein allostery. Overall, the study presented here provides a robust and systematic approach to estimate the contribution of individual residue internal motion to overall protein dynamics and allostery.

Dynamical Network of HIV-1 Protease Mutants Reveals the Mechanism of Drug Resistance and Unhindered Activity

HIV-1 protease variants resist drugs by active and non-active-site mutations. The active-site mutations, which are the primary or first set of mutations, hamper the stability of the enzyme and resist the drugs minimally. As a result, secondary mutations that not only increase protein stability for unhindered catalytic activity but also resist drugs very effectively arise. While the mechanism of drug resistance of the active-site mutations is through modulating the active-site pocket volume, the mechanism of drug resistance of the non-active-site mutations is unclear. Moreover, how these allosteric mutations, which are 8-21 Å distant, communicate to the active site for drug efflux is completely unexplored. Results from molecular dynamics simulations suggest that the primary mechanism of drug resistance of the secondary mutations involves opening of the flexible protease flaps. Results from both residue- and community-based network analyses reveal that this precise action of protease is accomplished by the presence of robust communication paths between the mutational sites and the functionally relevant regions: active site and flaps. While the communication is more direct in the wild type, it traverses across multiple intermediate residues in mutants, leading to weak signaling and unregulated motions of flaps. The global integrity of the protease network is, however, maintained through the neighboring residues, which exhibit high degrees of conservation, consistent with clinical data and mutagenesis studies.

Long-Range Conformational Response of a PDZ Domain to Ligand Binding and Release: A Molecular Dynamics Study

The binding of a ligand to a protein may induce long-range structural or dynamical changes in the biomacromolecule even at sites physically well separated from the binding pocket. A system for which such behavior has been widely discussed is the PDZ2 domain of human tyrosine phosphatase 1E. Here, we present results from equilibrium trajectories of the PDZ2 domain in the free and ligand-bound state, as well as nonequilibrium simulations of the relaxation of PDZ2 after removal of its peptide ligand. The study reveals changes in inter-residue contacts, backbone dihedral angles, and C(α) positions upon ligand release. Our findings show a long-range conformational response of the PDZ2 domain to ligand release in the form of a collective shift of the secondary structure elements α2, β2, β3, α1-β4, and the C terminal loop relative to the rest of the protein away from the N-terminus, and a shift of the loops β2-β3 and β1-β2 in the opposite direction. The shifts lead to conformational changes in the backbone, especially in the β2-β3 loop but also in the β5-α2 and the α2-β6 loop, and are accompanied by changes of inter-residue contacts mainly within the β2-β3 loop as well as between the α2 helix and other segments. The residues showing substantial changes of inter-residue contacts, backbone conformations, or C(α) positions are considered "key residues" for the long-range conformational response of PDZ2. By comparing these residues with various sets of residues highlighted by previous studies of PDZ2, we investigate the statistical correlation of the various approaches. Interestingly, we find a considerable correlation of our findings with several works considering structural changes but no significant correlations with approaches considering energy flow or networks based on inter-residue energies.

Use of network model to explore dynamic and allosteric properties of three GPCR homodimers

We used an elastic network model and protein structure network to study three class A GPCR homodimers.

Molecular dynamics simulations and modelling of the residue interaction networks in the BRAF kinase complexes with small molecule inhibitors: probing the allosteric effects of ligand-induced kinase dimerization and paradoxical activation

The computational analysis of residue interaction networks dissects the allosteric effects of inhibitor-induced BRAF kinase dimerization and paradoxical activation.

Mechanisms of Differential Allosteric Modulation in Homologous Proteins: Insights from the Analysis of Internal Dynamics and Energetics of PDZ Domains

Allostery is a general phenomenon in proteins whereby a perturbation at one site reverberates into a functional change at another one, through modulation of its conformational dynamics. Herein, we address the problem of how the molecular signal encoded by a ligand is differentially transmitted through the structures of two homologous PDZ proteins: PDZ2, which responds to binding with structural and dynamical changes in regions distal from the ligand site, and PDZ3, which is characterized by less-intense dynamical variations. We use novel methods of analysis of MD simulations in the unbound and bound states to investigate the determinants of the differential allosteric behavior of the two proteins. The analysis of the correlations between the redistribution of stabilization energy and local fluctuation patterns highlights the nucleus of residues responsible for the stabilization of the 3D fold, the stability core, as the substructure that defines the difference in the allosteric response: in PDZ2, it undergoes a consistent dynamic and energetic reorganization, whereas in PDZ3, it remains largely unperturbed. Specifically, we observe for PDZ2 a significant anticorrelation between the motions of distal loops and residues of the stability core and differences in the correlation patterns between the bound and unbound states. Such variation is not observed in PDZ3, indicating that its energetics and internal dynamics are less affected by the presence/absence of the ligand. Finally, we propose a model with a direct link between the modulation of the structural, energetic and dynamic properties of a protein, and its allosteric response to a perturbation.

Dancing through Life: Molecular Dynamics Simulations and Network-Centric Modeling of Allosteric Mechanisms in Hsp70 and Hsp110 Chaperone Proteins

Hsp70 and Hsp110 chaperones play an important role in regulating cellular processes that involve protein folding and stabilization, which are essential for the integrity of signaling networks. Although many aspects of allosteric regulatory mechanisms in Hsp70 and Hsp110 chaperones have been extensively studied and significantly advanced in recent experimental studies, the atomistic picture of signal propagation and energetics of dynamics-based communication still remain unresolved. In this work, we have combined molecular dynamics simulations and protein stability analysis of the chaperone structures with the network modeling of residue interaction networks to characterize molecular determinants of allosteric mechanisms. We have shown that allosteric mechanisms of Hsp70 and Hsp110 chaperones may be primarily determined by nucleotide-induced redistribution of local conformational ensembles in the inter-domain regions and the substrate binding domain. Conformational dynamics and energetics of the peptide substrate binding with the Hsp70 structures has been analyzed using free energy calculations, revealing allosteric hotspots that control negative cooperativity between regulatory sites. The results have indicated that cooperative interactions may promote a population-shift mechanism in Hsp70, in which functional residues are organized in a broad and robust allosteric network that can link the nucleotide-binding site and the substrate-binding regions. A smaller allosteric network in Hsp110 structures may elicit an entropy-driven allostery that occurs in the absence of global structural changes. We have found that global mediating residues with high network centrality may be organized in stable local communities that are indispensable for structural stability and efficient allosteric communications. The network-centric analysis of allosteric interactions has also established that centrality of functional residues could correlate with their sensitivity to mutations across diverse chaperone functions. This study reconciles a wide spectrum of structural and functional experiments by demonstrating how integration of molecular simulations and network-centric modeling may explain thermodynamic and mechanistic aspects of allosteric regulation in chaperones.

Analysis of core-periphery organization in protein contact networks reveals groups of structurally and functionally critical residues

The representation of proteins as networks of interacting amino acids, referred to as protein contact networks (PCN), and their subsequent analyses using graph theoretic tools, can provide novel insights into the key functional roles of specific groups of residues. We have characterized the networks corresponding to the native states of 66 proteins (belonging to different families) in terms of their core-periphery organization. The resulting hierarchical classification of the amino acid constituents of a protein arranges the residues into successive layers - having higher core order - with increasing connection density, ranging from a sparsely linked periphery to a densely intra-connected core (distinct from the earlier concept of protein core defined in terms of the three-dimensional geometry of the native state, which has least solvent accessibility). Our results show that residues in the inner cores are more conserved than those at the periphery. Underlining the functional importance of the network core, we see that the receptor sites for known ligand molecules of most proteins occur in the innermost core. Furthermore, the association of residues with structural pockets and cavities in binding or active sites increases with the core order. From mutation sensitivity analysis, we show that the probability of deleterious or intolerant mutations also increases with the core order. We also show that stabilization centre residues are in the innermost cores, suggesting that the network core is critically important in maintaining the structural stability of the protein. A publicly available Web resource for performing core-periphery analysis of any protein whose native state is known has been made available by us at http://www.imsc.res.in/ ~sitabhra/proteinKcore/index.html.

Dynamics-Driven Allostery in Protein Kinases

Protein kinases have very dynamic structures and their functionality strongly depends on their dynamic state. Active kinases reveal a dynamic pattern with residues clustering into semirigid communities that move in μs-ms timescale. Previously detected hydrophobic spines serve as connectors between communities. Communities do not follow the traditional subdomain structure of the kinase core or its secondary structure elements. Instead they are organized around main functional units. Integration of the communities depends on the assembly of the hydrophobic spine and phosphorylation of the activation loop. Single mutations can significantly disrupt the dynamic infrastructure and thereby interfere with long-distance allosteric signaling that propagates throughout the whole molecule. Dynamics is proposed to be the underlying mechanism for allosteric regulation in protein kinases.

Characteristics Analyses and Comparisons of the Protein Structure Networks Constructed by Different Methods

Protein structure networks (PSNs) were widely used in analyses of protein structure and function. In this work, we analyzed and compared the characters of PSNs by different methods. The degrees of the different types of the nodes were found to be associated with the amino acid characters, including SAS, secondary structure, hydropathy and the volume of amino acids. It showed that PSNs by the methods of CA10, SC10 and AT5 inherited more amino acid characters and had higher correlations with the original protein structures. And PSNs by these three methods would be powerful tools in understanding the characters of protein structures.

Weighted protein residue networks based on joint recurrences between residues

BACKGROUND

Weighted and un-weighted protein residue networks can predict key functional residues in proteins based on the closeness centrality C and betweenness centrality B values for each residue. A static snapshot of the protein structure, and a cutoff distance, are used to define edges between the network nodes. In this work we apply the weighted network approach to study the β-Lactamase Inhibitory Protein (BLIP). Joint recurrences extracted from molecular dynamics MD trajectory positions of the protein residue carbon alpha atoms are used to define edge weights between nodes, and no cutoff distance is used. The results for B and C from our approach are compared with those extracted from an un-weighted network, and a weighted network that uses interatomic contacts to define edge weights between nodes, respectively.

RESULTS

The joint recurrence weighted network approach performs well in pointing out key protein residues. Furthermore, it seems to emphasize residues with medium to high relative solvent accessibility that lie in loop regions between secondary structure elements of the protein.

CONCLUSIONS

Protein residue networks that use joint recurrences extracted from molecular dynamics simulations of a solvated protein perform well in pointing to hotspot residues and hotspot clusters. This approach uses no distance cutoff threshold, and does not exclude any interactions between the residues, including water-mediated interactions.

Communities in the iron superoxide dismutase amino acid network

Amino acid networks (AANs) analysis is a new way to reveal the relationship between protein structure and function. We constructed six different types of AANs based on iron superoxide dismutase (Fe-SOD) three-dimensional structure information. These Fe-SOD AANs have clear community structures when they were modularized by different methods. Especially, detected communities are related to Fe-SOD secondary structures. Regular structures show better correlations with detected communities than irregular structures, and loops weaken these correlations, which suggest that secondary structure is the unit element in Fe-SOD folding process. In addition, a comparative analysis of mesophilic and thermophilic Fe-SOD AANs' communities revealed that thermostable Fe-SOD AANs had more highly associated community structures than mesophilic one. Thermophilic Fe-SOD AANs also had more high similarity between communities and secondary structures than mesophilic Fe-SOD AANs. The communities in Fe-SOD AANs show that dense interactions in modules can help to stabilize thermophilic Fe-SOD.

Identifying key residues for protein allostery through rigid residue scan

Allostery is a ubiquitous process for protein regulatory activity in which a binding event can change a protein's function carried out at a distal site. Despite intensive theoretical and experimental investigation of protein allostery in the past five decades, effective methods have yet to be developed that can systematically identify key residues involved in allosteric mechanisms. In this study, we propose the rigid residue scan as a systematic approach to identify important allosteric residues. The third PDZ domain (PDZ3) in the postsynaptic density 95 protein (PSD-95) is used as a model system, and each amino acid residue is treated as a single rigid body during independent molecular dynamics simulations. Various indices based on cross-correlation matrices are used, which allow for two groups of residues with different functions to be identified. The first group is proposed as "switches" that are needed to "turn on" the binding effect of protein allostery. The second group is proposed as "wire residues" that are needed to propagate energy or information from the binding site to distal locations within the same protein. Among the nine residues suggested as important for PDZ3 intramolecular communication in this study, eight have been reported as critical for allostery in PDZ3. Therefore, the rigid residue scan approach is demonstrated to be an effective method for systemically identifying key residues in protein intramolecular communication and allosteric mechanisms.

Structure-based network analysis of activation mechanisms in the ErbB family of receptor tyrosine kinases: the regulatory spine residues are global mediators of structural stability and allosteric interactions

The ErbB protein tyrosine kinases are among the most important cell signaling families and mutation-induced modulation of their activity is associated with diverse functions in biological networks and human disease. We have combined molecular dynamics simulations of the ErbB kinases with the protein structure network modeling to characterize the reorganization of the residue interaction networks during conformational equilibrium changes in the normal and oncogenic forms. Structural stability and network analyses have identified local communities integrated around high centrality sites that correspond to the regulatory spine residues. This analysis has provided a quantitative insight to the mechanism of mutation-induced “superacceptor” activity in oncogenic EGFR dimers. We have found that kinase activation may be determined by allosteric interactions between modules of structurally stable residues that synchronize the dynamics in the nucleotide binding site and the αC-helix with the collective motions of the integrating αF-helix and the substrate binding site. The results of this study have pointed to a central role of the conserved His-Arg-Asp (HRD) motif in the catalytic loop and the Asp-Phe-Gly (DFG) motif as key mediators of structural stability and allosteric communications in the ErbB kinases. We have determined that residues that are indispensable for kinase regulation and catalysis often corresponded to the high centrality nodes within the protein structure network and could be distinguished by their unique network signatures. The optimal communication pathways are also controlled by these nodes and may ensure efficient allosteric signaling in the functional kinase state. Structure-based network analysis has quantified subtle effects of ATP binding on conformational dynamics and stability of the EGFR structures. Consistent with the NMR studies, we have found that nucleotide-induced modulation of the residue interaction networks is not limited to the ATP site, and may enhance allosteric cooperativity with the substrate binding region by increasing communication capabilities of mediating residues.

Long-range conformational transition of a photoswitchable allosteric protein: molecular dynamics simulation study

A local perturbation of a protein may lead to functional changes at some distal site. An example is the PDZ2 domain of human tyrosine phosphatase 1E, which shows an allosteric transition upon binding to a peptide ligand. Recently Buchli et al. presented a time-resolved study of this transition by covalently linking an azobenzene photoswitch across the binding groove and using a femtosecond laser pulse that triggers the cis-trans photoisomerization of azobenzene. To aid the interpretation of these experiments, in this work seven microsecond runs of all-atom molecular dynamics simulations each for the wild-type PDZ2 in the ligand-bound and -free state, as well as the photoswitchable protein (PDZ2S) in the cis and trans states of the photoswitch, in explicit water were conducted. First the theoretical model is validated by recalculating the available NMR data from the simulations. By comparing the results for PDZ2 and PDZ2S, it is analyzed to what extent the photoswitch indeed mimics the free-bound transition. A detailed description of the conformational rearrangement following the cis-trans photoisomerization of PDZ2S reveals a series of photoinduced structural changes that propagate from the anchor residues of the photoswitch via intermediate secondary structure segments to the C-terminus of PDZ2S. The changes of the conformational distribution of the C-terminal region is considered as the distal response of the isolated allosteric protein.

Computational modeling of allosteric regulation in the hsp90 chaperones: a statistical ensemble analysis of protein structure networks and allosteric communications

A fundamental role of the Hsp90 chaperone in regulating functional activity of diverse protein clients is essential for the integrity of signaling networks. In this work we have combined biophysical simulations of the Hsp90 crystal structures with the protein structure network analysis to characterize the statistical ensemble of allosteric interaction networks and communication pathways in the Hsp90 chaperones. We have found that principal structurally stable communities could be preserved during dynamic changes in the conformational ensemble. The dominant contribution of the inter-domain rigidity to the interaction networks has emerged as a common factor responsible for the thermodynamic stability of the active chaperone form during the ATPase cycle. Structural stability analysis using force constant profiling of the inter-residue fluctuation distances has identified a network of conserved structurally rigid residues that could serve as global mediating sites of allosteric communication. Mapping of the conformational landscape with the network centrality parameters has demonstrated that stable communities and mediating residues may act concertedly with the shifts in the conformational equilibrium and could describe the majority of functionally significant chaperone residues. The network analysis has revealed a relationship between structural stability, global centrality and functional significance of hotspot residues involved in chaperone regulation. We have found that allosteric interactions in the Hsp90 chaperone may be mediated by modules of structurally stable residues that display high betweenness in the global interaction network. The results of this study have suggested that allosteric interactions in the Hsp90 chaperone may operate via a mechanism that combines rapid and efficient communication by a single optimal pathway of structurally rigid residues and more robust signal transmission using an ensemble of suboptimal multiple communication routes. This may be a universal requirement encoded in protein structures to balance the inherent tension between resilience and efficiency of the residue interaction networks.

Functional versatility and structural adaptability of the Hsp90 chaperones are regulated by allosteric interactions that allow for diverse functions including modulation of ATP hydrolysis and binding with cochaperones and client proteins. By integrating molecular simulations and network-based approaches we have characterized conformational dynamics and allosteric interactions in different functional forms of Hsp90. The network centrality analysis and structural mapping of allosteric communications have revealed a small-world organization of the interaction network that is mediated by functionally important residues of the Hsp90 activity. We have found that effective allosteric communications in the Hsp90 chaperone may be provided by structurally stable residues that exhibit high centrality properties. Nucleotide-specific rewiring of the network topology and assortative organization of functional residues may protect the active form of the chaperone from random perturbations and detrimental mutations. These results have confirmed that allosteric interactions in the Hsp90 chaperone may be determined by a small-world network of functional residues that can regulate the network efficiency and resiliency by modulating the statistical ensemble of communication pathways in response to functional requirements of the ATPase cycle.

PyInteraph: a framework for the analysis of interaction networks in structural ensembles of proteins

In the last years, a growing interest has been gathering around the ability of Molecular Dynamics (MD) to provide insight into the paths of long-range structural communication in biomolecules. The knowledge of the mechanisms related to structural communication helps in the rationalization in atomistic details of the effects induced by mutations, ligand binding, and the intrinsic dynamics of proteins. We here present PyInteraph, a tool for the analysis of structural ensembles inspired by graph theory. PyInteraph is a software suite designed to analyze MD and structural ensembles with attention to binary interactions between residues, such as hydrogen bonds, salt bridges, and hydrophobic interactions. PyInteraph also allows the different classes of intra- and intermolecular interactions to be represented, combined or alone, in the form of interaction graphs, along with performing network analysis on the resulting interaction graphs. The program also integrates the network description with a knowledge-based force field to estimate the interaction energies between side chains in the protein. It can be used alone or together with the recently developed xPyder PyMOL plugin through an xPyder-compatible format. The software capabilities and associated protocols are here illustrated by biologically relevant cases of study. The program is available free of charge as Open Source software via the GPL v3 license at http://linux.btbs.unimib.it/pyinteraph/.

Amino acid contact energy networks impact protein structure and evolution

One of the most challenging tasks in structural proteomics is to understand the relationship between protein structure, biological function, and evolution. An understanding of amino acid networks based on protein topology has an important role in the study of this relationship; however, the relationship between network parameters underlying protein topology with structural properties or evolutionary rate is still unknown. To investigate this further, we modeled the three dimensional structure of proteins as amino acid contact energy networks (AACENs) with nodes represented as amino acid residues and edges established according to environment-dependent residue-residue contact energies. Five other types of networks were also constructed to investigate their topological parameters and compare their effect on protein structure and evolution: (1) a random contact network (RCN), (2) a rewiring network with the same degree of distribution as AACEN (RNDD), (3) long-range contact energy networks with and without the backbone connectivity (LCEN_BBs and LCENs), and (4) short range contact energy networks (SCENs). The results indicated that the long-range link percentage and the network clustering coefficient showed a significantly positive and negative correlation, respectively, with protein secondary structure density. In addition, the long-range link percentage and network diameter had a significantly positive and negative correlation, respectively, with evolutionary rate. According to our knowledge, this is the first study to identify the potential role of long-range links and network diameter in protein evolution.

The construction of an amino acid network for understanding protein structure and function

Amino acid networks (AANs) are undirected networks consisting of amino acid residues and their interactions in three-dimensional protein structures. The analysis of AANs provides novel insight into protein science, and several common amino acid network properties have revealed diverse classes of proteins. In this review, we first summarize methods for the construction and characterization of AANs. We then compare software tools for the construction and analysis of AANs. Finally, we review the application of AANs for understanding protein structure and function, including the identification of functional residues, the prediction of protein folding, analyzing protein stability and protein-protein interactions, and for understanding communication within and between proteins.

Energetic pathway sampling in a protein interaction domain

The affinity and specificity of protein-ligand interactions are influenced by energetic crosstalk within the protein domain. However, the molecular details of such intradomain allostery are still unclear. Here, we have experimentally detected and computationally predicted interaction pathways in the postsynaptic density 95/discs large/zonula occludens 1 (PDZ)-peptide ligand model system using wild-type and circularly permuted PDZ proteins. The circular permutant introduced small perturbations in the tertiary structure and a concomitant rewiring of allosteric pathways, allowing us to describe how subtle changes may reshape energetic signaling. The results were analyzed in the context of other members of the PDZ family, which were found to contain distinct interaction pathways for different peptide ligands. The data reveal a fascinating scenario whereby several energetic pathways are sampled within one single domain and distinct pathways are activated by specific protein ligands.

Allosteric regulation of the Hsp90 dynamics and stability by client recruiter cochaperones: protein structure network modeling

The fundamental role of the Hsp90 chaperone in supporting functional activity of diverse protein clients is anchored by specific cochaperones. A family of immune sensing client proteins is delivered to the Hsp90 system with the aid of cochaperones Sgt1 and Rar1 that act cooperatively with Hsp90 to form allosterically regulated dynamic complexes. In this work, functional dynamics and protein structure network modeling are combined to dissect molecular mechanisms of Hsp90 regulation by the client recruiter cochaperones. Dynamic signatures of the Hsp90-cochaperone complexes are manifested in differential modulation of the conformational mobility in the Hsp90 lid motif. Consistent with the experiments, we have determined that targeted reorganization of the lid dynamics is a unifying characteristic of the client recruiter cochaperones. Protein network analysis of the essential conformational space of the Hsp90-cochaperone motions has identified structurally stable interaction communities, interfacial hubs and key mediating residues of allosteric communication pathways that act concertedly with the shifts in conformational equilibrium. The results have shown that client recruiter cochaperones can orchestrate global changes in the dynamics and stability of the interaction networks that could enhance the ATPase activity and assist in the client recruitment. The network analysis has recapitulated a broad range of structural and mutagenesis experiments, particularly clarifying the elusive role of Rar1 as a regulator of the Hsp90 interactions and a stability enhancer of the Hsp90-cochaperone complexes. Small-world organization of the interaction networks in the Hsp90 regulatory complexes gives rise to a strong correspondence between highly connected local interfacial hubs, global mediator residues of allosteric interactions and key functional hot spots of the Hsp90 activity. We have found that cochaperone-induced conformational changes in Hsp90 may be determined by specific interaction networks that can inhibit or promote progression of the ATPase cycle and thus control the recruitment of client proteins.

ASD v2.0: updated content and novel features focusing on allosteric regulation

Allostery is the most direct and efficient way for regulation of biological macromolecule function and is induced by the binding of a ligand at an allosteric site topographically distinct from the orthosteric site. AlloSteric Database (ASD, http://mdl.shsmu.edu.cn/ASD) has been developed to provide comprehensive information on allostery. Owing to the inherent high receptor selectivity and lower target-based toxicity, allosteric regulation is expected to assume a more prominent role in drug discovery and bioengineering, leading to the rapid growth of allosteric findings. In this updated version, ASD v2.0 has expanded to 1286 allosteric proteins, 565 allosteric diseases and 22 008 allosteric modulators. A total of 907 allosteric site-modulator structural complexes and >200 structural pairs of orthosteric/allosteric sites in the allosteric proteins were constructed for researchers to develop allosteric site and pathway tools in response to community demands. Up-to-date allosteric pathways were manually curated in the updated version. In addition, both the front-end and the back-end of ASD have been redesigned and enhanced to allow more efficient access. Taken together, these updates are useful for facilitating the investigation of allosteric mechanisms, allosteric target identification and allosteric drug discovery.

The E104D mutation increases the susceptibility of human triosephosphate isomerase to proteolysis. Asymmetric cleavage of the two monomers of the homodimeric enzyme

The deficiency of human triosephosphate isomerase (HsTIM) generates neurological alterations, cardiomyopathy and premature death. The mutation E104D is the most frequent cause of the disease. Although the wild type and mutant exhibit similar kinetic parameters, it has been shown that the E104D substitution induces perturbation of an interfacial water network that, in turn, reduces the association constant between subunits promoting enzyme inactivation. To gain further insight into the effects of the mutation on the structure, stability and function of the enzyme, we measured the sensitivity of recombinant E104D mutant and wild type HsTIM to limited proteolysis. The mutation increases the susceptibility to proteolysis as consequence of the loss of rigidity of its overall 3-D structure. Unexpectedly, it was observed that proteolysis of wild type HsTIM generated two different stable nicked dimers. One was formed in relatively short times of incubation with proteinase K; as shown by spectrometric and crystallographic data, it corresponded to a dimer containing a nicked monomer and an intact monomer. The formation of the other nicked species requires relatively long incubation times with proteinase K and corresponds to a dimer with two clipped subunits. The first species retains 50% of the original activity, whereas the second species is inactive. Collectively, we found that the E104D mutant is highly susceptible to proteolysis, which in all likelihood contributes to the pathogenesis of enzymopathy. In addition, the proteolysis data on wild type HsTIM illustrate an asymmetric conduct of the two monomers.

Differential modulation of functional dynamics and allosteric interactions in the Hsp90-cochaperone complexes with p23 and Aha1: a computational study

Allosteric interactions of the molecular chaperone Hsp90 with a large cohort of cochaperones and client proteins allow for molecular communication and event coupling in signal transduction networks. The integration of cochaperones into the Hsp90 system is driven by the regulatory mechanisms that modulate the progression of the ATPase cycle and control the recruitment of the Hsp90 clientele. In this work, we report the results of computational modeling of allosteric regulation in the Hsp90 complexes with the cochaperones p23 and Aha1. By integrating protein docking, biophysical simulations, modeling of allosteric communications, protein structure network analysis and the energy landscape theory we have investigated dynamics and stability of the Hsp90-p23 and Hsp90-Aha1 interactions in direct comparison with the extensive body of structural and functional experiments. The results have revealed that functional dynamics and allosteric interactions of Hsp90 can be selectively modulated by these cochaperones via specific targeting of the regulatory hinge regions that could restrict collective motions and stabilize specific chaperone conformations. The protein structure network parameters have quantified the effects of cochaperones on conformational stability of the Hsp90 complexes and identified dynamically stable communities of residues that can contribute to the strengthening of allosteric interactions. According to our results, p23-mediated changes in the Hsp90 interactions may provide "molecular brakes" that could slow down an efficient transmission of the inter-domain allosteric signals, consistent with the functional role of p23 in partially inhibiting the ATPase cycle. Unlike p23, Aha1-mediated acceleration of the Hsp90-ATPase cycle may be achieved via modulation of the equilibrium motions that facilitate allosteric changes favoring a closed dimerized form of Hsp90. The results of our study have shown that Aha1 and p23 can modulate the Hsp90-ATPase activity and direct the chaperone cycle by exerting the precise control over structural stability, global movements and allosteric communications in Hsp90.

Identification of functionally key residues in AMPA receptor with a thermodynamic method

AMPA receptor mediates the fast excitatory synaptic transmission in the central nervous system, and it is activated by the binding of glutamate that results in the opening of the transmembrane ion channel. In the present work, the thermodynamic method developed by our group was improved and then applied to identify the functionally key residues that regulate the glutamate-binding affinity of AMPA receptor. In our method, the key residues are identified as those whose perturbation largely changes the ligand binding free energy of the protein. It is found that besides the ligand binding sites, other residues distant from the binding cleft can also influence the glutamate binding affinity through a long-range allosteric regulation. These allosteric sites include the hinge region of the ligand binding cleft, the dimer interface of the ligand binding domain, the linkers between the ligand binding domain and the transmembrane domain, and the interface between the N-terminal domain and the ligand binding domain. Our calculation results are consistent with the available experimental data. The results are helpful for our understanding of the mechanism of long-range allosteric communication in the AMPA receptor and the mechanism of channel opening triggered by glutamate binding.

MCPath: Monte Carlo path generation approach to predict likely allosteric pathways and functional residues

Allosteric mechanism of proteins is essential in biomolecular signaling. An important aspect underlying this mechanism is the communication pathways connecting functional residues. Here, a Monte Carlo (MC) path generation approach is proposed and implemented to define likely allosteric pathways through generating an ensemble of maximum probability paths. The protein structure is considered as a network of amino acid residues, and inter-residue interactions are described by an atomistic potential function. PDZ domain structures are presented as case studies. The analysis for bovine rhodopsin and three myosin structures are also provided as supplementary case studies. The suggested pathways and the residues constituting the pathways are maximally probable and mostly agree with the previous studies. Overall, it is demonstrated that the communication pathways could be multiple and intrinsically disposed, and the MC path generation approach provides an effective tool for the prediction of key residues that mediate the allosteric communication in an ensemble of pathways and functionally plausible residues. The MCPath server is available at http://safir.prc.boun.edu.tr/clbet_server.

A Mixed Protein Structure Network and Elastic Network Model Approach to Predict the Structural Communication in Biomolecular Systems: The PDZ2 Domain from Tyrosine Phosphatase 1E As a Case Study

Graph theory is being increasingly used to study the structural communication in biomolecular systems. This requires incorporating information on the system's dynamics, which is time-consuming and not suitable for high-throughput investigations. We propose a mixed Protein Structure Network (PSN) and Elastic Network Model (ENM)-based strategy, i.e., PSN-ENM, for fast investigation of allosterism in biological systems. PSN analysis and ENM-Normal Mode Analysis (ENM-NMA) are implemented in the structural analysis software Wordom, freely available at http://wordom.sourceforge.net/ . The method performs a systematic search of the shortest communication pathways that traverse a protein structure. A number of strategies to compare the structure networks of a protein in different functional states and to get a global picture of communication pathways are presented as well. The approach was validated on the PDZ2 domain from tyrosine phosphatase 1E (PTP1E) in its free (APO) and peptide-bound states. PDZ domains are, indeed, the systems whose structural communication and allosteric features are best characterized both in vitro and in silico. The agreement between predictions by the PSN-ENM method and in vitro evidence is remarkable and comparable to or higher than that reached by more time-consuming computational approaches tested on the same biological system. Finally, the PSN-ENM method was able to reproduce the salient communication features of unbound and bound PTP1E inferred from molecular dynamics simulations. High speed makes this method suitable for high throughput investigation of the communication pathways in large sets of biomolecular systems in different functional states.

Small world network strategies for studying protein structures and binding

Small world network concepts provide many new opportunities to investigate the complex three dimensional structures of protein molecules. This mini-review explores the published literature on using small-world network approaches to study protein structure, with emphasis on the different combinations of descriptors that have been tested, on studies involving ligand binding in protein-ligand complexes, and on protein-protein complexes. The benefits and success of small world network approaches, which change the focus from specific interactions to the local environment, even to non-local phenomenon, are described. The purpose is to show the different ways that small world network concepts have been used for building new computational models for studying protein structure and function, and for extending and improving existing modelling approaches.