Mechanical properties of BiP protein determined by nano‐rheology

Immunoglobulin Binding Protein (BiP) is a chaperone and molecular motor belonging to the Hsp70 family, involved in the regulation of important biological processes such as synthesis, folding and translocation of proteins in the Endoplasmic Reticulum. BiP has two highly conserved domains: the N‐terminal Nucleotide‐Binding Domain (NBD), and the C‐terminal Substrate‐Binding Domain (SBD), connected by a hydrophobic linker. ATP binds and it is hydrolyzed to ADP in the NBD, and BiP's extended polypeptide substrates bind in the SBD. Like many molecular motors, BiP function depends on both structural and catalytic properties that may contribute to its performance. One novel approach to study the mechanical properties of BiP considers exploring the changes in the viscoelastic behavior upon ligand binding, using a technique called nano‐rheology. This technique is essentially a traditional rheology experiment, in which an oscillatory force is directly applied to the protein under study, and the resulting average deformation is measured. Our results show that the folded state of the protein behaves like a viscoelastic material, getting softer when it binds nucleotides‐ ATP, ADP, and AMP‐PNP‐, but stiffer when binding HTFPAVL peptide substrate. Also, we observed that peptide binding dramatically increases the affinity for ADP, decreasing it dissociation constant (K_D) around 1000 times, demonstrating allosteric coupling between SBD and NBD domains.     

Interaction of the chaperone BiP with an antibody domain: implications for the chaperone cycle

BiP is an Hsp70 homologue found in the endoplasmic reticulum of eukaryotic cells. Like other Hsp70 chaperones, BiP interacts with its substrate proteins in an ATP-dependent manner. The functional analysis has so far been performed mainly with short, synthetic peptides. Here, we present an experimental system that allows to study the partial reactions of the BiP chaperone cycle for a natural substrate protein domain in its soluble, stably unfolded conformation. This unfolded antibody domain forms a binary complex with BiP in the absence of ATP. The dissociation of the BiP dimer seems to be the rate-limiting step in this reaction. The BiP-C(H)3 complexes dissociate rapidly in the presence of ATP. The affinity for BiP-binding peptides and the non-native antibody domain was determined to be similar, suggesting that only the peptide binding site is involved in these interactions. Furthermore, these results imply that, also in the context of the antibody domain, an extended peptide sequence is recognized. However, the accessibility of the BiP-binding site in the non-native protein seems to influence the kinetics of complex formation.     

Epitope‐specificity of recombinant antibodies reveals promiscuous peptide‐binding properties

Protein–peptide interactions are a common occurrence and essential for numerous cellular processes, and frequently explored in broad applications within biology, medicine, and proteomics. Therefore, understanding the molecular mechanism(s) of protein–peptide recognition, specificity, and binding interactions will be essential. In this study, we report the first detailed analysis of antibody–peptide interaction characteristics, by combining large‐scale experimental peptide binding data with the structural analysis of eight human recombinant antibodies and numerous peptides, targeting tryptic mammalian and eukaryote proteomes. The results consistently revealed that promiscuous peptide‐binding interactions, that is, both specific and degenerate binding, were exhibited by all antibodies, and the discovery was corroborated by orthogonal data, indicating that this might be a general phenomenon for low‐affinity antibody–peptide interactions. The molecular mechanism for the degenerate peptide‐binding specificity appeared to be executed through the use of 2–3 semi‐conserved anchor residues in the C‐terminal part of the peptides, in analogue to the mechanism utilized by the major histocompatibility complex–peptide complexes. In the long‐term, this knowledge will be instrumental for advancing our fundamental understanding of protein–peptide interactions, as well as for designing, generating, and applying peptide specific antibodies, or peptide‐binding proteins in general, in various biotechnical and medical applications.     

Crystallographic structure of the tetratricopeptide repeat domain of Plasmodium falciparum FKBP35 and its molecular interaction with Hsp90 C‐terminal pentapeptide

