Experimentally Guided Structural Modeling and Dynamics Analysis of Hsp90–p53 Interactions: Allosteric Regulation of the Hsp90 Chaperone by a Client Protein

Amyloid aggregates induced by the p53-R280T mutation lead to loss of p53 function in nasopharyngeal carcinoma

Nasopharyngeal carcinoma (NPC) is a malignant tumor that is highly prevalent in Southeast Asia, especially in South China. The pathogenesis of NPC is complex, and genetic alterations of tumor suppressors and proto-oncogenes play important roles in NPC carcinogenesis. p53 is unexpectedly highly expressed in NPC and possesses an uncommon mutation of R280T, which is different from a high frequency of hotspot mutations or low expression in other tumors. However, the mechanism of p53 loss of function and its correlation with R280T in NPC are still unclear. In this study, p53 amyloid aggregates were found to be widespread in NPC and can be mainly induced by the R280T mutation. Aggregated p53-R280T impeded its entry into the nucleus and was unable to initiate the transcription of downstream target genes, resulting in decreased NPC cell cycle arrest and apoptosis. In addition, NPC cells with p53-R280T amyloid aggregates also contributed aggressively to tumor growth in vivo. Transcriptome analysis suggested that p53 amyloid aggregation dysregulated major signaling pathways associated with the cell cycle, proliferation, apoptosis, and unfolded protein response (UPR). Further studies revealed that Hsp90, as a key molecular chaperone in p53 folding, was upregulated in NPC cells with p53-R280T aggregation, and the upregulated Hsp90 facilitated p53 aggregation in turn, forming positive feedback. Therefore, Hsp90 inhibitors could dissociate p53-R280T aggregation and restore the suppressor function of p53 in vitro and in vivo. In conclusion, our study demonstrated that p53-R280T may misfold to form aggregates with the help of Hsp90, resulting in the inability of sequestered p53 to initiate the transcription of downstream target genes. These results revealed a new mechanism for the loss of p53 function in NPC and provided novel mechanistic insight into NPC pathogenesis.

Conformational Dynamics and Mechanisms of Client Protein Integration into the Hsp90 Chaperone Controlled by Allosteric Interactions of Regulatory Switches: Perturbation-Based Network Approach for Mutational Profiling of the Hsp90 Binding and Allostery

Understanding the allosteric mechanisms of the Hsp90 chaperone interactions with cochaperones and client protein clientele is fundamental to dissect activation and regulation of many proteins. In this work, atomistic simulations are combined with perturbation-based approaches and dynamic network modeling for a comparative mutational profiling of the Hsp90 binding and allosteric interaction networks in the three Hsp90 maturation complexes with FKBP51 and P23 cochaperones and the glucocorticoid receptor (GR) client. The conformational dynamics signatures of the Hsp90 complexes and dynamics fluctuation analysis revealed how the intrinsic plasticity of the Hsp90 dimer can be modulated by cochaperones and client proteins to stabilize the closed dimer state required at the maturation stage of the ATPase cycle. In silico deep mutational scanning of the protein residues characterized the hot spots of protein stability and binding affinity in the Hsp90 complexes, showing that binding hot spots may often coincide with the regulatory centers that modulate dynamic allostery in the Hsp90 dimer. We introduce a perturbation-based network approach for mutational scanning of allosteric residue potentials and characterize allosteric switch clusters that control mechanism of cochaperone-dependent client recognition and remodeling by the Hsp90 chaperone. The results revealed a conserved network of allosteric switches in the Hsp90 complexes that allow cochaperones and GR protein to become integrated into the Hsp90 system by anchoring to the conformational switch points in the functional Hsp90 regions. This study suggests that the Hsp90 binding and allostery may operate under a regulatory mechanism in which activation or repression of the Hsp90 activity can be pre-encoded in the allosterically regulated Hsp90 dimer motions. By binding directly to the conformational switch centers on the Hsp90, cochaperones and interacting proteins can efficiently modulate the allosteric interactions and long-range communications required for client remodeling and activation.

