Hsp90 is a molecular chaperone essential for protein folding and activation in normal homeostasis and stress response. ATP binding and hydrolysis facilitate Hsp90 conformational changes required for client activation. Hsp90 plays an important role in disease states, particularly in cancer, where chaperoning of the mutated and overexpressed oncoproteins is important for function. Recent studies have illuminated mechanisms related to the chaperone function. However, an atomic resolution view of Hsp90 conformational dynamics, determined by the presence of different binding partners, is critical to define communication pathways between remote residues in different domains intimately affecting the chaperone cycle. Here, we present a computational analysis of signal propagation and long-range communication pathways in Hsp90. We carried out molecular dynamics simulations of the full-length Hsp90 dimer, combined with essential dynamics, correlation analysis, and a signal propagation model. All-atom MD simulations with timescales of 70 ns have been performed for complexes with the natural substrates ATP and ADP and for the unliganded dimer. We elucidate the mechanisms of signal propagation and determine “hot spots” involved in interdomain communication pathways from the nucleotide-binding site to the C-terminal domain interface. A comprehensive computational analysis of the Hsp90 communication pathways and dynamics at atomic resolution has revealed the role of the nucleotide in effecting conformational changes, elucidating the mechanisms of signal propagation. Functionally important residues and secondary structure elements emerge as effective mediators of communication between the nucleotide-binding site and the C-terminal interface. Furthermore, we show that specific interdomain signal propagation pathways may be activated as a function of the ligand. Our results support a “conformational selection model” of the Hsp90 mechanism, whereby the protein may exist in a dynamic equilibrium between different conformational states available on the energy landscape and binding of a specific partner can bias the equilibrium toward functionally relevant complexes.

Summary Swe1 (Saccharomyces WEE1), the only “true” tyrosine kinase in budding yeast, is an Hsp90 client protein. Here we show that Swe1Wee1 phosphorylates a conserved tyrosine residue (Y24 in yeast Hsp90 and Y38 in human Hsp90α) in the N-domain of Hsp90. Phosphorylation is cell cycle-associated and modulates the ability of Hsp90 to chaperone a selected clientele, including v-Src and several other kinases. Non-phosphorylatable mutants have normal ATPase activity, support yeast viability, and productively chaperone the Hsp90 client glucocorticoid receptor. Deletion of SWE1 in yeast increases Hsp90 binding to its inhibitor geldanamycin, and pharmacologic inhibition/silencing of Wee1 sensitizes cancer cells to Hsp90 inhibitor-induced apoptosis. These findings demonstrate that Hsp90 chaperoning of distinct client proteins is differentially regulated by specific post-translational modification of a unique subcellular pool of the chaperone, and they provide a novel strategy to increase the cellular potency of Hsp90 inhibitors.

Understanding how local protein modifications, such as binding small-molecule ligands, can trigger and regulate large-scale motions of large protein domains is a major open issue in molecular biology. We address various aspects of this problem by analyzing and comparing atomistic simulations of Hsp90 family representatives for which crystal structures of the full length protein are available: mammalian Grp94, yeast Hsp90 and E.coli HtpG. These chaperones are studied in complex with the natural ligands ATP, ADP and in the Apo state. Common key aspects of their functional dynamics are elucidated with a novel multi-scale comparison of their internal dynamics. Starting from the atomic resolution investigation of internal fluctuations and geometric strain patterns, a novel analysis of domain dynamics is developed. The results reveal that the ligand-dependent structural modulations mostly consist of relative rigid-like movements of a limited number of quasi-rigid domains, shared by the three proteins. Two common primary hinges for such movements are identified. The first hinge, whose functional role has been demonstrated by several experimental approaches, is located at the boundary between the N-terminal and Middle-domains. The second hinge is located at the end of a three-helix bundle in the Middle-domain and unfolds/unpacks going from the ATP- to the ADP-state. This latter site could represent a promising novel druggable allosteric site common to all chaperones.

The study of allosteric functional modulation in dynamic proteins is attracting increasing attention. In particular, the discovery of new allosteric sites may generate novel opportunities and strategies for drug development, overcoming the limits of classical active-site oriented drug design. In this paper, we report on the results of a novel, ab initio, fully computational approach for the discovery of allosteric inhibitors based on the physical characterization of signal propagation mechanisms in proteins and apply it to the important molecular chaperone Hsp90. We first characterize the allosteric "hot spots" involved in interdomain communication pathways from the nucleotide-binding site in the N-domain to the distal C-domain. On this basis, we develop dynamic pharmacophore models to screen drug libraries in the search for small molecules with the functional and conformational properties necessary to bind these "hot spot" allosteric sites. Experimental tests show that the selected moelcules bind the Hsp90 C-domain, exhibit antiproliferative activity in different tumor cell lines, while not affecting proliferation of normal human cells, destabilize Hsp90 client proteins, and disrupt association with several cochaperones known to bind the N- and M-domains of Hsp90. These results prove that the hits alter Hsp90 function by affecting its conformational dynamics and recognition properties through an allosteric mechanism. These findings provide us with new insights on the discovery and development of new allosteric inhibitors that are active on important cellular pathways through computational biology. Though based on the specific case of Hsp90, our approach is general and can readily be extended to other target proteins and pathways.

