Gaussian process: an alternative approach for QSAM modeling of peptides

Different statistical modeling methods (SMMs) are used for nonlinear system classification and regression. On the basis of Bayesian probabilistic inference, Gaussian process (GP) is preliminarily used in the field of quantitative structure-activity relationship (QSAR) but has not yet been applied to quantitative sequence-activity model (QSAM) of biosystems. This paper proposes the application of GP as an alternative tool for the QSAM modeling of peptides. To investigate the modeling performance of GP, three classical peptide panels were used: Angiotensin-I converting enzyme inhibitory dipeptides, bradykinin-potentiating pentapeptides and cationic antimicrobial pentadecapeptides. On this basis, we made a comprehensive comparison between the GP and some widely used SMMs such as PLS, artificial neural network (ANN) and support vector machine (SVM), and gave the conclusions as follow: (1) for those of structurally complicated peptides, particularly the polypeptides, linear PLS was incapable of capturing all dependences hidden in the peptide systems, (2) even in assistance with the monitoring technique, ANN was inclined to be overtrained in the cases of insufficient number of peptide samples, (3) SVM and GP performed best for the three peptide panels. Moreover, since GP was able to correlate the linear and nonlinear-hybrid relationship, it was slightly superior to SVM at most peptide sets.     

Characterization of domain-peptide interaction interface: a case study on the amphiphysin-1 SH3 domain

Many important protein-protein interactions are mediated by peptide recognition modular domains, such as the Src homology 3 (SH3), SH2, PDZ, and WW domains. Characterizing the interaction interface of domain-peptide complexes and predicting binding specificity for modular domains are critical for deciphering protein-protein interaction networks. Here, we propose the use of an energetic decomposition analysis to characterize domain-peptide interactions and the molecular interaction energy components (MIECs), including van der Waals, electrostatic, and desolvation energy between residue pairs on the binding interface. We show a proof-of-concept study on the amphiphysin-1 SH3 domain interacting with its peptide ligands. The structures of the human amphiphysin-1 SH3 domain complexed with 884 peptides were first modeled using virtual mutagenesis and optimized by molecular mechanics (MM) minimization. Next, the MIECs between domain and peptide residues were computed using the MM/generalized Born decomposition analysis. We conducted two types of statistical analyses on the MIECs to demonstrate their usefulness for predicting binding affinities of peptides and for classifying peptides into binder and non-binder categories. First, combining partial least squares analysis and genetic algorithm, we fitted linear regression models between the MIECs and the peptide binding affinities on the training data set. These models were then used to predict binding affinities for peptides in the test data set; the predicted values have a correlation coefficient of 0.81 and an unsigned mean error of 0.39 compared with the experimentally measured ones. The partial least squares-genetic algorithm analysis on the MIECs revealed the critical interactions for the binding specificity of the amphiphysin-1 SH3 domain. Next, a support vector machine (SVM) was employed to build classification models based on the MIECs of peptides in the training set. A rigorous training-validation procedure was used to assess the performances of different kernel functions in SVM and different combinations of the MIECs. The best SVM classifier gave satisfactory predictions for the test set, indicated by average prediction accuracy rates of 78% and 91% for the binding and non-binding peptides, respectively. We also showed that the performance of our approach on both binding affinity prediction and binder/non-binder classification was superior to the performances of the conventional MM/Poisson-Boltzmann solvent-accessible surface area and MM/generalized Born solvent-accessible surface area calculations. Our study demonstrates that the analysis of the MIECs between peptides and the SH3 domain can successfully characterize the binding interface, and it provides a framework to derive integrated prediction models for different domain-peptide systems.     

Gaussian process: a promising approach for the modeling and prediction of Peptide binding affinity to MHC proteins