Plasmodium falciparum FK506‐binding protein 35 (PfFKBP35) that binds to FK506 contains a conserved tetratricopeptide repeat (TPR) domain. Several known TPR domains such as Hop, PPP5, CHIP, and FKBP52 are structurally conserved and are able to interact with molecular chaperones such as Hsp70/Hsp90. Here, we present the crystal structure of PfFKBP35‐TPR and demonstrate its interaction with Hsp90 C‐terminal pentapeptide (MEEVD) by surface plasmon resonance and nuclear magnetic resonance spectroscopy‐based binding studies. Our sequence and structural analyses reveal that PfFKBP35 is similar to Hop and PPP5 in possessing all the conserved residues which are important for carboxylate clamping with Hsp90. Mutational studies were carried out on positively charged clamp residues that are crucial for binding to carboxylate groups of aspartate, showing that all the mutated residues are important for Hsp90 binding. Molecular docking and electrostatic calculations demonstrated that the MEEVD peptide of Hsp90 can form aspartate clamp unlike FKBP52. Our results provide insightful information and structural basis about the molecular interaction between PfFKBP35‐TPR and Hsp90.     

Ab initio prediction of peptide-MHC binding geometry for diverse class I MHC allotypes

Since determining the crystallographic structure of all peptide-MHC complexes is infeasible, an accurate prediction of the conformation is a critical computational problem. These models can be useful for determining binding energetics, predicting the structures of specific ternary complexes with T-cell receptors, and designing new molecules interacting with these complexes. The main difficulties are (1) adequate sampling of the large number of conformational degrees of freedom for the flexible peptide, (2) predicting subtle changes in the MHC interface geometry upon binding, and (3) building models for numerous MHC allotypes without known structures. Whereas previous studies have approached the sampling problem by dividing the conformational variables into different sets and predicting them separately, we have refined the Biased-Probability Monte Carlo docking protocol in internal coordinates to optimize a physical energy function for all peptide variables simultaneously. We also imitated the induced fit by docking into a more permissive smooth grid representation of the MHC followed by refinement and reranking using an all-atom MHC model. Our method was tested by a comparison of the results of cross-docking 14 peptides into HLA-A*0201 and 9 peptides into H-2K(b) as well as docking peptides into homology models for five different HLA allotypes with a comprehensive set of experimental structures. The surprisingly accurate prediction (0.75 A backbone RMSD) for cross-docking of a highly flexible decapeptide, dissimilar to the original bound peptide, as well as docking predictions using homology models for two allotypes with low average backbone RMSDs of less than 1.0 A illustrate the method's effectiveness. Finally, energy terms calculated using the predicted structures were combined with supervised learning on a large data set to classify peptides as either HLA-A*0201 binders or nonbinders. In contrast with sequence-based prediction methods, this model was also able to predict the binding affinity for peptides to a different MHC allotype (H-2K(b)), not used for training, with comparable prediction accuracy.     

Specificity of peptide-induced depolymerization of the recombinant carboxy-terminal fragment of BiP/GRP78.

In the present study, we have used a non-denaturing gel electrophoresis assay to characterize the specificity of the peptide-induced depolymerization process of the isolated recombinant C-terminal domain (C30) of the molecular chaperone BiP, in the presence of specific synthetic peptides and with the neuropeptide Substance P. In the absence of peptidic ligand, C30 self-associates readily into multiple oligomeric species. Upon peptide addition, C30 oligomers convert into dimers, then into monomers. Our data indicate that the algorithm we previously developed to predict putative BiP binding sites in any protein sequence is also a good indicator as to whether a peptide can efficiently induce depolymerization of the C-terminal peptide binding domain and stimulate the ATPase activity of the full-length protein.     

Flexible docking of peptides to proteins using CABS‐dock

Molecular docking of peptides to proteins can be a useful tool in the exploration of the possible peptide binding sites and poses. CABS‐dock is a method for protein–peptide docking that features significant conformational flexibility of both the peptide and the protein molecules during the peptide search for a binding site. The CABS‐dock has been made available as a web server and a standalone package. The web server is an easy to use tool with a simple web interface. The standalone package is a command‐line program dedicated to professional users. It offers a number of advanced features, analysis tools and support for large‐sized systems. In this article, we outline the current status of the CABS‐dock method, its recent developments, applications, and challenges ahead.     

