ATP-induced conformational changes in Hsp70: molecular dynamics and experimental validation of an in silico predicted conformation

Structural insights into GrpEL1-mediated nucleotide and substrate release of human mitochondrial Hsp70

Maintenance of protein homeostasis is necessary for cell viability and depends on a complex network of chaperones and co-chaperones, including the heat-shock protein 70 (Hsp70) system. In human mitochondria, mitochondrial Hsp70 (mortalin) and the nucleotide exchange factor (GrpEL1) work synergistically to stabilize proteins, assemble protein complexes, and facilitate protein import. However, our understanding of the molecular mechanisms guiding these processes is hampered by limited structural information. To elucidate these mechanistic details, we used cryoEM to determine the first structures of full-length human mortalin-GrpEL1 complexes in previously unobserved states. Our structures and molecular dynamics simulations allow us to delineate specific roles for mortalin-GrpEL1 interfaces and to identify steps in GrpEL1-mediated nucleotide and substrate release by mortalin. Subsequent analyses reveal conserved mechanisms across bacteria and mammals and facilitate a complete understanding of sequential nucleotide and substrate release for the Hsp70 chaperone system.

Temperature Matters: Bacterial Response to Temperature Change

Temperature is one of the most important factors in all living organisms for survival. Being a unicellular organism, bacterium requires sensitive sensing and defense mechanisms to tolerate changes in temperature. During a temperature shift, the structure and composition of various cellular molecules including nucleic acids, proteins, and membranes are affected. In addition, numerous genes are induced during heat or cold shocks to overcome the cellular stresses, which are known as heat- and cold-shock proteins. In this review, we describe the cellular phenomena that occur with temperature change and bacterial responses from a molecular perspective, mainly in Escherichia coli.

Cytoplasmic molecular chaperones in Pseudomonas species

Pseudomonas is widespread in various environmental and host niches. To promote rejuvenation, cellular protein homeostasis must be finely tuned in response to diverse stresses, such as extremely high and low temperatures, oxidative stress, and desiccation, which can result in protein homeostasis imbalance. Molecular chaperones function as key components that aid protein folding and prevent protein denaturation. Pseudomonas, an ecologically important bacterial genus, includes human and plant pathogens as well as growth-promoting symbionts and species useful for bioremediation. In this review, we focus on protein quality control systems, particularly molecular chaperones, in ecologically diverse species of Pseudomonas, including the opportunistic human pathogen Pseudomonas aeruginosa, the plant pathogen Pseudomonas syringae, the soil species Pseudomonas putida, and the psychrophilic Pseudomonas antarctica.

The structural changes in the signaling mechanism of bacteriophytochromes in solution revealed by a multiscale computational investigation

Phytochromes are red-light sensing proteins, with important light-regulatory roles in different organisms, which are capturing an increasing interest in bioimaging and optogenetics. Upon absorption of light by the embedded bilin chromophore, they undergo structural changes that extend from the chromophore to the protein and finally drive the biological function. Up to now, the underlying mechanism still has to be characterized fully. Here we investigate the Pfr activated form of a bacterial phytochrome, by combining extensive molecular dynamics simulations with a polarizable QM/MM description of the spectroscopic properties, revealing a large structure relaxation in solution, compared to the crystal structure, both in the chromophore-binding pocket and in the overall structure of the phytochrome. Our results indicate that the final opening of the dimeric structure is preceded by an important internal reorganization of the phytochrome specific (PHY) domain involving a bend of the helical spine connecting the PHY domain with the chromophore-binding domain, opening the way to a new understanding of the activation pathway.

Expression, Rapid Purification and Functional Analysis of DnaK from Rhodococcus ruber

BACKGROUND

Heat shock proteins (HSPs) represent a group of important proteins which are produced by all kinds of organisms especially under stressful conditions. DnaK, an Hsp70 homolog in prokaryotes, has indispensable roles when microbes was confronted with stress conditions. However, few data on DnaK from Rhodococcus sp. were available in the literature. In a previous study, we reported that toluene and phenol stress gave rise to a 29.87-fold and 3.93-fold increase for the expression of DnaK from R. ruber SD3, respectively. Thus, we deduced DnaK was in correlation with the organic solvent tolerance of R. ruber SD3.

OBJECTIVE

To elucidate the role of DnaK in the organic solvent tolerance of R. ruber SD3, expression, purification and functional analysis of Dnak from R. ruber SD3 were performed in the present paper.

METHODS

In this article, DnaK from R. ruber SD3 was heterologously expressed in E. coli BL21(DE3) and purified by affinity chromatography. Functional analysis of DnaK was performed using determination of kinetics, docking, assay of chaperone activity and microbial growth.

RESULTS

The recombinant DnaK was rapidly purified by affinity chromatography with the purification fold of 1.9 and the recovery rate of 57.9%. Km, Vmax and Kcat for Dnak from R. ruber SD3 were 80.8 μM, 58.1 nmol/min and 374.3 S-1, respectively. The recombinant protein formed trimer in vitro, with the calculated molecular weight of 214 kDa. According to In-silico analysis, DnaK interacted with other molecular chaperones and some important proteins in the metabolism. The specific activity of catalase in the presence of recombinant DnaK was 1.85 times or 2.00 times that in the presence of BSA or Tris-HCl buffer after exposure to 54 °C for 1h. E. coli transformant with pET28-dnak showed higher growth than E. coli transformant with pET28 at 43°C and in the presence of phenol, respectively.

CONCLUSION

The biochemical properties and the interaction analysis of DnaK from R. ruber SD3 deepened our understanding of DnaK function. DnaK played an important role in microbial growth when R. ruber was subjected to various stress such as heating and organic solvent.

Large-Scale Conformational Changes and Protein Function: Breaking the in silico Barrier

