Accurate prediction of DnaK-peptide binding via homology modelling and experimental data

Mechanism of chaperone coordination during cotranslational protein folding in bacteria

Protein folding is assisted by molecular chaperones that bind nascent polypeptides during mRNA translation. Several structurally distinct classes of chaperones promote de novo folding, suggesting that their activities are coordinated at the ribosome. We used biochemical reconstitution and structural proteomics to explore the molecular basis for cotranslational chaperone action in bacteria. We found that chaperone binding is disfavored close to the ribosome, allowing folding to precede chaperone recruitment. Trigger factor recognizes compact folding intermediates that expose an extensive unfolded surface, and dictates DnaJ access to nascent chains. DnaJ uses a large surface to bind structurally diverse intermediates and recruits DnaK to sequence-diverse solvent-accessible sites. Neither Trigger factor, DnaJ, nor DnaK destabilize cotranslational folding intermediates. Instead, the chaperones collaborate to protect incipient structure in the nascent polypeptide well beyond the ribosome exit tunnel. Our findings show how the chaperone network selects and modulates cotranslational folding intermediates.

ANXA2 (annexin A2) is crucial to ATG7-mediated autophagy, leading to tumor aggressiveness in triple-negative breast cancer cells

Triple-negative breast cancer (TNBC) is associated with a poor prognosis and metastatic growth. TNBC cells frequently undergo macroautophagy/autophagy, contributing to tumor progression and chemotherapeutic resistance. ANXA2 (annexin A2), a potential therapeutic target for TNBC, has been reported to stimulate autophagy. In this study, we investigated the role of ANXA2 in autophagic processes in TNBC cells. TNBC patients exhibited high levels of ANXA2, which correlated with poor outcomes. ANXA2 increased LC3B-II levels following bafilomycin A1 treatment and enhanced autophagic flux in TNBC cells. Notably, ANXA2 upregulated the phosphorylation of HSF1 (heat shock transcription factor 1), resulting in the transcriptional activation of ATG7 (autophagy related 7). The mechanistic target of rapamycin kinase complex 2 (MTORC2) played an important role in ANXA2-mediated ATG7 transcription by HSF1. MTORC2 did not affect the mRNA level of ANXA2, but it was involved in the protein stability of ANXA2. HSPA (heat shock protein family A (Hsp70)) was a potential interacting protein with ANXA2, which may protect ANXA2 from lysosomal proteolysis. ANXA2 knockdown significantly increased sensitivity to doxorubicin, the first-line chemotherapeutic regimen for TNBC treatment, suggesting that the inhibition of autophagy by ANXA2 knockdown may overcome doxorubicin resistance. In a TNBC xenograft mouse model, we demonstrated that ANXA2 knockdown combined with doxorubicin administration significantly inhibited tumor growth compared to doxorubicin treatment alone, offering a promising avenue to enhance the effectiveness of chemotherapy. In summary, our study elucidated the molecular mechanism by which ANXA2 modulates autophagy, suggesting a potential therapeutic approach for TNBC treatment.Abbreviation: ATG: autophagy related; ChIP: chromatin-immunoprecipitation; HBSS: Hanks' balanced salt solution; HSF1: heat shock transcription factor 1; MTOR: mechanistic target of rapamycin kinase; TNBC: triple-negative breast cancer; TFEB: transcription factor EB; TFE3: transcription factor binding to IGHM enhancer 3.

Protein-folding chaperones predict structure-function relationships and cancer risk in BRCA1 mutation carriers

Predicting the risk of cancer mutations is critical for early detection and prevention, but differences in allelic severity of human carriers confound risk predictions. Here, we elucidate protein folding as a cellular mechanism driving differences in mutation severity of tumor suppressor BRCA1. Using a high-throughput protein-protein interaction assay, we show that protein-folding chaperone binding patterns predict the pathogenicity of variants in the BRCA1 C-terminal (BRCT) domain. HSP70 selectively binds 94% of pathogenic BRCA1-BRCT variants, most of which engage HSP70 more than HSP90. Remarkably, the magnitude of HSP70 binding linearly correlates with loss of folding and function. We identify a prevalent class of human hypomorphic BRCA1 variants that bind moderately to chaperones and retain partial folding and function. Furthermore, chaperone binding signifies greater mutation penetrance and earlier cancer onset in the clinic. Our findings demonstrate the utility of chaperones as quantitative cellular biosensors of variant folding, phenotypic severity, and cancer risk.

Human Hsp70 Substrate-Binding Domains Recognize Distinct Client Proteins

The 13 Hsp70 proteins in humans act on unique sets of substrates with diversity often being attributed to J-domain-containing protein (Hsp40 or JDP) cofactors. We were therefore surprised to find drastically different binding affinities for Hsp70-peptide substrates, leading us to probe substrate specificity among the 8 canonical Hsp70s from humans. We used peptide arrays to characterize Hsp70 binding and then mined these data using machine learning to develop an algorithm for isoform-specific prediction of Hsp70 binding sequences. The results of this algorithm revealed recognition patterns not predicted based on local sequence alignments. We then showed that none of the human isoforms can complement heat-shocked DnaK knockout Escherichia coli cells. However, chimeric Hsp70s consisting of the human nucleotide-binding domain and the substrate-binding domain of DnaK complement during heat shock, providing further evidence in vivo of the divergent function of the Hsp70 substrate-binding domains. We also demonstrated that the differences in heat shock complementation among the chimeras are not due to loss of DnaJ binding. Although we do not exclude JDPs as additional specificity factors, our data demonstrate substrate specificity among the Hsp70s, which has important implications for inhibitor development in cancer and neurodegeneration.

Tandem repeats of highly bioluminescent NanoLuc are refolded noncanonically by the Hsp70 machinery

Chaperones are a large family of proteins crucial for maintaining cellular protein homeostasis. One such chaperone is the 70 kDa heat shock protein (Hsp70), which plays a crucial role in protein (re)folding, stability, functionality, and translocation. While the key events in the Hsp70 chaperone cycle are well established, a relatively small number of distinct substrates were repetitively investigated. This is despite Hsp70 engaging with a plethora of cellular proteins of various structural properties and folding pathways. Here we analyzed novel Hsp70 substrates, based on tandem repeats of NanoLuc (Nluc), a small and highly bioluminescent protein with unique structural characteristics. In previous mechanical unfolding and refolding studies, we have identified interesting misfolding propensities of these Nluc-based tandem repeats. In this study, we further investigate these properties through in vitro bulk experiments. Similar to monomeric Nluc, engineered Nluc dyads and triads proved to be highly bioluminescent. Using the bioluminescence signal as the proxy for their structural integrity, we determined that heat-denatured Nluc dyads and triads can be efficiently refolded by the E. coli Hsp70 chaperone system, which comprises DnaK, DnaJ, and GrpE. In contrast to previous studies with other substrates, we observed that Nluc repeats can be efficiently refolded by DnaK and DnaJ, even in the absence of GrpE co-chaperone. Taken together, our study offers a new powerful substrate for chaperone research and raises intriguing questions about the Hsp70 mechanisms, particularly in the context of structurally diverse proteins.

Protein-Folding Chaperones Predict Structure-Function Relationships and Cancer Risk in BRCA1 Mutation Carriers

Identifying pathogenic mutations and predicting their impact on protein structure, function and phenotype remain major challenges in genome sciences. Protein-folding chaperones participate in structure-function relationships by facilitating the folding of protein variants encoded by mutant genes. Here, we utilize a high-throughput protein-protein interaction assay to test HSP70 and HSP90 chaperone interactions as predictors of pathogenicity for variants in the tumor suppressor BRCA1. Chaperones bind 77% of pathogenic BRCA1-BRCT variants, most of which engaged HSP70 more than HSP90. Remarkably, the magnitude of chaperone binding to variants is proportional to the degree of structural and phenotypic defect induced by BRCA1 mutation. Quantitative chaperone interactions identified BRCA1-BRCT separation-of-function variants and hypomorphic alleles missed by pathogenicity prediction algorithms. Furthermore, increased chaperone binding signified greater cancer risk in human BRCA1 carriers. Altogether, our study showcases the utility of chaperones as quantitative cellular biosensors of variant folding and phenotypic severity.