Effect of HSP90AB1 and CC domain interaction on Bcr-Abl protein cytoplasm localization and function in chronic myeloid leukemia cells

Background

The fusion oncoprotein Bcr-Abl is mostly located in the cytoplasm, which causes chronic myeloid leukemia (CML). After moving into the nucleus, the fusion protein can induce apoptosis of CML cells. The coiled-coil domain (CC domain) of Bcr-Abl protein plays a central role in the subcellular localization. However, how CC domain affects subcellular localization of Bcr-Abl remains unclear.

Methods

Herein, the key proteins interacting with the Bcr-Abl CC domain were screened by immunoprecipitation binding mass spectrometry. The specific site of Bcr-Abl CC domain binding to target protein was predicted by Deep Viewer. Immunoprecipitation assay was used to confirmed the specific sites of protein binding. IF and western blot were used to observe the subcellular localization of target protein. Western blot was used to examine the protein changes. CCK-8, clonal formation test and FCM cycle detection were used to observe the effect of inhibitor on the proliferation ability of CML cells. FCM apoptosis detection was used to observe the level of cells apoptosis.

Results

HSP90AB1 interacts with Bcr-Abl CC domain via N-terminal domain (NTD), preventing the transport of Bcr-Abl protein to the nucleus and maintaining the activation of Bcr-Abl tyrosine kinase. The nucleus-entrapped Bcr-Abl markedly inhibits the proliferation and induces apoptosis of CML cells by activating p73 and repressing the expression of cytoplasmic oncogenic signaling pathways mediated by Bcr-Abl. Moreover, the combination of 17AAG (Tanespimycin) with Leptomycin B (LMB) considerably decreased the proliferation of CML cells.

Conclusion

Our study provides evidence that it is feasible to transport Bcr-Abl into the nucleus as an alternative strategy for the treatment of CML, and targeting the NTD of HSP90AB1 to inhibit the interaction with Bcr-Abl is more accurate for the development and application of HSP90 inhibitor in the treatment of CML and other Bcr-Abl-addicted malignancies.

Dissecting Molecular Principles of the Hsp90 Chaperone Regulation by Allosteric Modulators Using a Hierarchical Simulation Approach and Network Modeling of Allosteric Interactions: Conformational Selection Dictates the Diversity of Protein Responses and Ligand-Specific Functional Mechanisms

Conformational plasticity of the Hsp90 molecular chaperones underlies the diversity of functional mechanisms that these versatile molecular machines employ to coordinate their vast protein clientele in the cellular environment. Despite a steady progress in studies of the Hsp90 machinery, a great deal remains unknown about molecular principles and ligand-specific functional mechanisms of the Hsp90 regulation by allosteric modulators that attracted significant attention because of their therapeutic potential. Due to structural complexity and dynamic nature of the Hsp90 responses to allosteric modulators, the atomistic details about the mode of action of these small molecules continue to be fairly scarce and controversial. In this work, we employ an integrative strategy that encompassed atomistic simulations of the Hsp90 proteins and hierarchical modeling of Hsp90-ligand binding with network analysis to explore functional mechanisms of the Hsp90 regulation by a panel of allosteric modulators (novobiocin, KU-135, KU-174, and KU-32) with different models of action. The results show that functional mechanisms of allosteric modulation in the Hsp90 proteins may be driven by conformational selection principles in which ligands elicit pre-existing states of the unbound chaperone to drive ligand-specific protein responses and distinct scenarios of Hsp90 regulation. We found that novobiocin can selectively sequester an ensemble of open chaperone conformations and inhibit the progression of the functional cycle through a cascade of cumulative dynamic changes. In contrast, KU-32 displayed unique preferences toward partially closed dynamic states, inducing robust allosteric signaling and stimulation of the ATPase cycle. The proposed model of the Hsp90 regulation by allosteric modulators reconciled diverse experimental data and showed that allosteric modulators may operate via targeted exploitation of dynamic landscapes eliciting vastly different protein responses and diverse mechanisms of action.

