Modeling Hsp70/Hsp40 interaction by multi-scale molecular simulations and coevolutionary sequence analysis

CISD3/MiNT is required for complex I function, mitochondrial integrity, and skeletal muscle maintenance

Mitochondria play a central role in muscle metabolism and function. A unique family of iron-sulfur proteins, termed CDGSH Iron Sulfur Domain-containing (CISD/NEET) proteins, support mitochondrial function in skeletal muscles. The abundance of these proteins declines during aging leading to muscle degeneration. Although the function of the outer mitochondrial CISD/NEET proteins, CISD1/mitoNEET and CISD2/NAF-1, has been defined in skeletal muscle cells, the role of the inner mitochondrial CISD protein, CISD3/MiNT, is currently unknown. Here, we show that CISD3 deficiency in mice results in muscle atrophy that shares proteomic features with Duchenne muscular dystrophy. We further reveal that CISD3 deficiency impairs the function and structure of skeletal muscles, as well as their mitochondria, and that CISD3 interacts with, and donates its [2Fe-2S] clusters to, complex I respiratory chain subunit NADH Ubiquinone Oxidoreductase Core Subunit V2 (NDUFV2). Using coevolutionary and structural computational tools, we model a CISD3-NDUFV2 complex with proximal coevolving residue interactions conducive of [2Fe-2S] cluster transfer reactions, placing the clusters of the two proteins 10 to 16 Å apart. Taken together, our findings reveal that CISD3/MiNT is important for supporting the biogenesis and function of complex I, essential for muscle maintenance and function. Interventions that target CISD3 could therefore impact different muscle degeneration syndromes, aging, and related conditions.

The Interplay between Heat Shock Proteins and Cancer Pathogenesis: A Novel Strategy for Cancer Therapeutics

Heat shock proteins (HSPs) are developmentally conserved families of protein found in both prokaryotic and eukaryotic organisms. HSPs are engaged in a diverse range of physiological processes, including molecular chaperone activity to assist the initial protein folding or promote the unfolding and refolding of misfolded intermediates to acquire the normal or native conformation and its translocation and prevent protein aggregation as well as in immunity, apoptosis, and autophagy. These molecular chaperonins are classified into various families according to their molecular size or weight, encompassing small HSPs (e.g., HSP10 and HSP27), HSP40, HSP60, HSP70, HSP90, and the category of large HSPs that include HSP100 and ClpB proteins. The overexpression of HSPs is induced to counteract cell stress at elevated levels in a variety of solid tumors, including anticancer chemotherapy, and is closely related to a worse prognosis and therapeutic resistance to cancer cells. HSPs are also involved in anti-apoptotic properties and are associated with processes of cancer progression and development, such as metastasis, invasion, and cell proliferation. This review outlines the previously mentioned HSPs and their significant involvement in diverse mechanisms of tumor advancement and metastasis, as well as their contribution to identifying potential targets for therapeutic interventions.

Bacterial J-Domains with C-Terminal Tags Contact the Substrate Binding Domain of DnaK and Sequester Chaperone Activity

Functional interactions between the molecular chaperone DnaK and cofactor J-proteins (DnaJs), as well as their homologs, are crucial to the maintenance of proteostasis across cell types. In the bacterial pathogen Mycobacterium tuberculosis, DnaK-DnaJ interactions are essential for cell growth and represent potential targets for antibiotic or adjuvant development. While the N-terminal J-domains of J-proteins are known to form important contacts with DnaK, C-terminal domains have varied roles. Here, we have studied the effect of adding C-terminal tags to N-terminal J-domain truncations of mycobacterial DnaJ1 and DnaJ2 to promote additional interactions with DnaK. We found that His6 tags uniquely promote binding to additional sites in the substrate binding domain at the C-terminus of DnaK. Other C-terminal tags attached to J-domains, even peptides known to interact with DnaK, do not produce the same effects. Expression of C-terminally modified DnaJ1 or DnaJ2 J-domains in mycobacterial cells suppresses chaperone activity following proteotoxic stress, which is exaggerated in the presence of a small-molecule DnaK inhibitor. Hence, this work uncovers genetically encodable J-protein variants that may be used to study chaperone-cofactor interactions in other organisms.

The Hsp70 and JDP proteins: Structure-function perspective on molecular chaperone activity

Proteins are the most structurally diverse cellular biomolecules that act as molecular machines driving essential activities of all living organisms. To be functional, most of the proteins need to fold into a specific three-dimensional structure, which on one hand should be stable enough to oppose disruptive conditions and on the other hand flexible enough to allow conformational dynamics necessary for their biological functions. This compromise between stability and dynamics makes proteins susceptible to stress-induced misfolding and aggregation. Moreover, the folding process itself is intrinsically prone to conformational errors. Molecular chaperones are proteins that mitigate folding defects and maintain the structural integrity of the cellular proteome. Promiscuous Hsp70 chaperones are central to these processes and their activity depends on the interaction with obligatory J-domain protein (JDP) partners. In this review, we discuss structural aspects of Hsp70s, JDPs, and their interaction in the context of biological activities.

Data-driven large-scale genomic analysis reveals an intricate phylogenetic and functional landscape in J-domain proteins

