Functional Role and Hierarchy of the Intermolecular Interactions in

Feb 12, 2018 - BiKi Life Sciences: A New Suite for Molecular Dynamics and Related Methods in Drug Discovery. Journal of Chemical Information and Model...
0 downloads 12 Views 6MB Size
Subscriber access provided by UNIV OF DURHAM

Article

Functional Role and Hierarchy of the Intermolecular Interactions in Binding of Protein Kinase Clients to the Hsp90-Cdc37 Chaperone: Structure-Based Network Modeling of Allosteric Regulation Gabrielle Stetz, and Gennady M. Verkhivker J. Chem. Inf. Model., Just Accepted Manuscript • DOI: 10.1021/acs.jcim.7b00638 • Publication Date (Web): 12 Feb 2018 Downloaded from http://pubs.acs.org on February 13, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Journal of Chemical Information and Modeling is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

Functional Role and Hierarchy of the Intermolecular Interactions in Binding of Protein Kinase Clients to the

Hsp90-Cdc37

Chaperone:

Structure-Based

Network Modeling of Allosteric Regulation Gabrielle Stetz1, Gennady M. Verkhivker1, 2 ‡ 1

Graduate Program in Computational and Data Sciences, Department of Computational Sciences,

Schmid College of Science and Technology, Chapman University, One University Drive, Orange, CA 92866, USA 2

Chapman University School of Pharmacy, Irvine, CA 92618, USA

‡corresponding author E-mail: [email protected]

1 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Abstract A fundamental role of the Hsp90-Cdc37 chaperone machinery in mediating conformational development and activation of diverse protein kinase clients is essential for signal transduction. Structural and biochemical studies have demonstrated that characterization of global conformational changes and allosteric interactions in the Hsp90-Cdc37-kinase complexes are central to our understanding of the mechanisms underlying kinase recruitment and processing by the Hsp90-Cdc37 chaperone. The recent cryo-electron microscopy structure of the Hsp90Cdc37-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 characterized functional role and hierarchy of the intermolecular interactions in binding of protein kinase clients to the Hsp90-Cdc37 system. The network analysis revealed important relationships between structural stability, global centrality and functional significance of regulatory hotspots in chaperone regulation and client recognition. A unique asymmetric topography of the intermolecular communities in the chaperone-kinase complex has quantified a central mediating role of Cdc37 in client recognition and allosteric regulation of the chaperone-kinase complex. Modeling of allosteric pathways in the chaperone complex has further clarified structural and energetic signatures of allosteric hotspots, particularly linking sites of post-translational modifications in Hsp90 with their role in allosteric interactions and client regulation. The results of this integrative computational study are compared with a wide range of structural, biochemical and cell-based experiments, offering a robust network-centric model of allosteric regulation and client kinase recognition by the Hsp90Cdc37 chaperone machine.

2 ACS Paragon Plus Environment

Page 2 of 65

Page 3 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

Introduction The 90 kDa heat-shock proteins Hsp90s are evolutionary conserved ubiquitous molecular chaperones that are vital for functioning of all eukaryotic organisms.1,2 The Hsp90 chaperones assist in folding of protein substrates, maintaining stability and activity of numerous signaling client proteins associated with cell growth, proliferation, differentiation, and cell death.3-7 The Hsp90 family consists of four different isoforms: cytosolic constitutive (β) and inducible (α) forms, a mitochondrial member (TRAP1) and an ER-resident protein (GRP94). The monomers in the Hsp90 homodimer feature functional domain organization.

The N-terminal domain

(Hsp90-NTD) is involved in the nucleotide binding and ATP hydrolysis, the middle domain (Hsp90-MD) is often implicated in binding of cochaperones and recognition of client proteins, and the

is involved in constitutive dimerization.8-16

C-terminal domain (Hsp90-CTD)

Crystallographic and biophysical studies have characterized conformational landscape and structural transformations of Hsp90 that are associated with the ATP binding and hydrolysis cycle.17-21 Hydrogen/deuterium exchange mass spectrometry (HX-MS) experiments have revealed conformational plasticity of eukaryotic Hsp90s that are considerably more dynamic than the HtpG chaperones. Importantly, these studies have demonstrated that conformational dynamics of the eukaryotic cytosolic Hsp90 is largely independent on the nucleotide binding and can be mainly driven by stochastic thermal fluctuations.22-24 The cytosolic Hsp90 proteins utilize assistance of cochaperones in managing protein clients in the course of the ATPase cycle.25-27 The Hsp90 interactions with client proteins can be also regulated through a multitude of posttranslational modification (PTM) sites that are broadly distributed across the Hsp90 domains.28-30 Cdc37 is a client-specific cochaperone that assists Hsp90 in managing conformational folding, maturation and activation of protein kinases.,31 3 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 65

The fundamental mechanisms and thermodynamics of the Hsp90-client interactions were first explored in a seminal cell-based interactome mapping study32 revealing that the Hsp90-Cdc37 machinery can recognize kinase clients by sensing their conformational instability, while rejecting stable native folds presented to the chaperone by nonclient kinases. HX-MS studies and functional assays of structurally similar kinase chimeras have revealed that the Hsp90-Cdc37 system can differentiate between weak c-Src kinase client and strong constitutive client v-Src kinase.33,34 According to these studies, the interaction strength of chaperone-client interactions correlates with the degree of unfolding cooperativity and solvent exposure of the regulatory regions in the kinase catalytic domain.33,34 The emerging mechanism of the Hsp90-Cdc37 client recognition implies that Cdc37 may function as a principal gate-keeper to the Hsp90 system by recognizing conformational instability of client kinases. Cdc37 plays a central role in chaperoning of kinase clients

since

impairment of Cdc37 function can dramatically affect

binding of oncogenic kinase clients with the Hsp90-Cdc37 chaperone system.35 NMR spectroscopy studies have shown that Cdc37 facilitates binding of client substrates by eliciting a cascade of

concerted conformational changes and orchestrating

an allosteric cross-talk

between multiple interaction sites of the interacting proteins.36,37 However, atomistic details and mechanics of client processing by the Hsp90-Cdc37 chaperone system have been poorly understood until very recently. The preliminary EM reconstruction of the Hsp90-Cdc37-Cdk4 kinase complex indicated an asymmetric intertwined topology of multiprotein assembly, but the atomistic details of the binding interfaces hypothetical,

and functional interactions remained largely

offering different and often conflicting explanations of the underlying

mechanism.38 The latest cryo–electron microscopy structure of the Hsp90-Cdc37-Cdk4 kinase complex has produced shock waves in the chaperone community, 4 ACS Paragon Plus Environment

unveiling a completely

Page 5 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

unexpected structural arrangement, in which the Hsp90 and Cdc37 trap a partially disordered kinase state, providing framework for understanding allosteric mechanisms of kinase recruitment to the chaperone.39,40

Allosteric regulation of the Hsp90 interactions

with

cochaperones and

binding partners

determines versatility of the chaperone system in coordinating signaling responses and timely processing of

multiple clients. Computational studies of allosteric signaling in single domain

proteins and complex multiprotein assemblies have gained a significant momentum in recent years41-49 reshaping the conventional view of allostery. The statistical thermodynamics-based model of allosteric regulation has clarified many aspects of population-shift and dynamic-based mechanisms of allostery50-52 suggesting that protein ensembles of preexisting conformational states and communication pathways can be modulated through allosteric interactions with substrates and binding partners.53 Recent investigations from our laboratory and several other groups have suggested that Hsp90 interactions with cochaperones and client proteins can be highly dynamic and cooperative in nature, leading to modulation of conformational stability in the clients and eliciting global conformational rearrangements in the chaperones.54-62 These studies indicated that mediating sites in the Hsp90 interaction networks may signify hotspots of chaperone regulation.

Despite significant insights in biochemical and structural characterizations of the Hsp90-Cdc37 interactions with the diverse kinase clientele, the mechanism by which allosteric cross-talk of the binding interfaces promotes integration of client kinases to the Hsp90 system remains elusive. The recent reevaluation of the role of the Hsp90 chaperones in the kinase lifecycle has pointed to 5 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

several conceptual issues concerning regulation of the Hsp90-Cdc37 interactions with client proteins.40 Some of these outstanding questions include a) detailing of functional specialization of Cdc37 N-terminal (Cdc37-NTD) and C-terminal domains (Cdc37-M/C) in recruitment of partially unfolded kinase clients; b) quantifying allosteric cross-talk and hierarchy of the intermolecular binding interfaces in the Hsp90-Cdc37-client complexes; c) clarifying regulatory role of functional sites in allosteric mechanisms of Hsp90-Cdc37 binding with kinase clients; c) characterizing allosteric communication pathways and mediating centers of allosteric signaling in the Hsp90-client complexes. In the current study, we combined molecular simulations and protein stability analysis with modeling of the interaction networks and communication pathways to characterize

allosteric

interactions in the Hsp90-Cdc37-kinase complex and

identify regulatory hotspots that control client recognition by the chaperone machine. The incorporation of both dynamic residue correlations and coevolutionary residue dependencies provided a novel theoretical framework for

reconstruction of allosteric interaction networks

and communication pathways in the Hsp90 structures. A community analysis reveals modularity of allosteric interaction networks and hierarchy of binding interfaces. By identifying principal stabilizing communities and communication pathways at the atomic details, we obtain novel molecular insights into regulatory role and coupling of conserved PTM sites that serve as allosteric switch points of the Hsp90-Cdc37 system. The results of this investigation provide a novel rationale to the latest experimental data on client recruitment by the Hsp90-Cdc37 system, detailing how Hsp90 and Cdc37 act cooperatively in protecting unstable kinase clients and restraining their undesirable oncogenic activity.

6 ACS Paragon Plus Environment

Page 6 of 65

Page 7 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

Materials and Methods Discrete Molecular Dynamics and All-Atom Reconstruction of Conformational Ensembles

The structures of the human Hsp90-Cdc37-Cdk4 complex (pdb id 5FWM, 5FWK, 5FWL) were obtained from the Protein Data Bank (www.rcsb.org). During structure preparation process, hydrogen atoms and missing residues were initially added and assigned according to the WHATIF program web interface.63 The unresolved segments in the Hsp90-Cdc37-Cdk4 structures were modeled and reconstructed using template-based loop prediction approach ArchPRED.64 We employed the formalism and implementation of the discrete molecular dynamics (DMD) simulations65,66 to rapidly generate conformational ensembles based on structures of the human Hsp90-Cdc37-Cdk4 complex. In the DMD approach, the protein structures were modeled as systems consisting of C α residue-based beads interacting through a discontinuous square well potential. The details of the DMD model implemented in this study were outlined in our earlier investigations.67

Conformational ensembles generated in DMD

