BiPPred: Combined sequence- and structure-based prediction of peptide binding to the Hsp70 chaperone BiP

Molecular Modelling in Bioactive Peptide Discovery and Characterisation

Molecular modelling is a vital tool in the discovery and characterisation of bioactive peptides, providing insights into their structural properties and interactions with biological targets. Many models predicting bioactive peptide function or structure rely on their intrinsic properties, including the influence of amino acid composition, sequence, and chain length, which impact stability, folding, aggregation, and target interaction. Homology modelling predicts peptide structures based on known templates. Peptide-protein interactions can be explored using molecular docking techniques, but there are challenges related to the inherent flexibility of peptides, which can be addressed by more computationally intensive approaches that consider their movement over time, called molecular dynamics (MD). Virtual screening of many peptides, usually against a single target, enables rapid identification of potential bioactive peptides from large libraries, typically using docking approaches. The integration of artificial intelligence (AI) has transformed peptide discovery by leveraging large amounts of data. AlphaFold is a general protein structure prediction tool based on deep learning that has greatly improved the predictions of peptide conformations and interactions, in addition to providing estimates of model accuracy at each residue which greatly guide interpretation. Peptide function and structure prediction are being further enhanced using Protein Language Models (PLMs), which are large deep-learning-derived statistical models that learn computer representations useful to identify fundamental patterns of proteins. Recent methodological developments are discussed in the context of canonical peptides, as well as those with modifications and cyclisations. In designing potential peptide therapeutics, the main outstanding challenge for these methods is the incorporation of diverse non-canonical amino acids and cyclisations.

The Hsc70 system maintains the synaptic SNARE protein SNAP-25 in an assembly-competent state and delays its aggregation

The complex mechanism of synaptic vesicle fusion with the plasma membrane for neurotransmitter release is initiated by the formation of the SNARE complex at the presynaptic terminal of the neuron. The SNARE complex is composed of four helices contributed by three proteins: one from syntaxin (localized at the plasma membrane), one from synaptobrevin (localized at the synaptic vesicle), and two from the intrinsically disordered and aggregation-prone synaptosomal-associated 25 kDa protein (SNAP-25), which is localized to the plasma membrane by virtue of palmitoylation of cysteine residues. The fusion process is tightly regulated and requires the constitutively expressed Hsp70 chaperone (Hsc70) and its J-protein co-chaperone CSPα. We hypothesize that Hsc70 and CSPα cooperate to chaperone SNAP-25, disfavoring its aggregation and keeping it in a folding state competent for SNARE complex formation. To test this hypothesis, we used a bottom-up approach and studied the interaction between Hsc70 and CSPα with SNAP-25 in vitro. We showed that the aggregation of SNAP-25 is delayed in the presence of Hsc70 and CSPα. Using a peptide array that spans the sequence of SNAP-25, we identified three potential Hsc70-interacting sequences and designed peptides containing these sequences to test binding in solution. We characterized the interaction of SNAP-25-derived peptides with Hsc70 and CSPα using a combination of biochemical and biophysical techniques, including native-PAGE, binding affinity by fluorescence anisotropy, ATPase-activity of Hsc70, and NMR. We have identified an Hsc70 binding site within SNAP-25 that is likely to represent the site used in the cell to facilitate SNARE complex formation.

Probing Electrostatic and Hydrophobic Associative Interactions in Cells

Weak nonspecific interactions between biomacromolecules determine the cytoplasmic organization. Despite their importance, it is challenging to determine these interactions in the intracellular dense and heterogeneous mixture of biomacromolecules. Here, we develop a method to indicate electrostatic and hydrophobic associative interactions and map these interactions. The method relies on a genetically encoded probe containing a sensing peptide and a circularly permuted green fluorescent protein that provides a ratiometric readout. Inside bacterial and mammalian cells, we see that the cytoplasmic components interact strongly with cationic and hydrophobic probes but not with neutral hydrophilic probes, which remain inert. The Escherichia coli cytoplasm interacts strongly with highly negatively charged hydrophilic probes, but the HEK293T cytoplasm does not. These associative interactions are modulated by ATP depletion. Hence, the nonspecific associative interaction profile in cells is condition- and species-dependent.

Chaperone BiP controls ER stress sensor Ire1 through interactions with its oligomers

