PeptideMine - A webserver for the design of peptides for protein-peptide binding studies derived from protein-protein interactomes

Background: Signal transduction events often involve transient, yet specific, interactions between structurally conserved protein domains and polypeptide sequences in target proteins. The identification and validation of these associating domains is crucial to understand signal transduction pathways that modulate different cellular or developmental processes. Bioinformatics strategies to extract and integrate information from diverse sources have been shown to facilitate the experimental design to understand complex biological events. These methods, primarily based on information from high-throughput experiments, have also led to the identification of new connections thus providing hypothetical models for cellular events. Such models, in turn, provide a framework for directing experimental efforts for validating the predicted molecular rationale for complex cellular processes. In this context, it is envisaged that the rational design of peptides for protein-peptide binding studies could substantially facilitate the experimental strategies to evaluate a predicted interaction. This rational design procedure involves the integration of protein-protein interaction data, gene ontology, physico-chemical calculations, domain-domain interaction data and information on functional sites or critical residues. Results: Here we describe an integrated approach called ``PeptideMine'' for the identification of peptides based on specific functional patterns present in the sequence of an interacting protein. This approach based on sequence searches in the interacting sequence space has been developed into a webserver, which can be used for the identification and analysis of peptides, peptide homologues or functional patterns from the interacting sequence space of a protein. To further facilitate experimental validation, the PeptideMine webserver also provides a list of physico-chemical parameters corresponding to the peptide to determine the feasibility of using the peptide for in vitro biochemical or biophysical studies. Conclusions: The strategy described here involves the integration of data and tools to identify potential interacting partners for a protein and design criteria for peptides based on desired biochemical properties. Alongside the search for interacting protein sequences using three different search programs, the server also provides the biochemical characteristics of candidate peptides to prune peptide sequences based on features that are most suited for a given experiment. The PeptideMine server is available at the URL: http://caps.ncbs.res.in/peptidemine

Dengue Envelope Domain III-Peptide Binding Analysis via Tryptophan Fluorescence Quenching Assay

In the efforts to find an anti-viral treatment for dengue, a simple tryptophan fluorescence-screening assay aimed at identifying dengue domain III envelope (EIII) protein inhibitors was developed. Residue Trp391 of EIII was used as an intrinsic probe to monitor the change in fluorescence of the tryptophan residue upon binding to a peptide. The analysis was based on the electron excitation at 280 nm and fluorescence emission at 300–400 nm of EIII, followed by quenching of fluorescence in the presence of potential peptidic inhibitors coded DS36wt, DS36opt, DN58wt and DN58opt. The present study found that the fluorescence of the recombinant EIII was quenched following the binding of DS36opt, DN58wt and DN58opt ina concentration-dependent manner. Since the λmax for emission remained unchanged, the effect was not dueto a change in the environment of the tryptophan side chain. In contrast, a minimal fluorescence-quenching effect of DS36wt at 20 and 40 µM suggested that the DS36wt does not have any binding ability to EIII. This was supported by a simple native-page gel retardation assay that showed a band shift of EIII domain whenincubated with DS36opt, DN58wt and DN58opt but not with DS36wt. We thus developed a low-cost and convenientspectrophotometric binding assay for the analysis of EIII–peptide interactions in a drug screening application.

Differential peptide binding to CD40 evokes counteractive responses.

The antigen-presenting cell-expressed CD40 is implied in the regulation of counteractive immune responses such as induction of pro-inflammatory and anti-inflammatory cytokines interleukin (IL)-12 and IL-10, respectively. The mechanism of this duality in CD40 function remains unknown. Here, we investigated whether such duality depends on ligand binding. Based on CD40 binding, we identifed two dodecameric peptides, peptide-7 and peptide-19, from the phage peptide library. Peptide-7 induces IL-10 and increases Leishmania donovani infection in macrophages, whereas peptide-19 induces IL-12 and reduces L. donovani infection. CD40-peptide interaction analyses by surface plasmon resonance and atomic force microscopy suggest that the functional differences are not associated with the studied interaction parameters. The molecular dynamic simulation of the CD40-peptides interaction suggests that these two peptides bind to two different places on CD40. Thus, we suggest for the first time that differential binding of the ligands imparts functional duality to CD40.