Large-scale conformational changes are essential to link protein structures with their function at the cell and organism scale, but have been elusive both experimentally and computationally. Over the past few years developments in cryo-electron microscopy and crystallography techniques have started to reveal multiple snapshots of increasingly large and flexible systems, deemed impossible only short time ago. As structural information accumulates, theoretical methods become central to understand how different conformers interconvert to mediate biological function. Here we briefly survey current in silico methods to tackle large conformational changes, reviewing recent examples of cross-validation of experiments and computational predictions, which show how the integration of different scale simulations with biological information is already starting to break the barriers between the in silico, in vitro, and in vivo worlds, shedding new light onto complex biological problems inaccessible so far.

Expression and molecular characterization of stress-responsive genes (hsp70 and Mn-sod) and evaluation of antioxidant enzymes (CAT and GPx) in heavy metal exposed freshwater ciliate, Tetmemena sp

Response of heavy metals namely cadmium (Cd) and copper (Cu) on the expression of stress responsive genes in the fresh water ciliate, Tetmemena sp. (single cell eukaryote) was studied. Stress responsive genes include heat shock protein genes and genes involved in antioxidant defence system. Quantitative real time PCR (qRT-PCR) was employed to evaluate the effects of Cd and Cu on the expression of cytosolic hsp70 and Mn-sod genes. Increase in the expression of these genes was observed after exposure with the heavy metals. The macronuclear cytosolic hsp70 and Mn-sod (SOD2) genes were also sequenced and characterized using various bioinformatics tools. In antioxidant defence system, the superoxide dismutase (SOD) family is a first line antioxidant enzyme group involved in catalysing reactive oxygen species (ROS) to hydrogen peroxide and molecular oxygen. Influence of Cd and Cu on the activity of SOD has already been reported by our group. Therefore, the enzymatic activities of antioxidant enzymes, catalase (CAT) and glutathione peroxidase (GPx) were studied in the presence of Cd and Cu and there was significant increase in activity of these enzymes in concentration dependent manner. This study suggests that cytosolic hsp70, Mn-sod and the antioxidant enzymes such as CAT and GPx can be used as effective molecular biomarkers for heavy metal toxicity and Tetmemena sp. can be used as potential model for understanding the molecular response to heavy metal contamination in aquatic ecosystems.

Discorhabdin N, a South African Natural Compound, for Hsp72 and Hsc70 Allosteric Modulation: Combined Study of Molecular Modeling and Dynamic Residue Network Analysis

The human heat shock proteins (Hsps), predominantly Hsp72 and Hsp90, have been strongly implicated in various critical stages of oncogenesis and progression of human cancers. While drug development has extensively focused on Hsp90 as a potential anticancer target, much less effort has been put against Hsp72. This work investigated the therapeutic potential of Hsp72 and its constitutive isoform, Hsc70, via in silico-based screening against the South African Natural Compounds Database (SANCDB). A comparative modeling approach was used to obtain nearly full-length 3D structures of the closed conformation of Hsp72 and Hsc70 proteins. Molecular docking of SANCDB compounds identified one potential allosteric modulator, Discorhabdin N, binding to the allosteric β substrate binding domain (SBDβ) back pocket, with good binding affinities in both cases. This allosteric region was identified in one of our previous studies. Subsequent all-atom molecular dynamics simulations and free energy calculations exhibited promising protein⁻ligand association characteristics, indicative of strong binding qualities. Further, we utilised dynamic residue network analysis (DRN) to highlight protein regions actively involved in cross-domain communication. Most residues identified agreed with known allosteric signal regulators from literature, and were further investigated for the purpose of deducing meaningful insights into the allosteric modulation properties of Discorhabdin N.

A Local Allosteric Network in Heat Shock Protein 70 (Hsp70) Links Inhibitor Binding to Enzyme Activity and Distal Protein-Protein Interactions

Allosteric inhibitors can be more difficult to optimize without an understanding of how their binding influences the conformational motions of the target. Here, we used an integrated computational and experimental approach to probe the molecular mechanism of an allosteric inhibitor of heat shock protein 70 (Hsp70). The anticancer compound, MKT-077, is known to bind a conserved site in members of the Hsp70 family, which favors the ADP-bound state and interferes with a protein-protein interaction (PPI) at long range. However, the binding site does not overlap with either the nucleotide-binding cleft or the PPI contact surface, so its mechanism is unclear. To this end, we modeled Hsp70's internal dynamics and studied how MKT-077 alters local sampling of its allosteric states. The results pointed to a set of concerted motions between five loops in Hsp70's nucleotide-binding domain (NBD), surrounding the MKT-077 binding site. To test this prediction, we mutated key residues and monitored chaperone activities in vitro. Together, the results indicate that MKT-077 interacts with loop222 to favor a pseudo-ADP bound conformer of Hsp70's NBD, even when ATP is present. We used this knowledge to synthesize an analog of MKT-077 that would better prevent motions of loop222 and confirmed that it had improved antiproliferative activity in breast cancer cells. These results provide an example of how to unlock and leverage the complex mechanisms of allosteric inhibitors.

Using steered molecular dynamics to study the interaction between ADP and the nucleotide-binding domain of yeast Hsp70 protein Ssa1