HIGHLIGHTS

Chaperones detect an abundance of pathogenic folding variants of BRCA1-BRCT.Degree of chaperone binding reflects severity of structural and phenotypic defect.Chaperones identify separation-of-function and hypomorphic variants. Chaperone interactions indicate penetrance and expressivity of BRCA1 alleles.

Incorporating physics to overcome data scarcity in predictive modeling of protein function: A case study of BK channels

Machine learning has played transformative roles in numerous chemical and biophysical problems such as protein folding where large amount of data exists. Nonetheless, many important problems remain challenging for data-driven machine learning approaches due to the limitation of data scarcity. One approach to overcome data scarcity is to incorporate physical principles such as through molecular modeling and simulation. Here, we focus on the big potassium (BK) channels that play important roles in cardiovascular and neural systems. Many mutants of BK channel are associated with various neurological and cardiovascular diseases, but the molecular effects are unknown. The voltage gating properties of BK channels have been characterized for 473 site-specific mutations experimentally over the last three decades; yet, these functional data by themselves remain far too sparse to derive a predictive model of BK channel voltage gating. Using physics-based modeling, we quantify the energetic effects of all single mutations on both open and closed states of the channel. Together with dynamic properties derived from atomistic simulations, these physical descriptors allow the training of random forest models that could reproduce unseen experimentally measured shifts in gating voltage, ∆V1/2, with a RMSE ~ 32 mV and correlation coefficient of R ~ 0.7. Importantly, the model appears capable of uncovering nontrivial physical principles underlying the gating of the channel, including a central role of hydrophobic gating. The model was further evaluated using four novel mutations of L235 and V236 on the S5 helix, mutations of which are predicted to have opposing effects on V1/2 and suggest a key role of S5 in mediating voltage sensor-pore coupling. The measured ∆V1/2 agree quantitatively with prediction for all four mutations, with a high correlation of R = 0.92 and RMSE = 18 mV. Therefore, the model can capture nontrivial voltage gating properties in regions where few mutations are known. The success of predictive modeling of BK voltage gating demonstrates the potential of combining physics and statistical learning for overcoming data scarcity in nontrivial protein function prediction.

Incorporating physics to overcome data scarcity in predictive modeling of protein function: a case study of BK channels

Machine learning has played transformative roles in numerous chemical and biophysical problems such as protein folding where large amount of data exists. Nonetheless, many important problems remain challenging for data-driven machine learning approaches due to the limitation of data scarcity. One approach to overcome data scarcity is to incorporate physical principles such as through molecular modeling and simulation. Here, we focus on the big potassium (BK) channels that play important roles in cardiovascular and neural systems. Many mutants of BK channel are associated with various neurological and cardiovascular diseases, but the molecular effects are unknown. The voltage gating properties of BK channels have been characterized for 473 site-specific mutations experimentally over the last three decades; yet, these functional data by themselves remain far too sparse to derive a predictive model of BK channel voltage gating. Using physics-based modeling, we quantify the energetic effects of all single mutations on both open and closed states of the channel. Together with dynamic properties derived from atomistic simulations, these physical descriptors allow the training of random forest models that could reproduce unseen experimentally measured shifts in gating voltage, ΔV 1/2 , with a RMSE ∼ 32 mV and correlation coefficient of R ∼ 0.7. Importantly, the model appears capable of uncovering nontrivial physical principles underlying the gating of the channel, including a central role of hydrophobic gating. The model was further evaluated using four novel mutations of L235 and V236 on the S5 helix, mutations of which are predicted to have opposing effects on V 1/2 and suggest a key role of S5 in mediating voltage sensor-pore coupling. The measured ΔV 1/2 agree quantitatively with prediction for all four mutations, with a high correlation of R = 0.92 and RMSE = 18 mV. Therefore, the model can capture nontrivial voltage gating properties in regions where few mutations are known. The success of predictive modeling of BK voltage gating demonstrates the potential of combining physics and statistical learning for overcoming data scarcity in nontrivial protein function prediction.

Author Summary

Deep machine learning has brought many exciting breakthroughs in chemistry, physics and biology. These models require large amount of training data and struggle when the data is scarce. The latter is true for predictive modeling of the function of complex proteins such as ion channels, where only hundreds of mutational data may be available. Using the big potassium (BK) channel as a biologically important model system, we demonstrate that a reliable predictive model of its voltage gating property could be derived from only 473 mutational data by incorporating physics-derived features, which include dynamic properties from molecular dynamics simulations and energetic quantities from Rosetta mutation calculations. We show that the final random forest model captures key trends and hotspots in mutational effects of BK voltage gating, such as the important role of pore hydrophobicity. A particularly curious prediction is that mutations of two adjacent residues on the S5 helix would always have opposite effects on the gating voltage, which was confirmed by experimental characterization of four novel mutations. The current work demonstrates the importance and effectiveness of incorporating physics in predictive modeling of protein function with scarce data.

Distribution and solvent exposure of Hsp70 chaperone binding sites across the Escherichia coli proteome

Many proteins must interact with molecular chaperones to achieve their native state in the cell. Yet, how chaperone binding-site characteristics affect the folding process is poorly understood. The ubiquitous Hsp70 chaperone system prevents client-protein aggregation by holding unfolded conformations and by unfolding misfolded states. Hsp70 binding sites of client proteins comprise a nonpolar core surrounded by positively charged residues. However, a detailed analysis of Hsp70 binding sites on a proteome-wide scale is still lacking. Further, it is not known whether proteins undergo some degree of folding while chaperone bound. Here, we begin to address the above questions by identifying Hsp70 binding sites in 2258 Escherichia coli (E. coli) proteins. We find that most proteins bear at least one Hsp70 binding site and that the number of Hsp70 binding sites is directly proportional to protein size. Aggregation propensity upon release from the ribosome correlates with number of Hsp70 binding sites only in the case of large proteins. Interestingly, Hsp70 binding sites are more solvent-exposed than other nonpolar sites, in protein native states. Our findings show that the majority of E. coli proteins are systematically enabled to interact with Hsp70 even if this interaction only takes place during a fraction of the protein lifetime. In addition, our data suggest that some conformational sampling may take place within Hsp70-bound states, due to the solvent exposure of some chaperone binding sites in native proteins. In all, we propose that Hsp70-chaperone-binding traits have evolved to favor Hsp70-assisted protein folding devoid of aggregation.

HSP70-binding motifs function as protein quality control degrons

Protein quality control (PQC) degrons are short protein segments that target misfolded proteins for proteasomal degradation, and thus protect cells against the accumulation of potentially toxic non-native proteins. Studies have shown that PQC degrons are hydrophobic and rarely contain negatively charged residues, features which are shared with chaperone-binding regions. Here we explore the notion that chaperone-binding regions may function as PQC degrons. When directly tested, we found that a canonical Hsp70-binding motif (the APPY peptide) functioned as a dose-dependent PQC degron both in yeast and in human cells. In yeast, Hsp70, Hsp110, Fes1, and the E3 Ubr1 target the APPY degron. Screening revealed that the sequence space within the chaperone-binding region of APPY that is compatible with degron function is vast. We find that the number of exposed Hsp70-binding sites in the yeast proteome correlates with a reduced protein abundance and half-life. Our results suggest that when protein folding fails, chaperone-binding sites may operate as PQC degrons, and that the sequence properties leading to PQC-linked degradation therefore overlap with those of chaperone binding.

Direct observation of chemo-mechanical coupling in DnaK by single-molecule force experiments