On the basis of Bayesian probabilistic inference, Gaussian process (GP) is a powerful machine learning method for nonlinear classification and regression, but has only very limited applications in the new areas of computational vaccinology and immunoinformatics. In the current work, we present a paradigmatic study of using GP regression technique to quantitatively model and predict the binding affinities of over 7000 immunodominant peptide epitopes to six types of human major histocompatibility complex (MHC) proteins. In this procedure, the sequence patterns of diverse peptides are characterized quantitatively and the resulting variables are then correlated with the experimentally measured affinities between different MHC and their peptide ligands, by using a linearity- and nonlinearity-hybrid GP approach. We also make systematical comparisons between the GP and two sophisticated modeling methods as partial least square (PLS) regression and support vector machine (SVM) with respect to their fitting ability, predictive power and generalization capability. The results suggest that GP could be a new and effective tool for the modeling and prediction of MHC-peptide interactions and would be promising in the field of computer-aided vaccine design (CAVD).     

Ligand-induced strain in hydrogen bonds of the c-Src SH3 domain detected by NMR

Changes in the molecular conformation of proteins can result from a variety of perturbations, and can play crucial roles in the regulation of biological activity. A new solution NMR method has been applied to monitor ligand-induced changes in hydrogen bond geometry in the chicken c-Src SH3 domain. The structural response of this domain to ligand binding has been investigated by measuring trans-hydrogen bond (15)N-(13)C' scalar couplings in the free state and when bound to the high affinity class I ligand RLP2, containing residues RALPPLPRY. A comparison between hydrogen bonds in high resolution X-ray structures of this domain and those observed via (h3)J(NC') couplings in solution shows remarkable agreement. Two backbone-to-side-chain hydrogen bonds are observed in solution, and each appears to play a role in stabilization of loop structure. Reproducible ligand-induced changes in trans-hydrogen bond scalar couplings are observed across the domain that translate into changes in hydrogen bond length ranging between 0.02 to 0.12 A. The observed changes can be rationalized by an induced fit mechanism in which hydrogen bonds across the protein participate in a compensatory response to forces imparted at the protein-ligand interface. Upon ligand binding, mutual intercalation of the two Leu-Pro segments of the ligand between three aromatic side-chains protruding from the SH3 surface wedges apart secondary structural elements within the SH3 domain. This disruption is transmitted in a domino-like effect across the domain through networks of hydrogen bonded peptide planes. The unprecedented resolution obtained demonstrates the ability to characterize subtle structural rearrangements within a protein upon perturbation, and represents a new step in the endeavor to understand how hydrogen bonds contribute to the stabilization and function of biological macromolecules.     

Acan125 binding to the SH3 domain of acanthamoeba myosin-IC

The domain organization of Acanthamoeba myosin-I, an oligomodular motor protein, includes a potentially important protein interaction module that is mostly uncharacterized. The Src homology 3, SH3, domain of myosin-I binds Acan125, a protein containing at least two consensus ligand binding domains: C-terminal SH3 binding motifs (PXXP) and N-terminal leucine-rich repeats. We report the first affinities determined for an SH3 domain of any myosin, namely, K(d) = 7 microM for a 21-residue synthetic peptide based on the PXXP domain sequence and K(d) = 0.15 microM for the PXXP domain included in the C-terminus of Acan125. These values are consistent with affinities reported for peptides and proteins that associate with SH3. By deletional analysis we show that only the PXXP domain is required for Acan125 to interact with the SH3 domain of Acanthamoeba myosin-IC (AmyoC(SH3)). The synthetic peptide described above at a concentration near the K(d) for SH3 binding blocked the interaction between native AmyoC and Acan125, mapping the interaction to the PXXP domain of Acan125 and the SH3 domain of myosin-I. These results are consistent with prototypical SH3 binding and suggest that a PXXP module is both necessary and sufficient to interact with an SH3 module of myosin-I.     

Investigation of the PDZ domain ligand binding site using chemically modified peptides