The 70-kD heat shock protein (Hsp70) chaperone system is a central hub of the proteostasis network that helps maintain protein homeostasis in all organisms. The recruitment of Hsp70 to perform different and specific cellular functions is regulated by the J-domain protein (JDP) co-chaperone family carrying the small namesake J-domain, required to interact and drive the ATPase cycle of Hsp70s. Besides the J-domain, prokaryotic and eukaryotic JDPs display a staggering diversity in domain architecture, function, and cellular localization. Very little is known about the overall JDP family, despite their essential role in cellular proteostasis, development, and its link to a broad range of human diseases. In this work, we leverage the exponentially increasing number of JDP gene sequences identified across all kingdoms owing to the advancements in sequencing technology and provide a broad overview of the JDP repertoire. Using an automated classification scheme based on artificial neural networks (ANNs), we demonstrate that the sequences of J-domains carry sufficient discriminatory information to reliably recover the phylogeny, localization, and domain composition of the corresponding full-length JDP. By harnessing the interpretability of the ANNs, we find that many of the discriminatory sequence positions match residues that form the interaction interface between the J-domain and Hsp70. This reveals that key residues within the J-domains have coevolved with their obligatory Hsp70 partners to build chaperone circuits for specific functions in cells.

The Plasmodium falciparum exported J domain proteins fine-tune human and malarial Hsp70s: pathological exploitation of proteostasis machinery

Cellular proteostasis requires a network of molecular chaperones and co-chaperones, which facilitate the correct folding and assembly of other proteins, or the degradation of proteins misfolded beyond repair. The function of the major chaperones, heat shock protein 70 (Hsp70) and heat shock protein 90 (Hsp90), is regulated by a cohort of co-chaperone proteins. The J domain protein (JDP) family is one of the most diverse co-chaperone families, playing an important role in functionalizing the Hsp70 chaperone system to form a powerful protein quality control network. The intracellular malaria parasite, Plasmodium falciparum, has evolved the capacity to invade and reboot mature human erythrocytes, turning them into a vehicles of pathology. This process appears to involve the harnessing of both the human and parasite chaperone machineries. It is well known that malaria parasite-infected erythrocytes are highly enriched in functional human Hsp70 (HsHsp70) and Hsp90 (HsHsp90), while recent proteomics studies have provided evidence that human JDPs (HsJDPs) may also be enriched, but at lower levels. Interestingly, P. falciparum JDPs (PfJDPs) are the most prominent and diverse family of proteins exported into the infected erythrocyte cytosol. We hypothesize that the exported PfJPDs may be an evolutionary consequence of the need to boost chaperone power for specific protein folding pathways that enable both survival and pathogenesis of the malaria parasite. The evidence suggests that there is an intricate network of PfJDP interactions with the exported malarial Hsp70 (PfHsp70-x) and HsHsp70, which appear to be important for the trafficking of key malarial virulence factors, and the proteostasis of protein complexes of human and parasite proteins associated with pathology. This review will critically evaluate the current understanding of the role of exported PfJDPs in pathological exploitation of the proteostasis machinery by fine-tuning the chaperone properties of both human and malarial Hsp70s.

Genome-Wide Identification and Analysis of the Hsp40/J-Protein Family Reveals Its Role in Soybean (Glycine max) Growth and Development

The J-protein family comprises molecular chaperones involved in plant growth, development, and stress responses. Little is known about this gene family in soybean. Hence, we characterized J-protein genes in soybean, with the most highly expressed and responsive during flower and seed development. We also revealed their phylogeny, structure, motif analysis, chromosome location, and expression. Based on their evolutionary links, we divided the 111 potential soybean J-proteins into 12 main clades (I-XII). Gene-structure estimation revealed that each clade had an exon-intron structure resembling or comparable to others. Most soybean J-protein genes lacked introns in Clades I, III, and XII. Moreover, transcriptome data obtained from a publicly accessible soybean database and RT-qPCR were used to examine the differential expression of DnaJ genes in various soybean tissues and organs. The expression level of DnaJ genes indicated that, among 14 tissues, at least one tissue expressed the 91 soybean genes. The findings suggest that J-protein genes could be involved in the soybean growth period and offer a baseline for further functional research into J-proteins' role in soybean. One important application is the identification of J-proteins that are highly expressed and responsive during flower and seed development in soybean. These genes likely play crucial roles in these processes, and their identification can contribute to breeding programs to improve soybean yield and quality.

Combining phylogeny and coevolution improves the inference of interaction partners among paralogous proteins

Predicting protein-protein interactions from sequences is an important goal of computational biology. Various sources of information can be used to this end. Starting from the sequences of two interacting protein families, one can use phylogeny or residue coevolution to infer which paralogs are specific interaction partners within each species. We show that these two signals can be combined to improve the performance of the inference of interaction partners among paralogs. For this, we first align the sequence-similarity graphs of the two families through simulated annealing, yielding a robust partial pairing. We next use this partial pairing to seed a coevolution-based iterative pairing algorithm. This combined method improves performance over either separate method. The improvement obtained is striking in the difficult cases where the average number of paralogs per species is large or where the total number of sequences is modest.

J-domain protein chaperone circuits in proteostasis and disease

The J-domain proteins (JDP) form the largest protein family among cellular chaperones. In cooperation with the Hsp70 chaperone system, these co-chaperones orchestrate a plethora of distinct functions, including those that help maintain cellular proteostasis and development. JDPs evolved largely through the fusion of a J-domain with other protein subdomains. The highly conserved J-domain facilitates the binding and activation of Hsp70s. How JDPs (re)wire Hsp70 chaperone circuits and promote functional diversity remains insufficiently explained. Here, we discuss recent advances in our understanding of the JDP family with a focus on the regulation built around J-domains to ensure correct pairing and assembly of JDP-Hsp70 machineries that operate on different clientele under various cellular growth conditions.

Protein Cofactor Mimics Disrupt Essential Chaperone Function in Stressed Mycobacteria