Insight into the bovine milk peptide LPcin‐YK3 selection in the proteolytic system of Lactobacillus species

Antimicrobial peptides are class of small, positively charged peptides known for their broad‐spectrum antimicrobial activity. Antimicrobial activities for most antimicrobial peptides have largely remained elusive, particularly in the lactic acid bacteria. However, recently our investigation using LPcin‐YK3, an antimicrobial peptide from bovine milk, suggests that in vitro antimicrobial activity was reduced over 100‐fold compared with pathogenic bacteria. Additionally, for the structural study of how antimicrobial peptide undergoes its reaction at the proteolytic pathway of lactic acid bacteria based on degradation assay and propidium iodide staining, we performed molecular docking for interaction between oligopeptide‐binding protein A and LPcin‐YK3 peptide. Given that degradation related to the LPcin‐YK3 peptide in lactic acid bacteria proteolytic system, the inhibitory inactivity of LPcin‐YK3 against beneficial lactic acid bacteria strains may be one of the primary pharmacological properties of recombinant peptide discovered in bovine milk. These results provide structural and functional insights into the proteolytic mechanism and possibility as a putative substrate of oligopeptide‐binding protein A in respect of LPcin‐YK3 peptide.     

Structure-based identification of MHC binding peptides: Benchmarking of prediction accuracy

Identification of MHC binding peptides is essential for understanding the molecular mechanism of immune response. However, most of the prediction methods use motifs/profiles derived from experimental peptide binding data for specific MHC alleles, thus limiting their applicability only to those alleles for which such data is available. In this work we have developed a structure-based method which does not require experimental peptide binding data for training. Our method models MHC-peptide complexes using crystal structures of 170 MHC-peptide complexes and evaluates the binding energies using two well known residue based statistical pair potentials, namely Betancourt-Thirumalai (BT) and Miyazawa-Jernigan (MJ) matrices. Extensive benchmarking of prediction accuracy on a data set of 1654 epitopes from class I and class II alleles available in the SYFPEITHI database indicate that BT pair-potential can predict more than 60% of the known binders in case of 14 MHC alleles with AUC values for ROC curves ranging from 0.6 to 0.9. Similar benchmarking on 29,522 class I and class II MHC binding peptides with known IC(50) values in the IEDB database showed AUC values higher than 0.6 for 10 class I alleles and 9 class II alleles in predictions involving classification of a peptide to be binder or non-binder. Comparison with recently available benchmarking studies indicated that, the prediction accuracy of our method for many of the class I and class II MHC alleles was comparable to the sequence based methods, even if it does not use any experimental data for training. It is also encouraging to note that the ranks of true binding peptides could further be improved, when high scoring peptides obtained from pair potential were re-ranked using all atom forcefield and MM/PBSA method.     

Sub‐angstrom modeling of complexes between flexible peptides and globular proteins

A wide range of regulatory processes in the cell are mediated by flexible peptides that fold upon binding to globular proteins. Computational efforts to model these interactions are hindered by the large number of rotatable bonds in flexible peptides relative to typical ligand molecules, and the fact that different peptides assume different backbone conformations within the same binding site. In this study, we present Rosetta FlexPepDock, a novel tool for refining coarse peptide–protein models that allows significant changes in both peptide backbone and side chains. We obtain high resolution models, often of sub‐angstrom backbone quality, over an extensive and general benchmark that is based on a large nonredundant dataset of 89 peptide–protein interactions. Importantly, side chains of known binding motifs are modeled particularly well, typically with atomic accuracy. In addition, our protocol has improved modeling quality for the important application of cross docking to PDZ domains. We anticipate that the ability to create high resolution models for a wide range of peptide–protein complexes will have significant impact on structure‐based functional characterization, controlled manipulation of peptide interactions, and on peptide‐based drug design. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Modulation of the ATPase cycle of BiP by peptides and proteins

