Energetics of nucleotide-induced DnaK conformational states

Single-molecule mechanical studies of chaperones and their clients

Single-molecule force spectroscopy provides access to the mechanics of biomolecules. Recently, magnetic and laser optical tweezers were applied in the studies of chaperones and their interaction with protein clients. Various aspects of the chaperone-client interactions can be revealed based on the mechanical probing strategies. First, when a chaperone is probed under load, one can examine the inner workings of the chaperone while it interacts with and works on the client protein. Second, when protein clients are probed under load, the action of chaperones on folding clients can be studied in great detail. Such client folding studies have given direct access to observing actions of chaperones in real-time, like foldase, unfoldase, and holdase activity. In this review, we introduce the various single molecule mechanical techniques and summarize recent single molecule mechanical studies on heat shock proteins, chaperone-mediated folding on the ribosome, SNARE folding, and studies of chaperones involved in the folding of membrane proteins. An outlook on significant future developments is given.

An allosteric inhibitor of bacterial Hsp70 chaperone potentiates antibiotics and mitigates resistance

DnaK is the bacterial homolog of Hsp70, an ATP-dependent chaperone that helps cofactor proteins to catalyze nascent protein folding and salvage misfolded proteins. In the pathogen Mycobacterium tuberculosis, the causative agent of tuberculosis (TB), DnaK and its cofactors are proposed antimycobacterial targets, yet few small-molecule inhibitors or probes exist for these families of proteins. Here, we describe the repurposing of a drug called telaprevir that is able to allosterically inhibit the ATPase activity of DnaK and to prevent chaperone function by mimicking peptide substrates. In mycobacterial cells, telaprevir disrupts DnaK- and cofactor-mediated cellular proteostasis, resulting in enhanced efficacy of aminoglycoside antibiotics and reduced resistance to the frontline TB drug rifampin. Hence, this work contributes to a small but growing collection of protein chaperone inhibitors, and it demonstrates that these molecules disrupt bacterial mechanisms of survival in the presence of different antibiotic classes.

Fine Tuning: Effects of Post-Translational Modification on Hsp70 Chaperones

The discovery of heat shock proteins shaped our view of protein folding in the cell. Since their initial discovery, chaperone proteins were identified in all domains of life, demonstrating their vital and conserved functional roles in protein homeostasis. Chaperone proteins maintain proper protein folding in the cell by utilizing a variety of distinct, characteristic mechanisms to prevent aberrant intermolecular interactions, prevent protein aggregation, and lower entropic costs to allow for protein refolding. Continued study has found that chaperones may exhibit alternative functions, including maintaining protein folding during endoplasmic reticulum (ER) import and chaperone-mediated degradation, among others. Alternative chaperone functions are frequently controlled by post-translational modification, in which a given chaperone can switch between functions through covalent modification. This review will focus on the Hsp70 class chaperones and their Hsp40 co-chaperones, specifically highlighting the importance of post-translational control of chaperones. These modifications may serve as a target for therapeutic intervention in the treatment of diseases of protein misfolding and aggregation.

Conserved conformational selection mechanism of Hsp70 chaperone-substrate interactions

Molecular recognition is integral to biological function and frequently involves preferred binding of a molecule to one of several exchanging ligand conformations in solution. In such a process the bound structure can be selected from the ensemble of interconverting ligands a priori (conformational selection, CS) or may form once the ligand is bound (induced fit, IF). Here we focus on the ubiquitous and conserved Hsp70 chaperone which oversees the integrity of the cellular proteome through its ATP-dependent interaction with client proteins. We directly quantify the flux along CS and IF pathways using solution NMR spectroscopy that exploits a methyl TROSY effect and selective isotope-labeling methodologies. Our measurements establish that both bacterial and human Hsp70 chaperones interact with clients by selecting the unfolded state from a pre-existing array of interconverting structures, suggesting a conserved mode of client recognition among Hsp70s and highlighting the importance of molecular dynamics in this recognition event.