Bacterial DnaK is an ATP-dependent molecular chaperone important for maintaining cellular proteostasis in concert with cofactor proteins. The cofactor DnaJ delivers non-native client proteins to DnaK and activates its ATPase activity, which is required for protein folding. In the bacterial pathogen Mycobacterium tuberculosis, DnaK is assisted by two DnaJs, DnaJ1 and DnaJ2. Functional protein-protein interactions (PPIs) between DnaK and at least one DnaJ are essential for survival of mycobacteria; hence, these PPIs represent untapped antibacterial targets. Here, we synthesize peptide-based mimetics of DnaJ1 and DnaJ2 N-terminal domains as rational inhibitors of DnaK-cofactor interactions. We find that covalently stabilized DnaJ mimetics are capable of disrupting DnaK-cofactor activity in vitro and prevent mycobacterial recovery from proteotoxic stress in vivo, leading to cell death. Since chaperones and cofactors are highly conserved, we anticipate these results will inform the design of other mimetics to modulate chaperone function across cell types.

Heat shock proteins and the calcineurin-crz1 signaling regulate stress responses in fungi

The heat shock proteins (Hsps) act as a molecular chaperone to stabilize client proteins involved in various cell functions in fungi. Hsps are classified into different families such as HSP90, HSP70, HSP60, HSP40, and small HSPs (sHsps). Hsp90, a well-studied member of the Hsp family proteins, plays a role in growth, cell survival, and pathogenicity in fungi. Hsp70 and sHsps are involved in the development, tolerance to stress conditions, and drug resistance in fungi. Hsp60 is a mitochondrial chaperone, and Hsp40 regulates fungal ATPase machinery. In this review, we describe the cell functions, regulation, and the molecular link of the Hsps with the calcineurin-crz1 calcium signaling pathway for their role in cell survival, growth, virulence, and drug resistance in fungi and related organisms.

Modeling the Dynamics of Protein–Protein Interfaces, How and Why?

Protein–protein assemblies act as a key component in numerous cellular processes. Their accurate modeling at the atomic level remains a challenge for structural biology. To address this challenge, several docking and a handful of deep learning methodologies focus on modeling protein–protein interfaces. Although the outcome of these methods has been assessed using static reference structures, more and more data point to the fact that the interaction stability and specificity is encoded in the dynamics of these interfaces. Therefore, this dynamics information must be taken into account when modeling and assessing protein interactions at the atomistic scale. Expanding on this, our review initially focuses on the recent computational strategies aiming at investigating protein–protein interfaces in a dynamic fashion using enhanced sampling, multi-scale modeling, and experimental data integration. Then, we discuss how interface dynamics report on the function of protein assemblies in globular complexes, in fuzzy complexes containing intrinsically disordered proteins, as well as in active complexes, where chemical reactions take place across the protein–protein interface.

Divergence in Dimerization and Activity of Primate APOBEC3C

The APOBEC3 (A3) family of single-stranded DNA cytidine deaminases are host restriction factors that inhibit lentiviruses, such as HIV-1, in the absence of the Vif protein that causes their degradation. Deamination of cytidine in HIV-1 (−)DNA forms uracil that causes inactivating mutations when uracil is used as a template for (+)DNA synthesis. For APOBEC3C (A3C), the chimpanzee and gorilla orthologues are more active than human A3C, and we determined that Old World Monkey A3C from rhesus macaque (rh) is not active against HIV-1. Biochemical, virological, and coevolutionary analyses combined with molecular dynamics simulations showed that the key amino acids needed to promote rhA3C antiviral activity, 44, 45, and 144, also promoted dimerization and changes to the dynamics of loop 1, near the enzyme active site. Although forced evolution of rhA3C resulted in a similar dimer interface with hominid A3C, the key amino acid contacts were different. Overall, our results determine the basis for why rhA3C is less active than human A3C and establish the amino acid network for dimerization and increased activity. Based on identification of the key amino acids determining Old World Monkey antiviral activity we predict that other Old World Monkey A3Cs did not impart anti-lentiviral activity, despite fixation of a key residue needed for hominid A3C activity. Overall, the coevolutionary analysis of the A3C dimerization interface presented also provides a basis from which to analyze dimerization interfaces of other A3 family members.

Reversible phosphorylation of cyclin T1 promotes assembly and stability of P-TEFb

The positive transcription elongation factor b (P-TEFb) is a critical coactivator for transcription of most cellular and viral genes, including those of HIV. While P-TEFb is regulated by 7SK snRNA in proliferating cells, P-TEFb is absent due to diminished levels of CycT1 in quiescent and terminally differentiated cells, which has remained unexplored. In these cells, we found that CycT1 not bound to CDK9 is rapidly degraded. Moreover, productive CycT1:CDK9 interactions are increased by PKC-mediated phosphorylation of CycT1 in human cells. Conversely, dephosphorylation of CycT1 by PP1 reverses this process. Thus, PKC inhibitors or removal of PKC by chronic activation results in P-TEFb disassembly and CycT1 degradation. This finding not only recapitulates P-TEFb depletion in resting CD4+ T cells but also in anergic T cells. Importantly, our studies reveal mechanisms of P-TEFb inactivation underlying T cell quiescence, anergy, and exhaustion as well as proviral latency and terminally differentiated cells.

Targeting Chaperone/Co-Chaperone Interactions with Small Molecules: A Novel Approach to Tackle Neurodegenerative Diseases