Protein allostery requires a communication channel for functional regulation between distal sites within a protein. In the molecular chaperone Hsp70, a two-domain enzyme, the ATP/ADP status of an N-terminal nucleotide-binding domain regulates the substrate affinity of a C-terminal substrate-binding domain. Recently available three-dimensional structures of Hsp70 in ATP/ADP states have provided deep insights into molecular pathways of allosteric signals. However, direct mechanical probing of long-range allosteric coupling between the ATP hydrolysis step and domain states is missing. Using laser optical tweezers, we examined the mechanical properties of a truncated two-domain DnaK(1-552ye) in apo/ADP/ATP- and peptide-bound states. We find that in the apo and ADP states, DnaK domains are mechanically stable and rigid. However, in the ATP state, substrate-binding domain (SBD)∗ye is mechanically destabilized as the result of interdomain docking followed by the unfolding of the α-helical lid. By observing the folding state of the SBD, we could observe the continuous ATP/ADP cycling of the enzyme in real time with a single molecule. The SBD lid closure is strictly coupled to the chemical steps of the ATP hydrolysis cycle even in the presence of peptide substrate.

Characterization of the interactome profiling of Mycoplasma fermentans DnaK in cancer cells reveals interference with key cellular pathways

Chaperone proteins are redundant in nature and, to achieve their function, they bind a large repertoire of client proteins. DnaK is a bacterial chaperone protein that recognizes misfolded and aggregated proteins and drives their folding and intracellular trafficking. Some Mycoplasmas are associated with cancers, and we demonstrated that infection with a strain of Mycoplasma fermentans isolated in our lab promoted lymphoma in a mouse model. Its DnaK is expressed intracellularly in infected cells, it interacts with key proteins to hamper essential pathways related to DNA repair and p53 functions and uninfected cells can take-up extracellular DnaK. We profile here for the first time the eukaryotic proteins interacting with DnaK transiently expressed in five cancer cell lines. A total of 520 eukaryotic proteins were isolated by immunoprecipitation and identified by Liquid Chromatography Mass Spectrometry (LC-MS) analysis. Among the cellular DnaK-binding partners, 49 were shared between the five analyzed cell lines, corroborating the specificity of the interaction of DnaK with these proteins. Enrichment analysis revealed multiple RNA biological processes, DNA repair, chromatin remodeling, DNA conformational changes, protein-DNA complex subunit organization, telomere organization and cell cycle as the most significant ontology terms. This is the first study to show that a bacterial chaperone protein interacts with key eukaryotic components thus suggesting DnaK could become a perturbing hub for the functions of important cellular pathways. Given the close interactions between bacteria and host cells in the local microenvironment, these results provide a foundation for future mechanistic studies on how bacteria interfere with essential cellular processes.

Comparative Study on Feature Selection in Protein Structure and Function Prediction

Many effective methods extract and fuse different protein features to study the relationship between protein sequence, structure, and function, but different methods have preferences in solving the research of protein structure and function, which requires selecting valuable and contributing features to design more effective prediction methods. This work mainly focused on the feature selection methods in the study of protein structure and function, and systematically compared and analyzed the efficiency of different feature selection methods in the prediction of protein structures, protein disorders, protein molecular chaperones, and protein solubility. The results show that the feature selection method based on nonlinear SVM performs best in protein structure prediction, protein solubility prediction, protein molecular chaperone prediction, and protein solubility prediction. After selection, the accuracy of features is improved by 13.16% ~71%, especially the Kmer features and PSSM features of proteins.

Chaperone requirements for de novo folding of Saccharomyces cerevisiae septins

Polymers of septin protein complexes play cytoskeletal roles in eukaryotic cells. The specific subunit composition within complexes controls functions and higher-order structural properties. All septins have globular GTPase domains. The other eukaryotic cytoskeletal NTPases strictly require assistance from molecular chaperones of the cytosol, particularly the cage-like chaperonins, to fold into oligomerization-competent conformations. We previously identified cytosolic chaperones that bind septins and influence the oligomerization ability of septins carrying mutations linked to human disease, but it was unknown to what extent wild-type septins require chaperone assistance for their native folding. Here we use a combination of in vivo and in vitro approaches to demonstrate chaperone requirements for de novo folding and complex assembly by budding yeast septins. Individually purified septins adopted nonnative conformations and formed nonnative homodimers. In chaperonin- or Hsp70-deficient cells, septins folded slower and were unable to assemble posttranslationally into native complexes. One septin, Cdc12, was so dependent on cotranslational chaperonin assistance that translation failed without it. Our findings point to distinct translation elongation rates for different septins as a possible mechanism to direct a stepwise, cotranslational assembly pathway in which general cytosolic chaperones act as key intermediaries.

Computationally-Aided Modeling of Hsp70-Client Interactions: Past, Present, and Future

Hsp70 molecular chaperones play central roles in maintaining a healthy cellular proteome. Hsp70s function by binding to short peptide sequences in incompletely folded client proteins, thus preventing them from misfolding and/or aggregating, and in many cases holding them in a state that is competent for subsequent processes like translocation across membranes. There is considerable interest in predicting the sites where Hsp70s may bind their clients, as the ability to do so sheds light on the cellular functions of the chaperone. In addition, the capacity of the Hsp70 chaperone family to bind to a broad array of clients and to identify accessible sequences that enable discrimination of those that are folded from those that are not fully folded, which is essential to their cellular roles, is a fascinating puzzle in molecular recognition. In this article we discuss efforts to harness computational modeling with input from experimental data to develop a predictive understanding of the promiscuous yet selective binding of Hsp70 molecular chaperones to accessible sequences within their client proteins. We trace how an increasing understanding of the complexities of Hsp70-client interactions has led computational modeling to new underlying assumptions and design features. We describe the trend from purely data-driven analysis toward increased reliance on physics-based modeling that deeply integrates structural information and sequence-based functional data with physics-based binding energies. Notably, new experimental insights are adding to our understanding of the molecular origins of "selective promiscuity" in substrate binding by Hsp70 chaperones and challenging the underlying assumptions and design used in earlier predictive models. Taking the new experimental findings together with exciting progress in computational modeling of protein structures leads us to foresee a bright future for a predictive understanding of selective-yet-promiscuous binding exploited by Hsp70 molecular chaperones; the resulting new insights will also apply to substrate binding by other chaperones and by signaling proteins.

The Impact of Hidden Structure on Aggregate Disassembly by Molecular Chaperones

Protein aggregation, or the uncontrolled self-assembly of partially folded proteins, is an ever-present danger for living organisms. Unimpeded, protein aggregation can result in severe cellular dysfunction and disease. A group of proteins known as molecular chaperones is responsible for dismantling protein aggregates. However, how protein aggregates are recognized and disassembled remains poorly understood. Here we employ a single particle fluorescence technique known as Burst Analysis Spectroscopy (BAS), in combination with two structurally distinct aggregate types grown from the same starting protein, to examine the mechanism of chaperone-mediated protein disaggregation. Using the core bi-chaperone disaggregase system from Escherichia coli as a model, we demonstrate that, in contrast to prevailing models, the overall size of an aggregate particle has, at most, a minor influence on the progression of aggregate disassembly. Rather, we show that changes in internal structure, which have no observable impact on aggregate particle size or molecular chaperone binding, can dramatically limit the ability of the bi-chaperone system to take aggregates apart. In addition, these structural alterations progress with surprising speed, rendering aggregates resistant to disassembly within minutes. Thus, while protein aggregate structure is generally poorly defined and is often obscured by heterogeneous and complex particle distributions, it can have a determinative impact on the ability of cellular quality control systems to process protein aggregates.

Multivalent protein-protein interactions are pivotal regulators of eukaryotic Hsp70 complexes

Heat shock protein 70 (Hsp70) is a molecular chaperone and central regulator of protein homeostasis (proteostasis). Paramount to this role is Hsp70's binding to client proteins and co-chaperones to produce distinct complexes, such that understanding the protein-protein interactions (PPIs) of Hsp70 is foundational to describing its function and dysfunction in disease. Mounting evidence suggests that these PPIs include both "canonical" interactions, which are universally conserved, and "non-canonical" (or "secondary") contacts that seem to have emerged in eukaryotes. These two categories of interactions involve discrete binding surfaces, such that some clients and co-chaperones engage Hsp70 with at least two points of contact. While the contributions of canonical interactions to chaperone function are becoming increasingly clear, it can be challenging to deconvolute the roles of secondary interactions. Here, we review what is known about non-canonical contacts and highlight examples where their contributions have been parsed, giving rise to a model in which Hsp70's secondary contacts are not simply sites of additional avidity but are necessary and sufficient to impart unique functions. From this perspective, we propose that further exploration of non-canonical contacts will generate important insights into the evolution of Hsp70 systems and inspire new approaches for developing small molecules that tune Hsp70-mediated proteostasis.