Roles of computational modelling in understanding p53 structure, biology, and its therapeutic targeting

The transcription factor p53 plays pivotal roles in numerous biological processes, including the suppression of tumours. The rich availability of biophysical data aimed at understanding its structure–function relationships since the 1990s has enabled the application of a variety of computational modelling techniques towards the establishment of mechanistic models. Together they have provided deep insights into the structure, mechanics, energetics, and dynamics of p53. In parallel, the observation that mutations in p53 or changes in its associated pathways characterize several human cancers has resulted in a race to develop therapeutic modulators of p53, some of which have entered clinical trials. This review describes how computational modelling has played key roles in understanding structural-dynamic aspects of p53, formulating hypotheses about domains that are beyond current experimental investigations, and the development of therapeutic molecules that target the p53 pathway.

Exploring Mechanisms of Communication Switching in the Hsp90-Cdc37 Regulatory Complexes with Client Kinases through Allosteric Coupling of Phosphorylation Sites: Perturbation-Based Modeling and Hierarchical Community Analysis of Residue Interaction Networks

Understanding molecular principles underlying chaperone-based modulation of kinase client activity is critically important to dissect functions and activation mechanisms of many oncogenic proteins. The recent experimental studies have suggested that phosphorylation sites in the Hsp90 and Cdc37 proteins can serve as conformational communication switches of chaperone regulation and kinase interactions. However, a mechanism of allosteric coupling between phosphorylation sites in the Hsp90 and Cdc37 during client binding is poorly understood, and the molecular signatures underpinning specific roles of phosphorylation sites in the Hsp90 regulation remain unknown. In this work, we employed a combination of evolutionary analysis, coarse-grained molecular simulations together with perturbation-based network modeling and scanning of the unbound and bound Hsp90 and Cdc37 structures to quantify allosteric effects of phosphorylation sites and identify unique signatures that are characteristic for communication switches of kinase-specific client binding. By using network-based metrics of the dynamic intercommunity bridgeness and community centrality, we characterize specific signatures of phosphorylation switches involved in allosteric regulation. Through perturbation-based analysis of the dynamic residue interaction networks, we show that mutations of kinase-specific phosphorylation switches can induce long-range effects and lead to a global rewiring of the allosteric network and signal transmission in the Hsp90-Cdc37-kinase complex. We determine a specific group of phosphorylation sites in the Hsp90 where mutations may have a strong detrimental effect on allosteric interaction network, providing insight into the mechanism of phosphorylation-induced communication switching. The results demonstrate that kinase-specific phosphorylation switches of communications in the Hsp90 may be partly predisposed for their regulatory role based on preexisting allosteric propensities.

Atomistic Modeling of the ABL Kinase Regulation by Allosteric Modulators Using Structural Perturbation Analysis and Community-Based Network Reconstruction of Allosteric Communications

In this work, we have examined the molecular mechanisms of allosteric regulation of the ABL tyrosine kinase at the atomic level. Atomistic modeling of the ABL complexes with a panel of allosteric modulators has been performed using a combination of molecular dynamics simulations, structural residue perturbation scanning, and a novel community analysis of the residue interaction networks. Our results have indicated that allosteric inhibitors and activators may exert a differential control on allosteric signaling between the kinase binding sites and functional regions. While the inhibitor binding can strengthen the closed ABL state and induce allosteric communications directed from the allosteric pocket to the ATP binding site, the DPH activator may induce a more dynamic open form and activate allosteric couplings between the ATP and substrate binding sites. By leveraging a network-centric theoretical framework, we have introduced a novel community analysis method and global topological parameters that have unveiled the hierarchical modularity and the intercommunity bridging sites in the residue interaction network. We have found that allosteric functional hotspots responsible for the kinase regulation may serve the intermodular bridges in the global interaction network. The central conclusion from this analysis is that the regulatory switch centers play a fundamental role in the modular network organization of ABL as the unique intercommunity bridges that connect the SH2 and SH3 domains with the catalytic core into a functional kinase assembly. The hierarchy of network organization in the ABL regulatory complexes may allow for the synergistic action of dense intercommunity links required for the robust signal transfer in the catalytic core and sparse network bridges acting as the regulatory control points that orchestrate allosteric transitions between the inhibited and active kinase forms.