The dysfunction of the proteostasis network is a molecular hallmark of neurodegenerative diseases such as Alzheimer’s disease, Parkinson’s disease, Huntington’s disease, and amyotrophic lateral sclerosis. Molecular chaperones are a major component of the proteostasis network and maintain cellular homeostasis by folding client proteins, assisting with intracellular transport, and interfering with protein aggregation or degradation. Heat shock protein 70 kDa (Hsp70) and 90 kDa (Hsp90) are two of the most important chaperones whose functions are dependent on ATP hydrolysis and collaboration with their co-chaperones. Numerous studies implicate Hsp70, Hsp90, and their co-chaperones in neurodegenerative diseases. Targeting the specific protein–protein interactions between chaperones and their particular partner co-chaperones with small molecules provides an opportunity to specifically modulate Hsp70 or Hsp90 function for neurodegenerative diseases. Here, we review the roles of co-chaperones in Hsp70 or Hsp90 chaperone cycles, the impacts of co-chaperones in neurodegenerative diseases, and the development of small molecules modulating chaperone/co-chaperone interactions. We also provide a future perspective of drug development targeting chaperone/co-chaperone interactions for neurodegenerative diseases.

The Hsp70-Chaperone Machines in Bacteria

The ATP-dependent Hsp70s are evolutionary conserved molecular chaperones that constitute central hubs of the cellular protein quality surveillance network. None of the other main chaperone families (Tig, GroELS, HtpG, IbpA/B, ClpB) have been assigned with a comparable range of functions. Through a multitude of functions Hsp70s are involved in many cellular control circuits for maintaining protein homeostasis and have been recognized as key factors for cell survival. Three mechanistic properties of Hsp70s are the basis for their high versatility. First, Hsp70s bind to short degenerate sequence motifs within their client proteins. Second, Hsp70 chaperones switch in a nucleotide-controlled manner between a state of low affinity for client proteins and a state of high affinity for clients. Third, Hsp70s are targeted to their clients by a large number of cochaperones of the J-domain protein (JDP) family and the lifetime of the Hsp70-client complex is regulated by nucleotide exchange factors (NEF). In this review I will discuss advances in the understanding of the molecular mechanism of the Hsp70 chaperone machinery focusing mostly on the bacterial Hsp70 DnaK and will compare the two other prokaryotic Hsp70s HscA and HscC with DnaK.

Large-scale discovery of protein interactions at residue resolution using co-evolution calculated from genomic sequences

Increasing numbers of protein interactions have been identified in high-throughput experiments, but only a small proportion have solved structures. Recently, sequence coevolution-based approaches have led to a breakthrough in predicting monomer protein structures and protein interaction interfaces. Here, we address the challenges of large-scale interaction prediction at residue resolution with a fast alignment concatenation method and a probabilistic score for the interaction of residues. Importantly, this method (EVcomplex2) is able to assess the likelihood of a protein interaction, as we show here applied to large-scale experimental datasets where the pairwise interactions are unknown. We predict 504 interactions de novo in the E. coli membrane proteome, including 243 that are newly discovered. While EVcomplex2 does not require available structures, coevolving residue pairs can be used to produce structural models of protein interactions, as done here for membrane complexes including the Flagellar Hook-Filament Junction and the Tol/Pal complex.

Our understanding of the residue-level details of protein interactions remains incomplete. Here, the authors show sequence coevolution can be used to infer interacting proteins with residue-level details, including predicting 467 interactions de novo in the Escherichia coli cell envelope proteome.

Structural-functional diversity of malaria parasite's PfHSP70-1 and PfHSP40 chaperone pair gives an edge over human orthologs in chaperone-assisted protein folding

Plasmodium falciparum, the human malaria parasite harbors a metastable proteome which is vulnerable to proteotoxic stress conditions encountered during its lifecycle. How parasite's chaperone machinery is able to maintain its aggregation-prone proteome in functional state, is poorly understood. As HSP70-40 system forms the central hub in cellular proteostasis, we investigated the protein folding capacity of PfHSP70-1 and PfHSP40 chaperone pair and compared it with human orthologs (HSPA1A and DNAJA1). Despite the structural similarity, we observed that parasite chaperones and their human orthologs exhibit striking differences in conformational dynamics. Comprehensive biochemical investigations revealed that PfHSP70-1 and PfHSP40 chaperone pair has better protein folding, aggregation inhibition, oligomer remodeling and disaggregase activities than their human orthologs. Chaperone-swapping experiments suggest that PfHSP40 can also efficiently cooperate with human HSP70 to facilitate the folding of client-substrate. SPR-derived kinetic parameters reveal that PfHSP40 has higher binding affinity towards unfolded substrate than DNAJA1. Interestingly, the observed slow dissociation rate of PfHSP40-substrate interaction allows PfHSP40 to maintain the substrate in folding-competent state to minimize its misfolding. Structural investigation through small angle x-ray scattering gave insights into the conformational architecture of PfHSP70-1 (monomer), PfHSP40 (dimer) and their complex. Overall, our data suggest that the parasite has evolved functionally diverged and efficient chaperone machinery which allows the human malaria parasite to survive in hostile conditions. The distinct allosteric landscapes and interaction kinetics of plasmodial chaperones open avenues for the exploration of small-molecule based antimalarial interventions.

Exploiting Folding and Degradation Machineries To Target Undruggable Proteins: What Can a Computational Approach Tell Us?