Universal protein misfolding intermediates can bypass the proteostasis network and remain soluble and less functional

Some misfolded protein conformations can bypass proteostasis machinery and remain soluble in vivo. This is an unexpected observation, as cellular quality control mechanisms should remove misfolded proteins. Three questions, then, are: how do long-lived, soluble, misfolded proteins bypass proteostasis? How widespread are such misfolded states? And how long do they persist? We address these questions using coarse-grain molecular dynamics simulations of the synthesis, termination, and post-translational dynamics of a representative set of cytosolic E. coli proteins. We predict that half of proteins exhibit misfolded subpopulations that bypass molecular chaperones, avoid aggregation, and will not be rapidly degraded, with some misfolded states persisting for months or longer. The surface properties of these misfolded states are native-like, suggesting they will remain soluble, while self-entanglements make them long-lived kinetic traps. In terms of function, we predict that one-third of proteins can misfold into soluble less-functional states. For the heavily entangled protein glycerol-3-phosphate dehydrogenase, limited-proteolysis mass spectrometry experiments interrogating misfolded conformations of the protein are consistent with the structural changes predicted by our simulations. These results therefore provide an explanation for how proteins can misfold into soluble conformations with reduced functionality that can bypass proteostasis, and indicate, unexpectedly, this may be a wide-spread phenomenon.

Heat shock protein Grp78/BiP/HspA5 binds directly to TDP-43 and mitigates toxicity associated with disease pathology

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease with no cure or effective treatment in which TAR DNA Binding Protein of 43 kDa (TDP-43) abnormally accumulates into misfolded protein aggregates in affected neurons. It is widely accepted that protein misfolding and aggregation promotes proteotoxic stress. The molecular chaperones are a primary line of defense against proteotoxic stress, and there has been long-standing interest in understanding the relationship between chaperones and aggregated protein in ALS. Of particular interest are the heat shock protein of 70 kDa (Hsp70) family of chaperones. However, defining which of the 13 human Hsp70 isoforms is critical for ALS has presented many challenges. To gain insight into the specific Hsp70 that modulates TDP-43, we investigated the relationship between TDP-43 and the Hsp70s using proximity-dependent biotin identification (BioID) and discovered several Hsp70 isoforms associated with TDP-43 in the nucleus, raising the possibility of an interaction with native TDP-43. We further found that HspA5 bound specifically to the RNA-binding domain of TDP-43 using recombinantly expressed proteins. Moreover, in a Drosophila strain that mimics ALS upon TDP-43 expression, the mRNA levels of the HspA5 homologue (Hsc70.3) were significantly increased. Similarly we observed upregulation of HspA5 in prefrontal cortex neurons from human ALS patients. Finally, overexpression of HspA5 in Drosophila rescued TDP-43-induced toxicity, suggesting that upregulation of HspA5 may have a compensatory role in ALS pathobiology.

A polypeptide model for toxic aberrant proteins induced by aminoglycoside antibiotics

Aminoglycoside antibiotics interfere with the selection of cognate tRNAs during translation, resulting in the synthesis of aberrant proteins that are the ultimate cause of cell death. However, the toxic potential of aberrant proteins and how they avoid degradation by the cell’s protein quality control (QC) machinery are not understood. Here we report that levels of the heat shock (HS) transcription factor σ32 increased sharply following exposure of Escherichia coli to the aminoglycoside kanamycin (Kan), suggesting that at least some of the aberrant proteins synthesized in these cells were recognized as substrates by DnaK, a molecular chaperone that regulates the HS response, the major protein QC pathway in bacteria. To further investigate aberrant protein toxic potential and interaction with cell QC factors, we studied an acutely toxic 48-residue polypeptide (ARF48) that is encoded by an alternate reading frame in a plant cDNA. As occurred in cells exposed to Kan, σ32 levels were strongly elevated following ARF48 expression, suggesting that ARF48 was recognized as a substrate by DnaK. Paradoxically, an internal 10-residue region that was tightly bound by DnaK in vitro also was required for the ARF48 toxic effect. Despite the increased levels of σ32, levels of several HS proteins were unchanged following ARF48 expression, suggesting that the HS response had been aborted. Nucleoids were condensed and cell permeability increased rapidly following ARF48 expression, together suggesting that ARF48 disrupts DNA-membrane interactions that could be required for efficient gene expression. Our results are consistent with earlier studies showing that aberrant proteins induced by aminoglycoside antibiotics disrupt cell membrane integrity. Insights into the mechanism for this effect could be gained by further study of the ARF48 model system.

Active unfolding of the glucocorticoid receptor by the Hsp70/Hsp40 chaperone system in single-molecule mechanical experiments

One of the key unresolved questions in the field of molecular chaperones is how they can actively unfold proteins. In this study, we discovered that the Hsp70/Hsp40 chaperone system completely unfolds a native soluble substrate protein, the ligand-binding domain of the glucocorticoid receptor, in a concerted action. Our high-resolution optical tweezers data show in real time how the substrate is attacked by the chaperone machinery. As soon as the hormone has left the binding pocket, up to five Hsp70/Hsp40 complexes bind and unfold the protein in a stepwise manner. This finding constitutes direct evidence that the chaperone machinery can bind to the folded core of the receptor, thus providing a mechanism for Hsp70-induced protein unfolding.

The glucocorticoid receptor (GR) is an important transcription factor and drug target linked to a variety of biological functions and diseases. It is one of the most stringent physiological clients of the Hsp90/Hsp70/Hsp40 chaperone system. In this study, we used single-molecule force spectroscopy by optical tweezers to observe the interaction of the GR’s ligand-binding domain (GR-LBD) with the Hsp70/Hsp40 chaperone system (Hsp70/40). We show in real time that Hsp70/40 can unfold the complete GR-LBD in a stepwise manner. Each unfolding step involves binding of an Hsp70 to the GR-LBD and subsequent adenosine triphosphate (ATP) hydrolysis, stimulated by Hsp40. The kinetics of chaperone-mediated unfolding depend on chaperone concentrations as well as the presence of the nucleotide exchange factor BAG1. We find that Hsp70/40 can stabilize new unfolding intermediates, showing that Hsp70/40 can directly interact with the folded core of the protein when working as an unfoldase. Our results support an unfolding mechanism where Hsp70 can directly bind to folded protein structures and unfold them upon ATP hydrolysis. These results provide important insights into the regulation of GR by Hsp70/40.

The E. coli DnaK chaperone stimulates the α-complementation of β-galactosidase

pUC18 and pUC19 are well-known high copy-number plasmid vectors routinely used for DNA cloning purposes. We show here that, in E. coli transformed by native pUC18, the α-complementation of β-galactosidase (i.e., mediated by the peptide LacZα18) is intrinsically weak and slow, but is greatly stimulated by the DnaK/DnaJ/GrpE chaperone system. In contrast, the α-complementation mediated by the peptide LacZα19 (in E. coli transformed by the native pUC19) is much more efficient, and therefore does not require the assistance of the DnaK chaperone machinery. The marked difference between these two LacZα peptides is reproduced in cell-free protein expression system coupled with α-complementation. We conclude that: (i) α-complementation of β-galactosidase is DnaK-mediated depending upon the LacZα peptide donor. (ii) DnaK, sensu stricto, is not necessary for α-complementation, but can enhance it to a great extent. (iii) this observation could be used to establish an easy and inexpensive method for screening small molecules libraries in search of DnaK inhibitors and also for deciphering the DnaK-mediated protein quality control mechanism. This article is protected by copyright. All rights reserved.

Physics-based modeling provides predictive understanding of selectively promiscuous substrate binding by Hsp70 chaperones