Statistical deconvolution of enthalpic energetic contributions to MHC-peptide binding affinity

Background: MHC Class I molecules present antigenic peptides to cytotoxic T cells, which forms an integral part of the adaptive immune response. Peptides are bound within a groove formed by the MHC heavy chain. Previous approaches to MHC Class I-peptide binding prediction have largely concentrated on the peptide anchor residues located at the P2 and C-terminus positions. Results: A large dataset comprising MHC-peptide structural complexes was created by remodelling pre-determined x-ray crystallographic structures. Static energetic analysis, following energy minimisation, was performed on the dataset in order to characterise interactions between bound peptides and the MHC Class I molecule, partitioning the interactions within the groove into van der Waals, electrostatic and total non-bonded energy contributions. Conclusion: The QSAR techniques of Genetic Function Approximation (GFA) and Genetic Partial Least Squares (G/PLS) algorithms were used to identify key interactions between the two molecules by comparing the calculated energy values with experimentally-determined BL50 data. Although the peptide termini binding interactions help ensure the stability of the MHC Class I-peptide complex, the central region of the peptide is also important in defining the specificity of the interaction. As thermodynamic studies indicate that peptide association and dissociation may be driven entropically, it may be necessary to incorporate entropic contributions into future calculations.

Solution structure and peptide binding studies of the C-terminal Src homology 3-like domain of the diphtheria toxin repressor protein

The diphtheria toxin repressor (DtxR) is the best-characterized member of a family of homologous proteins that regulate iron uptake and virulence gene expression in the Gram-positive bacteria. DtxR contains two domains that are separated by a short, unstructured linker. The N-terminal domain is structurally well-defined and is responsible for Fe2+ binding, dimerization, and DNA binding. The C-terminal domain adopts a fold similar to eukaryotic Src homology 3 domains, but the functional role of the C-terminal domain in repressor activity is unknown. The solution structure of the C-terminal domain, consisting of residues N130-L226 plus a 13-residue N-terminal extension, has been determined by using NMR spectroscopy. Residues before A147 are highly mobile and adopt a random coil conformation, but residues A147-L226 form a single structured domain consisting of five β-strands and three helices arranged into a partially orthogonal, two-sheet β-barrel, similar to the structure observed in the crystalline Co2+ complex of full-length DtxR. Chemical shift perturbation studies demonstrate that a proline-rich peptide corresponding to residues R125-G139 of intact DtxR binds to the C-terminal domain in a pocket formed by residues in β-strands 2, 3, and 5, and helix 3. Binding of the proline-rich peptide by the C-terminal domain of DtxR presents an example of peptide binding by a prokaryotic Src homology 3-like protein. The results of this study, combined with previous x-ray studies of intact DtxR, provide insights into a possible biological function of the C-terminal domain in regulating repressor activity.

Allosteric crosstalk between peptide-binding, transport, and ATP hydrolysis of the ABC transporter TAP

The transporter associated with antigen processing (TAP) is essential for intracellular transport of protein fragments into the endoplasmic reticulum for loading of major histocompatibility complex (MHC) class I molecules. On the cell surface, these peptide–MHC complexes are monitored by cytotoxic T lymphocytes. To study the ATP hydrolysis of TAP, we developed an enrichment and reconstitution procedure, by which we fully restored TAP function in proteoliposomes. A TAP-specific ATPase activity was identified that could be stimulated by peptides and blocked by the herpes simplex virus protein ICP47. Strikingly, the peptide-binding motif of TAP directly correlates with the stimulation of the ATPase activity, demonstrating that the initial peptide-binding step is responsible for TAP selectivity. ATP hydrolysis follows Michaelis–Menten kinetics with a maximal velocity Vmax of 2 μmol/min per mg TAP, corresponding to a turnover number of approximately 5 ATP per second. This turnover rate is sufficient to account for the role of TAP in peptide loading of MHC molecules and the overall process of antigen presentation. Interestingly, sterically restricted peptides that bind but are not transported by TAP do not stimulate ATPase activity. These results point to coordinated dialogue between the peptide-binding site, the nucleotide-binding domain, and the translocation site via conformational changes within the TAP complex.