Advances in genomics and proteomics have unveiled an ever-growing number of key proteins and provided mechanistic insights into the genesis of pathologies. This wealth of data showed that changes in expression levels of specific proteins, mutations, and post-translational modifications can result in (often subtle) perturbations of functional protein-protein interaction networks, which ultimately determine disease phenotypes. Although many such validated pathogenic proteins have emerged as ideal drug targets, there are also several that escape traditional pharmacological regulation; these proteins have thus been labeled "undruggable". The challenges posed by undruggable targets call for new sorts of molecular intervention. One fascinating solution is to perturb a pathogenic protein's expression levels, rather than blocking its activities. In this Concept paper, we shall discuss chemical interventions aimed at recruiting undruggable proteins to the ubiquitin proteasome system, or aimed at disrupting protein-protein interactions in the chaperone-mediated cellular folding machinery: both kinds of intervention lead to a decrease in the amount of active pathogenic protein expressed. Specifically, we shall discuss the role of computational strategies in understanding the molecular determinants characterizing the function of synthetic molecules typically designed for either type of intervention. Finally, we shall provide our perspectives and views on the current limitations and possibilities to expand the scope of rational approaches to the design of chemical regulators of protein levels.

FilterDCA: Interpretable supervised contact prediction using inter-domain coevolution

Predicting three-dimensional protein structure and assembling protein complexes using sequence information belongs to the most prominent tasks in computational biology. Recently substantial progress has been obtained in the case of single proteins using a combination of unsupervised coevolutionary sequence analysis with structurally supervised deep learning. While reaching impressive accuracies in predicting residue-residue contacts, deep learning has a number of disadvantages. The need for large structural training sets limits the applicability to multi-protein complexes; and their deep architecture makes the interpretability of the convolutional neural networks intrinsically hard. Here we introduce FilterDCA, a simpler supervised predictor for inter-domain and inter-protein contacts. It is based on the fact that contact maps of proteins show typical contact patterns, which results from secondary structure and are reflected by patterns in coevolutionary analysis. We explicitly integrate averaged contacts patterns with coevolutionary scores derived by Direct Coupling Analysis, improving performance over standard coevolutionary analysis, while remaining fully transparent and interpretable. The FilterDCA code is available at http://gitlab.lcqb.upmc.fr/muscat/FilterDCA.

The de novo prediction of tertiary and quaternary protein structures has recently seen important advances, by combining unsupervised, purely sequence-based coevolutionary analyses with structure-based supervision using deep learning for contact-map prediction. While showing impressive performance, deep-learning methods require large training sets and pose severe obstacles for their interpretability. Here we construct a simple, transparent and therefore fully interpretable inter-domain contact predictor, which uses the results of coevolutionary Direct Coupling Analysis in combination with explicitly constructed filters reflecting typical contact patterns in a training set of known protein structures, and which improves the accuracy of predicted contacts significantly. Our approach thereby sheds light on the question how contact information is encoded in coevolutionary signals.

Post-translational modifications of Hsp70 family proteins: Expanding the chaperone code

Cells must be able to cope with the challenge of folding newly synthesized proteins and refolding those that have become misfolded in the context of a crowded cytosol. One such coping mechanism that has appeared during evolution is the expression of well-conserved molecular chaperones, such as those that are part of the heat shock protein 70 (Hsp70) family of proteins that bind and fold a large proportion of the proteome. Although Hsp70 family chaperones have been extensively examined for the last 50 years, most studies have focused on regulation of Hsp70 activities by altered transcription, co-chaperone "helper" proteins, and ATP binding and hydrolysis. The rise of modern proteomics has uncovered a vast array of post-translational modifications (PTMs) on Hsp70 family proteins that include phosphorylation, acetylation, ubiquitination, AMPylation, and ADP-ribosylation. Similarly to the pattern of histone modifications, the histone code, this complex pattern of chaperone PTMs is now known as the "chaperone code." In this review, we discuss the history of the Hsp70 chaperone code, its currently understood regulation and functions, and thoughts on what the future of research into the chaperone code may entail.

Two-step mechanism of J-domain action in driving Hsp70 function

J-domain proteins (JDPs), obligatory Hsp70 cochaperones, play critical roles in protein homeostasis. They promote key allosteric transitions that stabilize Hsp70 interaction with substrate polypeptides upon hydrolysis of its bound ATP. Although a recent crystal structure revealed the physical mode of interaction between a J-domain and an Hsp70, the structural and dynamic consequences of J-domain action once bound and how Hsp70s discriminate among its multiple JDP partners remain enigmatic. We combined free energy simulations, biochemical assays and evolutionary analyses to address these issues. Our results indicate that the invariant aspartate of the J-domain perturbs a conserved intramolecular Hsp70 network of contacts that crosses domains. This perturbation leads to destabilization of the domain-domain interface—thereby promoting the allosteric transition that triggers ATP hydrolysis. While this mechanistic step is driven by conserved residues, evolutionarily variable residues are key to initial JDP/Hsp70 recognition—via electrostatic interactions between oppositely charged surfaces. We speculate that these variable residues allow an Hsp70 to discriminate amongst JDP partners, as many of them have coevolved. Together, our data points to a two-step mode of J-domain action, a recognition stage followed by a mechanistic stage.

It is well appreciated that Hsp70-based systems are the most versatile among molecular chaperones—functioning in all cell types and in all subcellular compartments. Via cyclic binding to protein substrates, Hsp70s facilitate their folding, trafficking, degradation and ability to interact with other proteins. Hsp70 function, however, depends on transient interaction with J-domain protein cochaperones that not only deliver substrates, but also activate the structural changes needed for efficient Hsp70 binding to substrate. But how J-domain proteins mechanistically function to drive these changes and how an Hsp70 discriminates among multiple J-domain partners have remained challenging central questions. Here, by using a combination of computational, evolutionary and experimental approaches, we provide evidence for a two-step mechanism. The initial recognition step involves variable residues that allow fine tuning of both the specificity and strength of J-domain protein interaction with Hsp70. The second, that is the mechanistic step, involves conserved residues that act to disrupt a conserved network of intramolecular interactions within Hsp70, thus ensuring robust activation of the structural changes necessary for effective substrate binding. We suggest that our findings are likely applicable to most Hsp70 systems.