To help cells cope with protein misfolding and aggregation, Hsp70 molecular chaperones selectively bind a variety of sequences ("selective promiscuity"). Statistical analyses from substrate-derived peptide arrays reveal that DnaK, the E. coli Hsp70, binds to sequences containing three to five branched hydrophobic residues, although otherwise the specific amino acids can vary considerably. Several high-resolution structures of the substrate -binding domain (SBD) of DnaK bound to peptides reveal a highly conserved configuration of the bound substrate and further suggest that the substrate-binding cleft consists of five largely independent sites for interaction with five consecutive substrate residues. Importantly, both substrate backbone orientations (N- to C- and C- to N-) allow essentially the same backbone hydrogen-bonding and side-chain interactions with the chaperone. In order to rationalize these observations, we performed atomistic molecular dynamics simulations to sample the interactions of all 20 amino acid side chains in each of the five sites of the chaperone in the context of the conserved substrate backbone configurations. The resulting interaction energetics provide the basis set for deriving a predictive model that we call Paladin (Physics-based model of DnaK-Substrate Binding). Trained using available peptide array data, Paladin can distinguish binders and nonbinders of DnaK with accuracy comparable to existing predictors and further predicts the detailed configuration of the bound sequence. Tested using existing DnaK-peptide structures, Paladin correctly predicted the binding register in 10 out of 13 substrate sequences that bind in the N- to C- orientation, and the binding orientation in 16 out of 22 sequences. The physical basis of the Paladin model provides insight into the origins of how Hsp70s bind substrates with a balance of selectivity and promiscuity. The approach described here can be extended to other Hsp70s where extensive peptide array data is not available.

Computational algorithmic and molecular dynamics study of functional and structural impacts of non-synonymous single nucleotide polymorphisms in human DHFR gene

Human dihydrofolate reductase (DHFR) is a conserved enzyme that is central to folate metabolism and is widely targeted in pathogenic diseases as well as cancers. Although studies have reported the fact that genetic mutations in DHFR leads to a rare autosomal recessive inborn error of folate metabolism and drug resistance, there is a lack of an extensive study on how the deleterious non-synonymous SNPs (nsSNPs) disrupt its phenotypic effects. In this study, we aim at discovering the structural and functional consequences of nsSNPs in DHFR by employing a combined computational approach consisting of ten recently developed in silico tools for identification of damaging nsSNPs and molecular dynamics (MD) simulation for getting deeper insights into the magnitudes of damaging effects. Our study revealed the presence of 12 most deleterious nsSNPs affecting the native phenotypic effects, with three (R71T, G118D, Y122D) identified in the co-factor and ligand binding active sites. MD simulations also suggested that these three SNPs particularly Y122D, alter the overall structural flexibility and dynamics of the native DHFR protein which can provide more understandings into the crucial roles of these mutants in influencing the loss of DHFR function.

A dual-reporter system for investigating and optimizing protein translation and folding in E. coli

Strategies for investigating and optimizing the expression and folding of proteins for biotechnological and pharmaceutical purposes are in high demand. Here, we describe a dual-reporter biosensor system that simultaneously assesses in vivo protein translation and protein folding, thereby enabling rapid screening of mutant libraries. We have validated the dual-reporter system on five different proteins and find an excellent correlation between reporter signals and the levels of protein expression and solubility of the proteins. We further demonstrate the applicability of the dual-reporter system as a screening assay for deep mutational scanning experiments. The system enables high throughput selection of protein variants with high expression levels and altered protein stability. Next generation sequencing analysis of the resulting libraries of protein variants show a good correlation between computationally predicted and experimentally determined protein stabilities. We furthermore show that the mutational experimental data obtained using this system may be useful for protein structure calculations.

Selective promiscuity in the binding of E. coli Hsp70 to an unfolded protein

Heat shock protein 70 (Hsp70) chaperones bind many different sequences and discriminate between incompletely folded and folded clients. Most research into the origins of this "selective promiscuity" has relied on short peptides as substrates to dissect the binding, but much less is known about how Hsp70s bind full-length client proteins. Here, we connect detailed structural analyses of complexes between the Escherichia coli Hsp70 (DnaK) substrate-binding domain (SBD) and peptides encompassing five potential binding sites in the precursor to E. coli alkaline phosphatase (proPhoA) with SBD binding to full-length unfolded proPhoA. Analysis of SBD complexes with proPhoA peptides by a combination of X-ray crystallography, methyl-transverse relaxation optimized spectroscopy (methyl-TROSY), and paramagnetic relaxation enhancement (PRE) NMR and chemical cross-linking experiments provided detailed descriptions of their binding modes. Importantly, many sequences populate multiple SBD binding modes, including both the canonical N to C orientation and a C to N orientation. The favored peptide binding mode optimizes substrate residue side-chain compatibility with the SBD binding pockets independent of backbone orientation. Relating these results to the binding of the SBD to full-length proPhoA, we observe that multiple chaperones may bind to the protein substrate, and the binding sites, well separated in the proPhoA sequence, behave independently. The hierarchy of chaperone binding to sites on the protein was generally consistent with the apparent binding affinities observed for the peptides corresponding to these sites. Functionally, these results reveal that Hsp70s "read" sequences without regard to the backbone direction and that both binding orientations must be considered in current predictive algorithms.

Protein structure and aggregation: a marriage of necessity ruled by aggregation gatekeepers

Protein aggregation propensity is a pervasive and seemingly inescapable property of proteomes. Strikingly, a significant fraction of the proteome is supersaturated, meaning that, for these proteins, their native conformation is less stable than the aggregated state. Maintaining the integrity of a proteome under such conditions is precarious and requires energy-consuming proteostatic regulation. Why then is aggregation propensity maintained at such high levels over long evolutionary timescales? Here, we argue that the conformational stability of the native and aggregated states are correlated thermodynamically and that codon usage strengthens this correlation. As a result, the folding of stable proteins requires kinetic control to avoid aggregation, provided by aggregation gatekeepers. These unique residues are evolutionarily selected to kinetically favor native folding, either on their own or by coopting chaperones.

Using CHOU'S 5-Steps Rule to Predict O-Linked Serine Glycosylation Sites by Blending Position Relative Features and Statistical Moment

Glycosylation of proteins in eukaryote cells is an important and complicated post-translation modification due to its pivotal role and association with crucial physiological functions within most of the proteins. Identification of glycosylation sites in a polypeptide chain is not an easy task due to multiple impediments. Analytical identification of these sites is expensive and laborious. There is a dire need to develop a reliable computational method for precise determination of such sites which can help researchers to save time and effort. Herein, we propose a novel predictor namely iGlycoS-PseAAC by integrating the Chou's Pseudo Amino Acid Composition (PseAAC) and relative/absolute position-based features. The self-consistency results show that the accuracy revealed by the model using the benchmark dataset for prediction of O-linked glycosylation having serine sites is 98.8 percent. The overall accuracy of predictor achieved through 10-fold cross validation by combining the positive and negative results is 97.2 percent. The overall accuracy achieved through Jackknife test is 96.195 percent by aggregating of all the prediction results. Thus the proposed predictor can help in predicting the O-linked glycosylated serine sites in an efficient and accurate way. The overall results show that the accuracy of the iGlycoS-PseAAC is higher than the existing tools.

Synergistic stabilization of a double mutant in chymotrypsin inhibitor 2 from a library screen in E. coli

Most single point mutations destabilize folded proteins. Mutations that stabilize a protein typically only have a small effect and multiple mutations are often needed to substantially increase the stability. Multiple point mutations may act synergistically on the stability, and it is often not straightforward to predict their combined effect from the individual contributions. Here, we have applied an efficient in-cell assay in E. coli to select variants of the barley chymotrypsin inhibitor 2 with increased stability. We find two variants that are more than 3.8 kJ mol-1 more stable than the wild-type. In one case, the increased stability is the effect of the single substitution D55G. The other case is a double mutant, L49I/I57V, which is 5.1 kJ mol-1 more stable than the sum of the effects of the individual mutations. In addition to demonstrating the strength of our selection system for finding stabilizing mutations, our work also demonstrate how subtle conformational effects may modulate stability.

The Proteome Folding Problem and Cellular Proteostasis