Genetics experiments have identified six mutations located in the subdomain IA (A17V, R23H, G32D, G32S, R34K, V372I) of Ssa1 that influence propagation of the yeast [PSI⁺] prion. However, the underlining molecular mechanisms of these mutations are still unclear. The six mutation sites are present in the IA subdomain of the nucleotide-binding domain (NBD). The ATPase subdomain IA is a critical mediator of inter-domain allostery in Hsp70 molecular chaperones, so the mutation and changes in this subdomain may influence the function of the substrate-binding domain. In addition, ADP release is a rate-limiting step of the ATPase cycle and dysregulation of the ATPase cycle influences the propagation of the yeast [PSI⁺] prion. In this work, steered molecular dynamics (SMD) simulations were performed to explore the interaction between ADP and NBD. Results suggest that during the SMD simulations, hydrophobic interactions are predominant and variations in the binding state of ADP within the mutants is a potential reason for in vivo effects on yeast [PSI⁺] prion propagation. Additionally, we identify the primary residues in the ATPase domain that directly constitute the main hydrophobic interaction network and directly influence the ADP interaction state with the NBD of Ssa1. Furthermore, this in silico analysis reaffirms the importance of previously experimentally-determined residues in the Hsp70 ATPase domain involved in ADP binding and also identifies new residues potentially involved in this process.

Disrupted Hydrogen-Bond Network and Impaired ATPase Activity in an Hsc70 Cysteine Mutant

The ATPase domain of members of the 70 kDa heat shock protein (Hsp70) family shows a high degree of sequence, structural, and functional homology across species. A broadly conserved residue within the Hsp70 ATPase domain that captured our attention is an unpaired cysteine, positioned proximal to the site of nucleotide binding. Prior studies of several Hsp70 family members show this cysteine is not required for Hsp70 ATPase activity, yet select amino acid replacements of the cysteine can dramatically alter ATP hydrolysis. Moreover, post-translational modification of the cysteine has been reported to limit ATP hydrolysis for several Hsp70s. To better understand the underlying mechanism for how perturbation of this noncatalytic residue modulates Hsp70 function, we determined the structure for a cysteine-to-tryptophan mutation in the constitutively expressed, mammalian Hsp70 family member Hsc70. Our work reveals that the steric hindrance produced by a cysteine-to-tryptophan mutation disrupts the hydrogen-bond network within the active site, resulting in a loss of proper catalytic magnesium coordination. We propose that a similarly altered active site is likely observed upon post-translational oxidation. We speculate that the subtle changes we detect in the hydrogen-bonding network may relate to the previously reported observation that cysteine oxidation can influence Hsp70 interdomain communication.

Molecular cloning, prokaryotic expression, purification, structural studies and functional implications of Heat Shock Protein 70 (Hsp70) from Rutilus frisii kutum

A novel Hsp70 chaperone from Rutilus frisii kutum was identified, cloned, expressed, purified and its functional characteristics revealed. The 3D structure of Hsp70 from Rutilus kutum was constructed using the crystal structure of E. coli Hsp70 as the template, with 47% sequence identity. The in vitro ATPase activity assay after 60min, ATP hydrolysis of purified recombinant Hsp70 (8μM) was improved by binding to denatured thermally luciferase (3μM) about 2.5-fold compared with that of Hsp70 alone. Based on the results, it was found that the purified Hsp70 chaperone was able to considerably suppress heat-induced aggregation of luciferase by binding to DnaJ co-chaperone (5μM) more than 70% after 10min at 42°C. In addition, Hsp70 DnaJ complex improved the refolding of heat-shocked luciferase nearly 40% after 60min at 25°C. It was concluded that Hsp70 protein from Rutilus frisii kutum has the critical role in preventing heat-induced aggregation of luciferase and refolding of heat-denatured luciferase was strictly dependent on the activity of Hsp70, thus, this protein can potentially be used for improving the functional properties of luciferase in various applications.

Role and mechanism of the Hsp70 molecular chaperone machines in bacterial pathogens

Heat shock proteins are highly conserved, stress-inducible, ubiquitous proteins that maintain homeostasis in both eukaryotes and prokaryotes. Hsp70 proteins belong to the heat shock protein family and enhance bacterial survival in hostile environments. Hsp70, known as DnaK in prokaryotes, supports numerous processes such as the assembly and disassembly of protein complexes, the refolding of misfolded and clustered proteins, membrane translocation and the regulation of regulatory proteins. The chaperone-based activity of Hsp70 depends on dynamic interactions between its two domains, known as the ATPase domain and the substrate-binding domain. It also depends on interactions between these domains and other co-chaperone molecules such as the Hsp40 protein family member DnaJ and nucleotide exchange factors. DnaJ is the primary chaperone that interacts with nascent polypeptide chains and functions to prevent their premature release from the ribosome and misfolding before it is targeted by DnaK. Adhesion of bacteria to host cells is mediated by both host and bacterial Hsp70. Following infection of the host, bacterial Hsp70 (DnaK) is in a position to initiate bacterial survival processes and trigger an immune response by the host. Any mutations in the dnaK gene have been shown to decrease the viability of bacteria inside the host. This review will give insights into the structure and mechanism of Hsp70 and its role in regulating the protein activity that contributes to pathogenesis.

Structural insights into a unique Hsp70-Hsp40 interaction in the eukaryotic ribosome-associated complex

Cotranslational chaperones assist de novo folding of nascent polypeptides, prevent them from aggregating and modulate translation. The ribosome-associated complex (RAC) is unique in that the Hsp40 protein Zuo1 and the atypical Hsp70 chaperone Ssz1 form a stable heterodimer, which acts as a cochaperone for the Hsp70 chaperone Ssb. Here we present the structure of the Chaetomium thermophilum RAC core comprising Ssz1 and the Zuo1 N terminus. We show how the conserved allostery of Hsp70 proteins is abolished and this Hsp70-Hsp40 pair is molded into a functional unit. Zuo1 stabilizes Ssz1 in trans through interactions that in canonical Hsp70s occur in cis. Ssz1 is catalytically inert and cannot adopt the closed conformation, but the substrate binding domain β is completed by Zuo1. Our study offers insights into the coupling of a special Hsp70-Hsp40 pair, which evolved to link protein folding and translation.