simulations were subsequently subjected to all-atom reconstruction by PULCHRA method68 and CG2AA tool.69 The all-atom conformations were optimized using 3Drefine method70 that utilizes atomic-level minimization with a composite physics and knowledge-based force fields. Protein Stability Calculations A systematic alanine scanning of protein residues was performed using FoldX approach.71 We utilized a graphical user interface for the FoldX calculations72 that was implemented as a plugin for the YASARA molecular graphics suite package.73 To obtain ensemble-based estimates of protein stability changes and ensure efficiency of FoldX computations, we evaluated the

7 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 65

average ∆∆G values using representative samples (~200-300) of the equilibrium ensembles within a modified FoldX protocol.74 Protein Structure Networks and Modeling of Communication Pathways For network-based analysis, we employed a graph-based representation of protein structures in which residues as network nodes and inter-residue edges represent residue interactions.75-77 The details of graph construction using residue interaction cut-off strength ( Imin ) were presented in our previous studies of molecular chaperones.78-80. The edges in the residue interaction network are weighted based on coevolutionary mutual information80 and dynamic residue correlations couplings obtained from molecular simulations.81 in this model, weight wij is defined as the element of a matrix measuring the generalized correlation coefficient rMI ( x i , x j ) between residue fluctuations in structural and coevolutionary dimensions. The composite residue vector describes structural residue positions and respective proximity-based coevolutionary score80: wij = − log[rMI ( x i , x j )] (1)

Mutual Information (MI) analysis is used to estimate the extent of the mutual coevolutionary relationship between pairs of positions in the protein family.82,

83

MI was calculated using

MISTIC approach.83 The edge lengths in the residue interaction network are obtained using the generalized correlation coefficients rMI ( x i , x j ) associated with both dynamic correlation extracted from DMD trajectories and mutual information shared by each pair of residues. The length (i.e. weight) of the edge that connects nodes i and j wij = − log[rMI ( x i , x j )] is calculated from the generalized correlation coefficient between this nodes.80 The ensemble of shortest paths is determined from matrix of communication distances by the Floyd-Warshall algorithm84 that compares all possible paths between each pair of residue nodes. 8 ACS Paragon Plus Environment

Page 9 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

Network Centrality Analysis and Community Detection Using the constructed protein structure networks, we computed the residue-based betweenness parameter. The betweenness of residue i is defined to be the sum of the fraction of shortest paths between all pairs of residues that pass through residue i :

N

Cb (ni ) = ∑ j 2.0 kcal/moll in these positions signified a central functional role of this interfacial hotspot. Protein stability changes in the Cdk4 kinase revealed differential stabilization of the kinase lobes, which is manifested in a structurally stable C-lobe and a more flexible N-lobe. Indeed, the N-lobe residues showed a greater tolerance to alanine substitutions, while mutations of the C-lobe residues produced a more significant destabilization effect (Figure 3D). In the complex, energetic polarization of the kinase lobes was even more apparent, revealing only small changes in protein stability caused by mutations of the N-lobe residues. This effect was particularly pronounced in the partially disordered N-lobe regions: the αC-helix (residues 50-67), αC-β4 loop (residues 67-74) and a β4-β5 sheet (residues 75-95). In contrast, substitutions in the C-lobe regions resulted in the markedly larger energetic changes, which were comparable with the corresponding effects in the folded kinase form (Figure 3D). Hence, client recruitment into the Hsp90-Cdc37 chaperone may enhance energetic instability of the N-lobe, while maintaining rigidity of the C-lobe. The observed differential stabilization of the kinase lobes may reflect the intrinsic dynamic preferences of kinase clients that can be further amplified 17 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 65

upon integration into the chaperone system. These findings are consistent with several recent studies suggesting that the excessive conformational

mobility of the N-lobe regions

in

oncogenic kinases can shape up structural preferences of activation cancer mutations often populating these regions which allows for robust modulation of functional dynamics.99, 100 Structural mapping highlighted the fact that high stability sites in the Hsp90 dimer were often spatially proximal and formed tight interaction modules. Some of these hydrophobic clusters stability hotspots (I288-F304-L332 and P287-W289-F304-Y356) were immediately adjacent to the NTD-MD regions (L388, I390 and R392 residues), facilitating the inter-domain interactions and rigidifying the catalytic site environment (Figure 4A, B). Structural analysis of stability hotspots also illustrated the interconnectivity of these centers in the complex (Figure 4). Several highly stable residues are assembled into an interacting cluster (I404-V381-I408-L369-I370L401) that is directly linked to the interfacial lysine residues (K399, K402, K406, K411) making specific interactions with Cdc37-NTD (Figure 4 C). Another hotspot of protein stability is formed by the residue cluster F376-L343-F433-N436 that anchored the Src-loop at the binding interface with the kinase N-lobe (Figure 4D). These clusters include or directly connected with several regulatory centers. In particular, a significant number of conserved PTM sites (K402, K406, K411, K438) are strategically positioned to couple the energetic hotspot centers in Hsp90MD with Cdc37-NTD and Cdk4 N-lobe regions (Figure 4C,D). We argue that structural preferences of these PTM centers in Hsp90 may be related with their fundamental role as dynamic carriers of allosteric signaling and mediators of binding interactions in the chaperoneclient complexes. In general, we found that stable energetic centers in the Hsp90 chaperone are coupled with more flexible interfacial residues that become immobilized in the chaperone-client complex. Our findings corroborated with the evidence that regions of high and low structural 18 ACS Paragon Plus Environment

Page 19 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

stabilities could be often proximal that enables coupling between global conformational changes and binding to protein substrates. The observed partitioning of stable and flexible regions in the chaperone-kinase complex is consistent with the idea that protein ability to evolve new functions may be linked with structural modularity, whereby flexible recognition regions are linked to a highly ordered central scaffold.101,102 This pattern was observed in the Hsp90-Cdc37-Cdk4 complex, where the central Hsp90-MD core is enriched by highly stable energetic hotspots, while flanking recognition loops are flexible in the unbound form and become ordered in the chaperone-kinase complex.

19 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 4. Structural mapping of protein stability hotspot clusters. (A) The overview of the intermolecular interfaces and energetic hotspot centers in the Hsp90-Cdc37-Cdk4 complex. Hsp90-A monomer is shown in orange ribbons and Hsp90-B monomer B is in pink ribbons. Cdc37 is in red ribbons, and Cdk4 kinase is in yellow ribbons. The location of three major centers of energetic hotspots is indicated by rectangular boxes colored in red (with a close-up in panel B), blue (a close-up is presented in panel C), and green (a close-up is in panel D). (B) Structural mapping of the intramolecular cluster connecting energetic hotspots near the NTD-MD interface. For clarity only Hsp90-B monomer (colored in pink ribbons) is shown. The high stability clusters I288-F304-L332 and P287-W289-F304-Y356 are shown in domaincolored spheres. These residues are directly linked to L388, I390 and R392 (in sticks) that bridge NTD and MD regions near the catalytic site. Catalytic residue R392 and ATP are shown in atomcolored spheres. (C) A close-up of the protein stability hotspot formed by the hydrophobic cluster (I404-V381-I408-L369-I370-L401). The hotspot residues are shown in domain-colored pink spheres (Hsp90-B monomer of the dimer). This stability hotspot is immediately proximal to the Hsp90 binding residues K399, K402, R405, K406 (shown in pink sticks ) that form a network of salt bridges with Cdc37-NTD residues D14,D15, E16,D17,E18 ( in red sticks). Hsp90-B in pink ribbons, Cdc37 in red ribbons, and Cdk4 is in yellow ribbons. (D) The protein stability hotspot cluster F376-L343-F433-N436 (in pink spheres) anchors the Src-loop residues R338, F341, L343, F344 (in sticks) that make contacts with the Cdk4 N-lobe regions. The Cdc37-M/C domain is shown in red ribbons and Cdk4-N lobe is shown in yellow ribbons.

20 ACS Paragon Plus Environment

Page 20 of 65

Page 21 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

Network Analysis Identifies Mediating Centers of Allosteric Interactions To investigate allosteric coupling of the binding interfaces and identify global mediating centers that control kinase client recruitment, we characterized the residue interaction network in the Hsp90-Cdc37-Cdk4 complex. Using a graph-based representation of protein structures77, we integrated topological connectivity of protein residues and dynamic contact maps of residue cross-correlations with coevolutionary residue dependencies in the construction of the residue interaction network.80 A global centrality parameter, residue betweenness (also termed residue centrality in the text), was explored to identify mediating centers of allosteric interactions in the unbound and bound Hsp90 structures. The majority of high centrality residues were located in the Hsp90-MD regions and assembled in tight interaction clusters, likely to strengthen robustness of key mediating centers (Figure 5A, B). The notable peaks in the Hsp90-MD corresponded to W289, F341, L343, Y356, I370, Y373, F376, I400, F433, E515, and D518 residues. These sites coincided with major stabilization centers, suggesting that a small group of hotspot clusters may be involved in mediating global interaction network. To determine client-induced modulation of the global interaction network, we evaluated differences in the residue centrality of the Hsp90 residues in the unbound and client-bound forms of Hsp90 (Figure 5 C, D). Several peaks in the Hsp90-MD regions were induced upon complex formation, particularly a pronounced increase in the network centrality of the Src loop (residues 341-350) involved in the interface with the disordered kinase regions (Figure 5D). These high centrality sites also overlapped with the important stabilization clusters in this region (F341-L343-F344-F376 and F376-L343-F433N436). Another noticeable change in the network organization of the complex was a markedly increased centrality of the CTD helix 15 (residues 599-611).

21 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 5. Residue centrality profiles of the Hsp90 dimer. Residue-based centrality distributions of the Hsp90 monomers in the chaperone-kinase complex (A, B). The client-induced differential changes in the residue centrality of the Hsp90 residues (C, D) are obtained from differences between residue betweenness values in the complex and unbound Hsp90 dimer. The distributions are derived by averaging computations of network parameters over equilibrium ensembles extracted from DMD trajectories. The distributions are shown in domain-colored bars. The Hsp90-NTD residues are shown in red bars, Hsp90-MD residues are in blue bars, and Hsp90CTD residues are in green bars.

22 ACS Paragon Plus Environment

Page 22 of 65

Page 23 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

Hence, kinase recruitment may induce significant changes in the network centrality of the Hsp90 recognition regions that could serve as mediators of allosteric interactions. The residue centrality profiling of the Cdc37 and Cdk4 domains pointed to the key mediating role of the Cdc37-NTD in the allosteric interaction network. Two pronounced sharp peaks were seen in the Cdc37-NTD profile (Figure 6A). The first peak corresponded to cluster of residues (D14, D15, E16, D17, and E18) involved in the interface with Hsp90-B monomer. The second peak (residues H20, P21, N22, I23, F29, W31, and N32) corresponded to the binding interface formed by Cdc37-NTD with the Cdk4 C-lobe (residues 147-153). The importance of this interface was also reflected in the centrality distribution of Cdk4 client, featuring peak near C-lobe residues 147-153 (Figure 6B). These profiles reflected the importance of the Cdc37-NTD/Cdk4 N-lobe interface in mediating allosteric communications.