Stunning advances have been achieved in addressing the protein folding problem, providing deeper understanding of the mechanisms by which proteins navigate energy landscapes to reach their native states and enabling powerful algorithms to connect sequence to structure. However, the realities of the in vivo protein folding problem remain a challenge to reckon with. Here, we discuss the concept of the "proteome folding problem"-the problem of how organisms build and maintain a functional proteome-by admitting that folding energy landscapes are characterized by many misfolded states and that cells must deploy a network of chaperones and degradation enzymes to minimize deleterious impacts of these off-pathway species. The resulting proteostasis network is an inextricable part of in vivo protein folding and must be understood in detail if we are to solve the proteome folding problem. We discuss how the development of computational models for the proteostasis network's actions and the relationship to the biophysical properties of the proteome has begun to offer new insights and capabilities.

Hsp70-mediated quality control: should I stay or should I go?

Chaperones of the 70 kDa heat shock protein (Hsp70) superfamily are key components of the cellular proteostasis system. Together with its co-chaperones, Hsp70 forms proteostasis subsystems that antagonize protein damage during physiological and stress conditions. This function stems from highly regulated binding and release cycles of protein substrates, which results in a flow of unfolded, partially folded and misfolded species through the Hsp70 subsystem. Specific factors control how Hsp70 makes decisions regarding folding and degradation fates of the substrate proteins. In this review, we summarize how the flow of Hsp70 substrates is controlled in the cell with special emphasis on recent advances regarding substrate release mechanisms.

Function, Therapeutic Potential, and Inhibition of Hsp70 Chaperones

Hsp70s are among the most highly conserved proteins in all of biology. Through an iterative binding and release of exposed hydrophobic residues on client proteins, Hsp70s can prevent aggregation and promote folding to the native state of their client proteins. The human proteome contains eight canonical Hsp70s. Because Hsp70s are relatively promiscuous they play a role in folding a large proportion of the proteome. Hsp70s are implicated in disease through their ability to regulate protein homeostasis. In recent years, researchers have attempted to develop selective inhibitors of Hsp70 isoforms to better understand the role of individual isoforms in biology and as potential therapeutics. Selective inhibitors have come from rational design, forced localization, and serendipity, but the development of completely selective inhibitors remains elusive. In the present review, we discuss the Hsp70 structure and function, the known Hsp70 client proteins, the role of Hsp70s in disease, and current efforts to discover Hsp70 modulators.

The functions and regulation of heat shock proteins; key orchestrators of proteostasis and the heat shock response

Cells respond to protein-damaging (proteotoxic) stress by activation of the Heat Shock Response (HSR). The HSR provides cells with an enhanced ability to endure proteotoxic insults and plays a crucial role in determining subsequent cell death or survival. The HSR is, therefore, a critical factor that influences the toxicity of protein stress. While named for its vital role in the cellular response to heat stress, various components of the HSR system and the molecular chaperone network execute essential physiological functions as well as responses to other diverse toxic insults. The effector molecules of the HSR, the Heat Shock Factors (HSFs) and Heat Shock Proteins (HSPs), are also important regulatory targets in the progression of neurodegenerative diseases and cancers. Modulation of the HSR and/or its extended network have, therefore, become attractive treatment strategies for these diseases. Development of effective therapies will, however, require a detailed understanding of the HSR, important features of which continue to be uncovered and are yet to be completely understood. We review recently described and hallmark mechanistic principles of the HSR, the regulation and functions of HSPs, and contexts in which the HSR is activated and influences cell fate in response to various toxic conditions.

Mapping the degradation pathway of a disease-linked aspartoacylase variant

Canavan disease is a severe progressive neurodegenerative disorder that is characterized by swelling and spongy degeneration of brain white matter. The disease is genetically linked to polymorphisms in the aspartoacylase (ASPA) gene, including the substitution C152W. ASPA C152W is associated with greatly reduced protein levels in cells, yet biophysical experiments suggest a wild-type like thermal stability. Here, we use ASPA C152W as a model to investigate the degradation pathway of a disease-causing protein variant. When we expressed ASPA C152W in Saccharomyces cerevisiae, we found a decreased steady state compared to wild-type ASPA as a result of increased proteasomal degradation. However, molecular dynamics simulations of ASPA C152W did not substantially deviate from wild-type ASPA, indicating that the native state is structurally preserved. Instead, we suggest that the C152W substitution interferes with the de novo folding pathway resulting in increased proteasomal degradation before reaching its stable conformation. Systematic mapping of the protein quality control components acting on misfolded and aggregation-prone species of C152W, revealed that the degradation is highly dependent on the molecular chaperone Hsp70, its co-chaperone Hsp110 as well as several quality control E3 ubiquitin-protein ligases, including Ubr1. In addition, the disaggregase Hsp104 facilitated refolding of aggregated ASPA C152W, while Cdc48 mediated degradation of insoluble ASPA protein. In human cells, ASPA C152W displayed increased proteasomal turnover that was similarly dependent on Hsp70 and Hsp110. Our findings underscore the use of yeast to determine the protein quality control components involved in the degradation of human pathogenic variants in order to identify potential therapeutic targets.

Hsp70 Inhibits Aggregation of IAPP by Binding to the Heterogeneous Prenucleation Oligomers

Molecular chaperone Hsp70 plays important roles in the pathology of amyloid diseases by inhibiting aberrant aggregation of proteins. However, the biophysical mechanism of the interaction of Hsp70 with the intrinsically disordered proteins (IDPs) is unclear. Here, we report that Hsp70 inhibits aggregation of islet amyloid polypeptide (IAPP) at substoichiometric concentrations under diverse solution conditions, including in the absence of ATP. The inhibitory effect is strongest if Hsp70 is added in the beginning of aggregation but progressively less if added later, indicating a role for Hsp70 in preventing nucleation of IAPP. However, ensemble measurement of the binding affinity suggests poor interactions between Hsp70 and IAPP. Therefore, we hypothesize that the interaction must involve a rare species (e.g., the oligomeric intermediates of IAPP). Size exclusion chromatography and field flow fractionation are then used to fractionate the constituent species. Multiangle light scattering and fluorescence correlation spectroscopy measurements indicate that the dominant fraction in size exclusion chromatography contains a few nanomolar Hsp70-IAPP complexes amid several μmoles of free Hsp70. Using single-particle two-color coincidence detection measurements, we detected a minor fraction that exhibits fluorescence bursts arising from heterogeneous oligomeric complexes of IAPP and Hsp70. Taken together, our results indicate that Hsp70 interacts poorly with the monomers but strongly with oligomers of IAPP. This is likely a generic feature of the interactions of Hsp70 chaperones with the amyloidogenic IDPs. Whereas high-affinity interactions with the oligomers prevent aberrant aggregation, poor interaction with the monomers averts interference with the physiological functions of the IDPs.

In Silico Tools and Approaches for the Prediction of Functional and Structural Effects of Single-Nucleotide Polymorphisms on Proteins: An Expert Review

Single-nucleotide polymorphisms (SNPs) are single-base variants that contribute to human biological variation and pathogenesis of many human diseases. Among all SNP types, nonsynonymous single-nucleotide polymorphisms (nsSNPs) can alter many structural, biochemical, and functional features of a protein such as folding characteristics, charge distribution, stability, dynamics, and interactions with other proteins/nucleotides. These modifications in the protein structure can lead nsSNPs to be closely associated with many multifactorial diseases such as cancer, diabetes, and neurodegenerative diseases. Predicting structural and functional effects of nsSNPs with experimental approaches can be time-consuming and costly; hence, computational prediction tools and algorithms are being widely and increasingly utilized in biology and medical research. This expert review examines the in silico tools and algorithms for the prediction of functional or structural effects of SNP variants, in addition to the description of the phenotypic effects of nsSNPs on protein structure, association between pathogenicity of variants, and functional or structural features of disease-associated variants. Finally, case studies investigating the functional and structural effects of nsSNPs on selected protein structures are highlighted. We conclude that creating a consistent workflow with a combination of in silico approaches or tools should be considered to increase the performance, accuracy, and precision of the biological and clinical predictions made in silico.

Exposure of a cryptic Hsp70 binding site determines the cytotoxicity of the ALS-associated SOD1-mutant A4V