The complex multistep activation cascade of Ire1 involves changes in the Ire1 conformation and oligomeric state. Ire1 activation enhances ER folding capacity, in part by overexpressing the ER Hsp70 molecular chaperone BiP; in turn, BiP provides tight negative control of Ire1 activation. This study demonstrates that BiP regulates Ire1 activation through a direct interaction with Ire1 oligomers. Particularly, we demonstrated that the binding of Ire1 luminal domain (LD) to unfolded protein substrates not only trigger conformational changes in Ire1-LD that favour the formation of Ire1-LD oligomers but also exposes BiP binding motifs, enabling the molecular chaperone BiP to directly bind to Ire1-LD in an ATP-dependent manner. These transient interactions between BiP and two short motifs in the disordered region of Ire1-LD are reminiscent of interactions between clathrin and another Hsp70, cytoplasmic Hsc70. BiP binding to substrate-bound Ire1-LD oligomers enables unfolded protein substrates and BiP to synergistically and dynamically control Ire1-LD oligomerisation, helping to return Ire1 to its deactivated state when an ER stress response is no longer required.

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.

Sequence elements within the PEXEL motif and its downstream region modulate PTEX-dependent protein export in Plasmodium falciparum

The parasite Plasmodium falciparum causes the most severe form of malaria and to invade and replicate in red blood cells (RBCs), it exports hundreds of proteins across the encasing parasitophorous vacuole membrane (PVM) into this host cell. The exported proteins help modify the RBC to support rapid parasite growth and avoidance of the human immune system. Most exported proteins possess a conserved Plasmodium export element (PEXEL) motif with the consensus RxLxE/D/Q amino acid sequence, which acts as a proteolytic cleavage recognition site within the parasite's endoplasmic reticulum (ER). Cleavage occurs after the P1 L residue and is thought to help release the protein from the ER so it can be putatively escorted by the HSP101 chaperone to the parasitophorous vacuole space surrounding the intraerythrocytic parasite. HSP101 and its cargo are then thought to assemble with the rest of a Plasmodium translocon for exported proteins (PTEX) complex, that then recognises the xE/D/Q capped N-terminus of the exported protein and translocates it across the vacuole membrane into the RBC compartment. Here, we present evidence that supports a dual role for the PEXEL's conserved P2 ' position E/Q/D residue, first, for plasmepsin V cleavage in the ER, and second, for efficient PTEX mediated export across the PVM into the RBC. We also present evidence that the downstream 'spacer' region separating the PEXEL motif from the folded functional region of the exported protein controls cargo interaction with PTEX as well. The spacer must be of a sufficient length and permissive amino acid composition to engage the HSP101 unfoldase component of PTEX to be efficiently translocated into the RBC compartment.

Protein Biosynthesis and Maturation in the ER

The endoplasmic reticulum takes care of the folding, assembly, and quality control of thousands of proteins destined to the different compartments of the endomembrane system or to be secreted in the apoplast. Here we describe how these early events in the life of all these proteins can be followed biochemically by using velocity or isopycnic ultracentrifugation, metabolic labelling with radioactive amino acids, drug treatments, and immunoselection in various conditions and, in certain cases, predicted in silico by algorithms.

Effects of HSP70 chaperones Ssa1 and Ssa2 on Ste5 scaffold and the mating mitogen-activated protein kinase (MAPK) pathway in Saccharomyces cerevisiae

Ste5 is a prototype of scaffold proteins that regulate activation of mitogen-activated protein kinase (MAPK) cascades in all eukaryotes. Ste5 associates with many proteins including Gβγ (Ste4), Ste11 MAPKKK, Ste7 MAPKK, Fus3 and Kss1 MAPKs, Bem1, Cdc24. Here we show that Ste5 also associates with heat shock protein 70 chaperone (Hsp70) Ssa1 and that Ssa1 and its ortholog Ssa2 are together important for Ste5 function and efficient mating responses. The majority of purified overexpressed Ste5 associates with Ssa1. Loss of Ssa1 and Ssa2 has deleterious effects on Ste5 abundance, integrity, and localization particularly when Ste5 is expressed at native levels. The status of Ssa1 and Ssa2 influences Ste5 electrophoresis mobility and formation of high molecular weight species thought to be phosphorylated, ubiquitinylated and aggregated and lower molecular weight fragments. A Ste5 VWA domain mutant with greater propensity to form punctate foci has reduced predicted propensity to bind Ssa1 near the mutation sites and forms more punctate foci when Ssa1 Is overexpressed, supporting a dynamic protein quality control relationship between Ste5 and Ssa1. Loss of Ssa1 and Ssa2 reduces activation of Fus3 and Kss1 MAPKs and FUS1 gene expression and impairs mating shmoo morphogenesis. Surprisingly, ssa1, ssa2, ssa3 and ssa4 single, double and triple mutants can still mate, suggesting compensatory mechanisms exist for folding. Additional analysis suggests Ssa1 is the major Hsp70 chaperone for the mating and invasive growth pathways and reveals several Hsp70-Hsp90 chaperone-network proteins required for mating morphogenesis.

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.