The observed rewiring of the global interaction network near the Hsp90 intermolecular interfaces may be linked with dynamic exchanges between the Hsp90 and Cdk4 that manifested in stabilization of the Hsp90-CTD binding regions involved in protection of unstable kinase client. Our observations also suggested that the hierarchy of spatially separated binding interfaces may be important for allosteric mechanism of client recruitment. According to the network analysis, Cdc37-NTD/Cdk4 C-lobe and Hsp90-Cdk4 N-lobe interfacial regions are central for client binding, whereas binding interface formed between Cdc37-M/C and Cdk4 N-lobe regions may play a supporting role. In this mechanism, Cdc37-NTD and Hsp90 can act cooperatively during client recognition and loading.

23 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 6. Residue centrality profiles of Cdc37 and Cdk4 in the Hsp90-Cdc37-Cdk4 complex. Residue-based centrality distributions of Cdc37 (A) and Cdk4 (B) are shown in domain-colored bars. The Cdc37-NTD (residues 1-147) is in blue bars, Cdc37-M/C (residues 148-260) in red bars. The Cdk4-N lobe (residues 1-99) is in blue bars, and Cdk4 C-lobe (residues 100-295) is in red bars. (C) Structural mapping of high centrality residues in the Hsp90-Cdc37-Cdk4 complex. For clarity, only Hsp90-B (pink ribbons) is shown. Cdc37 in red ribbons, Cdk4 in yellow ribbons. (D) A close-up of high centrality residues at the Hsp90-B/Cdc37-NTD interface. The Cdc37-NTD residues (pS13, D14, D15, E16, E18, H20, P21, N22, I23, F29, W31, and R33) are highlighted in red spheres. The Hsp90-B interfacial residues (K399, K402, R405, and K406) and adjacent hydrophobic sites (I370, F376, and Y373) are shown in pink spheres. 24 ACS Paragon Plus Environment

Page 24 of 65

Page 25 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

Asymmetric Topography of Intermolecular Communities Reveals Functional Role of Cdc37 in Protecting Kinase Client in the Repressed State Recent experimental revelations indicated that in the absence of binding partners the closed Hsp90 state may be intrinsically asymmetric at the MD-CTD interface, whereas Cdc37 and client binding force Hsp90 to adopt a symmetric state.40 In this context, we examined organization and hierarchy of allosteric interaction networks in the chaperone-kinase complex. By using modular decomposition of the global interaction network, we quantified topography of allosteric interactions revealing an asymmetry of the community maps that may have functional relevance for the mechanism of kinase recruitment. Importantly, we found that high centrality residues may be assembled into local modules of strongly interacting nodes (Figure 7). Structural map of communities in the Hsp90-Cdc37-Cdk4 complex highlighted the prevalence of the intermolecular modules anchored by Cdc37-NTD residues and illustrated an asymmetric topography of mediating centers (Figure 7). The spatial distribution of the intermolecular communities displayed a clear bias towards the Hsp90-Cdc37-NTD and Cdc37-NTD/Cdk4-Clobe interfaces. We detected a densely populated network of the intermolecular communities at the Hsp90-MD interface with Cdc37-NTD (Figure 7B). Stable modules in this region are formed through an extensive network of electrostatic interactions and salt bridges: K399 (Hsp90)-D15 (Cdc37-NTD)-D18 (Cdc37-NTD), K402 (Hsp90)-R405 (Hsp90)-E16 (Cdc37-NTD), and N407 (Hsp90)-V403 (Hsp90)-V116 (Cdc37-NTD). The largest and most stable community in the Hsp90-Cdc37-NTD region is formed through interactions of the Hsp90 residues (L119, Q326, L327, and L395) and functional sites of Cdc37 (Y4, W7, I10) (Figure 7B). This broad community links phosphorylated sites on Cdc37 (Y4) with the energetic hotspots of Hsp90 stability and contributes decisively to the stability of the Hsp90-Cdc37 interface. 25 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 7. Structural map of major intermolecular communities. (A) An overview of the communities (in colored spheres), Hsp90-A monomer is shown in orange ribbons and Hsp90-B monomer is in pink ribbons. Cdc37 is in red, and Cdk4 is in yellow ribbons. The residues in major communities are shown in spheres colored according to domain annotation. (B) The communities linking Hsp90-B and Cdc37-NTD : K399(Hsp90)-D15(Cdc37)-D18(Cdc37) (in orange spheres), K402(Hsp90)-R405(Hsp90)-E16(Cdc37) (in blue spheres), N407(Hsp90)V403(Hsp90)-V116(Cdc37) (in salmon spheres), L327(Hsp90)-Q326(Hsp90)-L119(Cdc37)W7(Cdc37)-I400(Hsp90)-V403(Hsp90)-L3095(Hsp90)-P387(Hsp90)-Y4(Cdc37)-Q397(Hsp90)I10(Cdc37) (in hot pink spheres). The key residues in the largest community are annotated. (C) The communities linking Cdc37-NTD and Cdk4 C-lobe: H20(Cdc37)-I23(Cdc37)-L28(Cdc37)26 ACS Paragon Plus Environment

Page 26 of 65

Page 27 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

F130(Cdk4)-L131(Cdk4)-F159(Cdk4)-I136(Cdk4)-W31(Cdc37)-L128(Cdk4) (in blue spheres), and N22(Cdc37)-Q123(Cdk4)-V154(Cdk4)-V148(Cdk4)-L120(Cdk4) (in salmon spheres). (D) A single largest community at the Cdc37-M/C interface with the kinase N-lobe. The contributing kinase residues (Y17, V37, E56, V57, L60, and L74) are shown in yellow spheres, and Cdc37M/C residues (W168, D169, I241, K242) are shown in red spheres. Note contribution of conserved regulatory kinase sites E56, L60, and L74 from the regulatory αC-helix.

The analysis of the Hsp90-Cdc37 communities provided an interesting and novel rationale to various functional and biochemical experiments on regulation of the Hsp90-Cdc37 interactions. It has been recognized that both Hsp90 and Cdc37 are delicately regulated by phosphorylation. In particular, Cdc37 can be phosphorylated on Y4 and Y298 by Yes kinase103 and on S13 position by CKII kinase104 that are both Hsp90 clients. Phosphomimetic mutations of these Cdc37 sites can lead to impaired Cdc37-kinase binding and prevent kinase recruitment to the Hsp90-Cdc37. In our analysis, Cdc37-Y4 phosphorylation site is embedded into the largest intermolecular community. As a result, mutations in this single position can simultaneously disrupt multiple interactions with community neighbors and severely impair the stability of the complex. In the Hsp90–Cdc37–Cdk4 structure, Cdc37-S13 is phosphorylated and contributes to the intercommunity connectivity in the complex. We observed that the interaction cluster formed by pS13 with H33 and R36 of the Cdc37-NTD can link the Hsp90-Cdc37 and Cdc37-NTD-Cdk4 communities. A mediating role of S13 in connecting multiple communities in this region can explain a severe detrimental effect of phosphomimetic mutants in this position on regulation of the Hsp90-Cdc37 interactions. The network-centric interpretation of functional significance of 27 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

these sites for Hsp90-Cdc37 regulation is based on their integrating role in large intermolecular communities. In network terms, modifications of these sites would simultaneously weaken numerous interactions, compromising modularity and interconnectivity of the global network and consequently undermining efficiency of allosteric signaling. The asymmetry of the interfacial communities is manifested in a dense network of local modules is anchored by Cdc37-NTD residues, while Cdc37-M/C regions can be only weakly coupled to the Hsp90 and Cdk4 binding partners. The principal stabilizing community at the Cdc37-NTD interface with Cdk4 linked together Cdc37-NTD residues (H20, I23, L28, and W31) and Cdk4C-lobe residues (G127, L128, F130, L131, I136, H138, and F159). Strikingly, this community connected the recognition Cdc37-NTD motif (HPN) with the regulatory HRD and DFG kinase motifs (H138, F159). Notably, H138 and F159 residues from Cdk4 belong to the regulatory Rspine (L60, L74, H138, F159, and D196) whose cooperative assembly controls kinase activation.105 Through dense interactions formed in this community Cdc37-NTD can protect a partially unfolded kinase conformation, while also constraining the R-spine in the inactive arrangement (Figure 7C). These findings provide a compelling rationale to the experimental observations that the Hsp90 chaperone can force kinase to remain in the inactive repressed state and only after dissociation clients can be activated.106-108 Another important component of the global network was formed at the interface between Hsp90-MD/CTD and unfolded kinase segment of the N-lobe. The intermolecular communities W598(Hsp90)-M606(Hsp90)R612(Hsp90)-E94(Cdk4) and L343(Hsp90)-F344(Hsp90)-L91(Cdk4)-V89(Cdk4) protect the unfolded kinase segment. Several small communities can be formed by the Cdc37-M/C residues, including F341(Hsp90)-E235(Cdc37)-Q236(Cdc37), Q208(Cdc37)-R166(Cdc37)-V72(Cdk4) and Y17(Cdc37)-K35(Cdk4)-D169(Cdc37)-V20(Cdk4). The most interesting finding was the 28 ACS Paragon Plus Environment

Page 28 of 65

Page 29 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