The accumulation of toxic protein aggregates is thought to play a key role in a range of degenerative pathologies, but it remains unclear why aggregation of polypeptides into non-native assemblies is toxic and why cellular clearance pathways offer ineffective protection. We here study the A4V mutant of SOD1, which forms toxic aggregates in motor neurons of patients with familial amyotrophic lateral sclerosis (ALS). A comparison of the location of aggregation prone regions (APRs) and Hsp70 binding sites in the denatured state of SOD1 reveals that ALS-associated mutations promote exposure of the APRs more than the strongest Hsc/Hsp70 binding site that we could detect. Mutations designed to increase the exposure of this Hsp70 interaction site in the denatured state promote aggregation but also display an increased interaction with Hsp70 chaperones. Depending on the cell type, in vitro this resulted in cellular inclusion body formation or increased clearance, accompanied with a suppression of cytotoxicity. The latter was also observed in a zebrafish model in vivo. Our results suggest that the uncontrolled accumulation of toxic SOD1A4V aggregates results from insufficient detection by the cellular surveillance network.

Autonomous aggregation suppression by acidic residues explains why chaperones favour basic residues

Many chaperones favour binding to hydrophobic sequences that are flanked by basic residues while disfavouring acidic residues. However, the origin of this bias in protein quality control remains poorly understood. Here, we show that while acidic residues are the most efficient aggregation inhibitors, they are also less compatible with globular protein structure than basic amino acids. As a result, while acidic residues allow for chaperone-independent control of aggregation, their use is structurally limited. Conversely, we find that, while being more compatible with globular structure, basic residues are not sufficient to autonomously suppress protein aggregation. Using Hsp70, we show that chaperones with a bias towards basic residues are structurally adapted to prioritize aggregating sequences whose structural context forced the use of the less effective basic residues. The hypothesis that emerges from our analysis is that the bias of many chaperones for basic residues results from fundamental thermodynamic and kinetic constraints of globular structure. This also suggests the co-evolution of basic residues and chaperones allowed for an expansion of structural variety in the protein universe.

KLR-70: A Novel Cationic Inhibitor of the Bacterial Hsp70 Chaperone

The heat-shock factor Hsp70 and other molecular chaperones play a central role in nascent protein folding. Elucidating the task performed by individual chaperones within the complex cellular milieu, however, has been challenging. One strategy for addressing this goal has been to monitor protein biogenesis in the absence and presence of inhibitors of a specific chaperone, followed by analysis of folding outcomes under both conditions. In this way, the role of the chaperone of interest can be discerned. However, development of chaperone inhibitors, including well-known proline-rich antimicrobial peptides, has been fraught with undesirable side effects, including decreased protein expression yields. Here, we introduce KLR-70, a rationally designed cationic inhibitor of the Escherichia coli Hsp70 chaperone (also known as DnaK). KLR-70 is a 14-amino acid peptide bearing naturally occurring residues and engineered to interact with the DnaK substrate-binding domain. The interaction of KLR-70 with DnaK is enantioselective and is characterized by high affinity in a buffered solution. Importantly, KLR-70 does not significantly interact with the DnaJ and GroEL/ES chaperones, and it does not alter nascent protein biosynthesis yields across a wide concentration range. Some attenuation of the anti-DnaK activity of KLR-70, however, has been observed in the complex E. coli cell-free environment. Interestingly, the d enantiomer D-KLR-70, unlike its all-L KLR-70 counterpart, does not bind the DnaK and DnaJ chaperones, yet it strongly inhibits translation. This outcome suggests that the two enantiomers (KLR-70 and D-KLR-70) may serve as orthogonal inhibitors of chaperone binding and translation. In summary, KLR-70 is a novel chaperone inhibitor with high affinity and selectivity for bacterial Hsp70 and with considerable potential to help in parsing out the role of Hsp70 in nascent protein folding.

TTRMDB: A database for structural and functional analysis on the impact of SNPs over transthyretin (TTR) using bioinformatic tools

Hereditary Transthyretin-associated amyloidosis (ATTR) is an autosomal dominant protein-folding disorder with adult-onset caused by mutation of transthyretin (TTR). TTR is characterized by extracellular deposition of amyloid, leading to loss of autonomy and finally, death. More than 100 distinct mutations in TTR gene have been reported from variable age of onset, clinical expression and penetrance data. Besides, the cure for the disease remains still obscure. Further, the prioritizing of mutations concerning the characteristic features governing the stability and pathogenicity of TTR mutant proteins remains unanswered, to date and thus, a complex state of study for researchers. Herein, we provide a full report encompassing the effects of every reported mutant model of TTR protein about the stability, functionality and pathogenicity using various computational tools. In addition, the results obtained from our study were used to create TTRMDB (Transthyretin mutant database), which could be easy access to researchers at http://vit.ac.in/ttrmdb.

A Systematic and Comprehensive Review on Disease-Causing Genes in Amyotrophic Lateral Sclerosis

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disorder and is characterized by degeneration and axon loss from the upper motor neuron, that descends from the lower motor neuron in the brain. Over the period, assorted outcomes from medical findings, molecular pathogenesis, and structural and biophysical studies have abetted in providing thoughtful insights underlying the importance of disease-causing genes in ALS. Consequently, numerous mechanisms were proposed for the pathogenesis of ALS, considering protein mutations, aggregation, and misfolding. Besides, the answers to the majority of ALS cases that happen to be sporadic still remain obscure. The application in discovering susceptibility factors in ALS contemplating the genetic factors is to be further dissevered in the future years with innovation in research studies. Hence, this review targets in revisiting the breakthroughs on the disease-causing genes related with ALS.

Selective Binding of HSC70 and its Co-Chaperones to Structural Hotspots on CFTR

Mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) channel cause cystic fibrosis. Chaperones, including HSC70, DNAJA1 and DNAJA2, play key roles in both the folding and degradation of wild-type and mutant CFTR at multiple cellular locations. DNAJA1 and HSC70 promote the folding of newly synthesized CFTR at the endoplasmic reticulum (ER), but are required for the rapid turnover of misfolded channel at the plasma membrane (PM). DNAJA2 and HSC70 are also involved in the ER-associated degradation (ERAD) of misfolded CFTR, while they assist the refolding of destabilized channel at the PM. These outcomes may depend on the binding of chaperones to specific sites within CFTR, which would be exposed in non-native states. A CFTR peptide library was used to identify binding sites for HSC70, DNAJA1 and DNAJA2, validated by competition and functional assays. Each chaperone had a distinct binding pattern, and sites were distributed between the surfaces of the CFTR cytosolic domains, and domain interfaces known to be important for channel assembly. The accessibility of sites to chaperones will depend on the degree of CFTR folding or unfolding. Different folded states may be recognized by unique combinations of HSC70, DNAJA1 and DNAJA2, leading to divergent biological effects.

Mutagenesis-Based Characterization and Improvement of a Novel Inclusion Body Tag

Whereas, bacterial inclusion bodies (IBs) for long were regarded as undesirable aggregates emerging during recombinant protein production, they currently receive attention as promising nanoparticulate biomaterials with diverse applications in biotechnology and biomedicine. We previously identified ssTorA, a signal sequence that normally directs protein export via the Tat pathway in E. coli, as a tag that induces the accumulation of fused proteins into IBs under overexpression conditions. Here, we used targeted mutagenesis to identify features and motifs being either critical or dispensable for IB formation. We found that IB formation is neither related to the function of ssTorA as a Tat-signal sequence nor is it a general feature of this family of signal sequences. IB formation was inhibited by co-overexpression of ssTorA binding chaperones TorD and DnaK and by amino acid substitutions that affect the propensity of ssTorA to form an α-helix. Systematic deletion experiments identified a minimal region of ssTorA required for IB formation in the center of the signal sequence. Unbiased genetic screening of a library of randomly mutagenized ssTorA sequences for reduced aggregation properties allowed us to pinpoint residues that are critical to sustain insoluble expression. Together, the data point to possible mechanisms for the aggregation of ssTorA fusions. Additionally, they led to the design of a tag with superior IB-formation properties compared to the original ssTorA sequence.