Several chemically modified analogues to a tightly binding ligand for the second PDZ domain of MAGI-3 were synthesized and evaluated for their ability to compete with native peptide ligands. N-methyl scanning of the ligand backbone amides revealed the energetically important hydrogen bonds between the ligand backbone and the PDZ domain. Analogues to the ligand's conserved threonine/serine(-2) residue, involved in a side chain to side chain hydrogen bond with a conserved histidine in the PDZ domain, revealed that the interaction is highly sensitive to the steric structure around the hydroxyl group of this residue. Analogues of the ligand carboxy terminus revealed that the full hydrogen bond network of the GLGF loop is important in ligand binding.     

The identification of conserved interactions within the SH3 domain by alignment of sequences and structures

The SH3 domain, comprised of approximately 60 residues, is found within a wide variety of proteins, and is a mediator of protein-protein interactions. Due to the large number of SH3 domain sequences and structures in the databases, this domain provides one of the best available systems for the examination of sequence and structural conservation within a protein family. In this study, a large and diverse alignment of SH3 domain sequences was constructed, and the pattern of conservation within this alignment was compared to conserved structural features, as deduced from analysis of eighteen different SH3 domain structures. Seventeen SH3 domain structures solved in the presence of bound peptide were also examined to identify positions that are consistently most important in mediating the peptide-binding function of this domain. Although residues at the two most conserved positions in the alignment are directly involved in peptide binding, residues at most other conserved positions play structural roles, such as stabilizing turns or comprising the hydrophobic core. Surprisingly, several highly conserved side-chain to main-chain hydrogen bonds were observed in the functionally crucial RT-Src loop between residues with little direct involvement in peptide binding. These hydrogen bonds may be important for maintaining this region in the precise conformation necessary for specific peptide recognition. In addition, a previously unrecognized yet highly conserved beta-bulge was identified in the second beta-strand of the domain, which appears to provide a necessary kink in this strand, allowing it to hydrogen bond to both sheets comprising the fold.     

Distinct peptide binding specificities of Src homology 3 (SH3) protein domains can be determined by modulation of local energetics across the binding interface

The yeast Nbp2p SH3 and Bem1p SH3b domains bind certain target peptides with similar high affinities, yet display vastly different affinities for other targets. To investigate this unusual behavior, we have solved the structure of the Nbp2p SH3-Ste20 peptide complex and compared it with the previously determined structure of the Bem1p SH3b bound to the same peptide. Although the Ste20 peptide interacts with both domains in a structurally similar manner, extensive in vitro studies with domain and peptide mutants revealed large variations in interaction strength across the binding interface of the two complexes. Whereas the Nbp2p SH3 made stronger contacts with the peptide core RXXPXXP motif, the Bem1p SH3b domain made stronger contacts with residues flanking the core motif. Remarkably, this modulation of local binding energetics can explain the distinct and highly nuanced binding specificities of these two domains.     

Evolution of High-Affinity Peptide Probes to Detect the SH3 Domain of Cancer Biomarker BCR–ABL

The BCR–ABL fusion protein is closely associated with the pathological progression of chronic myelogenous leukaemia and some other myeloproliferative diseases, which has long been recognized as one of the most important cancer biomarkers in the tumor diagnosis community. The SH3 domain of BCR–ABL is a small, conserved protein module that specifically recognizes and binds proline-rich peptide fragments. In the current study, we used a synthetic strategy to discover new peptide probes with high affinity binding to the BCR–ABL SH3 domain. In the procedure, a sequence-based machine learning predictor was developed based on a set of affinity-known SH3 binders, and the predictor was then used to guide the evolutional optimization of numerous virtual peptides to enrich high binding potency for the SH3 domain. Subsequently, a evolved peptide population was generated, from which ten peptides with the highest affinity scores were selected and their interaction free energies with SH3 domain were characterized systematically using a combination of molecular dynamics simulation and binding free energy analysis. Consequently, four peptides were suggested as promising candidates and their affinities toward SH3 domain were assayed; two peptides, APTYTPPPPP and APTYAPPPPP, were identified to have potent binding capability with dissociation constants K _d of 3 and 8 μM, respectively. Further, the structural basis and energetic property of SH3 domain in complex with APTYTPPPPP were examined in detail, revealing a non-specific interaction in SH3–peptide recognition that should render a broad ligand spectrum for the domain.     