Statistical physics of interacting proteins: Impact of dataset size and quality assessed in synthetic sequences

Identifying protein-protein interactions is crucial for a systems-level understanding of the cell. Recently, algorithms based on inverse statistical physics, e.g., direct coupling analysis (DCA), have allowed to use evolutionarily related sequences to address two conceptually related inference tasks: finding pairs of interacting proteins and identifying pairs of residues which form contacts between interacting proteins. Here we address two underlying questions: How are the performances of both inference tasks related? How does performance depend on dataset size and the quality? To this end, we formalize both tasks using Ising models defined over stochastic block models, with individual blocks representing single proteins and interblock couplings protein-protein interactions; controlled synthetic sequence data are generated by Monte Carlo simulations. We show that DCA is able to address both inference tasks accurately when sufficiently large training sets of known interaction partners are available and that an iterative pairing algorithm allows to make predictions even without a training set. Noise in the training data deteriorates performance. In both tasks we find a quadratic scaling relating dataset quality and size that is consistent with noise adding in square-root fashion and signal adding linearly when increasing the dataset. This implies that it is generally good to incorporate more data even if their quality are imperfect, thereby shedding light on the empirically observed performance of DCA applied to natural protein sequences.

Kinetics of the conformational cycle of Hsp70 reveals the importance of the dynamic and heterogeneous nature of Hsp70 for its function

Heat shock protein 70 kDa (Hsp70) plays a central role in maintaining protein homeostasis. It cooperates with cochaperone Hsp40, which stimulates Hsp70 ATPase activity and presents protein substrates to Hsp70 to assist refolding. The mechanism by which Hsp40 regulates the intramolecular and intermolecular changes of Hsp70 is still largely unknown. Here, by bulk and single-molecule FRET, we report the conformational dynamics of Hsp70 and its regulation by Hsp40 as well as the kinetics of the multistep Hsp70–Hsp40 functional cycle. We show that Hsp40 modulates the conformations of ATP-bound Hsp70 to a domain-undocked ATPase-stimulated state, and facilitates the formation of a heterotetrameric Hsp70–Hsp40 complex. Our findings provide insights into the functional mechanism of this core chaperone machinery.

Hsp70 is a conserved molecular chaperone that plays an indispensable role in regulating protein folding, translocation, and degradation. The conformational dynamics of Hsp70 and its regulation by cochaperones are vital to its function. Using bulk and single-molecule fluorescence resonance energy transfer (smFRET) techniques, we studied the interdomain conformational distribution of human stress-inducible Hsp70A1 and the kinetics of conformational changes induced by nucleotide and the Hsp40 cochaperone Hdj1. We found that the conformations between and within the nucleotide- and substrate-binding domains show heterogeneity. The conformational distribution in the ATP-bound state can be induced by Hdj1 to form an “ADP-like” undocked conformation, which is an ATPase-stimulated state. Kinetic measurements indicate that Hdj1 binds to monomeric Hsp70 as the first step, then induces undocking of the two domains and closing of the substrate-binding cleft. Dimeric Hdj1 then facilitates dimerization of Hsp70 and formation of a heterotetrameric Hsp70–Hsp40 complex. Our results provide a kinetic view of the conformational cycle of Hsp70 and reveal the importance of the dynamic nature of Hsp70 for its function.

Coevolutionary Analysis of Protein Subfamilies by Sequence Reweighting

Extracting structural information from sequence co-variation has become a common computational biology practice in the recent years, mainly due to the availability of large sequence alignments of protein families. However, identifying features that are specific to sub-classes and not shared by all members of the family using sequence-based approaches has remained an elusive problem. We here present a coevolutionary-based method to differentially analyze subfamily specific structural features by a continuous sequence reweighting (SR) approach. We introduce the underlying principles and test its predictive capabilities on the Response Regulator family, whose subfamilies have been previously shown to display distinct, specific homo-dimerization patterns. Our results show that this reweighting scheme is effective in assigning structural features known a priori to subfamilies, even when sequence data is relatively scarce. Furthermore, sequence reweighting allows assessing if individual structural contacts pertain to specific subfamilies and it thus paves the way for the identification specificity-determining contacts from sequence variation data.

Engineering repressors with coevolutionary cues facilitates toggle switches with a master reset

Engineering allosteric transcriptional repressors containing an environmental sensing module (ESM) and a DNA recognition module (DRM) has the potential to unlock a combinatorial set of rationally designed biological responses. We demonstrated that constructing hybrid repressors by fusing distinct ESMs and DRMs provides a means to flexibly rewire genetic networks for complex signal processing. We have used coevolutionary traits among LacI homologs to develop a model for predicting compatibility between ESMs and DRMs. Our predictions accurately agree with the performance of 40 engineered repressors. We have harnessed this framework to develop a system of multiple toggle switches with a master OFF signal that produces a unique behavior: each engineered biological activity is switched to a stable ON state by different chemicals and returned to OFF in response to a common signal. One promising application of this design is to develop living diagnostics for monitoring multiple parameters in complex physiological environments and it represents one of many circuit topologies that can be explored with modular repressors designed with coevolutionary information.

Coevolutionary Couplings Unravel PAM-Proximal Constraints of CRISPR-SpCas9