Folding of heterologous proteins in bacterial cell factories: Cellular mechanisms and engineering strategies

The expression of correctly folded and functional heterologous proteins is important in many biotechnological production processes, whether it is enzymes, biopharmaceuticals or biosynthetic pathways for production of sustainable chemicals. For industrial applications, bacterial platform organisms, such as E. coli, are still broadly used due to the availability of tools and proven suitability at industrial scale. However, expression of heterologous proteins in these organisms can result in protein aggregation and low amounts of functional protein. This review provides an overview of the cellular mechanisms that can influence protein folding and expression, such as co-translational folding and assembly, chaperone binding, as well as protein quality control, across different model organisms. The knowledge of these mechanisms is then linked to different experimental methods that have been applied in order to improve functional heterologous protein folding, such as codon optimization, fusion tagging, chaperone co-production, as well as strain and protein engineering strategies.

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

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

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.

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.

Electrostatics drive the molecular chaperone BiP to preferentially bind oligomerized states of a client protein

Hsp70 chaperones bind short monomeric peptides with a weak characteristic affinity in the low micromolar range, but can also bind some aggregates, fibrils, and amyloids, with low nanomolar affinity. While this differential affinity enables Hsp70 to preferentially target potentially toxic aggregates, it is unknown how a chaperone can differentiate between monomeric and aggregated states of a client protein and why preferential binding is only observed for some aggregated clients but not others. Here we examine the interaction of BiP (the Hsp70 paralog in the endoplasmic reticulum) with the client proIGF2, the pro-protein form of IGF2 that includes a long and mostly disordered E-peptide region that promotes proIGF2 oligomerization. By dissecting the mechanism by which BiP targets proIGF2 and E-peptide oligomers we discover that electrostatic attraction is a powerful driving force for oligomer recognition. We identify the specific BiP binding sites on proIGF2 and as monomers they bind BiP with characteristically weak affinity in the low micromolar range, but electrostatic attraction to E-peptide oligomers boosts the affinity to the low nanomolar level. The dominant role of electrostatics is manifested kinetically as a steering force that accelerates the binding of BiP to E-peptide oligomers by approximately two orders of magnitude as compared against monomeric peptides. Electrostatic targeting of Hsp70 provides an explanation for why preferential binding has been observed for some aggregated clients but not others, as all the currently-documented cases in which Hsp70 binds aggregates with high-affinity involve clients that have an opposite charge to Hsp70.

Z-α1-antitrypsin polymers impose molecular filtration in the endoplasmic reticulum after undergoing phase transition to a solid state

Misfolding of secretory proteins in the endoplasmic reticulum (ER) features in many human diseases. In α1-antitrypsin deficiency, the pathogenic Z variant aberrantly assembles into polymers in the hepatocyte ER, leading to cirrhosis. We show that α1-antitrypsin polymers undergo a liquid:solid phase transition, forming a protein matrix that retards mobility of ER proteins by size-dependent molecular filtration. The Z-α1-antitrypsin phase transition is promoted during ER stress by an ATF6-mediated unfolded protein response. Furthermore, the ER chaperone calreticulin promotes Z-α1-antitrypsin solidification and increases protein matrix stiffness. Single-particle tracking reveals that solidification initiates in cells with normal ER morphology, previously assumed to represent a healthy pool. We show that Z-α1-antitrypsin-induced hypersensitivity to ER stress can be explained by immobilization of ER chaperones within the polymer matrix. This previously unidentified mechanism of ER dysfunction provides a template for understanding a diverse group of related proteinopathies and identifies ER chaperones as potential therapeutic targets.

Structure of Hsp90-Hsp70-Hop-GR reveals the Hsp90 client-loading mechanism