Dynamical Structures of Hsp70 and Hsp70-Hsp40 Complexes

Protein misfolding and aggregation are pathological events that place a significant amount of stress on the maintenance of protein homeostasis (proteostasis). For prevention and repair of protein misfolding and aggregation, cells are equipped with robust mechanisms that mainly rely on molecular chaperones. Two classes of molecular chaperones, heat shock protein 70 kDa (Hsp70) and Hsp40, recognize and bind to misfolded proteins, preventing their toxic biomolecular aggregation and enabling refolding or targeted degradation. Here, we review the current state of structural biology of Hsp70 and Hsp40-Hsp70 complexes and examine the link between their structures, dynamics, and functions. We highlight the power of nuclear magnetic resonance spectroscopy to untangle complex relationships behind molecular chaperones and their mechanism(s) of action.

Human Inducible Hsp70: Structures, Dynamics, and Interdomain Communication from All-Atom Molecular Dynamics Simulations

The 70 kDa human heat shock protein is a major molecular chaperone involved in de novo folding of proteins in vivo and refolding of proteins under stress conditions. Hsp70 is related to several "misfolding diseases" and other major pathologies, such as cancer, and is a target for new therapies. Hsp70 is comprised of two main domains: an N-terminal nucleotide binding domain (NBD) and a C-terminal substrate protein binding domain (SBD). The chaperone function of Hsp70 is based on an allosteric mechanism. Binding of ATP in NBD decreases the affinity of the substrate for SBD, and hydrolysis of ATP is promoted by binding of polypeptide segments in the SBD. No complete structure of human Hsp70 is known. Here, we report two models of human Hsp70, constructed by homology with Saccharomyces cerevisiae cochaperone protein Hsp110 (open model) and with Escherichia coli 70 kDa DnaK (closed model) and relaxed for several tens to hundreds of nanoseconds by using all-atom molecular dynamics simulations in explicit solvent. We obtain two stable states, Hsp70 with SBD open and SBD closed, which agree with experimental and structural information for ATP-Hsp70 and ADP-Hsp70, respectively. The dynamics of the transition from the open to closed states is investigated with a coarse-grained model and normal-mode analysis. The results show that the conformational change between the two states can be represented by a relatively small number of collective modes which involved major conformational changes in the two domains. These modes provide a mechanistic representation of the communication between NBD and SBD and allow us to identify subdomains and residues that appear to have a critical role in the conformational change mechanism that guides the chaperoning cycle of Hsp70.

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

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

BiP and its nucleotide exchange factors Grp170 and Sil1: mechanisms of action and biological functions

BiP (immunoglobulin heavy-chain binding protein) is the endoplasmic reticulum (ER) orthologue of the Hsp70 family of molecular chaperones and is intricately involved in most functions of this organelle through its interactions with a variety of substrates and regulatory proteins. Like all Hsp70 family members, the ability of BiP to bind and release unfolded proteins is tightly regulated by a cycle of ATP binding, hydrolysis, and nucleotide exchange. As a characteristic of the Hsp70 family, multiple DnaJ-like co-factors can target substrates to BiP and stimulate its ATPase activity to stabilize the binding of BiP to substrates. However, only in the past decade have nucleotide exchange factors for BiP been identified, which has shed light not only on the mechanism of BiP-assisted folding in the ER but also on Hsp70 family members that reside throughout the cell. We will review the current understanding of the ATPase cycle of BiP in the unique environment of the ER and how it is regulated by the nucleotide exchange factors, Grp170 (glucose-regulated protein of 170kDa) and Sil1, both of which perform unanticipated roles in various biological functions and disease states.

Common functionally important motions of the nucleotide-binding domain of Hsp70

The 70 kDa heat shock proteins (Hsp70) are a family of molecular chaperones involved in protein folding, aggregate prevention, and protein disaggregation. They consist of the substrate-binding domain (SBD) that binds client substrates, and the nucleotide-binding domain (NBD), whose cycles of nucleotide hydrolysis and exchange underpin the activity of the chaperone. To characterize the structure-function relationships that link the binding state of the NBD to its conformational behavior, we analyzed the dynamics of the NBD of the Hsp70 chaperone from Bos taurus (PDB 3C7N:B) by all-atom canonical molecular dynamics simulations. It was found that essential motions within the NBD fall into three major classes: the mutual class, reflecting tendencies common to all binding states, and the ADP- and ATP-unique classes, which reflect conformational trends that are unique to either the ADP- or ATP-bound states, respectively. "Mutual" class motions generally describe "in-plane" and/or "out-of-plane" (scissor-like) rotation of the subdomains within the NBD. This result is consistent with experimental nuclear magnetic resonance data on the NBD. The "unique" class motions target specific regions on the NBD, usually surface loops or sites involved in nucleotide binding and are, therefore, expected to be involved in allostery and signal transmission. For all classes, and especially for those of the "unique" type, regions of enhanced mobility can be identified; these are termed "hot spots," and their locations generally parallel those found by NMR spectroscopy. The presence of magnesium and potassium cations in the nucleotide-binding pocket was also found to influence the dynamics of the NBD significantly.

Structural mechanisms of chaperone mediated protein disaggregation