In silico identification and analysis of the binding site for aminocoumarin type inhibitors in the C-terminal domain of Hsp90

Hsp90 contains two Nucleotide Binding Sites (NBS): one each in its N-terminal domain (NTD) and C-terminal domain (CTD), respectively. Previously we used computational techniques to locate a nucleotide-binding site in the CTD. Nucleotide binding at this site stabilized the structurally labile region within this domain, thus providing a rationale for increased resistance to thermal denaturation and proteolysis. A scan for ligand-binding sites in CTD revealed four potential sites with the requisite volume to accommodate aminocoumarins and -derived inhibitors. Only one of these reproducibly formed docked complexes with inhibitors and showed excellent interactions with residues lining the site. Fortuitously, it was identical to the aforementioned nucleotide-binding site thus providing an explanation for the reported direct competition between inhibitors and nucleotides. Further studies with carefully chosen inhibitors and some inactive analogues provided an explanation for the known Structure-Activity Relationships (SAR) of aminocoumarin and -derived inhibitors. We also performed similar studies of the NTD to discern the reason(s) for its inability to bind aminocoumarins, given the family resemblance to prokaryotic Top-IV and Gyr-B. Our studies permitted the identification of the putative inhibitor binding site in the CTD, an explanation for increased resistance to thermal denaturation and proteolysis upon inhibitor binding as well as direct competition with ATP.

Allosteric Modulation of Human Hsp90α Conformational Dynamics

Central to Hsp90's biological function is its ability to interconvert between various conformational states. Drug targeting of Hsp90's regulatory mechanisms, including its modulation by cochaperone association, presents as an attractive therapeutic strategy for Hsp90 associated pathologies. In this study, we utilized homology modeling techniques to calculate full-length structures of human Hsp90α in closed and partially open conformations and used these structures as a basis for several molecular dynamics based analyses aimed at elucidating allosteric mechanisms and modulation sites in human Hsp90α. Atomistic simulations demonstrated that bound adenosine triphosphate (ATP) stabilizes the dimer by "tensing" each protomer, while adenosine diphosphate (ADP) and apo configurations "relax" the complex by increasing global flexibility, the former case resulting in a fully open "v-like" conformation. Dynamic residue network analysis revealed regions of the protein involved in intraprotein communication and identified several key communication hubs that correlate with known functional sites. Pairwise comparison of betweenness centrality, shortest path, and residue fluctuations revealed that a proportional relationship exists between the latter two measurables and an inverse relationship between these two and betweenness centrality. This analysis showed how protein flexibility, degree of compactness, and the distance cutoff used for network construction influence the correlations between these metrics. These findings are novel and suggest shortest path and betweenness centrality to be more relevant quantities to follow for detecting functional residues in proteins compared to residue fluctuations. Perturbation response scanning analysis identified several potential residue sites capable of modulating conformational change in favor of interstate conversion. For the ATP-bound open conformation, these sites were found to overlap with known Aha1 and client binding sites, demonstrating how naturally occurring forces associated with cofactor binding could allosterically modulate conformational dynamics.

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.

Atomistic simulations and network-based modeling of the Hsp90-Cdc37 chaperone binding with Cdk4 client protein: A mechanism of chaperoning kinase clients by exploiting weak spots of intrinsically dynamic kinase domains