Proteins are the workhorses of a cell and are involved in almost all biological processes. Newly made proteins need to ‘fold’ into precise three-dimensional shapes in order to carry out their roles. However, proteins sometimes fold incorrectly or unfold. These protein forms are not able to work effectively and in some cases may even cause diseases.

Chaperone proteins help other proteins to fold correctly and are found in living organisms ranging in complexity from bacteria to humans. There are many different types of chaperones that play different roles inside cells. One, called Hsp70, binds to proteins that are incorrectly folded to help them to mature into their correct structures. However, it was not clear whether Hsp70 can also associate with the mature, correctly folded form of the proteins.

A technique called Nuclear Magnetic Resonance (NMR) spectroscopy can distinguish between mature, unfolded and chaperone-bound forms of the same protein. Sekhar et al. therefore used NMR to investigate which forms of a protein Hsp70 binds to. This revealed that both the bacterial and human versions of the Hsp70 chaperone interact only with unfolded proteins.

The results presented by Sekhar et al. also explain why Hsp70 does not disrupt the routine workings of the cell: because it does not bind to mature forms of proteins. These observations extend our understanding of how chaperones assist in folding proteins, and fit into a broader research theme exploring how proteins recognize one another. It will now be interesting to see whether the same mechanism holds for more complex forms of proteins, such as aggregates, or larger protein structures with regions of both folded and unfolded elements.

Allosteric fine-tuning of the conformational equilibrium poises the chaperone BiP for post-translational regulation

BiP is the only Hsp70 chaperone in the endoplasmic reticulum (ER) and similar to other Hsp70s, its activity relies on nucleotide- and substrate-controllable docking and undocking of its nucleotide-binding domain (NBD) and substrate-binding domain (SBD). However, little is known of specific features of the BiP conformational landscape that tune BiP to its unique tasks and the ER environment. We present methyl NMR analysis of the BiP chaperone cycle that reveals surprising conformational heterogeneity of ATP-bound BiP that distinguishes BiP from its bacterial homologue DnaK. This unusual poise enables gradual post-translational regulation of the BiP chaperone cycle and its chaperone activity by subtle local perturbations at SBD allosteric 'hotspots'. In particular, BiP inactivation by AMPylation of its SBD does not disturb Hsp70 inter-domain allostery and preserves BiP structure. Instead it relies on a redistribution of the BiP conformational ensemble and stabilization the domain-docked conformation in presence of ADP and ATP.

The intervening domain from MeCP2 enhances the DNA affinity of the methyl binding domain and provides an independent DNA interaction site

Methyl-CpG binding protein 2 (MeCP2) preferentially interacts with methylated DNA and it is involved in epigenetic regulation and chromatin remodelling. Mutations in MeCP2 are linked to Rett syndrome, the leading cause of intellectual retardation in girls and causing mental, motor and growth impairment. Unstructured regions in MeCP2 provide the plasticity for establishing interactions with multiple binding partners. We present a biophysical characterization of the methyl binding domain (MBD) from MeCP2 reporting the contribution of flanking domains to its structural stability and dsDNA interaction. The flanking disordered intervening domain (ID) increased the structural stability of MBD, modified its dsDNA binding profile from an entropically-driven moderate-affinity binding to an overwhelmingly enthalpically-driven high-affinity binding. Additionally, ID provided an additional site for simultaneously and autonomously binding an independent dsDNA molecule, which is a key feature linked to the chromatin remodelling and looping activity of MeCP2, as well as its ability to interact with nucleosomes replacing histone H1. The dsDNA interaction is characterized by an unusually large heat capacity linked to a cluster of water molecules trapped within the binding interface. The dynamics of disordered regions together with extrinsic factors are key determinants of MeCP2 global structural properties and functional capabilities.

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

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

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

A Spatial Control for Correct Timing of Gene Expression during the Escherichia coli Cell Cycle