The ClpB/Hsp104 and Hsp70 classes of molecular chaperones use ATP hydrolysis to dissociate protein aggregates and complexes, and to move proteins through membranes. ClpB/Hsp104 are members of the AAA+ family of proteins which form ring-shaped hexamers. Loops lining the pore in the ring engage substrate proteins as extended polypeptides. Interdomain rotations and conformational changes in these loops coupled to ATP hydrolysis unfold and pull proteins through the pore. This provides a mechanism that progressively disrupts local secondary and tertiary structure in substrates, allowing these chaperones to dissociate stable aggregates such as β-sheet rich prions or coiled coil SNARE complexes. While the ClpB/Hsp104 mechanism appears to embody a true power-stroke in which an ATP powered conformational change in one protein is directly coupled to movement or structural change in another, the mechanism of force generation by Hsp70s is distinct and less well understood. Both active power-stroke and purely passive mechanisms in which Hsp70 captures spontaneous fluctuations in a substrate have been proposed, while a third proposed mechanism-entropic pulling-may be able to generate forces larger than seen in ATP-driven molecular motors without the conformational coupling required for a power-stroke. The disaggregase activity of these chaperones is required for thermotolerance, but unrestrained protein complex/aggregate dissociation is potentially detrimental. Disaggregating chaperones are strongly auto-repressed, and are regulated by co-chaperones which recruit them to protein substrates and activate the disaggregases via mechanisms involving either sequential transfer of substrate from one chaperone to another and/or simultaneous interaction of substrate with multiple chaperones. By effectively subjecting substrates to multiple levels of selection by multiple chaperones, this may insure that these potent disaggregases are only activated in the appropriate context.

Prediction of allosteric sites on protein surfaces with an elastic-network-model-based thermodynamic method

Allostery is a rapid and efficient way in many biological processes to regulate protein functions, where binding of an effector at the allosteric site alters the activity and function at a distant active site. Allosteric regulation of protein biological functions provides a promising strategy for novel drug design. However, how to effectively identify the allosteric sites remains one of the major challenges for allosteric drug design. In the present work, a thermodynamic method based on the elastic network model was proposed to predict the allosteric sites on the protein surface. In our method, the thermodynamic coupling between the allosteric and active sites was considered, and then the allosteric sites were identified as those where the binding of an effector molecule induces a large change in the binding free energy of the protein with its ligand. Using the proposed method, two proteins, i.e., the 70 kD heat shock protein (Hsp70) and GluA2 alpha-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid (AMPA) receptor, were studied and the allosteric sites on the protein surface were successfully identified. The predicted results are consistent with the available experimental data, which indicates that our method is a simple yet effective approach for the identification of allosteric sites on proteins.

ATPase subdomain IA is a mediator of interdomain allostery in Hsp70 molecular chaperones

The versatile functions of the heat shock protein 70 (Hsp70) family of molecular chaperones rely on allosteric interactions between their nucleotide-binding and substrate-binding domains, NBD and SBD. Understanding the mechanism of interdomain allostery is essential to rational design of Hsp70 modulators. Yet, despite significant progress in recent years, how the two Hsp70 domains regulate each other's activity remains elusive. Covariance data from experiments and computations emerged in recent years as valuable sources of information towards gaining insights into the molecular events that mediate allostery. In the present study, conservation and covariance properties derived from both sequence and structural dynamics data are integrated with results from Perturbation Response Scanning and in vivo functional assays, so as to establish the dynamical basis of interdomain signal transduction in Hsp70s. Our study highlights the critical roles of SBD residues D481 and T417 in mediating the coupled motions of the two domains, as well as that of G506 in enabling the movements of the α-helical lid with respect to the β-sandwich. It also draws attention to the distinctive role of the NBD subdomains: Subdomain IA acts as a key mediator of signal transduction between the ATP- and substrate-binding sites, this function being achieved by a cascade of interactions predominantly involving conserved residues such as V139, D148, R167 and K155. Subdomain IIA, on the other hand, is distinguished by strong coevolutionary signals (with the SBD) exhibited by a series of residues (D211, E217, L219, T383) implicated in DnaJ recognition. The occurrence of coevolving residues at the DnaJ recognition region parallels the behavior recently observed at the nucleotide-exchange-factor recognition region of subdomain IIB. These findings suggest that Hsp70 tends to adapt to co-chaperone recognition and activity via coevolving residues, whereas interdomain allostery, critical to chaperoning, is robustly enabled by conserved interactions.

Fungal invasion of epithelial cells

Interaction between host cells and invasive Candida plays a large role in the pathogenicity of Candida species. Fungal-induced endocytosis and active penetration are the two distinct, yet complementary invasion mechanisms of invasive candidiasis. Induced endocytosis is a microorganism-triggered, epithelial-driven, clathrin-mediated and actin-dependent process. During the fundamental pathological process of induced endocytosis, invasins (Als3 and Ssa1), which mediate the binding of host epithelial surface proteins, are expressed by Candida species on the hyphal surface. Sequentially, the interaction between invasins and host epithelial surface proteins stimulates the recruitment of clathrin, dynamin and cortactin to the sites where Candida enters epithelial cells, which in turn induce the actin cytoskeleton reorganization. Actin cytoskeleton provides the force required for fungal internalization. Parallely, active penetration of Candida can directly pass through epithelial cells possibly due to progressive elongation of hyphae and physical forces. Several molecules, such as secreted hydrolases and Als3, can affect the protective barrier of the epithelium and make Candida actively penetrate into epithelial cells through intercellular gaps of epithelial layers.

Heat shock protein 70 inhibitors. 1. 2,5'-thiodipyrimidine and 5-(phenylthio)pyrimidine acrylamides as irreversible binders to an allosteric site on heat shock protein 70

Heat shock protein 70 (Hsp70) is an important emerging cancer target whose inhibition may affect multiple cancer-associated signaling pathways and, moreover, result in significant cancer cell apoptosis. Despite considerable interest from both academia and pharmaceutical companies in the discovery and development of druglike Hsp70 inhibitors, little success has been reported so far. Here we describe structure–activity relationship studies in the first rationally designed Hsp70 inhibitor class that binds to a novel allosteric pocket located in the N-terminal domain of the protein. These 2,5′-thiodipyrimidine and 5-(phenylthio)pyrimidine acrylamides take advantage of an active cysteine embedded in the allosteric pocket to act as covalent protein modifiers upon binding. The study identifies derivatives 17a and 20a, which selectively bind to Hsp70 in cancer cells. Addition of high nanomolar to low micromolar concentrations of these inhibitors to cancer cells leads to a reduction in the steady-state levels of Hsp70-sheltered oncoproteins, an effect associated with inhibition of cancer cell growth and apoptosis. In summary, the described scaffolds represent a viable starting point for the development of druglike Hsp70 inhibitors as novel anticancer therapeutics.