A fundamental role of the Hsp90 and Cdc37 chaperones in mediating conformational development and activation of diverse protein kinase clients is essential in signal transduction. There has been increasing evidence that the Hsp90-Cdc37 system executes its chaperoning duties by recognizing conformational instability of kinase clients and modulating their folding landscapes. The recent cryo-electron microscopy structure of the Hsp90-Cdc37-Cdk4 kinase complex has provided a framework for dissecting regulatory principles underlying differentiation and recruitment of protein kinase clients to the chaperone machinery. In this work, we have combined atomistic simulations with protein stability and network-based rigidity decomposition analyses to characterize dynamic factors underlying allosteric mechanism of the chaperone-kinase cycle and identify regulatory hotspots that control client recognition. Through comprehensive characterization of conformational dynamics and systematic identification of stabilization centers in the unbound and client- bound Hsp90 forms, we have simulated key stages of the allosteric mechanism, in which Hsp90 binding can induce instability and partial unfolding of Cdk4 client. Conformational landscapes of the Hsp90 and Cdk4 structures suggested that client binding can trigger coordinated dynamic changes and induce global rigidification of the Hsp90 inter-domain regions that is coupled with a concomitant increase in conformational flexibility of the kinase client. This process is allosteric in nature and can involve reciprocal dynamic exchanges that exert global effect on stability of the Hsp90 dimer, while promoting client instability. The network-based rigidity analysis and emulation of thermal unfolding of the Cdk4-cyclin D complex and Hsp90-Cdc37-Cdk4 complex revealed weak spots of kinase instability that are present in the native Cdk4 structure and are targeted by the chaperone during client recruitment. Our findings suggested that this mechanism may be exploited by the Hsp90-Cdc37 chaperone to recruit and protect intrinsically dynamic kinase clients from degradation. The results of this investigation are discussed and interpreted in the context of diverse experimental data, offering new insights into mechanisms of chaperone regulation and binding.

Inhibition of heat shock protein 90 rescues glucocorticoid-induced bone loss through enhancing bone formation

Endogenous glucocorticoids (GCs) support normal bone development and bone mass maintenance, whereas long-term exposure to pharmacological dosages of GCs uncouples bone formation and resorption, resulting in GC-induced osteoporosis (GIOP). Heat shock protein 90 (HSP90) chaperoning glucocorticoid receptor (GR) signaling prompts us to speculate that HSP90 plays critical roles in GC-mediated bone formation and GIOP. In the present study, inhibition of HSP90 activity by 17-Demethoxy-17-allyaminogeldanmycin (17-AAG) or knockdown of HSP90 expression by siRNAs attenuated dexamethasone(Dex)-induced GR nuclear accumulation and transcriptional output of GR signaling, whereas overexpression of HSP90α or HSP90β enhanced GR transactivity in C3H10T1/2 cells. Though 17-AAG itself enhanced osteoblastic differentiation, it restored the Dex(10-8M)-induced and Dex(10-6M)-negated osteoblastic differentiation in C3H10T1/2 cells and primary calvarial osteoblasts. Moreover, systemic administration of 17-AAG to mice induced not only osteoclastogenesis but also osteoblastogenesis, whereas bone formation possibly exceeded bone resorption, eventually leading to the increased bone masses. Likewise, systemic administration of 17-AAG to mice restored GC-negated osteoblastogenesis and enhanced GC-induced osteoclastogenesis, similarly, 17-AAG-induced bone formation possibly exceeded both 17-AAG- and GC-induced bone resorption, eventually resulting in rescue of GIOP. Together, the present study has revealed that inhibition of HSP90 restores GIOP through enhancing bone formation, and our findings may help to shed light on the pathogenesis of GIOP and provide targets for the therapeutic intervention of the disease.

In silico identification and computational analysis of the nucleotide binding site in the C-terminal domain of Hsp90