Most proteins must fold to a well-defined structure with a minimal stability to perform their function. Here we use a simple, molecular dynamics-based, energy decomposition approach to map the principal energetic interactions in a set of proteins representative of different folds. This work involves the all-atom simulation and analysis of the native structures and mutants of five different proteins representative of an all-alpha (yACPB, Protein A), all-beta (SH3), and a mixed alpha/beta fold (Proteins G and L). Given a certain structure, a native sequence and a set of mutants, we show that our model discriminates the ability of a mutation to yield a more or less stable protein, in agreement with experimental data, catching the principal energetic determinants of protein stabilization. Our approach identifies the interaction determinants responsible to define a fold and shows that mutations can either modulate the strength of pair-wise coupling between residues important for folding, or modify the profile of the principal interactions. Furthermore, we address the question of how to evaluate the fitness of a sequence to a given structure by comparing the information contained in the energy map, which recapitulates the chemistry of the sequence, to that contained in the contact map, which recapitulates the fold topology. The results show that the better fit between the energetic properties of the sequence and the fold topology corresponds to a higher stabilization of the protein. We discuss the relevance of these observations to the analysis of protein designability and to the rational evolution of new sequences.

An increasing number of functional studies of proteins have shown that sequence and structural similarities alone may not be sufficient for reliable prediction of their interaction properties. This is particularly true for proteins recognizing specific antibodies, where the prediction of antibody-binding sites, called epitopes, has proven challenging. The antibody-binding properties of an antigen depend on its structure and related dynamics. Aiming to predict the antibody-binding regions of a protein, we investigate a new approach based on the integrated analysis of the dynamical and energetic properties of antigens, to identify nonoptimized, low-intensity energetic interaction networks in the protein structure isolated in solution. The method is based on the idea that recognition sites may correspond to localized regions with low-intensity energetic couplings with the rest of the protein, which allows them to undergo conformational changes, to be recognized by a binding partner, and to tolerate mutations with minimal energetic expense. Upon analyzing the results on isolated proteins and benchmarking against antibody complexes, it is found that the method successfully identifies binding sites located on the protein surface that are accessible to putative binding partners. The combination of dynamics and energetics can thus discriminate between epitopes and other substructures based only on physical properties. We discuss implications for vaccine design.

The absolute values of the one-electron redox potentials of the two quinones (Q(A) and Q(B)) in bacterial photosynthetic reaction centers from Rhodobacter sphaeroides were calculated by evaluating the electrostatic energies from the solution of the linearized Poisson-Boltzmann equation at pH 7.0. The redox potential for Q(A) was calculated to be between -173 and -160 mV, which is close to the lowest measured values that are assumed to refer to nonequilibrated protonation patterns in the redox state Q(A)(-). The redox potential of quinone Q(B) is found to be about 160-220 mV larger for the light-exposed than for the dark-adapted structure. These values of the redox potentials are obtained if Asp-L213 is nearly protonated (probability 0.75-1.0) before and after electron transfer from Q(A) to Q(B), while Glu-L212 is partially protonated (probability 0.6) in the initial state Q(A)(-)Q(B)(0) and fully protonated in the final state Q(A)(0)Q(B)(-). Conversely, if the charge state of the quinones is varied from Q(A)(-)Q(B)(0) to Q(A)(0)Q(B)(-) corresponding to the electron transfer from Q(A) to Q(B), Asp-L213 remains protonated, while Glu-L212 changes its protonation state from 0.15 H(+) to fully protonated. In agreement with results from FTIR spectra, there is proton uptake at Glu-L212 going along with the electron transfer, whereas Asp-L213 does not change its protonation state. However, in our simulations Asp-L213 is considered to be protonated rather than ionized as deduced from FTIR spectra. The calculated redox potential of Q(A) shows little dependence on the charge state of Asp-L213, which is due to a strong coupling with the protonation state of Asp-M17 but increases by 50 mV if Glu-L212 changes from the ionized to the protonated charge state. Both are in agreement with fluorescence measurements observing the decay of SP(+)Q(A)(-) in a wide pH regime. The computed difference in redox potential of Q(B) in the light-exposed and dark-adapted structure was traced back to the hydrogen bond of Q(B) with His-L190 that is lost in the dark-adapted structure and the charge of the non-heme iron atom, which is closer to Q(B) in the light-exposed than in the dark-adapted structure.

BackgroundThe conversion of the cellular prion protein (PrPC) into the infectious form (PrPSc) is the key event in prion induced neurodegenerations. This process is believed to involve a multi-step conformational transition from an α-helical (PrPC) form to a β-sheet-rich (PrPSc) state. In addition to the conformational difference, PrPSc exhibits as covalent signature the sulfoxidation of M213. To investigate whether such modification may play a role in the misfolding process we have studied the impact of methionine oxidation on the dynamics and energetics of the HuPrP(125–229) α-fold.Methodology/Principal FindingsUsing molecular dynamics simulation, essential dynamics, correlated motions and signal propagation analysis, we have found that substitution of the sulfur atom of M213 by a sulfoxide group impacts on the stability of the native state increasing the flexibility of regions preceding the site of the modification and perturbing the network of stabilizing interactions. Together, these changes favor the population of alternative states which maybe essential in the productive pathway of the pathogenic conversion. These changes are also observed when the sulfoxidation is placed at M206 and at both, M206 and M213.Conclusions/SignificanceOur results suggest that the sulfoxidation of Helix-3 methionines might be the switch for triggering the initial α-fold destabilization required for the productive pathogenic conversion.