Maintaining a healthy proteome is fundamental for the survival of all organisms¹. Integral to this are Hsp90 and Hsp70, molecular chaperones that together facilitate the folding, remodelling and maturation of the many 'client proteins' of Hsp90². The glucocorticoid receptor (GR) is a model client protein that is strictly dependent on Hsp90 and Hsp70 for activity^3-7. Chaperoning GR involves a cycle of inactivation by Hsp70; formation of an inactive GR-Hsp90-Hsp70-Hop 'loading' complex; conversion to an active GR-Hsp90-p23 'maturation' complex; and subsequent GR release⁸. However, to our knowledge, a molecular understanding of this intricate chaperone cycle is lacking for any client protein. Here we report the cryo-electron microscopy structure of the GR-loading complex, in which Hsp70 loads GR onto Hsp90, uncovering the molecular basis of direct coordination by Hsp90 and Hsp70. The structure reveals two Hsp70 proteins, one of which delivers GR and the other scaffolds the Hop cochaperone. Hop interacts with all components of the complex, including GR, and poises Hsp90 for subsequent ATP hydrolysis. GR is partially unfolded and recognized through an extended binding pocket composed of Hsp90, Hsp70 and Hop, revealing the mechanism of GR loading and inactivation. Together with the GR-maturation complex structure⁹, we present a complete molecular mechanism of chaperone-dependent client remodelling, and establish general principles of client recognition, inhibition, transfer and activation.

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.

Two γ-zeins induce the unfolded protein response

The rapid, massive synthesis of storage proteins that occurs during seed development stresses endoplasmic reticulum (ER) homeostasis, which activates the ER unfolded protein response (UPR). However, how different storage proteins contribute to UPR is not clear. We analyzed vegetative tissues of transgenic Arabidopsis (Arabidopsis thaliana) plants constitutively expressing the common bean (Phaseolus vulgaris) soluble vacuolar storage protein PHASEOLIN (PHSL) or maize (Zea mays) prolamins (27-kDa γ-zein or 16-kDa γ-zein) that participate in forming insoluble protein bodies in the ER. We show that 16-kDa γ-zein significantly activates the INOSITOL REQUIRING ENZYME1/BASIC LEUCINE ZIPPER 60 (bZIP60) UPR branch-but not the bZIP28 branch or autophagy-leading to induction of major UPR-controlled genes that encode folding helpers that function inside the ER. Protein blot analysis of IMMUNOGLOBULIN-BINDING PROTEIN (BIP) 1 and 2, BIP3, GLUCOSE REGULATED PROTEIN 94 (GRP94), and ER-localized DNAJ family 3A (ERDJ3A) polypeptides confirmed their higher accumulation in the plant expressing 16-kDa γ-zein. Expression of 27-kDa γ-zein significantly induced only BIP3 and ERDJ3A transcription even though an increase in GRP94 and BIP1/2 polypeptides also occurred in this plant. These results indicate a significant but weaker effect of 27-kDa γ-zein compared to 16-kDa γ-zein, which corresponds with the higher availability of 16-kDa γ-zein for BIP binding, and indicates subtle protein-specific modulations of plant UPR. None of the analyzed genes was significantly induced by PHSL or by a mutated, soluble form of 27-kDa γ-zein that traffics along the secretory pathway. Such variability in UPR induction may have influenced the evolution of storage proteins with different tissue and subcellular localization.

Influenza A viruses balance ER stress with host protein synthesis shutoff

Excessive production of viral glycoproteins during infections poses a tremendous stress potential on the endoplasmic reticulum (ER) protein folding machinery of the host cell. The host cell balances this by providing more ER resident chaperones and reducing translation. For viruses, this unfolded protein response (UPR) offers the potential to fold more glycoproteins. We postulated that viruses could have developed means to limit the inevitable ER stress to a beneficial level for viral replication. Using a relevant human pathogen, influenza A virus (IAV), we first established the determinant for ER stress and UPR induction during infection. In contrast to a panel of previous reports, we identified neuraminidase to be the determinant for ER stress induction, and not hemagglutinin. IAV relieves ER stress by expression of its nonstructural protein 1 (NS1). NS1 interferes with the host messenger RNA processing factor CPSF30 and suppresses ER stress response factors, such as XBP1. In vivo viral replication is increased when NS1 antagonizes ER stress induction. Our results reveal how IAV optimizes glycoprotein expression by balancing folding capacity.

Molecular Dynamics Scoring of Protein-Peptide Models Derived from Coarse-Grained Docking

One of the major challenges in the computational prediction of protein-peptide complexes is the scoring of predicted models. Usually, it is very difficult to find the most accurate solutions out of the vast number of sometimes very different and potentially plausible predictions. In this work, we tested the protocol for Molecular Dynamics (MD)-based scoring of protein-peptide complex models obtained from coarse-grained (CG) docking simulations. In the first step of the scoring procedure, all models generated by CABS-dock were reconstructed starting from their original C-alpha trace representations to all-atom (AA) structures. The second step included geometry optimization of the reconstructed complexes followed by model scoring based on receptor-ligand interaction energy estimated from short MD simulations in explicit water. We used two well-known AA MD force fields, CHARMM and AMBER, and a CG MARTINI force field. Scoring results for 66 different protein-peptide complexes show that the proposed MD-based scoring approach can be used to identify protein-peptide models of high accuracy. The results also indicate that the scoring accuracy may be significantly affected by the quality of the reconstructed protein receptor structures.