Molecular details of Itk activation by prolyl isomerization and phospholigand binding: the NMR structure of the Itk SH2 domain bound to a phosphopeptide

The Src homology 2 (SH2) domain of interleukin-2 tyrosine kinase (Itk) is a critical component of the regulatory apparatus controlling the activity of this immunologically important enzyme. To gain insight into the structural features associated with the activated form of Itk, we have solved the NMR structure of the SH2 domain bound to a phosphotyrosine-containing peptide (pY) and analyzed changes in trans-hydrogen bond scalar couplings ((3h)J(NC')) that result from pY binding. Isomerization of a single prolyl imide bond in this domain is responsible for simultaneous existence of two distinct SH2 conformers. Prolyl isomerization directs ligand recognition: the trans conformer preferentially binds pY. The structure of the SH2/pY complex provides insight into the ligand specificity; the BG loop in the ligand-free trans SH2 conformer is pre-arranged for optimal contacts with the pY+3 residue of the ligand. Analysis of (3h)J(NC') couplings arising from hydrogen bonds has revealed propagation of structural changes from the pY binding pocket to the CD loop containing conformationally heterogeneous proline as well as to the alphaB helix, on the opposite site of the domain. These findings offer a structural framework for understanding the roles of prolyl isomerization and pY binding in Itk regulation.     

Novel Src homology 3 domain-binding motifs identified from proteomic screen of a Pro-rich region

The Src homology (SH) 3 domain has been shown recently to bind peptide sequences that lack the canonical PXXP motif. The diverse specificity in ligand recognition for a group of 15 SH3 domains has now been investigated using arrays of peptides derived from the proline-rich region of the SH2 domain-containing leukocyte protein of 76 kDa (SLP-76). A screen of the peptide arrays using individual or mixed SH3 domains has allowed the identification of a number of candidate SH3-binding peptides. Although some peptides contain the conventional PXXP motif, most are devoid of such a motif and are instead enriched in basic residues. Fluorescent polarization measurements using soluble peptides and purified SH3 domains demonstrated that several SH3 domains, including those from growth factor receptor-bound protein 2 (Grb2), NCK, and phospholipase C (PLC)-gamma1, bound with moderate affinities (10-100 microm) to a group of non-conventional peptides. Of particular interest, the PLC-gamma1 SH3 domain was found to associate with SLP-76 through at least three distinct sites, two of which bore a novel KKPP motif and the other contained the classic PXXP sequence. Intriguingly mutation of critical residues for the three sites not only affected binding of SLP-76 to the PLC-gamma1 SH3 domain but also to the Grb2 C-terminal SH3 domain, indicating that the binding sites in SLP-76 for the two SH3 domains are overlapped. Our studies suggest that the SH3 domain is an inherently promiscuous interaction module capable of binding to peptides that may or may not contain a PXXP motif. Furthermore the identification of numerous non-conventional SH3-binding peptides in SLP-76 implies that the global ligand pool for SH3 domains in a mammalian proteome may be significantly greater than previously acknowledged.     

Prediction of binding affinities between the human amphiphysin-1 SH3 domain and its peptide ligands using homology modeling, molecular dynamics and molecular field analysis

The SH3 domain of the human protein amphiphysin-1, which plays important roles in clathrin-mediated endocytosis, actin function and signaling transduction, can recognize peptide motif PXRPXR (X is any amino acid) with high affinity and specificity. We have constructed a complex structure of the amphiphysin-1 SH3 domain and a high-affinity peptide ligand PLPRRPPRA using homology modeling and molecular docking, which was optimized by molecular dynamics (MD). Three-dimensional quantitative structure-affinity relationship (3D-QSAR) analyses on the 200 peptides with known binding affinities to the amphiphysin-1 SH3 domain was then performed using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). The best CoMSIA model showed promising predictive power, giving good predictions for about 95% of the peptides in the test set (absolute prediction errors less than 1.0). It was used to validate peptide-SH3 binding structure and provide insight into the structural requirements for binding of peptides to SH3 domains. Finally, MD simulations were performed to analyze the interaction between the SH3 domain and another peptide GFPRRPPPRG that contains with the PXRPXsR (s represents residues with small side chains) motif. MD simulations demonstrated that the binding conformation of GFPRRPPPRG is quite different from that of PLPRRPPRAA especially the four residues at the C terminal, which may explain why the CoMSIA model cannot give good predictions on the peptides of the PXRPXsR motif. Because of its efficiency and predictive power, the 3D-QSAR model can be used as a scoring filter for predicting peptide sequences bound to SH3 domains.     