identification of a single largest community at the Cdc37-M/C interface with the kinase N-lobe that linked kinase residues (Y17, V37, E56, V57, L60, L74, A79) with Cdc37-M/C residues (W168, D169, I241, K242). Of particular significance is evidence that conserved regulatory kinase sites, such as catalytic residue E56 (from catalytic salt bridge pair K35-E56) and R-spine residues L60, L74 from the regulatory αC-helix (Figure 7D). As a result, Cdc37-M/C domain can aid in protecting the inactive R-spine arrangement of the kinase client. An important result of this analysis is that Hsp90-Cdc37 interfacial regions are coupled through communities with the key regulatory residues of Cdk4 that are evolutionary conserved and functionally essential across the spectrum of kinase clients. These kinase residues belong to the regulatory αC-helix (L60, L74, H138, and F159), the catalytic HRD motif (H138) and DFG motif (F159). Notably, these regulatory residues are known to control a dynamic kinase.105 Consequently, Hsp90 and Cdc37 can act cooperatively at the binding interfaces to shield unstable kinase client and simultaneously control positions of the regulatory kinase residues responsible for Cdk4 activation. We argue that these interactions may underlie a mechanism by which how Hsp90 and Cdc37 can restrain undesirable oncogenic activity of the kinase clients. Our results suggested that preservation of these binding interactions and corresponding communities may have evolutionary benefit for efficient recognition and processing of kinase clientele by the chaperone system. To summarize, the asymmetric topography of the intermolecular communities in the chaperone-kinase complex supported the indispensable mediating role of Cdc37 in allosteric regulation, particularly revealing specific contributions of Cdc37 domains in protecting unstable kinase client in the inactive dormant state. The observed modularity and asymmetry of the intermolecular communities may be an important global signature of the Hsp90-client interactions that enables efficient allosteric signaling in the complex required for productive kinase processing. 29 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Modeling of Allosteric Communication Pathways: Communication Hotspots Define Regulatory Switch Points of the Hsp90-Cdc37 Interactions with Kinase Client Using the results of community decomposition, we explored edge betweenness (or edge centrality) in the global interaction network as a proxy for modeling of allosteric communication pathways. This parameter is defined as the ratio of all the shortest paths passing through a particular edge to the total number of shortest paths in the network. Accordingly, edges that interconnect many local communities (inter-modular bridges) would tend to have the higher edge betweenness values as compared to the intra-modular edges. We evaluated the edge centrality distribution based on all short inter-residue paths in the Hsp90-Cdc37-Cdk4 complex (Figure 8A). The distribution featured a sharp decline and long tail indicative of a small-world organization in which only a small number of edges have a significantly higher betweenness values as compared to the average. By considering the ensemble of short paths that connected the Cdc37-NTD/Cdk4 C-lobe interface with the ATP binding site, we detected a pronounced shift in the distribution towards high edge centrality values (Figure 8A). Only a moderate shift in the edge centrality was observed for the ensemble of short pathways connecting the kinase lobes (Figure 8A). Accordingly, the ensemble of signaling routes that pass through Cdc37-NTD interface may be important for coordinating allosteric communications in the chaperone complex. By considering short paths going through the high centrality edges, we also found that the Hsp90 and Cdc37-NTD regions harbor most of these communication hubs (Figure 8B). Strikingly, the inter-community bridges are often formed by phosphorylation and lysine acetylation sites, suggesting the previously unrecognized role of PTMs as potential transmitters of allosteric communications in the chaperone-client complexes.

30 ACS Paragon Plus Environment

Page 30 of 65

Page 31 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