The ER Chaperones BiP and Grp94 Regulate the Formation of Insulin-Like Growth Factor 2 (IGF2) Oligomers

While cytosolic Hsp90 chaperones have been extensively studied, less is known about how the ER Hsp90 paralog Grp94 recognizes clients and influences client folding. Here, we examine how Grp94 and the ER Hsp70 paralog, BiP, influence the folding of insulin-like growth factor 2 (IGF2), an established client protein of Grp94. ProIGF2 is composed of a disulfide-bonded insulin-like hormone and a C-terminal E-peptide that has sequence characteristics of an intrinsically disordered region. BiP and Grp94 have a minimal influence on folding whereby both chaperones slow proIGF2 folding and do not substantially alter the disulfide-bonded folding intermediates, suggesting that BiP and Grp94 may have an additional influence unrelated to proIGF2 folding. Indeed, we made the unexpected discovery that the E-peptide region allows proIGF2 to form dynamic oligomers. ProIGF2 oligomers can transition from a dynamic state that is capable of exchanging monomers to an irreversibly aggregated state, providing a plausible role for BiP and Grp94 in regulating proIGF2 oligomerization. In contrast to the modest influence on folding, BiP and Grp94 have a stronger influence on proIGF2 oligomerization and these chaperones exert counteracting effects. BiP suppresses proIGF2 oligomerization while Grp94 can enhance proIGF2 oligomerization in a nucleotide-dependent manner. We propose that BiP and Grp94 regulate the assembly and dynamic behavior of proIGF2 oligomers, although the biological role of proIGF2 oligomerization is not yet known.

Structure based peptide design, molecular dynamics and MM-PBSA studies for targeting C terminal dimerization of NFAT5 DNA binding domain

NFAT5 as a transcription factor with an established role in osmotic stress response, has also been revealed to be active under numerous settings, including pathological conditions such as diabetic microvascular complications, chronic arthritis and cancer. Despite these links, current strategies for downregulating NFAT5 activity only relies on indirect modulators, not directly targeting NFAT5, itself. With this study, through using a computational approach, an original peptide was explored to directly target C terminal dimerization of NFAT5 RHR, located in its DNA binding domain. At first, homodimeric NFAT5 RHR bound to its consensus DNA was used for prediction of a preliminary peptide sequence. Possible amino acid replacements for this preliminary peptide were predicted for optimization, which was followed by addition of a cell penetrating peptide sequence. These attempts yielded a small peptide library, which was further investigated for peptide affinities towards C terminal of NFAT5 RHR through molecular docking, 50 ns and 250 ns molecular dynamics simulations, followed by estimation of MM-PBSA based relative binding free energies. Results indicated that after receiving mutations on the preliminary peptide sequence for optimization, a unique peptide could target C terminal dimerization region of NFAT5 RHR through using its cell penetrating peptide sequence. In conclusion, this is the first study presenting computational evidence on identification of a novel peptide capable of directly targeting NFAT5 dimerization. Besides, future implications of these observations were also discussed in terms of methodology and possible applications.

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.

Combining Three-Dimensional Modeling with Artificial Intelligence to Increase Specificity and Precision in Peptide-MHC Binding Predictions

The reliable prediction of the affinity of candidate peptides for the MHC is important for predicting their potential antigenicity and thus influences medical applications, such as decisions on their inclusion in T cell-based vaccines. In this study, we present a rapid, predictive computational approach that combines a popular, sequence-based artificial neural network method, NetMHCpan 4.0, with three-dimensional structural modeling. We find that the ensembles of bound peptide conformations generated by the programs MODELLER and Rosetta FlexPepDock are less variable in geometry for strong binders than for low-affinity peptides. In tests on 1271 peptide sequences for which the experimental dissociation constants of binding to the well-characterized murine MHC allele H-2Db are known, by applying thresholds for geometric fluctuations the structure-based approach in a standalone manner drastically improves the statistical specificity, reducing the number of false positives. Furthermore, filtering candidates generated with NetMHCpan 4.0 with the structure-based predictor led to an increase in the positive predictive value (PPV) of the peptides correctly predicted to bind very strongly (i.e., K d < 100 nM) from 40 to 52% (p = 0.027). The combined method also significantly improved the PPV when tested on five human alleles, including some with limited data for training. Overall, an average increase of 10% in the PPV was found over the standalone sequence-based method. The combined method should be useful in the rapid design of effective T cell-based vaccines.

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.