Binding of the Src SH2 domain to phosphopeptides is determined by residues in both the SH2 domain and the phosphopeptides.

Src homology 2 (SH2) domains are found in a variety of signaling proteins and bind phosphotyrosine-containing peptide sequences. To explore the binding properties of the SH2 domain of the Src protein kinase, we used immobilized phosphopeptides to bind purified glutathione S-transferase-Src SH2 fusion proteins. With this assay, as well as a free-peptide competition assay, we have estimated the affinities of the Src SH2 domain for various phosphopeptides relative to a Src SH2-phosphopeptide interaction whose Kd has been determined previously (YEEI-P; Kd = 4 nM). Two Src-derived phosphopeptides, one containing the regulatory C-terminal Tyr-527 and another containing the autophosphorylation site Tyr-416, bind the Src SH2 domain in a specific though low-affinity manner (with about 10(4)-lower affinity than the YEEI-P peptide). A platelet-derived growth factor receptor (PDGF-R) phosphopeptide containing Tyr-857 does not bind appreciably to the Src SH2 domain, suggesting it is not the PDGF-R binding site for Src as previously reported. However, another PDGF-R-derived phosphopeptide containing Tyr-751 does bind the Src SH2 domain (with an affinity approximately 2 orders of magnitude lower than that of YEEI-P). All of the phosphopeptides which bind to the Src SH2 domain contain a glutamic acid at position -3 or -4 with respect to phosphotyrosine; changing this residue to alanine greatly diminishes binding. We have also tested Src SH2 mutants for their binding properties and have interpreted our results in light of the recent crystal structure solution for the Src SH2 domain. Mutations in various conserved and nonconserved residues (R155A, R155K, N198E, H201R, and H201L) cause slight reductions in binding, while two mutations cause severe reductions. The W148E mutant domain, which alters the invariant tryptophan that marks the N-terminal border of the SH2 domain, binds poorly to phosphopeptides. Inclusion of the SH3 domain in the fusion protein partially restores the binding by the W148E mutant. A change in the invariant arginine that coordinates twice with phosphotyrosine in the peptide (R175L) results in a nearly complete loss of binding. The R175L mutant does display high affinity for the PDGF-R peptide containing Tyr-751, via an interaction that is at least partly phosphotyrosine independent. We have used this interaction to show that the R175L mutation also disrupts the intramolecular interaction between the Src SH2 domain and the phosphorylated C terminus within the context of the entire Src protein; thus, the binding properties observed for mutant domains in an in vitro assay appear to mimic those that occur in vivo.     

Solution structure of the first SH3 domain of human vinexin and its interaction with vinculin peptides

Solution structure of the first Src homology (SH) 3 domain of human vinexin (V_SH3_1) was determined using nuclear magnetic resonance (NMR) method and revealed that it was a canonical SH3 domain, which has a typical beta-beta-beta-beta-alpha-beta fold. Using chemical shift perturbation and surface plasmon resonance experiments, we studied the binding properties of the SH3 domain with two different peptides from vinculin hinge regions: P856 and P868. The observations illustrated slightly different affinities of the two peptides binding to V_SH3_1. The interaction between P868 and V_SH3_1 belonged to intermediate exchange with a modest binding affinity, while the interaction between P856 and V_SH3_1 had a low binding affinity. The structure and ligand-binding interface of V_SH3_1 provide a structural basis for the further functional study of this important molecule.