BiP, the Hsp70 homologue of the endoplasmic reticulum, interacts with its non-native substrate proteins in an ATP-dependent manner. This interaction is coupled to the ATPase cycle of the chaperone. Binding of short, synthetic peptides stimulate the ATPase activity of BiP. In previous work, we showed that a stably unfolded antibody domain forms a binary complex with BiP. In this study we made use of this complex to analyse the effect of substrate proteins on the ATPase cycle of BiP. Kinetic constants of the partial reactions of the ATPase cycle were determined without substrate, in the presence of a short binding peptide and in the presence of the antibody domain. We show that, in contrast to smaller peptides, the non-native protein domain decelerates the rate limiting hydrolysis step of the ATPase cycle.     

Peptide array‐based characterization and design of ZnO‐high affinity peptides

Peptides with both an affinity for ZnO and the ability to generate ZnO nanoparticles have attracted attention for the self‐assembly and templating of nanoscale building blocks under ambient conditions with compositional uniformity. In this study, we have analyzed the specific binding sites of the ZnO‐binding peptide, EAHVMHKVAPRP, which was identified using a phage display peptide library. The peptide binding assay against ZnO nanoparticles was performed using peptides synthesized on a cellulose membrane using the spot method. Using randomized rotation of amino acids in the ZnO‐binding peptide, 125 spot‐synthesized peptides were assayed. The peptide binding activity against ZnO nanoparticles varied greatly. This indicates that ZnO binding does not depend on total hydrophobicity or other physical parameters of these peptides, but rather that ZnO recognizes the specific amino acid alignment of these peptides. In addition, several peptides were found to show higher binding ability compared with that of the original peptides. Identification of important binding sites in the EAHVMHKVAPRP peptide was investigated by shortened, stepwise sequence from both termini. Interestingly, two ZnO‐binding sites were found as 6‐mer peptides: HVMHKV and HKVAPR. The peptides identified by amino acid substitution of HKVAPR were found to show high affinity and specificity for ZnO nanoparticles. Biotechnol. Bioeng. 2010;106: 845–851. © 2010 Wiley Periodicals, Inc.     

Dynamic Aha1 co‐chaperone binding to human Hsp90

Hsp90 is an essential chaperone that requires large allosteric changes to determine its ATPase activity and client binding. The co‐chaperone Aha1, which is the major ATPase stimulator in eukaryotes, is important for regulation of Hsp90's allosteric timing. Little is known, however, about the structure of the Hsp90/Aha1 complex. Here, we characterize the solution structure of unmodified human Hsp90/Aha1 complex using NMR spectroscopy. We show that the 214‐kDa complex forms by a two‐step binding mechanism and adopts multiple conformations in the absence of nucleotide. Aha1 induces structural changes near Hsp90's nucleotide‐binding site, providing a basis for its ATPase‐enhancing activity. Our data reveal important aspects of this pivotal chaperone/co‐chaperone interaction and emphasize the relevance of characterizing dynamic chaperone structures in solution.     

DynaDock: A new molecular dynamics‐based algorithm for protein–peptide docking including receptor flexibility

Molecular docking programs play an important role in drug development and many well‐established methods exist. However, there are two situations for which the performance of most approaches is still not satisfactory, namely inclusion of receptor flexibility and docking of large, flexible ligands like peptides. In this publication a new approach is presented for docking peptides into flexible receptors. For this purpose a two step procedure was developed: first, the protein–peptide conformational space is scanned and approximate ligand poses are identified and second, the identified ligand poses are refined by a new molecular dynamics‐based method, optimized potential molecular dynamics (OPMD). The OPMD approach uses soft‐core potentials for the protein–peptide interactions and applies a new optimization scheme to the soft‐core potential. Comparison with refinement results obtained by conventional molecular dynamics and a soft‐core scaling approach shows significant improvements in the sampling capability for the OPMD method. Thus, the number of starting poses needed for successful refinement is much lower than for the other methods. The algorithm was evaluated on 15 protein–peptide complexes with 2–16mer peptides. Docking poses with peptide RMSD values <2.10 Å from the equilibrated experimental structures were obtained in all cases. For four systems docking into the unbound receptor structures was performed, leading to peptide RMSD values <2.12 Å. Using a specifically fitted scoring function in 11 of 15 cases the best scoring poses featured a peptide RMSD ≤2.10 Å. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Predicting binding modes of reversible peptide‐based inhibitors of falcipain‐2 consistent with structure–activity relationships