Temporal transcriptions of genes are achieved by different mechanisms such as dynamic interaction of activator and repressor proteins with promoters, and accumulation and/or degradation of key regulators as a function of cell cycle. We find that the TorR protein localizes to the old poles of the Escherichia coli cells, forming a functional focus. The TorR focus co-localizes with the nucleoid in a cell-cycle-dependent manner, and consequently regulates transcription of a number of genes. Formation of one TorR focus at the old poles of cells requires interaction with the MreB and DnaK proteins, and ATP, suggesting that TorR delivery requires cytoskeleton organization and ATP. Further, absence of the protein–protein interactions and ATP leads to loss in function of TorR as a transcription factor. We propose a mechanism for timing of cell-cycle-dependent gene transcription, where a transcription factor interacts with its target genes during a specific period of the cell cycle by limiting its own spatial distribution.

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.

Substrate-binding domain conformational dynamics mediate Hsp70 allostery

Binding of ATP to the N-terminal nucleotide-binding domain (NBD) of heat shock protein 70 (Hsp70) molecular chaperones reduces the affinity of their C-terminal substrate-binding domain (SBD) for unfolded protein substrates. ATP binding to the NBD leads to docking between NBD and βSBD and releasing of the α-helical lid that covers the substrate-binding cleft in the SBD. However, these structural changes alone do not fully account for the allosteric mechanism of modulation of substrate affinity and binding kinetics. Through a multipronged study of the Escherichia coli Hsp70 DnaK, we found that changes in conformational dynamics within the βSBD play a central role in interdomain allosteric communication in the Hsp70 DnaK. ATP-mediated NBD conformational changes favor formation of NBD contacts with lynchpin sites on the βSBD and force disengagement of SBD strand β8 from strand β7, which leads to repacking of a βSBD hydrophobic cluster and disruption of the hydrophobic arch over the substrate-binding cleft. In turn, these structural rearrangements drastically enhance conformational dynamics throughout the entire βSBD and particularly around the substrate-binding site. This negative, entropically driven allostery between two functional sites of the βSBD-the NBD binding interface and the substrate-binding site-confers upon the SBD the plasticity needed to bind to a wide range of chaperone clients without compromising precise control of thermodynamics and kinetics of chaperone-client interactions.

Modulation of the chaperone DnaK allosterism by the nucleotide exchange factor GrpE

Hsp70 chaperones comprise two domains, the nucleotide-binding domain (Hsp70NBD), responsible for structural and functional changes in the chaperone, and the substrate-binding domain (Hsp70SBD), involved in substrate interaction. Substrate binding and release in Hsp70 is controlled by the nucleotide state of DnaKNBD, with ATP inducing the open, substrate-receptive DnaKSBD conformation, whereas ADP forces its closure. DnaK cycles between the two conformations through interaction with two cofactors, the Hsp40 co-chaperones (DnaJ in Escherichia coli) induce the ADP state, and the nucleotide exchange factors (GrpE in E. coli) induce the ATP state. X-ray crystallography showed that the GrpE dimer is a nucleotide exchange factor that works by interaction of one of its monomers with DnaKNBD. DnaKSBD location in this complex is debated; there is evidence that it interacts with the GrpE N-terminal disordered region, far from DnaKNBD. Although we confirmed this interaction using biochemical and biophysical techniques, our EM-based three-dimensional reconstruction of the DnaK-GrpE complex located DnaKSBD near DnaKNBD. This apparent discrepancy between the functional and structural results is explained by our finding that the tail region of the GrpE dimer in the DnaK-GrpE complex bends and its tip contacts DnaKSBD, whereas the DnaKNBD-DnaKSBD linker contacts the GrpE helical region. We suggest that these interactions define a more complex role for GrpE in the control of DnaK function.

Crowding activates ClpB and enhances its association with DnaK for efficient protein aggregate reactivation