Bacterial Hsp70 resolves misfolded states and accelerates productive folding of a multi-domain protein

The ATP-dependent Hsp70 chaperones (DnaK in E. coli) mediate protein folding in cooperation with J proteins and nucleotide exchange factors (E. coli DnaJ and GrpE, respectively). The Hsp70 system prevents protein aggregation and increases folding yields. Whether it also enhances the rate of folding remains unclear. Here we show that DnaK/DnaJ/GrpE accelerate the folding of the multi-domain protein firefly luciferase (FLuc) ~20-fold over the rate of spontaneous folding measured in the absence of aggregation. Analysis by single-pair FRET and hydrogen/deuterium exchange identified inter-domain misfolding as the cause of slow folding. DnaK binding expands the misfolded region and thereby resolves the kinetically-trapped intermediates, with folding occurring upon GrpE-mediated release. In each round of release DnaK commits a fraction of FLuc to fast folding, circumventing misfolding. We suggest that by resolving misfolding and accelerating productive folding, the bacterial Hsp70 system can maintain proteins in their native states under otherwise denaturing stress conditions.

The Impact of Protein Architecture on Adaptive Evolution

Adaptive mutations play an important role in molecular evolution. However, the frequency and nature of these mutations at the intramolecular level are poorly understood. To address this, we analyzed the impact of protein architecture on the rate of adaptive substitutions, aiming to understand how protein biophysics influences fitness and adaptation. Using Drosophila melanogaster and Arabidopsis thaliana population genomics data, we fitted models of distribution of fitness effects and estimated the rate of adaptive amino-acid substitutions both at the protein and amino-acid residue level. We performed a comprehensive analysis covering genome, gene, and protein structure, by exploring a multitude of factors with a plausible impact on the rate of adaptive evolution, such as intron number, protein length, secondary structure, relative solvent accessibility, intrinsic protein disorder, chaperone affinity, gene expression, protein function, and protein-protein interactions. We found that the relative solvent accessibility is a major determinant of adaptive evolution, with most adaptive mutations occurring at the surface of proteins. Moreover, we observe that the rate of adaptive substitutions differs between protein functional classes, with genes encoding for protein biosynthesis and degradation signaling exhibiting the fastest rates of protein adaptation. Overall, our results suggest that adaptive evolution in proteins is mainly driven by intermolecular interactions, with host-pathogen coevolution likely playing a major role.

Successes and challenges in simulating the folding of large proteins

Computational simulations of protein folding can be used to interpret experimental folding results, to design new folding experiments, and to test the effects of mutations and small molecules on folding. However, whereas major experimental and computational progress has been made in understanding how small proteins fold, research on larger, multidomain proteins, which comprise the majority of proteins, is less advanced. Specifically, large proteins often fold via long-lived partially folded intermediates, whose structures, potentially toxic oligomerization, and interactions with cellular chaperones remain poorly understood. Molecular dynamics based folding simulations that rely on knowledge of the native structure can provide critical, detailed information on folding free energy landscapes, intermediates, and pathways. Further, increases in computational power and methodological advances have made folding simulations of large proteins practical and valuable. Here, using serpins that inhibit proteases as an example, we review native-centric methods for simulating the folding of large proteins. These synergistic approaches range from Gō and related structure-based models that can predict the effects of the native structure on folding to all-atom-based methods that include side-chain chemistry and can predict how disease-associated mutations may impact folding. The application of these computational approaches to serpins and other large proteins highlights the successes and limitations of current computational methods and underscores how computational results can be used to inform experiments. These powerful simulation approaches in combination with experiments can provide unique insights into how large proteins fold and misfold, expanding our ability to predict and manipulate protein folding.

Efficient conversion of chemical energy into mechanical work by Hsp70 chaperones

Hsp70 molecular chaperones are abundant ATP-dependent nanomachines that actively reshape non-native, misfolded proteins and assist a wide variety of essential cellular processes. Here, we combine complementary theoretical approaches to elucidate the structural and thermodynamic details of the chaperone-induced expansion of a substrate protein, with a particular emphasis on the critical role played by ATP hydrolysis. We first determine the conformational free-energy cost of the substrate expansion due to the binding of multiple chaperones using coarse-grained molecular simulations. We then exploit this result to implement a non-equilibrium rate model which estimates the degree of expansion as a function of the free energy provided by ATP hydrolysis. Our results are in quantitative agreement with recent single-molecule FRET experiments and highlight the stark non-equilibrium nature of the process, showing that Hsp70s are optimized to effectively convert chemical energy into mechanical work close to physiological conditions.

Heat shock-induced chaperoning by Hsp70 is enabled in-cell

Recent work has shown that weak protein-protein interactions are susceptible to the cellular milieu. One case in point is the binding of heat shock proteins (Hsps) to substrate proteins in cells under stress. Upregulation of the Hsp70 chaperone machinery at elevated temperature was discovered in the 1960s, and more recent studies have shown that ATPase activity in one Hsp70 domain is essential for control of substrate binding by the other Hsp70 domain. Although there are several denaturant-based assays of Hsp70 activity, reports of ATP-dependent binding of Hsp70 to a globular protein substrate under heat shock are scarce. Here we show that binding of heat-inducible Hsp70 to phosphoglycerate kinase (PGK) is remarkably different in vitro compared to in-cell. We use fluorescent-labeled mHsp70 and ePGK, and begin by showing that mHsp70 passes the standard β-galactosidase assay, and that it does not self-aggregate until 50°C in presence of ATP. Yet during denaturant refolding or during in vitro heat shock, mHsp70 shows only ATP-independent non-specific sticking to ePGK, as evidenced by nearly identical results with an ATPase activity-deficient K71M mutant of Hsp70 as a control. Addition of Hsp40 (co-factor) or Ficoll (crowder) does not reduce non-specific sticking, but cell lysate does. Therefore, Hsp70 does not act as an ATP-dependent chaperone on its substrate PGK in vitro. In contrast, we observe only specific ATP-dependent binding of mHsp70 to ePGK in mammalian cells, when compared to the inactive Hsp70 K71M mutant. We hypothesize that enhanced in-cell activity is not due to an unknown co-factor, but simply to a favorable shift in binding equilibrium caused by the combination of crowding and osmolyte/macromolecular interactions present in the cell. One candidate mechanism for such a favorable shift in binding equilibrium is the proven ability of Hsp70 to bind near-native states of substrate proteins in vitro. We show evidence for early onset of binding in-cell. Our results suggest that Hsp70 binds PGK preemptively, prior to its full unfolding transition, thus stabilizing it against further unfolding. We propose a "preemptive holdase" mechanism for Hsp70-substrate binding. Given our result for PGK, more proteins than one might think based on in vitro assays may be chaperoned by Hsp70 in vivo. The cellular environment thus plays an important role in maintaining proper Hsp70 function.

mRNA Engineering for the Efficient Chaperone-Mediated Co-Translational Folding of Recombinant Proteins in Escherichia coli

The production of soluble, functional recombinant proteins by engineered bacterial hosts is challenging. Natural molecular chaperone systems have been used to solubilize various recombinant proteins with limited success. Here, we attempted to facilitate chaperone-mediated folding by directing the molecular chaperones to their protein substrates before the co-translational folding process completed. To achieve this, we either anchored the bacterial chaperone DnaJ to the 3ʹ untranslated region of a target mRNA by fusing with an RNA-binding domain in the chaperone-recruiting mRNA scaffold (CRAS) system, or coupled the expression of DnaJ and a target recombinant protein using the overlapping stop-start codons 5ʹ-TAATG-3ʹ between the two genes in a chaperone-substrate co-localized expression (CLEX) system. By engineering the untranslated and intergenic sequences of the mRNA transcript, bacterial molecular chaperones are spatially constrained to the location of protein translation, expressing selected aggregation-prone proteins in their functionally active, soluble form. Our mRNA engineering methods surpassed the in-vivo solubilization efficiency of the simple DnaJ chaperone co-overexpression method, thus providing more effective tools for producing soluble therapeutic proteins and enzymes.