Falcipain‐2 (FP‐2) is a major hemoglobinase of Plasmodium falciparum, considered an important drug target for the development of antimalarials. A previous study reported a novel series of 20 reversible peptide‐based inhibitors of FP‐2. However, the lack of tridimensional structures of the complexes hinders further optimization strategies to enhance the inhibitory activity of the compounds. Here we report the prediction of the binding modes of the aforementioned inhibitors to FP‐2. A computational approach combining previous knowledge on the determinants of binding to the enzyme, docking, and postdocking refinement steps, is employed. The latter steps comprise molecular dynamics simulations and free energy calculations. Remarkably, this approach leads to the identification of near‐native ligand conformations when applied to a validation set of protein‐ligand structures. Overall, we proposed substrate‐like binding modes of the studied compounds fulfilling the structural requirements for FP‐2 binding and yielding free energy values that correlated well with the experimental data. Proteins 2017; 85:1666–1683. © 2017 Wiley Periodicals, Inc.     

Analysis of molecular interactions in heat-induced aggregation of a non-inhibitory serpin ovalbumin using a molecular chaperone

Aggregation occurs through hydrophobic interactions when a polypeptide chain refolds in non-native states or when genetic variants of biologically active proteins assume inappropriate conformations, as observed in the case of dysfunctional serpins. Here, using the molecular chaperone BiP from bovine liver microsomes, we characterized the hydrophobic nature of the peptide segment which is considered to be a site required for aggregation among a non-inhibitory serpin ovalbumin in a heat-denatured state. Screening of the peptide scan for binding of BiP showed that BiP-binding sites are mostly buried in the folded ovalbumin. When ovalbumin was heat-denatured, the denatured protein was recognized by the antibody that reacts with the hydrophobic surface of the amino-terminal segment of ovalbumin. This antibody significantly suppressed the binding of BiP to denatured ovalbumin. BiP also bound the immobilized peptide in an ATP-dependent manner and the peptide stimulated the ATPase activity of BiP with a Km of 165 microM and a Vmax of 0.4 nmol/min per milligram. Measurement of surface plasmon resonance showed that the peptide had a Kd of 0.52 microM by BiP, lower than that for RCMLA (Kd = 1.1 microM) and even lower than that of the peptide P10K, PLSRTLSVAAKK, (Kd = 21 microM). These results demonstrate that the aggregation-prone site on heat-denatured ovalbumin has almost the same hydrophobic nature of interacting with the molecular chaperone BiP as the conventionally known peptides that bind to the Escherichia coli chaperone DnaK.     

Integrity of N‐ and C‐termini is important for E. coli Hsp31 chaperone activity

Hsp31 is a stress‐inducible molecular chaperone involved in the management of protein misfolding at high temperatures and in the development of acid resistance in starved E. coli. Each subunit of the Hsp31 homodimer consists of two structural domains connected by a flexible linker that sits atop a continuous tract of nonpolar residues adjacent to a hydrophobic bowl defined by the dimerization interface. Previously, we proposed that while the bowl serves as a binding site for partially folded species at physiological temperatures, chaperone function under heat shock conditions requires that folding intermediates further anneal to high‐affinity binding sites that become uncovered upon thermally induced motion of the linker. In support of a mechanism requiring that client proteins first bind to the bowl, we show here that fusion of a 20‐residue‐long hexahistidine tag to the N‐termini of Hsp31 abolishes chaperone activity at all temperatures by inducing reversible structural changes that interfere with substrate binding. We further demonstrate that extending the C‐termini of Hsp31 with short His tags selectively suppresses chaperone function at high temperatures by interfering with linker movement. The structural and functional sensitivity of Hsp31 to lengthening is consistent with the high degree of conservation of class I Hsp31 orthologs and will serve as a cautionary tale on the implications of affinity tagging.     