Reactivation of intracellular protein aggregates after a severe stress is mandatory for cell survival. In bacteria, this activity depends on the collaboration between the DnaK system and ClpB, which in vivo occurs in a highly crowded environment. The reactivation reaction includes two steps: extraction of unfolded monomers from the aggregate and their subsequent refolding into the native conformation. Both steps might be compromised by excluded volume conditions that would favor aggregation of unstable protein folding intermediates. Here, we have investigated whether ClpB and the DnaK system are able to compensate this unproductive effect and efficiently reactivate aggregates of three different substrate proteins under crowding conditions. To this aim, we have compared the association equilibrium, biochemical properties, stability, and chaperone activity of the disaggregase ClpB in the absence and presence of an inert macromolecular crowding agent. Our data show that crowding i), increases three to four orders of magnitude the association constant of the functional hexamer; ii), shifts the conformational equilibrium of the protein monomer toward a compact state; iii), stimulates its ATPase activity; and iv), favors association of the chaperone with substrate proteins and with aggregate-bound DnaK. These effects strongly enhance protein aggregate reactivation by the DnaK-ClpB network, highlighting the importance of volume exclusion in complex processes in which several proteins have to work in a sequential manner.

Specialized Hsp70 chaperone (HscA) binds preferentially to the disordered form, whereas J-protein (HscB) binds preferentially to the structured form of the iron-sulfur cluster scaffold protein (IscU)

The Escherichia coli protein IscU serves as the scaffold for Fe-S cluster assembly and the vehicle for Fe-S cluster transfer to acceptor proteins, such as apoferredoxin. IscU populates two conformational states in solution, a structured conformation (S) that resembles the conformation of the holoprotein IscU-[2Fe-2S] and a dynamically disordered conformation (D) that does not bind metal ions. NMR spectroscopic results presented here show that the specialized Hsp70 chaperone (HscA), alone or as the HscA-ADP complex, preferentially binds to and stabilizes the D-state of IscU. IscU is released when HscA binds ATP. By contrast, the J-protein HscB binds preferentially to the S-state of IscU. Consistent with these findings, we propose a mechanism in which cluster transfer is coupled to hydrolysis of ATP bound to HscA, conversion of IscU to the D-state, and release of HscB.

Applications of isothermal titration calorimetry in pure and applied research--survey of the literature from 2010

Isothermal titration calorimetry (ITC) is a biophysical technique for measuring the formation and dissociation of molecular complexes and has become an invaluable tool in many branches of science from cell biology to food chemistry. By measuring the heat absorbed or released during bond formation, ITC provides accurate, rapid, and label-free measurement of the thermodynamics of molecular interactions. In this review, we survey the recent literature reporting the use of ITC and have highlighted a number of interesting studies that provide a flavour of the diverse systems to which ITC can be applied. These include measurements of protein-protein and protein-membrane interactions required for macromolecular assembly, analysis of enzyme kinetics, experimental validation of molecular dynamics simulations, and even in manufacturing applications such as food science. Some highlights include studies of the biological complex formed by Staphylococcus aureus enterotoxin C3 and the murine T-cell receptor, the mechanism of membrane association of the Parkinson's disease-associated protein α-synuclein, and the role of non-specific tannin-protein interactions in the quality of different beverages. Recent developments in automation are overcoming limitations on throughput imposed by previous manual procedures and promise to greatly extend usefulness of ITC in the future. We also attempt to impart some practical advice for getting the most out of ITC data for those researchers less familiar with the method.

Isothermal titration calorimetry for measuring macromolecule-ligand affinity

Isothermal titration calorimetry (ITC) is a useful tool for understanding the complete thermodynamic picture of a binding reaction. In biological sciences, macromolecular interactions are essential in understanding the machinery of the cell. Experimental conditions, such as buffer and temperature, can be tailored to the particular binding system being studied. However, careful planning is needed since certain ligand and macromolecule concentration ranges are necessary to obtain useful data. Concentrations of the macromolecule and ligand need to be accurately determined for reliable results. Care also needs to be taken when preparing the samples as impurities can significantly affect the experiment. When ITC experiments, along with controls, are performed properly, useful binding information, such as the stoichiometry, affinity and enthalpy, are obtained. By running additional experiments under different buffer or temperature conditions, more detailed information can be obtained about the system. A protocol for the basic setup of an ITC experiment is given.

Fine tuning of a biological machine: DnaK gains improved chaperone activity by altered allosteric communication and substrate binding