The clustered regularly interspaced short palindromic repeats (CRISPR) system, an immune system analog found in prokaryotes, allows a single-guide RNA to direct a CRISPR-associated protein (Cas) with combined helicase and nuclease activity to DNA. The presence of a specific protospacer adjacent motif (PAM) next to the DNA target site plays a crucial role in determining both efficacy and specificity of gene editing. Herein, we introduce a coevolutionary framework to computationally unveil nonobvious molecular interactions in CRISPR systems and experimentally probe their functional role. Specifically, we use direct coupling analysis, a statistical inference framework used to infer direct coevolutionary couplings, in the context of protein/nucleic acid interactions. Applied to Streptococcus pyogenes Cas9, a Hamiltonian metric obtained from coevolutionary relationships reveals, to our knowledge, novel PAM-proximal nucleotide preferences at the seventh position of S. pyogenes Cas9 PAM (5'-NGRNNNT-3'), which was experimentally confirmed by in vitro and functional assays in human cells. We show that coevolved and conserved interactions point to specific clues toward rationally engineering new generations of Cas9 systems and may eventually help decipher the diversity of this family of proteins.

Structural and functional analysis of the Hsp70/Hsp40 chaperone system

As one of the most abundant and highly conserved molecular chaperones, the 70-kDa heat shock proteins (Hsp70s) play a key role in maintaining cellular protein homeostasis (proteostasis), one of the most fundamental tasks for every living organism. In this role, Hsp70s are inextricably linked to many human diseases, most notably cancers and neurodegenerative diseases, and are increasingly recognized as important drug targets for developing novel therapeutics for these diseases. Hsp40s are a class of essential and universal partners for Hsp70s in almost all aspects of proteostasis. Thus, Hsp70s and Hsp40s together constitute one of the most important chaperone systems across all kingdoms of life. In recent years, we have witnessed significant progress in understanding the molecular mechanism of this chaperone system through structural and functional analysis. This review will focus on this recent progress, mainly from a structural perspective.

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.

The role of coevolutionary signatures in protein interaction dynamics, complex inference, molecular recognition, and mutational landscapes

Evolution imposes constraints at the interface of interacting biomolecules in order to preserve function or maintain fitness. This pressure may have a direct effect on the sequence composition of interacting biomolecules. As a result, statistical patterns of amino acid or nucleotide covariance that encode for physical and functional interactions are observed in sequences of extant organisms. In recent years, global pairwise models of amino acid and nucleotide coevolution from multiple sequence alignments have been developed and utilized to study molecular interactions in structural biology. In proteins, for which the energy landscape is funneled and minimally frustrated, a direct connection between the physical and sequence space landscapes can be established. Estimating coevolutionary information from sequences of interacting molecules has a broad impact in molecular biology. Applications include the accurate determination of 3D structures of molecular complexes, inference of protein interaction partners, models of protein-protein interaction specificity, the elucidation, and design of protein-nucleic acid recognition as well as the discovery of genome-wide epistatic effects. The current state of the art of coevolutionary analysis includes biomedical applications ranging from mutational landscapes and drug-design to vaccine development.

Coevolutionary Analysis of Protein Sequences for Molecular Modeling

Thanks to the explosion of genomic sequencing, coevolutionary analysis of protein sequences has gained great and ever-increasing popularity in the last decade, and it is currently an important and well-established tool in structural bioinformatics and computational biology. This chapter concisely introduces the theoretical foundation and the practical aspects of coevolutionary analysis, as well as discusses the molecular modeling strategies to exploit its results in the study of protein structure, dynamics, and interactions. We present here a complete pipeline from sequence extraction to contact prediction through two examples, focusing on the predictions of inter-residue contacts in a single protein domain and on the analysis of a multi-domain protein that undergoes functional, large-scale conformational transitions.

Function, evolution, and structure of J-domain proteins

Hsp70 chaperone systems are very versatile machines present in nearly all living organisms and in nearly all intracellular compartments. They function in many fundamental processes through their facilitation of protein (re)folding, trafficking, remodeling, disaggregation, and degradation. Hsp70 machines are regulated by co-chaperones. J-domain containing proteins (JDPs) are the largest family of Hsp70 co-chaperones and play a determining role functionally specifying and directing Hsp70 functions. Many features of JDPs are not understood; however, a number of JDP experts gathered at a recent CSSI-sponsored workshop in Gdansk (Poland) to discuss various aspects of J-domain protein function, evolution, and structure. In this report, we present the main findings and the consensus reached to help direct future developments in the field of Hsp70 research.

Mitotic phosphorylation regulates Hsp72 spindle localization by uncoupling ATP binding from substrate release

Hsp72 is a member of the 70-kDa heat shock family of molecular chaperones (Hsp70s) that comprise a nucleotide-binding domain (NBD) and a substrate-binding domain (SBD) connected by a linker that couples the exchange of adenosine diphosphate (ADP) for adenosine triphosphate (ATP) with the release of the protein substrate. Mitotic phosphorylation of Hsp72 by the kinase NEK6 at Thr⁶⁶ located in the NBD promotes the localization of Hsp72 to the mitotic spindle and is required for efficient spindle assembly and chromosome congression and segregation. We determined the crystal structure of the Hsp72 NBD containing a genetically encoded phosphoserine at position 66. This revealed structural changes that stabilized interactions between subdomains within the NBD. ATP binding to the NBD of unmodified Hsp72 resulted in the release of substrate from the SBD, but phosphorylated Hsp72 retained substrate in the presence of ATP. Mutations that prevented phosphorylation-dependent subdomain interactions restored the connection between ATP binding and substrate release. Thus, phosphorylation of Thr⁶⁶ is a reversible mechanism that decouples the allosteric connection between nucleotide binding and substrate release, providing further insight into the regulation of the Hsp70 family. We propose that phosphorylation of Hsp72 on Thr⁶⁶ by NEK6 during mitosis promotes its localization to the spindle by stabilizing its interactions with components of the mitotic spindle.