A Peptidic Unconjugated GRP78/BiP Ligand Modulates the Unfolded Protein Response and Induces Prostate Cancer Cell Death

The molecular chaperone GRP78/BiP is a key regulator of protein folding in the endoplasmic reticulum, and it plays a pivotal role in cancer cell survival and chemoresistance. Inhibition of its function has therefore been an important strategy for inhibiting tumor cell growth in cancer therapy. Previous efforts to achieve this goal have used peptides that bind to GRP78/BiP conjugated to pro-drugs or cell-death-inducing sequences. Here, we describe a peptide that induces prostate tumor cell death without the need of any conjugating sequences. This peptide is a sequence derived from the cochaperone Bag-1. We have shown that this sequence interacts with and inhibits the refolding activity of GRP78/BiP. Furthermore, we have demonstrated that it modulates the unfolded protein response in ER stress resulting in PARP and caspase-4 cleavage. Prostate cancer cells stably expressing this peptide showed reduced growth and increased apoptosis in in vivo xenograft tumor models. Amino acid substitutions that destroyed binding of the Bag-1 peptide to GRP78/BiP or downregulation of the expression of GRP78 compromised the inhibitory effect of this peptide. This sequence therefore represents a candidate lead peptide for anti-tumor therapy.     

The Large Hsp70 Grp170 Binds to Unfolded Protein Substrates in Vivo with a Regulation Distinct from Conventional Hsp70s

The Hsp70 superfamily is a ubiquitous chaperone class that includes conventional and large Hsp70s. BiP is the only conventional Hsp70 in the endoplasmic reticulum (ER) whose functions include: assisting protein folding, targeting misfolded proteins for degradation, and regulating the transducers of the unfolded protein response. The ER also possesses a single large Hsp70, the glucose-regulated protein of 170 kDa (Grp170). Like BiP it is an essential protein, but its cellular functions are not well understood. Here we show that Grp170 can bind directly to a variety of incompletely folded protein substrates in the ER, and as expected for a bona fide chaperone, it does not interact with folded secretory proteins. Our data demonstrate that Grp170 and BiP associate with similar molecular forms of two substrate proteins, but while BiP is released from unfolded substrates in the presence of ATP, Grp170 remains bound. In comparison to conventional Hsp70s, the large Hsp70s possess two unique structural features: an extended C-terminal α-helical domain and an unstructured loop in the putative substrate binding domain with an unknown function. We find that in the absence of the α-helical domain the interaction of Grp170 with substrates is reduced. In striking contrast, deletion of the unstructured loop results in increased binding to substrates, suggesting the presence of unique intramolecular mechanisms of control for the chaperone functions of large Hsp70s.     

Tumor suppressor Tsc1 is a new Hsp90 co‐chaperone that facilitates folding of kinase and non‐kinase clients

The tumor suppressors Tsc1 and Tsc2 form the tuberous sclerosis complex (TSC), a regulator of mTOR activity. Tsc1 stabilizes Tsc2; however, the precise mechanism involved remains elusive. The molecular chaperone heat‐shock protein 90 (Hsp90) is an essential component of the cellular homeostatic machinery in eukaryotes. Here, we show that Tsc1 is a new co‐chaperone for Hsp90 that inhibits its ATPase activity. The C‐terminal domain of Tsc1 (998–1,164 aa) forms a homodimer and binds to both protomers of the Hsp90 middle domain. This ensures inhibition of both subunits of the Hsp90 dimer and prevents the activating co‐chaperone Aha1 from binding the middle domain of Hsp90. Conversely, phosphorylation of Aha1‐Y223 increases its affinity for Hsp90 and displaces Tsc1, thereby providing a mechanism for equilibrium between binding of these two co‐chaperones to Hsp90. Our findings establish an active role for Tsc1 as a facilitator of Hsp90‐mediated folding of kinase and non‐kinase clients—including Tsc2—thereby preventing their ubiquitination and proteasomal degradation.