DnaK is a member of the Hsp70 family of molecular chaperones. This molecular machine couples the binding and hydrolysis of ATP to binding and release of substrate proteins. The switches that are involved in allosteric communication within this multidomain protein are mostly unknown. Previous insights were largely obtained by mutants, which displayed either wild-type activity or reduced folding assistance of substrate proteins. With a directed evolution approach for improved folding assistance we selected a DnaK variant characterized by a glycine to alanine substitution at position 384 (G384A); this resulted in a 2.5-fold higher chaperone activity in an in vitro DnaK-assisted firefly luciferase refolding assay. Quantitative biochemical characterization revealed several changes of key kinetic parameters compared to the wild type. Most pronounced is a 13-fold reduced rate constant for substrate release in the ATP-bound state, which we assume, in conjunction with the resulting increase in substrate affinity, to be related to improved chaperone activity. As the underlying mechanistic reason for this change we propose an altered interface of allosteric communication of mutant G384A, which is notably located at a hinge position between nucleotide and substrate binding domain.

Allosteric communication between the nucleotide binding domains of caseinolytic peptidase B

ClpB is a hexameric chaperone that solubilizes and reactivates protein aggregates in cooperation with the Hsp70/DnaK chaperone system. Each of the identical protein monomers contains two nucleotide binding domains (NBD), whose ATPase activity must be coupled to exert on the substrate the mechanical work required for its reactivation. However, how communication between these sites occurs is at present poorly understood. We have studied herein the affinity of each of the NBDs for nucleotides in WT ClpB and protein variants in which one or both sites are mutated to selectively impair nucleotide binding or hydrolysis. Our data show that the affinity of NBD2 for nucleotides (K(d) = 3-7 μm) is significantly higher than that of NBD1. Interestingly, the affinity of NBD1 depends on nucleotide binding to NBD2. Binding of ATP, but not ADP, to NBD2 increases the affinity of NBD1 (the K(d) decreases from ≈160-300 to 50-60 μm) for the corresponding nucleotide. Moreover, filling of the NBD2 ring with ATP allows the cooperative binding of this nucleotide and substrates to the NBD1 ring. Data also suggest that a minimum of four subunits cooperate to bind and reactivate two different aggregated protein substrates.

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.

Conformational sampling and nucleotide-dependent transitions of the GroEL subunit probed by unbiased molecular dynamics simulations

GroEL is an ATP dependent molecular chaperone that promotes the folding of a large number of substrate proteins in E. coli. Large-scale conformational transitions occurring during the reaction cycle have been characterized from extensive crystallographic studies. However, the link between the observed conformations and the mechanisms involved in the allosteric response to ATP and the nucleotide-driven reaction cycle are not completely established. Here we describe extensive (in total long) unbiased molecular dynamics (MD) simulations that probe the response of GroEL subunits to ATP binding. We observe nucleotide dependent conformational transitions, and show with multiple 100 ns long simulations that the ligand-induced shift in the conformational populations are intrinsically coded in the structure-dynamics relationship of the protein subunit. Thus, these simulations reveal a stabilization of the equatorial domain upon nucleotide binding and a concomitant “opening” of the subunit, which reaches a conformation close to that observed in the crystal structure of the subunits within the ADP-bound oligomer. Moreover, we identify changes in a set of unique intrasubunit interactions potentially important for the conformational transition.

Molecular machines convert chemical energy to mechanical work in the process of carrying out their specific tasks. Often these proteins are fueled by ATP binding and hydrolysis, enabling switching between different conformations. The ATP-dependent chaperone GroEL is a molecular machine that opens and closes its barrel-like structure in order to provide a folding cage for unfolded proteins. The quest to fully understand and control GroEL and other molecular machines is enhanced by complementing experimental work with computational approaches. Here, we provide a description of the molecular basis for the conformational changes in the GroEL subunit by performing extensive molecular dynamics simulations. The simulations sample the conformational population for the different nucleotide-free and bound states in the isolated subunit. The results reveal that the conformations of the subunit when isolated resemble those of the subunit integrated in the GroEL complex. Moreover, the molecular dynamics simulations allow following detailed changes in individual interatomic interactions brought about by ATP-binding.