The extent of Ssa1/Ssa2 Hsp70 chaperone involvement in nuclear protein quality control degradation varies with the substrate

Protein misfolding is a recurring phenomenon that cells must manage; otherwise misfolded proteins can aggregate and become toxic should they persist. To counter this burden, cells have evolved protein quality control (PQC) mechanisms that manage misfolded proteins. Two classes of systems that function in PQC are chaperones that aid in protein folding and ubiquitin-protein ligases that ubiquitinate misfolded proteins for proteasomal degradation. How folding and degradative PQC systems interact and coordinate their respective functions is not yet fully understood. Previous studies of PQC degradation pathways in the endoplasmic reticulum and cytosol have led to the prevailing idea that these pathways require the activity of Hsp70 chaperones. Here, we find that involvement of the budding yeast Hsp70 chaperones Ssa1 and Ssa2 in nuclear PQC degradation varies with the substrate. In particular, nuclear PQC degradation mediated by the yeast ubiquitin-protein ligase San1 often involves Ssa1/Ssa2, but San1 substrate recognition and ubiquitination can proceed without these Hsp70 chaperone functions in vivo and in vitro. Our studies provide new insights into the variability of Hsp70 chaperone involvement with a nuclear PQC degradation pathway.

Actin stabilizing compounds show specific biological effects due to their binding mode

Actin binding compounds are widely used tools in cell biology. We compare the biological and biochemical effects of miuraenamide A and jasplakinolide, a structurally related prototypic actin stabilizer. Though both compounds have similar effects on cytoskeletal morphology and proliferation, they affect migration and transcription in a distinctive manner, as shown by a transcriptome approach in endothelial cells. In vitro, miuraenamide A acts as an actin nucleating, F-actin polymerizing and stabilizing compound, just like described for jasplakinolide. However, in contrast to jasplakinolide, miuraenamide A competes with cofilin, but not gelsolin or Arp2/3 for binding to F-actin. We propose a binding mode of miuraenamide A, explaining both its similarities and its differences to jasplakinolide. Molecular dynamics simulations suggest that the bromophenol group of miurenamide A interacts with residues Tyr133, Tyr143, and Phe352 of actin. This shifts the D-loop of the neighboring actin, creating tighter packing of the monomers, and occluding the binding site of cofilin. Since relatively small changes in the molecular structure give rise to this selectivity, actin binding compounds surprisingly are promising scaffolds for creating actin binders with specific functionality instead of just "stabilizers".

Predicting the bioactive conformations of macrocycles: a molecular dynamics-based docking procedure with DynaDock

Macrocyclic compounds are of growing interest as a new class of therapeutics, especially as inhibitors binding to protein-protein interfaces. As molecular modeling is a well-established complimentary tool in modern drug design, the number of attempts to develop reliable docking strategies and algorithms to accurately predict the binding mode of macrocycles is rising continuously. Standard molecular docking approaches need to be adapted to this application, as a comprehensive yet efficient sampling of all ring conformations of the macrocycle is necessary. To overcome this issue, we designed a molecular dynamics-based docking protocol for macrocycles, in which the challenging sampling step is addressed by conventional molecular dynamics (750 ns) simulations performed at moderately high temperature (370 K). Consecutive flexible docking with the DynaDock approach based on multiple, pre-sampled ring conformations yields highly accurate poses with ligand RMSD values lower than 1.8 Å. We further investigated the value of molecular dynamics-based complex stability estimations for pose selection and discuss its applicability in combination with standard binding free energy estimations for assessing the quality of poses in future blind docking studies.

Protein Quality Control in the Endoplasmic Reticulum

The site of protein folding and maturation for the majority of proteins that are secreted, localized to the plasma membrane or targeted to endomembrane compartments is the endoplasmic reticulum (ER). It is essential that proteins targeted to the ER are properly folded in order to carry out their function, as well as maintain protein homeostasis, as accumulation of misfolded proteins could lead to the formation of cytotoxic aggregates. Because protein folding is an error-prone process, the ER contains protein quality control networks that act to optimize proper folding and trafficking of client proteins. If a protein is unable to reach its native state, it is targeted for ER retention and subsequent degradation. The protein quality control networks of the ER that oversee this evaluation or interrogation process that decides the fate of maturing nascent chains is comprised of three general types of families: the classical chaperones, the carbohydrate-dependent system, and the thiol-dependent system. The cooperative action of these families promotes protein quality control and protein homeostasis in the ER. This review will describe the families of the ER protein quality control network and discuss the functions of individual members.