Figure 8. The distribution of edge centrality in the Hsp90-Cdc37-Cdk4 complex. (A)The edge centrality distribution is based on all short inter-residue paths in the Hsp90-Cdc37-Cdk4 complex (red bars). The edge centrality distribution for the ensemble of short pathways connecting the kinase lobes blue bars). The edge centrality distribution for the ensemble of short paths connecting the Cdc37-NTD/Cdk4 C-lobe interfacial residues with the ATP binding site (in green bars). (B) The domain-based distribution of high centrality inter-community edges (red bars). (C) Structural mapping of the major communities in the Hsp90 dimer (shown in green spheres for Hsp90-A monomer and cyan spheres for Hsp90-B monomer. The positions of important conserved PTM sites in the human Hsp90β are highlighted in red spheres and annotated. Note that these PTM sites can often assume role of the inter-community bridges. 31 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

To highlight this finding, we mapped major communities in the unbound Hsp90 dimer along with positions of conserved PTM sites (Figure 8C). Notably, the PTM sites are distributed along the community edges or are located at the inter-community boundaries. These observations suggested that inter-domain communication pathways in the chaperone-kinase complex may exploit these regulatory sites for efficient modulation of allosteric signals.

Using a community-hopping model of allosteric pathways80 we reconstructed ensemble of allosteric communication pathways that pass through high centrality edges (Figure 9A). The topography of this map revealed major highways absorbing most of the communication ‘traffic’ between the ATP binding sites of Hsp90 and kinase client (Figure 9B). By “zooming in” on the residue pairs with the highest edge betweenness, we identified key interactions that enable allosteric communications in the Hsp90-Cdc37-Cdk4 complex. The inter-modular edges of highest centrality included M363-S365, K350-S365, K399 (Hsp90-MD)-D15 (Cdc37-NTD), and K399 (Hsp90-MD)-E18 (Cdc37-NTD) connections (Figure 9B). Among communication hotspots that coordinate cross-talk between binding interfaces are Cdc37-NTD residues S13, F29, and W31. These sites are involved in high centrality bridges with Hsp90-MD (S13) and the regulatory DFG motif of Cdk4 (F29, W31). Our model suggested that these communication hotspots may function as important network bridges or switches that coordinate allosteric interactions and signal transmission in the chaperone complex. Strikingly, many high centrality edges connected experimentally known points of allosteric regulation in the Hsp90. Furthermore, some of the key allosteric hubs in the complex corresponded to conserved PTM sites of the Hsp90 chaperone.

32 ACS Paragon Plus Environment

Page 32 of 65

Page 33 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

Figure 9. Structural map of major communication pathways passing through high centrality edges in the Hsp90-Cdc37-Cdk4 complex. (A) A general overview and topography of communication pathways. Residues along these pathways are shown in green and cyan spheres for Hsp90 A and B monomers respectively. The positions of PTM sites in the Hsp90 chaperone are highlighted by red spheres. (B) A close-up of high centrality inter-modular connections is shown in domain-colored spheres. The functional sites of post-translational modifications in this communication map are highlighted in blue spheres and annotated.

33 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

It has been recognized that phosphorylation of the Hsp90 permits regulation of the conformational cycle by targeting spatially remote switch points of the inter-domain communication and these changes affect client maturation and the interaction with cochaperones.30 Our results provided a rationale of the experimentally observed effects by showing that global communication traffic may depend on the inter-modular bridges formed by PTM sites. In particular, phosphorylation of S13 in Cdc37-NTD is necessary for Cdc37 interaction with Hsp90 and also for chaperoning of numerous kinase clients including Cdk4. Consistent with these experiments, our results showed that the interactions formed by pS13 with H33 and R36 of the Cdc37-NTD form high centrality bridges connecting Hsp90-Cdc37 and Cdc37-Cdk4 communities. The emergence of high centrality edges near Cdc37-NTD binding interfaces further strengthened our arguments that Cdc37-NTD may be central in kinase recognition. Recent experiments identified residue S365 as an isoform-specific phosphorylation site in human Hsp90β chaperone.109 Mutations in this position may drastically affect client maturation and is believed to impair progression of the Hsp90-kinase cycle by altering conformational equilibrium and disfavoring the closed form of the Hsp90 dimer. Mutational analysis of K294 in human Hsp90α (K286 in human Hsp90β) demonstrated that its acetylation status is a strong determinant of client and cochaperone binding, since Hsp90 variants that cannot be acetylated at this position could impair chaperone function.110 These experimental studies attributed regulatory role of this site, located at the juncture between the charged linker and Hsp90-MD, to the local dynamic changes and loss of the intramolecular interactions caused by proximity to residues T285 and P287 from a hydrophobic cluster at the NTD-MD interface.

34 ACS Paragon Plus Environment

Page 34 of 65

Page 35 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

Our results showed that the two key inter-modular bridges involve S365 site (M363-S365 and K350-S365) and may serve as central mediators of long-range communications in the complex (Figure 9B). Several other PTM sites, including K286, K319, K350 and K399, also contribute to the highest centrality edges (Figure 9B). These findings indicated that regulatory PTM sites can be interconnected along main communication pathways, assuming role of primary control points and gate-keepers of allosteric signaling in the chaperone-kinase complex. According to our results, the detrimental functional effect of mutations in these PTM sites can arise from global changes in the interaction network and alteration of primary communication routes in the system. In this novel interpretation of the experimental data, mutations of PTM sites may not necessarily signal structural changes or disfavor the closed form of Hsp90 dimer, but rather lead to rewiring of major communication pathways and consequently compromise efficiency of allosteric signaling in the system. Based on our results, we propose that several other phosphorylation and lysine acetylation PTM sites (K350, S391, K399, S452) emerged as major control points of longrange communications (Figure 9B) could also serve as potential allosteric switches of the Hsp90 interactions with kinase clients.

To determine specific role of communication hotspots in an allosteric cross-talk between binding interfaces, we generated ensembles of allosteric communication pathways connecting spatially separated functional sites in the Hsp90 and Cdk4 client. We applied a community-hopping model of allosteric pathways80 in which a pair of residues may be viewed as ‘source’ and ‘sink’ of the pathway. 100 suboptimal paths were computed for each of the studied residue connections. We generated the ensemble of short length pathways that connected the catalytic site (R392) in the Hsp90-NTD with functional sites in Cdk4 client kinase. One of the optimal pathways linked 35 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Hsp90-R392 site with phosphorylation site T172 in C-lobe of Cdk4 client (Figure 10A). The optimal pathway proceeded from the catalytic R392 residue through (L395-Q397-I400-V403) and (L388-I362-R392) communities, reaching to the inter-modular edge formed by M363 and S365 residues, and connecting to the interfacial hotspot K399. This route traversed through the Hsp90/Cdc37-NTD interfacial bridge to reach to the interfacial hotspots of Cdc37-NTD (W31, F29, D14, and E18), and then linked via the regulatory DFG motif of Cdk4 (F159) to the phosphorylation site T172 of Cdk4 (Figure 10A). We also modeled communication pathways that link R392 site in the Hsp90-NTD with the Cdk4 N-lobe. In this case, optimal routes passed through Hsp90 community (R392-L388-I362) and navigated via high centrality edges connecting I362, M363, S365, K350 and K348 residues (Figure 10A). By traversing through mediating centers of Hsp90-MD, the paths reached the Hsp90 Src-loop residues (E344 and L343) and connected to the unfolded region of the kinase N-lobe (I87, V89, and L91). Significantly, in both cases, the ensembles of communication pathways showed a strong tendency to proceed through evolutionary conserved regulatory PTM sites.

In another series of computational experiments, we generated the ensemble of communication pathways that connect functional regions in Cdc37-NTD with the Cdc37 N-lobe (Figure 9B). Two equally probable alternative routes emerged, where the first path connected Cdc37-NTD residues (F29, W31) with the DFG motif in the C-lobe of Cdk4 and interfacial hotspots of Hsp90-CTD (W598), reaching out to the regulatory αC-helix of Cdk4 (T53, E67) (Figure 9B). The alternative route proceeded high centrality edges that connect communities near another regulatory center in the Hsp90-MD (F344, F433). Hence, both optimal pathways utilized structurally stable allosteric hotspots to efficiently transmit signal between binding interfaces. 36 ACS Paragon Plus Environment

Page 36 of 65

Page 37 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

Figure 10. Modeling of allosteric communication pathways in the Hsp90-Cdc37-Cdk4 complex. (A) Representative optimal short paths from the ensemble of communication pathways connecting catalytic residue R392 (Hsp90-B NTD) with functional sites in Cdk4 client kinase. (B) The most probable pathways from the ensemble of short inter-residue routes connecting Cdc37-NTD functional sites with Cdk4 N-lobe regions. The most probable paths are shown and the respective residues are shown in domain-colored spheres and annotated.

37 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Atomistic reconstruction of allosteric pathways revealed a well-defined topography of communication traffic in the complex that enables a cross-talk between regulatory residues and binding interfaces, where Cdc37-NTD regions play a central integral role. We found that Cdc37NTD residues form major mediating bridges between communities and dominate allosteric communications in the complex. Interestingly, optimal routes that interconnect spatially separated binding interfaces do not include the interfacial Cdc37-M/C residues, supporting the notion that C-terminal domain of Cdc37 may not be play a secondary role in an allosteric crosstalk during client recruitment. This model of client recognition and processing is reminiscent of a molecular mechanism underlying client binding by other chaperones studied by NMR experiments.111 our analysis also suggested that the regulatory role of allosteric control points in Hsp90 and Cdc37 may be determined by their role in the inter-community bridges of the global interaction network. The results rationalize the NMR experiments of the Hsp90-Cdc37 binding with protein kinase clients36,

37

and provide novel insights by revealing atomistic view of

allosteric signaling through a privileged group of functional centers.

Conclusions

In this study, we attempted to examine several outstanding questions concerning regulation of the Hsp90-Cdc37 interactions with client proteins by using an array of computational methods. The methodological framework of our study is based on integration of dynamic and networkcentric approaches into a unified and robust computational strategy for quantifying mechanisms of allosteric regulation. Using these approaches, we characterized allosteric interaction networks and hierarchy of the intermolecular binding interfaces in the Hsp90-Cdc37-client complexes Our 38 ACS Paragon Plus Environment

Page 38 of 65

Page 39 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

results demonstrated that allosteric regulation of the Hs90 binding may couple the intra-domain and inter-domain interactions in a complex and dynamic interaction network. Through systematic identification of stabilization centers and centrality sites in the interaction networks, we characterized allosteric cross-talk and hierarchy of the intermolecular binding interfaces in the Hsp90-Cdc37-client complex. According to our findings, Cdc37-NTD can mediate kinase recognition by orchestrating a cascade of allosteric interactions with the Hsp90 chaperone and Cdk4 client. The results also quantified functional specialization of the Hsp90 and Cdc37 domains in protecting unstable kinase client in its repressed inactive form. Modeling of allosteric pathways in the chaperone complex has further clarified structural and energetic signatures of allosteric hotspots, particularly linking PTM sites with their important role in allosteric interactions. By identifying principal communities and reconstructing communication pathways, we obtained novel molecular insights about specific functional roles of key PTM sites serving as control points of allosteric signaling. The emerged interconnectivity and coupling of regulatory sites in the Hsp90 chaperone may be useful in expanding array of therapeutic strategies for design of novel allosteric inhibitors targeting allosteric switch points that control Hsp90 interactions with oncogenic kinase clients.

39 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

AUTHOR INFORMATION * Corresponding Author Phone: 714-516-4586 Fax: 714-532-6048 E-mail: [email protected]

Acknowledgment This work was partly supported by institutional funding from Chapman University.

ABBREVIATIONS 90-kilodalton Heat shock protein, Hsp90; human protein Hsp90, N-terminal domain of Hsp90, Hsp90-NTD; middle domain of Hsp90, Hsp90-MD; C-terminal domain of Hsp90, Hsp90-CTD; N-terminal domain of Cdc37, Cdc37-NTD; middle and C-terminal domain of Cdc37, Cdc37M/C; hydrogen/deuterium exchange mass spectrometry, HX-MS; molecular dynamics, MD; discrete molecular dynamics, DMD; flexibility-rigidity index, FRI; mutual Information, MI; multiple sequence alignment, MSA; mutual information server to infer coevolution, MISTIC.

40 ACS Paragon Plus Environment

Page 40 of 65

Page 41 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

References

(1) Pearl, L. H.; Prodromou, C. Structure, Function, and Mechanism of the Hsp90 Molecular Chaperone. Adv. Protein Chem. 2001, 59, 157-186. (2) Richter, K.; Buchner, J. Hsp90: Chaperoning Signal Transduction. J. Cell. Physiol. 2001, 188, 281-290. (3) Picard, D. Heat-shock Protein 90, a Chaperone for Folding and Regulation. Cell. Mol. Life Sci. 2002, 59, 1640-1648. (4) Young, J. C.; Agashe, V. R.; Siegers, K.; Hartl, F. U. Pathways of Chaperone-Mediated Protein Folding in the Cytosol. Nat. Rev. Mol. Cell. Biol. 2004, 5, 781-791. (5) McClellan, A. J.; Xia, Y.; Deutschbauer, A. M.; Davis, R. W.; Gerstein, M.; Frydman, J. Diverse Cellular Functions of the Hsp90 Molecular Chaperone Uncovered Using Systems Approaches. Cell 2007, 131, 121-135. (6) Whitesell, L.; Lindquist, S. L. HSP90 and the Chaperoning of Cancer. Nat. Rev. Cancer

2005, 5, 761-72. (7) Taipale, M.; Jarosz, D. F.; Lindquist, S. L. HSP90 at the Hub of Protein Homeostasis: Emerging Mechanistic Insights. Nat. Rev. Mol. Cell. Biol. 2010, 11, 515-528. (8) Pearl, L. H.; Prodromou, C. Structure and Mechanism of the Hsp90 Molecular Chaperone Machinery. Annu. Rev. Biochem. 2006, 75, 271-294. (9) Pearl, L. H.; Prodromou, C.; Workman, P. The Hsp90 Molecular Chaperone: An Open and Shut Case for Treatment. Biochem. J. 2008, 410, 439-453. (10) Krukenberg, K. A.; Street, T. O.; Lavery, L. A.; Agard, D. A. Conformational Dynamics of the Molecular Chaperone Hsp90. Q. Rev. Biophys. 2011, 44, 229-255.

41 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(11) Didenko, T.; Duarte, A. M.; Karagoz, G. E.; Rudiger, S. G. Hsp90 Structure and Function Studied by NMR Spectroscopy. Biochim. Biophys. Acta 2012, 1823, 636-647. (12) Jackson, S. E. Hsp90: Structure and Function. Top. Curr. Chem. 2013, 328, 155-240. (13) Li, J.; Buchner, J. Structure, Function and Regulation of the Hsp90 Machinery. Biomed. J.

2013, 36, 106-117. (14) Prodromou, C. Mechanisms of Hsp90 Regulation. Biochem. J. 2016, 473, 2439-2452. (15) Pearl, L. H. Review: The HSP90 Molecular Chaperone - An Enigmatic ATPase. Biopolymers 2016, 105, 594-607. (16) Schopf, F. H.; Biebl, M. M.; Buchner, J. The HSP90 Chaperone Machinery. Nat. Rev. Mol. Cell. Biol. 2017, 18, 345-360. (17) Ali, M. M.; Roe, S. M.; Vaughan, C. K.; Meyer, P.; Panaretou, B.; Piper, P. W.; Prodromou, C.; Pearl, L. H. Crystal Structure of an Hsp90-nucleotide-p23/Sba1 Closed Chaperone Complex. Nature 2006, 440, 1013-1017. (18) Shiau, A. K.; Harris, S. F.; Southworth, D. R.; Agard, D. A. Structural Analysis of E. coli Hsp90 Reveals Dramatic Nucleotide-Dependent Conformational Rearrangements. Cell 2006, 127, 329-340. (19) Krukenberg, K. A.; Forster, F.; Rice, L. M.; Sali, A.; Agard, D. A. Multiple Conformations of E. coli Hsp90 in Solution: Insights into the Conformational Dynamics of Hsp90. Structure

2008, 16, 755-765. (20) Southworth, D. R.; Agard, D. A. Species-Dependent Ensembles of Conserved Conformational States Define the Hsp90 Chaperone ATPase Cycle. Mol. Cell 2008, 32, 631-640.

42 ACS Paragon Plus Environment

Page 42 of 65

Page 43 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

(21) Lavery, L. A.; Partridge, J. R.; Ramelot, T. A.; Elnatan, D.; Kennedy, M. A.; Agard, D. A. Structural Asymmetry in the Closed State of Mitochondrial Hsp90 (TRAP1) Supports a TwoStep ATP Hydrolysis Mechanism. Mol. Cell 2014, 53, 330-343. (22) Mickler, M.; Hessling, M.; Ratzke, C.; Buchner, J.; Hugel, T. The Large Conformational Changes of Hsp90 are Only weakly Coupled to ATP Hydrolysis. Nat. Struct. Mol. Biol. 2009, 16, 281-286. (23) Ratzke, C.; Berkemeier, F.; Hugel, T. Heat shock Protein 90's Mechanochemical Cycle is Dominated by Thermal Fluctuations. Proc. Natl. Acad. Sci. U. S. A. 2012, 109, 161-166. (24) Graf, C.; Lee, C. T.; Eva Meier-Andrejszki, L.; Nguyen, M. T.; Mayer, M. P. Differences in Conformational Dynamics within the Hsp90 Chaperone Family Reveal Mechanistic Insights. Front. Mol. Biosci. 2014, 1, 4. (25) Li, J.; Soroka, J.; Buchner, J. The Hsp90 Chaperone Machinery: Conformational Dynamics and Regulation by Co-chaperones. Biochim. Biophys. Acta 2012, 1823, 624-635. (26) Rohl, A.; Rohrberg, J.; Buchner, J. The Chaperone Hsp90: Changing Partners for Demanding Clients. Trends Biochem. Sci. 2013, 38, 253-262. (27) Karagoz, G. E.; Rudiger, S. G. Hsp90 Interaction with Clients. Trends Biochem. Sci. 2015, 40, 117-125. (28) Mollapour, M.; Neckers, L. Post-translational Modifications of Hsp90 and their Contributions to Chaperone Regulation. Biochim. Biophys. Acta 2012, 1823, 648-655. (29) Cloutier, P.; Coulombe, B. Regulation of Molecular Chaperones Through Post-translational Modifications: Decrypting the Chaperone Code. Biochim. Biophys. Acta 2013, 1829, 443-454. 43 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(30) Soroka, J.; Wandinger, S. K.; Mausbacher, N.; Schreiber, T.; Richter, K.; Daub, H.; Buchner, J. Conformational Switching of the Molecular Chaperone Hsp90 via Regulated Phosphorylation. Mol. Cell 2012, 45, 517-528. (31) Mandal, A. K.; Lee, P.; Chen, J. A.; Nillegoda, N.; Heller, A.; DiStasio, S.; Oen, H.; Victor, J.; Nair, D. M.; Brodsky, J. L.; Caplan, A. J. Cdc37 has Distinct Roles in Protein Kinase Quality Control that Protect Nascent Chains from Degradation and Promote Posttranslational Maturation. J. Cell. Biol. 2007, 176, 319-328. (32) Taipale, M.; Krykbaeva, I.; Koeva, M.; Kayatekin, C.; Westover, K. D.; Karras, G. I.; Lindquist, S. Quantitative Analysis of HSP90-Client Interactions Reveals Principles of Substrate Recognition. Cell 2012, 150, 987-1001. (33) Boczek, E. E.; Reefschlager, L. G.; Dehling, M.; Struller, T. J.; Hausler, E.; Seidl, A.; Kaila, V. R.; Buchner, J. Conformational Processing of Oncogenic v-Src Kinase by the Molecular Chaperone Hsp90. Proc. Natl. Acad. Sci. U. S. A. 2015, 112, E3189-E3198. (34) Luo, Q.; Boczek, E. E.; Wang, Q.; Buchner, J.; Kaila, V. R. Hsp90 Dependence of a Kinase is Determined by its Conformational Landscape. Sci. Rep. 2017, 7, 43996. (35) Smith, J. R.; Clarke, P. A.; de Billy, E.; Workman, P. Silencing the Cochaperone CDC37 Destabilizes Kinase Clients and Sensitizes Cancer Cells to HSP90 Inhibitors. Oncogene 2009, 28, 157-169. (36) Eckl, J. M.; Scherr, M. J.; Freiburger, L.; Daake, M. A.; Sattler, M.; Richter, K. Hsp90.Cdc37 Complexes with Protein Kinases Form Cooperatively with Multiple Distinct Interaction Sites. J. Biol. Chem. 2015, 290, 30843-30854.

44 ACS Paragon Plus Environment

Page 44 of 65

Page 45 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

(37) Keramisanou, D.; Aboalroub, A.; Zhang, Z.; Liu, W.; Marshall, D.; Diviney, A.; Larsen, R. W.; Landgraf, R.; Gelis, I. Molecular Mechanism of Protein Kinase Recognition and Sorting by the Hsp90 Kinome-Specific Cochaperone Cdc37. Mol. Cell 2016, 62, 260-271. (38) Vaughan, C. K.; Gohlke, U.; Sobott, F.; Good, V. M.; Ali, M. M.; Prodromou, C.; Robinson, C. V.; Saibil, H. R.; Pearl, L. H. Structure of an Hsp90-Cdc37-Cdk4 Complex. Mol. Cell 2006, 23, 697-707. (39) Verba, K. A.; Wang, R. Y.; Arakawa, A.; Liu, Y.; Shirouzu, M.; Yokoyama, S.; Agard, D. A. Atomic structure of Hsp90-Cdc37-Cdk4 Reveals that Hsp90 Traps and Stabilizes an Unfolded Kinase. Science 2016, 352, 1542-1547. (40) Verba, K. A.; Agard, D. A. How Hsp90 and Cdc37 Lubricate Kinase Molecular Switches. Trends Biochem. Sci. 2017, 42, 799-811. (41) Ghosh, A.; Sakaguchi, R.; Liu, C.; Vishveshwara, S.; Hou, Y. M. Allosteric Communication in Cysteinyl tRNA Synthetase: A Network of Direct and Indirect Readout. J. Biol. Chem. 2011, 286, 37721-37731. (42) Sethi, A.; Tian, J.; Derdeyn, C. A.; Korber, B.; Gnanakaran, S. A Mechanistic Understanding of Allosteric Immune Escape Pathways in the HIV-1 Envelope Glycoprotein. PLoS Comput. Biol. 2013, 9, e1003046. (43) Rivalta, I.; Sultan, M. M.; Lee, N. S.; Manley, G. A.; Loria, J. P.; Batista, V. S. Allosteric Pathways in Imidazole Glycerol Phosphate Synthase. Proc. Natl. Acad. Sci. U. S. A. 2012, 109, E1428-E1436. (44) Ricci, C. G.; Silveira, R. L.; Rivalta, I.; Batista, V. S.; Skaf, M. S. Allosteric Pathways in the PPARgamma-RXRalpha Nuclear Receptor Complex. Sci. Rep. 2016, 6, 19940.

45 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(45) Bhattacharya, S.; Vaidehi, N. Differences in Allosteric Communication Pipelines in the Inactive and Active states of a GPCR. Biophys. J. 2014, 107, 422-434.

(46) General, I. J.; Liu, Y.; Blackburn, M. E.; Mao, W.; Gierasch, L. M.; Bahar, I. ATPase Subdomain IA is a Mediator of Interdomain Allostery in Hsp70 Molecular Chaperones. PLoS Comput. Biol. 2014, 10, e1003624.

(47) Farabella, I.; Pham, T.; Henderson, N. S.; Geibel, S.; Phan, G.; Thanassi, D. G.; Delcour, A. H.; Waksman, G.; Topf, M. Allosteric Signalling in the Outer Membrane Translocation Domain of PapC Usher. Elife 2014, 3, e03532.

(48) Gasper, P. M.; Fuglestad, B.; Komives, E. A.; Markwick, P. R.; McCammon, J. A. Allosteric Networks in Thrombin Distinguish Procoagulant vs. Anticoagulant Activities. Proc. Natl. Acad. Sci. U. S. A. 2012, 109, 21216-21222.

(49) Hertig, S.; Latorraca, N. R.; Dror, R. O. Revealing Atomic-Level Mechanisms of Protein Allostery with Molecular Dynamics Simulations. PLoS Comput. Biol. 2016, 12, e1004746.

(50) Popovych, N.; Sun, S.; Ebright, R. H.; Kalodimos, C. G. Dynamically Driven Protein Allostery. Nat. Struct. Mol. Biol. 2006, 13, 831-838.

(51) Kalodimos, C. G. Protein Function and Allostery: A Dynamic Relationship. Ann. N. Y. Acad. Sci. 2012, 1260, 81-86.

(52) Motlagh, H. N.; Wrabl, J. O.; Li, J.; Hilser, V. J. The Ensemble Nature of Allostery. Nature

2014, 508, 331-339.

46 ACS Paragon Plus Environment

Page 46 of 65

Page 47 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

(53) Tsai, C. J.; Nussinov, R. A Unified View of "How Allostery Works". PLoS Comput. Biol.

2014, 10, e1003394.

(54) Morra, G.; Verkhivker, G.; Colombo, G. Modeling Signal Propagation Mechanisms and Ligand-Based Conformational Dynamics of the Hsp90 Molecular Chaperone Full Length Dimer. PLoS Comput. Biol. 2009, 5, e1000323.

(55) Dixit, A.; Verkhivker, G.M. Probing Molecular Mechanisms of the Hsp90 Chaperone: Biophysical Modeling Identifies Key Regulators of Functional Dynamics. PLoS One 2012, 7, e37605. (56) Morra, G.; Potestio, R.; Micheletti, C.; Colombo, G. Corresponding Functional Dynamics Across the Hsp90 Chaperone Family: Insights from a Multiscale Analysis of MD Simulations. PLoS Comput. Biol. 2012, 8, e1002433. (57) Blacklock, K.; Verkhivker, G.M. Differential Modulation of Functional Dynamics and Allosteric Interactions in the Hsp90-Cochaperone Complexes with p23 and Aha1: A computational Study. PLoS One 2013, 8, e71936. (58) Blacklock, K.; Verkhivker, G. M. Experimentally Guided Structural Modeling and Dynamics Analysis of Hsp90-p53 Interactions: Allosteric Regulation of the Hsp90 Chaperone by a Client Protein. J. Chem. Inf. Model. 2013, 53, 2962-2978. (59) Blacklock, K.; Verkhivker, G. M. Allosteric regulation of the Hsp90 dynamics and stability by client recruiter cochaperones: protein structure network modeling. PLoS One 2014, 9, e86547. (60) Blacklock, K.; Verkhivker, G. M. Computational modeling of allosteric regulation in the hsp90 chaperones: a statistical ensemble analysis of protein structure networks and allosteric communications. PLoS Comput. Biol. 2014, 10, e1003679. 47 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 48 of 65

(61) Rehn, A.; Moroni, E.; Zierer, B. K.; Tippel, F.; Morra, G.; John, C.; Richter, K.; Colombo, G.; Buchner, J. Allosteric Regulation Points Control the Conformational Dynamics of the Molecular Chaperone Hsp90. J. Mol. Biol. 2016, 428, 4559-4571. (62) Vettoretti, G.; Moroni, E.; Sattin, S.; Tao, J.; Agard, D. A.; Bernardi, A.; Colombo, G. Molecular Dynamics Simulations Reveal the Mechanisms of Allosteric Activation of Hsp90 by Designed Ligands. Sci. Rep. 2016, 6, 23830. (63) Hekkelman, M. L.; Te Beek, T. A.; Pettifer, S. R.; Thorne, D.; Attwood, T. K.; Vriend, G. WIWS: A Protein Structure Bioinformatics Web Service Collection. Nucleic Acids Res. 2010, 38, W719-W723. (64) Fernandez-Fuentes, N.; Zhai, J.; Fiser, A. ArchPRED: A Template Based Loop Structure Prediction Server. Nucleic Acids Res. 2006, 34, W173-W176. (65) Shirvanyants, D.; Ding, F.; Tsao, D.; Ramachandran, S.; Dokholyan, N. V. Discrete Molecular Dynamics: An Efficient and Versatile Simulation Method for Fine Protein Characterization. J. Phys. Chem. B. 2012, 116, 8375-8382. (66) Emperador, A.; Meyer, T.; Orozco, M. Protein Flexibility from Discrete Molecular Dynamics Simulations Using Quasi-Physical Potentials. Proteins 2010, 78, 83-94. (67) Stetz, G.; Tse, A.; Verkhivker, G. M., Ensemble-Based Modeling and Rigidity Decomposition of Allosteric Interaction Networks and Communication Pathways in CyclinDependent Kinases: Differentiating Kinase Clients of the Hsp90-Cdc37 Chaperone. PLoS One

2017, 12, e0186089. (68) Rotkiewicz, P.; Skolnick, J. Fast Procedure for Reconstruction of Full-Atom Protein Models from

Reduced

Representations.

J.

Comput.

Chem.

48 ACS Paragon Plus Environment

2008,

29,

1460-1465.

Page 49 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

(69) Lombardi, L. E.; Marti, M. A.; Capece, L. CG2AA: Backmapping Protein Coarse-Grained Structures. Bioinformatics 2016, 32, 1235-1237. (70) Bhattacharya, D.; Nowotny, J.; Cao, R.; Cheng, J. 3Drefine: An Interactive Web Server for Efficient Protein Structure Refinement. Nucleic Acids Res. 2016, 44, W406-W409. (71) Guerois, R.; Nielsen, J.E.; Serrano, L. Predicting Changes in the Stability of Proteins and Protein Complexes: A Study of More than 1000 Mutations. J. Mol. Biol. 2002, 320, 369-387. (72) Schymkowitz, J.; Borg, J.; Stricher, F.; Nys, R.; Rousseau, F.; Serrano L. The FoldX Web Server: An Online Force Field. Nucleic Acids Res. 2005, 33, W382-W388. (73) Van Durme, J.; Delgado, J.; Stricher, F.; Serrano, L.; Schymkowitz, J.; Rousseau, F. A Graphical Interface for the FoldX Force Field. Bioinformatics 2011, 27, 1711-1712. (74) Christensen, N.J.; Kepp, K.P. Stability Mechanisms of Laccase Isoforms Using a Modified FoldX Protocol Applicable to Widely Different Proteins. J. Chem. Theory Comput. 2013, 9, 3210-3223. (75) del Sol, A.; Fujihashi, H.; Amoros, D.; Nussinov, R. Residue Centrality, Functionally Important Residues, and Active Site Shape: Analysis of Enzyme and Non-Enzyme Families. Protein Sci. 2006, 15, 2120-2128. (76) del Sol, A.; Fujihashi, H.; Amoros, D.; Nussinov, R. Residues Crucial for Maintaining Short Paths in Network Communication Mediate Signaling in Proteins. Mol. Syst. Biol. 2006, 2, 2006.0019. (77) Vijayabaskar, M.S.; Vishveshwara, S. Interaction Energy Based Protein Structure Networks. Biophys. J. 2010, 99, 3704-3715.

49 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(78) Stetz, G.; Verkhivker, G. M. Dancing through Life: Molecular Dynamics Simulations and Network-Centric Modeling of Allosteric Mechanisms in Hsp70 and Hsp110 Chaperone Proteins. PLoS One 2015, 10, e0143752.61. (79) Stetz, G.; Verkhivker, G. M. Probing Allosteric Inhibition Mechanisms of the Hsp70 Chaperone Proteins Using Molecular Dynamics Simulations and Analysis of the Residue Interaction Networks. J, Chem. Inf. Model. 2016, 56, 1490-1517. (80) Stetz, G.; Verkhivker, G. M. Computational Analysis of Residue Interaction Networks and Coevolutionary Relationships in the Hsp70 Chaperones: A Community-Hopping Model of Allosteric Regulation and Communication. PLoS Comput. Biol. 2017, 13, e1005299. (81) Sethi, A.; Eargle, J.; Black, A.A.; Luthey-Schulten, Z. Dynamical Networks in tRNA: Protein Complexes. Proc. Natl. Acad. Sci. U.S.A. 2009, 106, 6620-6625. (82) Marino Buslje, C.; Teppa, E.; Di Domenico, T.; Delfino, J. M.; Nielsen, M. Networks of High Mutual Information Define the Structural Proximity of Catalytic Sites: Implications for Catalytic Residue Identification. PLoS Comput. Biol. 2010, 6, e1000978. (83) Simonetti, F. L.; Teppa, E.; Chernomoretz, A.; Nielsen, M.; Marino Buslje, C. MISTIC: Mutual Information Server to Infer Coevolution. Nucleic Acids Res. 2013, 41, W8-W14. (84) Floyd, R.W. Algorithm 97: Shortest Path. Commun. A.C.M. 1962, 5, 345. (85) Palla, G.; Derényi, I.; Farkas, I.; Vicsek, T. Uncovering the Overlapping Community Structure of Complex Networks in Nature and Society. Nature 2005, 435, 814-818. (86) Adamcsek, B.; Palla, G.; Farkas, I.J., Derényi, I.; Vicsek, T. CFinder: Locating Cliques and Overlapping Modules in Biological Networks. Bioinformatics 2006, 22, 1021-1023. (87) Girvan, M.; Newman, M. E. Community Structure in Social and Biological Networks. Proc. Natl. Acad. Sci. U. S. A. 2002, 99, 7821-7826. 50 ACS Paragon Plus Environment

Page 50 of 65

Page 51 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

(88) Newman, M. E.; Girvan, M. Finding and Evaluating Community Structure in Networks. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 2004, 69, 026113. (89) Chen, B.; Piel, W. H.; Gui, L.; Bruford, E.; Monteiro, A. The HSP90 Family of Genes in the Human Genome: Insights into their Divergence and Evolution. Genomics 2005, 86, 627-637. (90) Chen, B.; Zhong, D.; Monteiro, A. Comparative Genomics and Evolution of the HSP90 Family of Genes Across all Kingdoms of Organisms. BMC Genomics 2006, 7, 156. (91) Opron, K.; Xia, K.; Wei, G. W. Fast and Anisotropic Flexibility-Rigidity Index for Protein Flexibility and Fluctuation Analysis. J. Chem. Phys. 2014, 140, 234105. (92) Nguyen, D. D.; Xia, K.; Wei, G. W. Generalized Flexibility-Rigidity Index. J. Chem. Phys.

2016, 144, 234106. (93) Opron, K.; Xia, K.; Burton, Z.; Wei, G. W. Flexibility-Rigidity Index for Protein-Nucleic Acid Flexibility and Fluctuation Analysis. J. Comput. Chem. 2016, 37, 1283-1295. (94) Czemeres, J.; Buse, K.; Verkhivker, G. M. 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. PLoS One 2017, 12, e0190267. (95) Synoradzki, K.; Bieganowski, P. Middle Domain of Human Hsp90 Isoforms Differentially Binds Aha1 in Human Cells and Alters Hsp90 Activity in Yeast. Biochim. Biophys. Acta 2015, 1853, 445-452. (96) Hawle, P.; Siepmann, M.; Harst, A.; Siderius, M.; Reusch, H. P.; Obermann, W. M. The Middle Domain of Hsp90 Acts as a Discriminator Between Different Types of Client Proteins. Mol. Cell. Biol. 2006, 26, 8385-8395.

51 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(97) Cunningham, C. N.; Krukenberg, K. A.; Agard, D. A. Intra- and Intermonomer Interactions are Required to Synergistically Facilitate ATP hydrolysis in Hsp90. J. Biol. Chem. 2008, 283, 21170-21178. (98) Lee, C. C.; Lin, T. W.; Ko, T. P.; Wang, A. H. The Hexameric Structures of Human Heat Shock Protein 90. PLoS One 2011, 6, e19961. (99) Paladino, A.; Morra, G.; Colombo, G., Structural Stability and Flexibility Direct the Selection of Activating Mutations in Epidermal Growth Factor Receptor Kinase. J.Chem Inf. Model. 2015, 55, 1377-1387. (100) Kiel, C.; Benisty, H.; Llorens-Rico, V.; Serrano, L. The Yin-Yang of Kinase Activation and Unfolding Explains the Peculiarity of Val600 in the Activation Segment of BRAF. Elife

2016, 5, e12814. (101) Dellus-Gur, E.; Toth-Petroczy, A.; Elias, M.; Tawfik, D. S. What Makes a Protein Fold Amenable to Functional Innovation? Fold Polarity and Stability Trade-Offs. J. Mol. Biol. 2013, 425, 2609-2621. (102) Toth-Petroczy, A.; Tawfik, D. S. The Robustness and Innovability of Protein Folds. Curr. Opin. Struct. Biol. 2014, 26, 131-138. (103) Xu, W.; Mollapour, M.; Prodromou, C.; Wang, S.; Scroggins, B. T.; Palchick, Z.; Beebe, K.; Siderius, M.; Lee, M. J.; Couvillon, A.; Trepel, J. B.; Miyata, Y.; Matts, R.; Neckers, L. Dynamic Tyrosine Phosphorylation Modulates Cycling of the HSP90-P50(CDC37)-AHA1 Chaperone Machine. Mol. Cell 2012, 47, 434-443. (104) Miyata, Y. Protein Kinase CK2 in Health and Disease: CK2: The Kinase Controlling the Hsp90 Chaperone Machinery. Cell. Mol. Life Sci. 2009, 66, 1840-1849.

52 ACS Paragon Plus Environment

Page 52 of 65

Page 53 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

(105) Kornev, A. P.; Taylor, S. S. Dynamics-Driven Allostery in Protein Kinases. Trends Biochem. Sci. 2015, 40, 628-647. (106) Xu, W.; Yuan, X.; Beebe, K.; Xiang, Z.; Neckers, L. Loss of Hsp90 Association Upregulates Src-Dependent ErbB2 Activity. Mol. Cell. Biol. 2007, 27, 220-228. (107) Citri, A.; Gan, J.; Mosesson, Y.; Vereb, G.; Szollosi, J.; Yarden, Y. Hsp90 Restrains ErbB2/HER2 Signalling by Limiting Heterodimer Formation. EMBO Rep. 2004, 5, 1165-1170. (108) Ota, A.; Wang, Y. Cdc37/Hsp90 Protein-Mediated Regulation of IRE1alpha Protein Activity in Endoplasmic Reticulum Stress Response and Insulin Synthesis in INS-1 Cells. J. Biol. Chem. 2012, 287, 6266-6274. (109) Nguyen, M. T. N.; Kniess, R. A.; Daturpalli, S.; Le Breton, L.; Ke, X.; Chen, X.; Mayer, M. P., Isoform-Specific Phosphorylation in Human Hsp90beta Affects Interaction with Clients and the Cochaperone Cdc37. J. Mol. Biol. 2017, 429, 732-752. (110) Scroggins, B. T.; Robzyk, K.; Wang, D.; Marcu, M. G.; Tsutsumi, S.; Beebe, K.; Cotter, R. J.; Felts, S.; Toft, D.; Karnitz, L.; Rosen, N.; Neckers, L. An Acetylation Site in the Middle Domain of Hsp90 Regulates Chaperone Function. Mol. Cell 2007, 25, 151-159. (111) He, L.; Sharpe, T.; Mazur, A.; Hiller, S. A Molecular Mechanism of Chaperone-Client Recognition. Sci. Adv. 2016, 2, e1601625.

53 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Functional Role and Hierarchy of the Binding Interfaces in Recruitment of Protein Kinase Clients to the Hsp90-Cdc37 Chaperone Machine: Structure-Based Network Analysis of Allosteric Interactions and Communication Pathways

Gabrielle Stetz, Gennady M. Verkhivker

For Table of Contents Use Only

54 ACS Paragon Plus Environment

Page 54 of 65

Page 55 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

Figure 1. Structural mapping and analysis of the binding interfaces in the Hsp90-Cdc37-Cdk4 complex. The overview of the global topology and intermolecular interfaces in the Hsp90-Cdc37-Cdk4 complex is shown in ribbons (A) and surface representations (B). Hsp90-A monomer is shown in orange ribbons and Hsp90-B monomer B is in pink ribbons. Cdc37 is in red ribbons, and Cdk4 kinase is in yellow ribbons. Each monomer in the full length human Hsp90β protein is divided into three domains: NTD (residues 1-215), MD (residues 216-552), and CTD (residues 553-690). The domain annotation is based on comprehensive sequence analysis of the Hsp90 family.89,90 Cdc37 domains are annotated as follows: Cdc37-NTD (residues 1-147), and Cdc37-M/C (residues 148-260). The Cdk4 N-lobe (residues 1-99) and Cdk4 C-lobe (residues 100-295) annotations are consistent with the original study.39 (C) The Hsp90-Cdc37-NTD intermolecular cluster is formed by phosphorylated S13 (Cdc37-NTD) interacting with Cdc37-NTD residues (H33, R36) and Hsp90K406. Another binding interface is mediated by a network of salt bridges formed by the Cdc37-NTD residues (D14, D15, E16, D17, and E18) and Hsp90-MD residues (K399, K402, R405, and K406). (D) The Hsp90-Cdk4 interface involves the Src-loop region in the Hsp90-MD (residues F341, L343, F344, E345) and Hsp90-CTD residues L611, R612 and M602. The residues are shown in domain-colored spheres. 133x99mm (300 x 300 DPI)

ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 2. FRI profiles of the Hsp90 dimer. Residue-based FRI distributions of the Hsp90 dimer (A), Cdc37 (B) and Cdk4 kinase client (C) in the unbound form (shown in green lines) and in the Hsp90-Cdc37-Cdk4 complex (shown in maroon lines). (D) Structural mapping of the inter-domain cluster of rigid residues from Hsp90-MD (P287, W289, Y356, V360 and Hsp90-NTD (F208). Hsp90-A residues are shown in orange spheres and Hsp90-B monomer residues are in pink spheres. The cluster is linked to catalytic residue R392 and ATP (shown in atom-colored spheres). (E) Structural map of stable hydrophobic residues (shown in pink spheres) forming an interacting cluster in the Hsp90-MD (L369-I370-V381-L401-I404). These hydrophobic residues are proximal to a group of lysine residues K399, K402, R405, K406, and K411 (shown in hot pink spheres) that form specific interactions with Cdc37-NTD residues (shown in red spheres). (F) Structural map of stable residues at the interface of Hsp90-CTD (R338, L343, F344, N346, and M602 shown in pink spheres) and unstructured segment of the kinase N-lobe (F92, V93, E94, V96 shown in yellow spheres). These residues are flexible in the unbound proteins but become rigidified in the chaperone-kinase complex. 133x99mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 56 of 65

Page 57 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

Figure 3. Alanine scanning and protein stability analysis of the Hsp90-Cdc37-Cdk4 complex. Protein stability changes are computed for the Hsp90 dimer (A, B), Cdc37 (C) and Cdk4 (D) using a systematic alanine scanning of the protein residues to alanine and computing the effect of each mutation on protein stability with the FoldX approach. The profiles computed for the unbound Hsp90 dimer, Cdc37 and native Cdk4 kinase structure are shown in brown filled bars. The stability changes for protein residues in the Hsp90Cdc37-Cdk4 complex are shown in domain-colored bars. The Hsp90-NTD residues are shown in red bars, Hsp90-MD residues are in blue bars, and Hsp90-CTD residues are in green bars. The Cdc37-NTD residues are in red bars, Cdc37-M/C residues are in blue bars. The Cdk4-N lobe residues are in red bars, and Cdk4 C-lobe residues are in blue bars. 133x99mm (300 x 300 DPI)

ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 4. Structural mapping of protein stability hotspot clusters. (A) The overview of the intermolecular interfaces and energetic hotspot centers in the Hsp90-Cdc37-Cdk4 complex. Hsp90-A monomer is shown in orange ribbons and Hsp90-B monomer B is in pink ribbons. Cdc37 is in red ribbons, and Cdk4 kinase is in yellow ribbons. The location of three major centers of energetic hotspots is indicated by rectangular boxes colored in red (with a close-up in panel B), blue (a close-up is presented in panel C), and green (a close-up is in panel D). (B) Structural mapping of the intramolecular cluster connecting energetic hotspots near the NTD-MD interface. For clarity only Hsp90-B monomer (colored in pink ribbons) is shown. The high stability clusters I288-F304-L332 and P287-W289-F304-Y356 are shown in domain-colored spheres. These residues are directly linked to L388, I390 and R392 (in sticks) that bridge NTD and MD regions near the catalytic site. Catalytic residue R392 and ATP are shown in atom-colored spheres. (C) A close-up of the protein stability hotspot formed by the hydrophobic cluster (I404-V381-I408-L369-I370-L401). The hotspot residues are shown in domain-colored pink spheres (Hsp90-B monomer of the dimer). This stability hotspot is immediately proximal to the Hsp90 binding residues K399, K402, R405, K406 (shown in pink sticks ) that form a network of salt bridges with Cdc37-NTD residues D14,D15, E16,D17,E18 ( in red sticks). Hsp90-B in pink ribbons, Cdc37 in red ribbons, and Cdk4 is in yellow ribbons. (D) The protein stability hotspot cluster F376-L343-F433-N436 (in pink spheres) anchors the Src-loop residues R338, F341, L343, F344 (in sticks) that make contacts with the Cdk4 N-lobe regions. The Cdc37-M/C domain is shown in red ribbons and Cdk4N lobe is shown in yellow ribbons.

133x99mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 58 of 65

Page 59 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

Figure 5. Residue centrality profiles of the Hsp90 dimer. Residue-based centrality distributions of the Hsp90 monomers in the chaperone-kinase complex (A, B). The client-induced differential changes in the residue centrality of the Hsp90 residues (C, D) are obtained from differences between residue betweenness values in the complex and unbound Hsp90 dimer. The distributions are derived by averaging computations of network parameters over equilibrium ensembles extracted from DMD trajectories. The distributions are shown in domain-colored bars. The Hsp90-NTD residues are shown in red bars, Hsp90-MD residues are in blue bars, and Hsp90-CTD residues are in green bars. 133x99mm (300 x 300 DPI)

ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 6. Residue centrality profiles of Cdc37 and Cdk4 in the Hsp90-Cdc37-Cdk4 complex. Residue-based centrality distributions of Cdc37 (A) and Cdk4 (B) are shown in domain-colored bars. The Cdc37-NTD (residues 1-147) is in blue bars, Cdc37-M/C (residues 148-260) in red bars. The Cdk4-N lobe (residues 199) is in blue bars, and Cdk4 C-lobe (residues 100-295) is in red bars. (C) Structural mapping of high centrality residues in the Hsp90-Cdc37-Cdk4 complex. For clarity, only Hsp90-B (pink ribbons) is shown. Cdc37 in red ribbons, Cdk4 in yellow ribbons. (D) A close-up of high centrality residues at the Hsp90B/Cdc37-NTD interface. The Cdc37-NTD residues (pS13, D14, D15, E16, E18, H20, P21, N22, I23, F29, W31, and R33) are highlighted in red spheres. The Hsp90-B interfacial residues (K399, K402, R405, and K406) and adjacent hydrophobic sites (I370, F376, and Y373) are shown in pink spheres. 133x99mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 60 of 65

Page 61 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

Figure 7. Structural map of major intermolecular communities. (A) An overview of the communities (in colored spheres), Hsp90-A monomer is shown in orange ribbons and Hsp90-B monomer is in pink ribbons. Cdc37 is in red, and Cdk4 is in yellow ribbons. The residues in major communities are shown in spheres colored according to domain annotation. (B) The communities linking Hsp90-B and Cdc37-NTD : K399(Hsp90)-D15(Cdc37)-D18(Cdc37) (in orange spheres), K402(Hsp90)-R405(Hsp90)-E16(Cdc37) (in blue spheres), N407(Hsp90)-V403(Hsp90)-V116(Cdc37) (in salmon spheres), L327(Hsp90)-Q326(Hsp90)L119(Cdc37)-W7(Cdc37)-I400(Hsp90)-V403(Hsp90)-L3095(Hsp90)-P387(Hsp90)-Y4(Cdc37)-Q397(Hsp90)I10(Cdc37) (in hot pink spheres). The key residues in the largest community are annotated. (C) The communities linking Cdc37-NTD and Cdk4 C-lobe: H20(Cdc37)-I23(Cdc37)-L28(Cdc37)-F130(Cdk4)L131(Cdk4)-F159(Cdk4)-I136(Cdk4)-W31(Cdc37)-L128(Cdk4) (in blue spheres), and N22(Cdc37)Q123(Cdk4)-V154(Cdk4)-V148(Cdk4)-L120(Cdk4) (in salmon spheres). (D) A single largest community at the Cdc37-M/C interface with the kinase N-lobe. The contributing kinase residues (Y17, V37, E56, V57, L60, and L74) are shown in yellow spheres, and Cdc37-M/C residues (W168, D169, I241, K242) are shown in red spheres. Note contribution of conserved regulatory kinase sites E56, L60, and L74 from the regulatory αChelix. 133x99mm (300 x 300 DPI)

ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 8. The distribution of edge centrality in the Hsp90-Cdc37-Cdk4 complex. (A)The edge centrality distribution is based on all short inter-residue paths in the Hsp90-Cdc37-Cdk4 complex (red bars). The edge centrality distribution for the ensemble of short pathways connecting the kinase lobes blue bars). The edge centrality distribution for the ensemble of short paths connecting the Cdc37-NTD/Cdk4 C-lobe interfacial residues with the ATP binding site (in green bars). (B) The domain-based distribution of high centrality intercommunity edges (red bars). (C) Structural mapping of the major communities in the Hsp90 dimer (shown in green spheres for Hsp90-A monomer and cyan spheres for Hsp90-B monomer. The positions of important conserved PTM sites in the human Hsp90β are highlighted in red spheres and annotated. Note that these PTM sites can often assume role of the inter-community bridges. 133x99mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 62 of 65

Page 63 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

Figure 9. Structural map of major communication pathways passing through high centrality edges in the Hsp90-Cdc37-Cdk4 complex. (A) A general overview and topography of communication pathways. Residues along these pathways are shown in green and cyan spheres for Hsp90 A and B monomers respectively. The positions of PTM sites in the Hsp90 chaperone are highlighted by red spheres. (B) A close-up of high centrality inter-modular connections is shown in domain-colored spheres. The functional sites of posttranslational modifications in this communication map are highlighted in blue spheres and annotated. 133x99mm (300 x 300 DPI)

ACS Paragon Plus Environment

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 10. Modeling of allosteric communication pathways in the Hsp90-Cdc37-Cdk4 complex. (A) Representative optimal short paths from the ensemble of communication pathways connecting catalytic residue R392 (Hsp90-B NTD) with functional sites in Cdk4 client kinase. (B) The most probable pathways from the ensemble of short inter-residue routes connecting Cdc37-NTD functional sites with Cdk4 N-lobe regions. The most probable paths are shown and the respective residues are shown in domain-colored spheres and annotated.

133x99mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 64 of 65

Page 65 of 65 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

Functional Role and Hierarchy of the Binding Interfaces in Recruitment of Protein Kinase Clients to the Hsp90-Cdc37 Chaperone Machine: Structure-Based Network Analysis of Allosteric Interactions and Communication Pathways

Gabrielle Stetz, Gennady M. Verkhivker

For Table of Contents Use Only

88x35mm (300 x 300 DPI)

ACS Paragon Plus Environment