Emerging fields in chaperone proteins: A French workshop

The "Bioénergétique et Ingénierie des Protéines (BIP)" laboratory, CNRS (France), organized its first French workshop on molecular chaperone proteins and protein folding in November 2017. The goal of this workshop was to gather scientists working in France on chaperone proteins and protein folding. This initiative was a great success with excellent talks and fruitful discussions. The highlights were on the description of unexpected functions and post-translational regulation of known molecular chaperones (such as Hsp90, Hsp33, SecB, GroEL) and on state-of-the-art methods to tackle questions related to this theme, including Cryo-electron microscopy, Nuclear Magnetic Resonance (NMR), Electron Paramagnetic Resonance (EPR), simulation and modeling. We expect to organize a second workshop in two years that will include more scientists working in France in the chaperone field.

Inter-residue, inter-protein and inter-family coevolution: bridging the scales

Interacting proteins coevolve at multiple but interconnected scales, from the residue-residue over the protein-protein up to the family-family level. The recent accumulation of enormous amounts of sequence data allows for the development of novel, data-driven computational approaches. Notably, these approaches can bridge scales within a single statistical framework. Although being currently applied mostly to isolated problems on single scales, their immense potential for an evolutionary informed, structural systems biology is steadily emerging.

Characterization of C-ring component assembly in flagellar motors from amino acid coevolution

Bacterial flagellar motility, an important virulence factor, is energized by a rotary motor localized within the flagellar basal body. The rotor module consists of a large framework (the C-ring), composed of the FliG, FliM and FliN proteins. FliN and FliM contacts the FliG torque ring to control the direction of flagellar rotation. We report that structure-based models constrained only by residue coevolution can recover the binding interface of atomic X-ray dimer complexes with remarkable accuracy (approx. 1 Å RMSD). We propose a model for FliM-FliN heterodimerization, which agrees accurately with homologous interfaces as well as in situ cross-linking experiments, and hence supports a proposed architecture for the lower portion of the C-ring. Furthermore, this approach allowed the identification of two discrete and interchangeable homodimerization interfaces between FliM middle domains that agree with experimental measurements and might be associated with C-ring directional switching dynamics triggered upon binding of CheY signal protein. Our findings provide structural details of complex formation at the C-ring that have been difficult to obtain with previous methodologies and clarify the architectural principle that underpins the ultra-sensitive allostery exhibited by this ring assembly that controls the clockwise or counterclockwise rotation of flagella.

Molecular Mechanism of J-Domain-Triggered ATP Hydrolysis by Hsp70 Chaperones

Efficient targeting of Hsp70 chaperones to substrate proteins depends on J-domain cochaperones, which in synergism with substrates trigger ATP hydrolysis in Hsp70s and concomitant substrate trapping. We present the crystal structure of the J-domain of Escherichia coli DnaJ in complex with the E. coli Hsp70 DnaK. The J-domain interacts not only with DnaK's nucleotide-binding domain (NBD) but also with its substrate-binding domain (SBD) and packs against the highly conserved interdomain linker. Mutational replacement of contacts between J-domain and SBD strongly reduces the ability of substrates to stimulate ATP hydrolysis in the presence of DnaJ and compromises viability at heat shock temperatures. Our data demonstrate that the J-domain and the substrate do not deliver completely independent signals for ATP hydrolysis, but the J-domain, in addition to its direct influence on Hsp70s catalytic center, makes Hsp70 more responsive for the hydrolysis-inducing signal of the substrate, resulting in efficient substrate trapping.

Biomolecular coevolution and its applications: Going from structure prediction toward signaling, epistasis, and function

Evolution leads to considerable changes in the sequence of biomolecules, while their overall structure and function remain quite conserved. The wealth of genomic sequences, the ‘Biological Big Data’, modern sequencing techniques provide allows us to investigate biomolecular evolution with unprecedented detail. Sophisticated statistical models can infer residue pair mutations resulting from spatial proximity. The introduction of predicted spatial adjacencies as constraints in biomolecular structure prediction workflows has transformed the field of protein and RNA structure prediction toward accuracies approaching the experimental resolution limit. Going beyond structure prediction, the same mathematical framework allows mimicking evolutionary fitness landscapes to infer signaling interactions, epistasis, or mutational landscapes.

Electroacupuncture regulates the stress-injury-repair chain of events after cerebral ischemia/reperfusion injury

Inflammation after stroke is the main cause of cerebral ischemia/reperfusion injury. Cascading events after injury can lead to cell death. Heat shock protein 70 and other endogenous injury-signaling molecules are released by damaged cells, which can lead to systemic stress reactions. Protecting the brain through repair begins with the stress-injury-repair signaling chain. This study aimed to verify whether acupuncture acts through this chain to facilitate effective treatment of ischemic stroke. Rat models of cerebral ischemia/reperfusion injury were established by Zea Longa's method, and injury sites were identified by assessing neurological function, 2,3,5-triphenyltetrazolium chloride staining, and hematoxylin-eosin staining. Electroacupuncture at acupoints Baihui (DU20) and Zusanli (ST36) was performed in the model rats with dilatational waves, delivered for 20 minutes a day at 2–100 Hz and an amplitude of 2 mA. We analyzed the blood serum from the rats and found that inflammatory cytokines affected the levels of adrenotrophin and heat shock protein 70, each of which followed a similar bimodal curve. Specifically, electroacupuncture lowered the peak levels of adrenocorticotrophic hormone and heat shock protein 70. Thus, electroacupuncture was able to inhibit excessive stress, reduce inflammation, and promote the repair of neurons, which facilitated healing of ischemic stroke.