Identification of key hinge residues important for nucleotide-dependent allostery in E. coli Hsp70/DnaK

DnaK is a molecular chaperone that has important roles in protein folding. The hydrolysis of ATP is essential to this activity, and the effects of nucleotides on the structure and function of DnaK have been extensively studied. However, the key residues that govern the conformational motions that define the apo, ATP-bound, and ADP-bound states are not entirely clear. Here, we used molecular dynamics simulations, mutagenesis, and enzymatic assays to explore the molecular basis of this process. Simulations of DnaK's nucleotide-binding domain (NBD) in the apo, ATP-bound, and ADP/P_i-bound states suggested that each state has a distinct conformation, consistent with available biochemical and structural information. The simulations further suggested that large shearing motions between subdomains I-A and II-A dominated the conversion between these conformations. We found that several evolutionally conserved residues, especially G228 and G229, appeared to function as a hinge for these motions, because they predominantly populated two distinct states depending on whether ATP or ADP/P_i was bound. Consistent with the importance of these “hinge” residues, alanine point mutations caused DnaK to have reduced chaperone activities in vitro and in vivo. Together, these results clarify how sub-domain motions communicate allostery in DnaK.

DnaK belongs to the highly conserved heat shock protein 70 (Hsp70) family, a group of ATP-dependent molecular chaperones that regulates proteostasis. Studies have suggested that global movements of the subdomains in the nucleotide-binding domain (NBD) of DnaK regulate its catalytic activity. However, there is less known about the key residues involved in these subdomain motions and whether these residues might also regulate inter-domain allostery with the substrate-binding domain (SBD). To examine the motions in the NBD, dynamics simulations of DnaK's NBD in the apo, ATP-bound, and ADP/P_i-bound states were performed. Through essential dynamics and torsion angle analyses, we identified motions and highly conserved hinge residues between subdomains IIA and IIB that are likely to be important for nucleotide cycling and for communicating the nucleotide state to the SBD. Supporting this model, mutating these conserved hinge residues affected ATPase activity and chaperone functions in vitro and in bacteria, suggesting their importance in the nucleotide-dependent motions in DnaK.

Identification of an allosteric pocket on human hsp70 reveals a mode of inhibition of this therapeutically important protein

Hsp70s are important cancer chaperones that act upstream of Hsp90 and exhibit independent anti-apoptotic activities. To develop chemical tools for the study of human Hsp70, we developed a homology model that unveils a previously unknown allosteric site located in the nucleotide binding domain of Hsp70. Combining structure-based design and phenotypic testing, we discovered a previously unknown inhibitor of this site, YK5. In cancer cells, this compound is a potent and selective binder of the cytosolic but not the organellar human Hsp70s and has biological activity partly by interfering with the formation of active oncogenic Hsp70/Hsp90/client protein complexes. YK5 is a small molecule inhibitor rationally designed to interact with an allosteric pocket of Hsp70 and represents a previously unknown chemical tool to investigate cellular mechanisms associated with Hsp70.

Tracking the interplay between bound peptide and the lid domain of DnaK, using molecular dynamics

Hsp70 chaperones consist of two functional domains: the 44 kDa Nucleotide Binding Domain (NBD), that binds and hydrolyses ATP, and the 26 kDa Substrate Binding Domain (SBD), which binds unfolded proteins and reactivates them, utilizing energy obtained from nucleotide hydrolysis. The structure of the SBD of the bacterial Hsp70, DnaK, consists of two sub-domains: A β-sandwich part containing the hydrophobic cavity to which the hepta-peptide NRLLLTG (NR) is bound, and a segment made of 5 α-helices, called the “lid” that caps the top of the β-sandwich domain. In the present study we used the Escherichia coli Hsp70, DnaK, as a model for Hsp70 proteins, focusing on its SBD domain, examining the changes in the lid conformation. We deliberately decoupled the NBD from the SBD, limiting the study to the structure of the SBD section, with an emphasis on the interaction between the charges of the peptide with the residues located in the lid. Molecular dynamics simulations of the complex revealed significant mobility within the lid structure; as the structure was released from the forces operating during the crystallization process, the two terminal helices established a contact with the positive charge at the tip of the peptide. This contact is manifested only in the presence of electrostatic attraction. The observed internal motions within the lid provide a molecular role for the function of this sub-domain during the reaction cycle of Hsp 70 chaperones.

Real-time observation of the conformational dynamics of mitochondrial Hsp70 by spFRET

The numerous functions of the important class of molecular chaperones, heat shock proteins 70 (Hsp70), rely on cycles of intricate conformational changes driven by ATP-hydrolysis and regulated by cochaperones and substrates. Here, we used Förster resonance energy transfer to study the conformational dynamics of individual molecules of Ssc1, a mitochondrial Hsp70, in real time. The intrinsic dynamics of the substrate-binding domain of Ssc1 was observed to be uncoupled from the dynamic interactions between substrate- and nucleotide-binding domains. Analysis of the fluctuations in the interdomain separation revealed frequent transitions to a nucleotide-free state. The nucleotide-exchange factor Mge1 did not induce ADP release, as expected, but rather facilitated binding of ATP. These results indicate that the conformational cycle of Ssc1 is more elaborate than previously thought and provide insight into how the Hsp70s can perform a wide variety of functions.