Functional principles and regulation of molecular chaperones

To be able to perform their biological function, a protein needs to be correctly folded into its three dimensional structure. The protein folding process is spontaneous and does not require the input of energy. However, in the crowded cellular environment where there is high risk of inter-molecular interactions that may lead to protein molecules sticking to each other, hence forming aggregates, protein folding is assisted. Cells have evolved robust machinery called molecular chaperones to deal with the protein folding problem and to maintain proteins in their functional state. Molecular chaperones promote efficient folding of newly synthesized proteins, prevent their aggregation and ensure protein homeostasis in cells. There are different classes of molecular chaperones functioning in a complex interplay. In this review, we discuss the principal characteristics of different classes of molecular chaperones, their structure-function relationships, their mode of regulation and their involvement in human disorders.

Recent advances in the structural and mechanistic aspects of Hsp70 molecular chaperones

Hsp70 chaperones are central hubs of the protein quality control network and collaborate with co-chaperones having a J-domain (an ∼70-residue-long helical hairpin with a flexible loop and a conserved His-Pro-Asp motif required for ATP hydrolysis by Hsp70s) and also with nucleotide exchange factors to facilitate many protein-folding processes that (re)establish protein homeostasis. The Hsp70s are highly dynamic nanomachines that modulate the conformation of their substrate polypeptides by transiently binding to short, mostly hydrophobic stretches. This interaction is regulated by an intricate allosteric mechanism. The J-domain co-chaperones target Hsp70 to their polypeptide substrates, and the nucleotide exchange factors regulate the lifetime of the Hsp70-substrate complexes. Significant advances in recent years are beginning to unravel the molecular mechanism of this chaperone machine and how they treat their substrate proteins.

Ca2+ binding induced sequential allosteric activation of sortase A: An example for ion-triggered conformational selection

The allosteric activation of the intrinsically disordered enzyme Staphylococcus aureus sortase A is initiated via binding of a Ca²⁺ ion. Although Ca²⁺ binding was shown to initiate structural changes inducing disorder-to-order transitions, the details of the allosteric activation mechanism remain elusive. We performed long-term molecular dynamics simulations of sortase A without (3 simulations of 1.6 μs) and with bound Ca²⁺ (simulations of 1.6 μs, 1.8 μs, and 2.5 μs). Our results show that Ca²⁺ binding causes not only ordering of the disordered β6/β7 loop of the protein, but also modulates hinge motions in the dynamic β7/β8 loop, which is important for the catalytic activity of the enzyme. Cation binding triggers signal transmission from the Ca²⁺ binding site to the dynamic β7/β8 loop via the repetitive folding/unfolding of short helical stretches of the disordered β6/β7 loop. These correlated structural rearrangements lead to several distinct conformational states of the binding groove, which show optimal binding features for the sorting signal motif and feature binding energies up to 20 kcal/mol more favorable than observed for the sortase A without Ca²⁺. The presented results indicate a highly correlated, conformational selection-based activation mechanism of the enzyme triggered by cation binding. They also demonstrate the importance of the dynamics of intrinsically disordered regions for allosteric regulation.

Cytosolic Hsp70 and Hsp40 chaperones enable the biogenesis of mitochondrial β-barrel proteins

Mitochondrial β-barrel proteins are imported from the cytosol into the organelle. Jores et al. provide new insights into the early events of this process by describing an array of cytosolic chaperones and cochaperones that associate with newly synthesized β-barrel proteins and assure their optimal biogenesis.

Mitochondrial β-barrel proteins are encoded in the nucleus, translated by cytosolic ribosomes, and then imported into the organelle. Recently, a detailed understanding of the intramitochondrial import pathway of β-barrel proteins was obtained. In contrast, it is still completely unclear how newly synthesized β-barrel proteins reach the mitochondrial surface in an import-competent conformation. In this study, we show that cytosolic Hsp70 chaperones and their Hsp40 cochaperones Ydj1 and Sis1 interact with newly synthesized β-barrel proteins. These interactions are highly relevant for proper biogenesis, as inhibiting the activity of the cytosolic Hsp70, preventing its docking to the mitochondrial receptor Tom70, or depleting both Ydj1 and Sis1 resulted in a significant reduction in the import of such substrates into mitochondria. Further experiments demonstrate that the interactions between β-barrel proteins and Hsp70 chaperones and their importance are conserved also in mammalian cells. Collectively, this study outlines a novel mechanism in the early events of the biogenesis of mitochondrial outer membrane β-barrel proteins.