Use of chimeric proteins to investigate the role of transporter associated with antigen processing (TAP) structural domains in peptide binding and translocation

The transporter associated with antigen processing (TAP) comprises two subunits, TAP1 and TAP2, each containing a hydrophobic membrane-spanning region (MSR) and a nucleotide binding domain (NBD). The TAP1/TAP2 complex is required for peptide translocation across the endoplasmic reticulum membrane. To understand the role of each structural unit of the TAP1/TAP2 complex, we generated two chimeras containing TAP1 MSR and TAP2 NBD (T1MT2C) or TAP2 MSR and TAP1 NBD (T2MT1C). We show that TAP1/T2MT1C, TAP2/T1MT2C, and T1MT2C/T2MT1C complexes bind peptide with an affinity comparable to wild-type complexes. By contrast, TAP1/T1MT2C and TAP2/T2MT1C complexes, although observed, are impaired for peptide binding. Thus, the MSRs of both TAP1 and TAP2 are required for binding peptide. However, neither NBD contains unique determinants required for peptide binding. The NBD-switched complexes, T1MT2C/T2MT1C, TAP1/T2MT1C, and TAP2/T1MT2C, all translocate peptides, but with progressively reduced efficiencies relative to the TAP1/TAP2 complex. These results indicate that both nucleotide binding sites are catalytically active and support an alternating catalytic sites model for the TAP transport cycle, similar to that proposed for P-glycoprotein. The enhanced translocation efficiency of TAP1/T2MT1C relative to TAP2/T1MT2C complexes correlates with enhanced binding of the TAP1 NBD-containing constructs to ATP-agarose beads. Preferential ATP interaction with TAP1, if occurring in vivo, might polarize the transport cycle such that ATP binding to TAP1 initiates the cycle. However, our observations that TAP complexes containing two identical TAP NBDs can mediate translocation indicate that distinct properties of the nucleotide binding site per se are not essential for the TAP catalytic cycle.

Mutations in the C-terminal fragment of DnaK affecting peptide binding.

Escherichia coli DnaK acts as a molecular chaperone through its ATP-regulated binding and release of polypeptide substrates. Overexpressing a C-terminal fragment (CTF) of DnaK (Gly-384 to Lys-638) containing the polypeptide substrate binding domain is lethal in wild-type E. coli. This dominant-negative phenotype may result from the nonproductive binding of CTF to cellular polypeptide targets of DnaK. Mutations affecting DnaK substrate binding were identified by selecting noncytotoxic CTF mutants followed by in vitro screening. The clustering of such mutations in the three-dimensional structure of CTF suggests the model that loops L1,2 and L4,5 form a rigid core structure critical for interactions with substrate.

Critical contribution of beta chain residue 57 in peptide binding ability of both HLA-DR and -DQ molecules.

Position 57 in the beta chain of HLA class II molecules maintains an Asp/non-Asp dimorphism that has been conserved through evolution and is implicated in susceptibility to some autoimmune diseases. The latter effect may be due to the influence of this residue on the ability of class II alleles to bind specific pathogenic peptides. We utilized highly homologous pairs of both DR and DQ alleles that varied at residue 57 to investigate the impact of this dimorphism on binding of model peptides. Using a direct binding assay of biotinylated peptides on whole cells expressing the desired alleles, we report several peptides that bind differentially to the allele pairs depending on the presence or absence of Asp at position 57. Peptides with negatively charged residues at anchor position 9 bind well to alleles not containing Asp at position 57 in the beta chain but cannot bind well to homologous Asp-positive alleles. By changing the peptides at the single residue predicted to interact with this position 57, we demonstrate a drastically altered or reversed pattern of binding. Ala analog peptides confirm these interactions and identify a limited set of interaction sites between the bound peptides and the class II molecules. Clarification of the impact of specific class II polymorphisms on generating unique allele-specific peptide binding "repertoires" will aid in our understanding of the development of specific immune responses and HLA-associated diseases.