Molecular mechanism of allosteric communication in Hsp70 revealed by molecular dynamics simulations

Investigating ligand-regulated allosteric coupling between protein domains is fundamental to understand cell-life regulation. The Hsp70 family of chaperones represents an example of proteins in which ATP binding and hydrolysis at the Nucleotide Binding Domain (NBD) modulate substrate recognition at the Substrate Binding Domain (SBD). Herein, a comparative analysis of an allosteric (Hsp70-DnaK) and a non-allosteric structural homolog (Hsp110-Sse1) of the Hsp70 family is carried out through molecular dynamics simulations, starting from different conformations and ligand-states. Analysis of ligand-dependent modulation of internal fluctuations and local deformation patterns highlights the structural and dynamical changes occurring at residue level upon ATP-ADP exchange, which are connected to the conformational transition between closed and open structures. By identifying the dynamically responsive protein regions and specific cross-domain hydrogen-bonding patterns that differentiate Hsp70 from Hsp110 as a function of the nucleotide, we propose a molecular mechanism for the allosteric signal propagation of the ATP-encoded conformational signal.

Allostery, or the capability of proteins to respond to ligand binding events with a variation in structure or dynamics at a distant site, is a common feature for biomolecular function and regulation in a large number of proteins. Intra-protein connections and inter-residue coordinations underlie allosteric mechanisms and react to binding primarily through a finely tuned modulation of motions and structures at the microscopic scale. Hence, all-atom molecular dynamics simulations are suitable to investigate the molecular basis of allostery. Moreover, understanding intra-protein communication pathways at atomistic resolutions offers unique opportunities in rational drug design. Proteins of the Hsp70 family are allosteric molecular chaperones involved in maintaining cellular protein homeostasis. These proteins are involved in several types of cancer, neurodegenerative diseases, aging and infections and are therefore pharmaceutically relevant targets. In this work we have analyzed, by multiple molecular dynamics simulations, the long-range dynamical and conformational effects of ligands bound to Hsp70, and found relevant differences in comparison to the known non-allosteric structural homolog Hsp110. The resulting model of the mechanism of allosteric propagation offers the opportunity of identifying on-pathway allosteric druggable sites, which we propose could guide rational drug-design efforts targeting Hsp70.

Simulation of the opening and closing of Hsp70 chaperones by coarse-grained molecular dynamics

Heat-shock proteins 70 (Hsp70s) are key molecular chaperones which assist in the folding and refolding/disaggregation of proteins. Hsp70s, which consist of a nucleotide-binding domain (NBD, consisting of NBD-I and NBD-II subdomains) and a substrate-binding domain [SBD, further split into the β-sheet (SBD-β) and α-helical (SBD-α) subdomains], occur in two major conformations having (a) a closed SBD, in which the SBD and NBD domains do not interact, (b) an open SBD, in which SBD-α interacts with NBD-I and SBD-β interacts with the top parts of NBD-I and NBD-II. In the SBD-closed conformation, SBD is bound to a substrate protein, with release occurring after transition to the open conformation. While the transition from the closed to the open conformation is triggered efficiently by binding of adenosine triphosphate (ATP) to the NBD, it also occurs, although less frequently, in the absence of ATP. The reverse transition occurs after ATP hydrolysis. Here, we report canonical and multiplexed replica exchange simulations of the conformational dynamics of Hsp70s using a coarse-grained molecular dynamics approach with the UNRES force field. The simulations were run in the following three modes: (i) with the two halves of the NBD unrestrained relative to each other, (ii) with the two halves of the NBD restrained in an "open" geometry as in the SBD-closed form of DnaK (2KHO), and (iii) the two halves of NBD restrained in a "closed" geometry as in known experimental structures of ATP-bound NBD forms of Hsp70. Open conformations, in which the SBD interacted strongly with the NBD, formed spontaneously during all simulations; the number of transitions was largest in simulations carried out with the "closed" NBD domain, and smallest in those carried out with the "open" NBD domain; this observation is in agreement with the experimentally-observed influence of ATP-binding on the transition of Hsp70's from the SBD-closed to the SBD-open form. Two kinds of open conformations were observed: one in which SBD-α interacts with NBD-I and SBD-β interacts with the top parts of NBD-I and NBD-II (as observed in the structures of nucleotide exchange factors), and another one in which this interaction pattern is swapped. A third type of motion, in which SBD-α binds to NBD without dissociating from SBD-β was also observed. It was found that the first stage of interdomain communication (approach of SBD-β, to NBD) is coupled with the rotation of the long axes of NBD-I and NBD-II towards each other. To the best of our knowledge, this is the first successful simulation of the full transition of an Hsp70 from the SBD-closed to the SBD-open conformation.

Chemical screens against a reconstituted multiprotein complex: myricetin blocks DnaJ regulation of DnaK through an allosteric mechanism

DnaK is a molecular chaperone responsible for multiple aspects of bacterial proteostasis. The intrinsically slow ATPase activity of DnaK is stimulated by its co-chaperone, DnaJ, and these proteins often work in concert. To identify inhibitors we screened plant-derived extracts against a reconstituted mixture of DnaK and DnaJ. This approach resulted in the identification of flavonoids, including myricetin, which inhibited activity by up to 75%. Interestingly, myricetin prevented DnaJ-mediated stimulation of ATPase activity, with minimal impact on either DnaK's intrinsic turnover rate or its stimulation by another co-chaperone, GrpE. Using NMR, we found that myricetin binds DnaK at an unanticipated site between the IB and IIB subdomains and that it allosterically blocked binding of DnaK to DnaJ. Together, these results highlight a "gray box" screening approach, which might facilitate the identification of inhibitors of other protein-protein interactions.