Protein-peptide docking: opportunities and challenges

Peptides have recently attracted much attention as promising drug candidates. Rational design of peptide-derived therapeutics usually requires structural characterization of the underlying protein-peptide interaction. Given that experimental characterization can be difficult, reliable computational tools are needed. In recent years, a variety of approaches have been developed for 'protein-peptide docking', that is, predicting the structure of the protein-peptide complex, starting from the protein structure and the peptide sequence, including variable degrees of information about the peptide binding site and/or conformation. In this review, we provide an overview of protein-peptide docking methods and outline their capabilities, limitations, and applications in structure-based drug design. Key challenges are also briefly discussed, such as modeling of large-scale conformational changes upon binding, scoring of predicted models, and optimal inclusion of varied types of experimental data and theoretical predictions into an integrative modeling process.

Recent Developments and Applications of the MMPBSA Method

The Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approach has been widely applied as an efficient and reliable free energy simulation method to model molecular recognition, such as for protein-ligand binding interactions. In this review, we focus on recent developments and applications of the MMPBSA method. The methodology review covers solvation terms, the entropy term, extensions to membrane proteins and high-speed screening, and new automation toolkits. Recent applications in various important biomedical and chemical fields are also reviewed. We conclude with a few future directions aimed at making MMPBSA a more robust and efficient method.

Bap (Sil1) regulates the molecular chaperone BiP by coupling release of nucleotide and substrate

BiP is the endoplasmic member of the Hsp70 family. BiP is regulated by several co-chaperones including the nucleotide-exchange factor (NEF) Bap (Sil1 in yeast). Bap is a two-domain protein. The interaction of the Bap C-terminal domain with the BiP ATPase domain is sufficient for its weak NEF activity. However, stimulation of the BiP ATPase activity requires full-length Bap, suggesting a complex interplay of these two factors. Here, single-molecule FRET experiments with mammalian proteins reveal that Bap affects the conformation of both BiP domains, including the lid subdomain, which is important for substrate binding. The largely unstructured Bap N-terminal domain promotes the substrate release from BiP. Thus, Bap is a conformational regulator affecting both nucleotide and substrate interactions. The preferential interaction with BiP in its ADP state places Bap at a late stage of the chaperone cycle, in which it coordinates release of substrate and ADP, thereby resetting BiP for ATP and substrate binding.

Protein Folding Mediated by Trigger Factor and Hsp70: New Insights from Single-Molecule Approaches

Chaperones assist in protein folding, but what this common phrase means in concrete terms has remained surprisingly poorly understood. We can readily measure chaperone binding to unfolded proteins, but how they bind and affect proteins along folding trajectories has remained obscure. Here we review recent efforts by our labs and others that are beginning to pry into this issue, with a focus on the chaperones trigger factor and Hsp70. Single-molecule methods are central, as they allow the stepwise process of folding to be followed directly. First results have already revealed contrasts with long-standing paradigms: rather than acting only "early" by stabilizing unfolded chain segments, these chaperones can bind and stabilize partially folded structures as they grow to their native state. The findings suggest a fundamental redefinition of the protein folding problem and a more extensive functional repertoire of chaperones than previously assumed.

Members of the Hsp70 Family Recognize Distinct Types of Sequences to Execute ER Quality Control

Protein maturation in the endoplasmic reticulum is controlled by multiple chaperones, but how they recognize and determine the fate of their clients remains unclear. We developed an in vivo peptide library covering substrates of the ER Hsp70 system: BiP, Grp170, and three of BiP's DnaJ-family co-factors (ERdj3, ERdj4, and ERdj5). In vivo binding studies revealed that sites for pro-folding chaperones BiP and ERdj3 were frequent and dispersed throughout the clients, whereas Grp170, ERdj4, and ERdj5 specifically recognized a distinct type of rarer sequence with a high predicted aggregation potential. Mutational analyses provided insights into sequence recognition characteristics for these pro-degradation chaperones, which could be readily introduced or disrupted, allowing the consequences for client fates to be determined. Our data reveal unanticipated diversity in recognition sequences for chaperones; establish a sequence-encoded interplay between protein folding, aggregation, and degradation; and highlight the ability of clients to co-evolve with chaperones, ensuring quality control.