The peptide binding site of the substance P (NK-1) receptor localized by a photoreactive analogue of substance P: presence of a disulfide bond.

Substance P (SP) is a neuropeptide that mediates multiple physiological responses including transmission of painful stimuli and inflammation via an interaction with a receptor of known primary sequence. To identify the regions of the SP receptor, also termed the NK-1 receptor, involved in peptide recognition, we are using analogues of SP containing the photoreactive amino acid p-benzoyl-L-phenylalanine (Bpa). In the present study, we used radioiodinated Bpa8-SP to covalently label with high efficiency the rat SP receptor expressed in a transfected mammalian cell line. To identify the amino acid residue that serves as the site of covalent attachment, a membrane preparation of labeled receptor was subjected to partial enzymatic cleavage by trypsin. A major digestion product of 22 kDa was identified. Upon reduction with 2-mercaptoethanol the mass of this product decreased to 14 kDa. The 22-kDa tryptic fragment was purified in excellent yield by preparative SDS/PAGE under nonreducing conditions. Subcleavage with Staphylococcus aureus V8 protease and endoproteinase ArgC yielded fragments of 8.2 and 9.0 kDa, respectively. Upon reductive cleavage, the V8 protease fragment decreased to 3.0 kDa while the endoproteinase ArgC fragment decreased to 3.2 kDa. Taking into consideration enzyme specificity, molecular size, determination of the presence or absence of N-glycosylation sites, and recognition by antibodies to specific sequences of the SP receptor, the V8 protease fragment is Thr-173 to Glu-183, while the endoproteinase ArgC fragment is Val-178 to Arg-190. These two fragments share the common sequence Val-Val-Cys-Met-Ile-Glu (residues 178-183). The site of covalent attachment of radioiodinated Bpa8-SP is thus restricted to a residue within this overlap sequence. The data presented here also establish that the cysteine residue in this sequence Cys-180, which is positioned in the middle of the second extracellular loop, participates in a disulfide bond that links the first and second extracellular loops of the receptor.

Structural features of the invariant chain fragment CLIP controlling rapid release from HLA-DR molecules and inhibition of peptide binding.

The invariant chain (Ii) prevents binding of ligands to major histocompatibility complex (MHC) class II molecules in the endoplasmic reticulum and during intracellular transport. Stepwise removal of the Ii in a trans-Golgi compartment renders MHC class II molecules accessible for peptide loading, with CLIP (class II-associated Ii peptides) as the final fragment to be released. Here we show that CLIP can be subdivided into distinct functional regions. The C-terminal segment (residues 92-105) of the CLIP-(81-105) fragment mediates inhibition of self- and antigenic peptide binding to HLA-DR2 molecules. In contrast, the N-terminal segment CLIP-(81-98) binds to the Staphylococcus aureus enterotoxin B contact site outside the peptide-binding groove on the alpha 1 domain and does not interfere with peptide binding. Its functional significance appears to lie in the contribution to CLIP removal: the dissociation of CLIP-(81-105) is characterized by a fast off-rate, which is accelerated at endosomal pH, whereas in the absence of the N-terminal CLIP-(81-91), the off-rate of C-terminal CLIP-(92-105) is slow and remains unaltered at low pH. Mechanistically, the N-terminal segment of CLIP seems to prevent tight interactions of CLIP side chains with specificity pockets in the peptide-binding groove that normally occurs during maturation of long-lived class II-peptide complexes.

Extreme learning machine for predicting HLA-peptide binding

Machine learning techniques have been recognized as powerful tools for learning from data. One of the most popular learning techniques, the Back-Propagation (BP) Artificial Neural Networks, can be used as a computer model to predict peptides binding to the Human Leukocyte Antigens (HLA). The major advantage of computational screening is that it reduces the number of wet-lab experiments that need to be performed, significantly reducing the cost and time. A recently developed method, Extreme Learning Machine (ELM), which has superior properties over BP has been investigated to accomplish such tasks. In our work, we found that the ELM is as good as, if not better than, the BP in term of time complexity, accuracy deviations across experiments, and most importantly - prevention from over-fitting for prediction of peptide binding to HLA.