Allosteric signal transmission in the nucleotide-binding domain of 70-kDa heat shock protein (Hsp70) molecular chaperones

The 70-kDa heat shock protein (Hsp70) chaperones perform a wide array of cellular functions that all derive from the ability of their N-terminal nucleotide-binding domains (NBDs) to allosterically regulate the substrate affinity of their C-terminal substrate-binding domains in a nucleotide-dependent mechanism. To explore the structural origins of Hsp70 allostery, we performed NMR analysis on the NBD of DnaK, the Escherichia coli Hsp70, in six different states (ligand-bound or apo) and in two constructs, one that retains the conserved and functionally crucial portion of the interdomain linker (residues ) and another that lacks the linker. Chemical-shift perturbation patterns identify residues at subdomain interfaces that constitute allosteric networks and enable the NBD to act as a nucleotide-modulated switch. Nucleotide binding results in changes in subdomain orientations and long-range perturbations along subdomain interfaces. In particular, our findings provide structural details for a key mechanism of Hsp70 allostery, by which information is conveyed from the nucleotide-binding site to the interdomain linker. In the presence of ATP, the linker binds to the edge of the IIA β-sheet, which structurally connects the linker and the nucleotide-binding site. Thus, a pathway of allosteric communication leads from the NBD nucleotide-binding site to the substrate-binding domain via the interdomain linker.

Heat, pH induced aggregation and surface hydrophobicity of S. cerevesiae Ssa1 protein

Heat shock protein 70 is a conserved protein among organisms. Hsp70 helps substrate proteins to fold correctly. Unfolded substrate proteins increase the probability of the aggregate formation. High level recombinant protein expression in biotechnology often leads insoluble inclusion bodies. To prevent aggregation and to obtain high levels of soluble proteins, Hsp co-expression with desired recombinant protein in yeast becomes a popular method. For this purpose, S. cerevesiae cytosolic Hsp70 (Ssa1) biochemical properties were characterized. Alteration of Ssa1 structure between ATP- and ADP-bound states regulates its function. Therefore, conformation-dependent Ssa1 hydrophobicity and as a result aggregation may also play a key role in Ssa1 function. Therefore, a combination of FTIR, acrylamide quenching, and ANS was used to investigate the effect of nucleotide binding on the structure of Ssa1. Ssa1 secondary structure alterations and hydrophobic properties in aqueous solutions with differing ionic strengths and temperature were also studied.

Substrate discrimination of the chaperone BiP by autonomous and cochaperone-regulated conformational transitions

The endoplasmic reticulum is the site of folding, assembly and quality control for proteins of the secretory pathway. The ATP-regulated Hsp70 chaperone BiP (heavy chain-binding protein), together with cochaperones, has important roles in all of these processes. The functional cycle of Hsp70s is determined by conformational transitions that are required for substrate binding and release. Here, we used the intrinsically disordered C(H)1 domain of antibodies as an authentic substrate protein and analyzed the conformational cycle of BiP by single-molecule and ensemble Förster resonance energy transfer (FRET) measurements. Nucleotide binding resulted in concerted domain movements of BiP. Conformational transitions of the lid domain allowed BiP to discriminate between peptide and protein substrates. A major BiP cochaperone in antibody folding, ERdj3, modulated the conformational space of BiP in a nucleotide-dependent manner, placing the lid subdomain in an open, protein-accepting state.

The conformational dynamics of the mitochondrial Hsp70 chaperone

Heat shock proteins 70 (Hsp70) represent a ubiquitous and conserved family of molecular chaperones involved in a plethora of cellular processes. The dynamics of their ATP hydrolysis-driven and cochaperone-regulated conformational cycle are poorly understood. We used fluorescence spectroscopy to analyze, in real time and at single-molecule resolution, the effects of nucleotides and cochaperones on the conformation of Ssc1, a mitochondrial member of the family. We report that the conformation of its ADP state is unexpectedly heterogeneous, in contrast to a uniform ATP state. Substrates are actively involved in determining the conformation of Ssc1. The J protein Mdj1 does not interact transiently with the chaperone, as generally believed, but rather is released slowly upon ATP hydrolysis. Analysis of the major bacterial Hsp70 revealed important differences between highly homologous members of the family, possibly explaining tuning of Hsp70 chaperones to meet specific functions in different organisms and cellular compartments.

Mutagenesis reveals the complex relationships between ATPase rate and the chaperone activities of Escherichia coli heat shock protein 70 (Hsp70/DnaK)

The Escherichia coli 70-kDa heat shock protein, DnaK, is a molecular chaperone that engages in a variety of cellular activities, including the folding of proteins. During this process, DnaK binds its substrates in coordination with a catalytic ATPase cycle. Both the ATPase and protein folding activities of DnaK are stimulated by its co-chaperones, DnaJ and GrpE. However, it is not yet clear how changes in the stimulated ATPase rate of DnaK impact the folding process. In this study, we performed mutagenesis throughout the nucleotide-binding domain of DnaK to generate a collection of mutants in which the stimulated ATPase rates varied from 0.7 to 13.6 pmol/microg/min(-1). We found that this range was largely established by differences in the ability of the mutants to be stimulated by one or both of the co-chaperones. Next, we explored how changes in ATPase rate might impact refolding of denatured luciferase in vitro and found that the two activities were poorly correlated. Unexpectedly, we found several mutants that refold luciferase normally in the absence of significant ATP turnover, presumably by increasing the flexibility of DnaK. Finally, we tested whether DnaK mutants could complement growth of DeltadnaK E. coli cells under heat shock and found that the ability to refold luciferase was more predictive of in vivo activity than ATPase rate. This study provides insights into how flexibility and co-chaperone interactions affect DnaK-mediated ATP turnover and protein folding.