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

Comparison of allosteric signaling in DnaK and BiP using mutual information between simulated residue conformations

The heat shock protein 70 kDa (Hsp70) chaperone system serves as a critical component of protein quality control across a wide range of prokaryotic and eukaryotic organisms. Divergent evolution and specialization to particular organelles have produced numerous Hsp70 variants which share similarities in structure and general function, but differ substantially in regulatory aspects, including conformational dynamics and activity modulation by cochaperones. The human Hsp70 variant BiP (also known as GRP78 or HSPA5) is of therapeutic interest in the context of cancer, neurodegenerative diseases, and viral infection, including for treatment of the pandemic virus SARS-CoV-2. Due to the complex conformational rearrangements and high sequential variance within the Hsp70 protein family, it is in many cases poorly understood which amino acid mutations are responsible for biochemical differences between protein variants. In this study, we predicted residues associated with conformational regulation of human BiP and Escherichia coli DnaK. Based on protein structure networks obtained from molecular dynamics simulations, we analyzed the shared information between interaction timelines to highlight residue positions with strong conformational coupling to their environment. Our predictions, which focus on the binding processes of the chaperone's substrate and cochaperones, indicate residues filling potential signaling roles specific to either DnaK or BiP. By combining predictions of individual residues into conformationally coupled chains connecting ligand binding sites, we predict a BiP specific secondary signaling pathway associated with substrate binding. Our study sheds light on mechanistic differences in signaling and regulation between Hsp70 variants, which provide insights relevant to therapeutic applications of these proteins.

Nucleotide Exchange Factors for Hsp70 Molecular Chaperones: GrpE, Hsp110/Grp170, HspBP1/Sil1, and BAG Domain Proteins

Molecular chaperones of the Hsp70 family are key components of the cellular protein-folding machinery. Substrate folding is accomplished by iterative cycles of ATP binding, hydrolysis, and release. The ATPase activity of Hsp70 is regulated by two main classes of cochaperones: J-domain proteins stimulate ATPase hydrolysis by Hsp70, while nucleotide exchange factors (NEFs) facilitate the conversion from the ADP-bound to the ATP-bound state, thus closing the chaperone folding cycle. NEF function can additionally be antagonized by ADP dissociation inhibitors. Beginning with the discovery of the prototypical bacterial NEF, GrpE, a large diversity of nucleotide exchange factors for Hsp70 have been identified, connecting it to a multitude of cellular processes in the eukaryotic cell. Here we review recent advances toward structure and function of nucleotide exchange factors from the Hsp110/Grp170, HspBP1/Sil1, and BAG domain protein families and discuss how these cochaperones connect protein folding with cellular quality control and degradation pathways.

Hsp70 molecular chaperones: multifunctional allosteric holding and unfolding machines

The Hsp70 family of chaperones works with its co-chaperones, the nucleotide exchange factors and J-domain proteins, to facilitate a multitude of cellular functions. Central players in protein homeostasis, these jacks-of-many-trades are utilized in a variety of ways because of their ability to bind with selective promiscuity to regions of their client proteins that are exposed when the client is unfolded, either fully or partially, or visits a conformational state that exposes the binding region in a regulated manner. The key to Hsp70 functions is that their substrate binding is transient and allosterically cycles in a nucleotide-dependent fashion between high- and low-affinity states. In the past few years, structural insights into the molecular mechanism of this allosterically regulated binding have emerged and provided deep insight into the deceptively simple Hsp70 molecular machine that is so widely harnessed by nature for diverse cellular functions. In this review, these structural insights are discussed to give a picture of the current understanding of how Hsp70 chaperones work.

Perturbation-Response Scanning Reveals Key Residues for Allosteric Control in Hsp70

Hsp70 molecular chaperones play an important role in maintaining cellular homeostasis, and are implicated in a wide array of cellular processes, including protein recovery from aggregates, cross-membrane protein translocation, and protein biogenesis. Hsp70 consists of two domains, a nucleotide binding domain (NBD) and a substrate binding domain (SBD), each of which communicates via an allosteric mechanism such that the protein interconverts between two functional states, an ATP-bound open conformation and an ADP-bound closed conformation. The exact mechanism for interstate conversion is not as yet fully understood. However, the ligand-bound states of the NBD and SBD as well as interactions with cochaperones such as DnaJ and nucleotide exchange factor are thought to play crucial regulatory roles. In this study, we apply the perturbation-response scanning (PRS) method in combination with molecular dynamics simulations as a computational tool for the identification of allosteric hot residues in the large multidomain Hsp70 protein. We find evidence in support of the hypothesis that substrate binding triggers ATP hydrolysis and that the ADP-substrate complex favors interstate conversion to the closed state. Furthermore, our data are in agreement with the proposal that there is an allosterically active intermediate state between the open and closed states and vice versa, as we find evidence that ATP binding to the closed structure and peptide binding to the open structure allosterically "activate" the respective complexes. We conclude our analysis by showing how our PRS data fit the current opinion on the Hsp70 conformational cycle and present several allosteric hot residues that may provide a platform for further studies to gain additional insight into Hsp70 allostery.

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.

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.

The nucleotide exchange factors of Hsp70 molecular chaperones

Molecular chaperones of the Hsp70 family form an important hub in the cellular protein folding networks in bacteria and eukaryotes, connecting translation with the downstream machineries of protein folding and degradation. The Hsp70 folding cycle is driven by two types of cochaperones: J-domain proteins stimulate ATP hydrolysis by Hsp70, while nucleotide exchange factors (NEFs) promote replacement of Hsp70-bound ADP with ATP. Bacteria and organelles of bacterial origin have only one known NEF type for Hsp70, GrpE. In contrast, a large diversity of Hsp70 NEFs has been discovered in the eukaryotic cell. These NEFs belong to the Hsp110/Grp170, HspBP1/Sil1, and BAG domain protein families. In this short review we compare the structures and molecular mechanisms of nucleotide exchange factors for Hsp70 and discuss how these cochaperones contribute to protein folding and quality control in the cell.

Substrate Interaction Networks of the Escherichia coli Chaperones: Trigger Factor, DnaK and GroEL

In the dense cellular environment, protein misfolding and inter-molecular protein aggregation compete with protein folding. Chaperones associate with proteins to prevent misfolding and to assist in folding to the native state. In Escherichia coli, the chaperones trigger factor, DnaK/DnaJ/GrpE, and GroEL/ES are the major chaperones responsible for insuring proper de novo protein folding. With multitudes of proteins produced by the bacterium, the chaperones have to be selective for their substrates. Yet, chaperone selectivity cannot be too specific. Recent biochemical and high-throughput studies have provided important insights highlighting the strategies used by chaperones in maintaining proteostasis in the cell. Here, we discuss the substrate networks and cooperation among these protein folding chaperones.

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.