Hsp90 contains two distinct Nucleotide Binding Sites (NBS), in its N-terminal domain (NTD) and C-terminal domain (CTD), respectively. The NTD site belongs to the GHKL super-family of ATPases and has been the subject of extensive characterization. However, a structure of the nucleotide-bound form of CTD is still unavailable. In this study molecular modeling was employed to incorporate experimental data using partial constructs of the CTD, from work published by many research groups, onto existing structural models of its apo- form. Our attempts to locate potential nucleotide ligand-binding sites or cavities yielded one major candidate-a structurally unconventional site-exhibiting the requisite shape and volume for accommodation of tri-phosphate nucleotides. Its structure was refined by molecular dynamics (MD)-based techniques. We reproducibly docked the Mg²⁺ complexed form of ATP, GTP, CTP, TTP and UTP to this putative NBS. These docking simulations and calculated ligand-binding scores are in general agreement with published data about experimentally measured binding to the CTD. The overall pattern of interactions between residues lining the site and docked nucleotides is conserved and broadly similar to that of other nucleotide-binding sites. Our docking simulations suggest that nucleotide binding stabilizes the only structurally labile region, thereby providing a rationale for the increased resistance to thermal denaturation and proteolysis. The docked nucleotides do not intrude onto the surface of residues involved in dimerization or chaperoning. Our molecular modeling permitted recognition of larger structural changes in the nucleotide-bound CTD dimer, including stabilization of helix-2 in both chains and intra- and inter- chain interactions between three residues (I613, Q617, R620).

Protein-Protein Interface and Disease: Perspective from Biomolecular Networks

Protein-protein interactions are involved in many important biological processes and molecular mechanisms of disease association. Structural studies of interfacial residues in protein complexes provide information on protein-protein interactions. Characterizing protein-protein interfaces, including binding sites and allosteric changes, thus pose an imminent challenge. With special focus on protein complexes, approaches based on network theory are proposed to meet this challenge. In this review we pay attention to protein-protein interfaces from the perspective of biomolecular networks and their roles in disease. We first describe the different roles of protein complexes in disease through several structural aspects of interfaces. We then discuss some recent advances in predicting hot spots and communication pathway analysis in terms of amino acid networks. Finally, we highlight possible future aspects of this area with respect to both methodology development and applications for disease treatment.

A phase II trial of ganetespib, a heat shock protein 90 Hsp90) inhibitor, in patients with docetaxel-pretreated metastatic castrate-resistant prostate cancer (CRPC)-a prostate cancer clinical trials consortium (PCCTC) study

INTRODUCTION

Heat shock protein 90 (Hsp90) has been studied as a therapeutic target in many cancers. In preclinical trials, the Hsp90 ATPase inhibitor ganetespib demonstrated potent inhibition of solid tumor growth, with superior potency than prior Hsp90 inhibitors. Given the promising preclinical outcome and favorable pharmacologic properties of ganetespib, we conducted a phase II trial of single-agent ganetespib in patients with metastatic, castrate-resistant prostate cancer (mCRPC). The primary objective of the study was to determine the 6-month progression-free survival (PFS) rate.

METHODS

Patients with mCRPC who had been previously treated with docetaxel were enrolled after meeting eligibility criteria. All patients received ganetespib at 200 mg/m(2) on days 1, 8, and 15 of every 28 days (one cycle). Subjects who tolerated therapy were continued on ganetespib until disease progression. Considering that Hsp90 acetylation may confer insensitivity to Hsp90 inhibitors and maspin inhibits protein deacetylation, maspin-associated molecular markers were evaluated.

RESULTS

Eighteen patients were recruited into the trial; most were Caucasian, had performance status 1, had received prior docetaxel, and were heavily pretreated. Of the 17 patients who were treated, none attained 6-month PFS. Only 2 patients achieved PFS > 4 months. The median PFS was 1.9 months. As per the study design, the trial was terminated after the interim analysis. The most frequent types of Grade 3 toxicity were dehydration, diarrhea, and fatigue. Molecular markers provided little additional insight regarding drug activity.

CONCLUSIONS

Ganetespib demonstrated minimal clinical activity in men with mCRPC. The true 6-month PFS rate was, at most, 0.20. Possible reasons for this include selection of a heavily pretreated patient population and lack of agent potency in patients with mCRPC.

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.

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.