The multiple-specificity landscape of modular peptide recognition domains

Anoctamin 9 determines Ca2+ signals during activation of T-lymphocytes

Background

Activation of T-cells is initiated by an increase in intracellular Ca2+, which underlies positive and negative regulation. Because the phospholipid scramblase and ion channel ANO9 (TMEM16J) was shown previously to regulated Ca2+ signals in renal epithelial cells, we asked whether ANO9 demonstrates a similar regulation in T-cells.

Methods

We used measurements of the intracellular Ca2+ concentration to examine the effects of ANO9 on intracellular Ca2+ signaling and demonstrated expression of ANO9 and its effects on cellular and molecular parameters.

Results

ANO9 was found to be expressed in human lymphocytes, including the Jurkat T-lymphocyte cell line and mouse lymphocytes. ANO9 has been shown to affect intracellular Ca2+ signals in renal epithelial cells. Here we demonstrate the essential role of ANO9 during initiation of intracellular Ca2+ signals in Jurkat T-cells and isolated mouse lymphocytes. ANO9 is essential for the initial rise in intracellular Ca2+ due to influx of extracellular Ca2+ through store-operated ORAI1 Ca2+ entry channels. ANO9 is indispensable for T-cell function, independent on whether cells are activated by stimulation of the T-cell receptor with CD3-antibody or by PMA/phytohemagglutinin.

Conclusions

Upon activation of T-cells and formation of the immunological synapse, ANO9 recruits the Ca2+-ATPase (PMCA) to the plasma membrane, which is supported by the scaffolding protein discs large 1 (DLG1). PMCAs maintain low Ca2+ levels near ORAI1 channels thereby suppressing Ca2+-inhibition of ORAI1 and thus retaining store-operated Ca2+ entry (SOCE). It is suggested that ANO9 has a role in interorganelle communication and regulation of cellular protein trafficking, which probably requires its phospholipid scramblase function.

A quantitative intracellular peptide-binding assay reveals recognition determinants and context dependence of short linear motifs

Transient protein-protein interactions play key roles in controlling dynamic cellular responses. Many examples involve globular protein domains that bind to peptide sequences known as short linear motifs (SLiMs), which are enriched in intrinsically disordered regions of proteins. Here we describe a novel functional assay for measuring SLiM binding, called systematic intracellular motif-binding analysis (SIMBA). In this method, binding of a foreign globular domain to its cognate SLiM peptide allows yeast cells to proliferate by blocking a growth arrest signal. A high-throughput application of the SIMBA method involving competitive growth and deep sequencing provides rapid quantification of the relative binding strength for thousands of SLiM sequence variants and a comprehensive interrogation of SLiM sequence features that control their recognition and potency. We show that multiple distinct classes of SLiM-binding domains can be analyzed by this method and that the relative binding strength of peptides in vivo correlates with their biochemical affinities measured in vitro. Deep mutational scanning provides high-resolution definitions of motif recognition determinants and reveals how sequence variations at noncore positions can modulate binding strength. Furthermore, mutational scanning of multiple parent peptides that bind human tankyrase ARC or YAP WW domains identifies distinct binding modes and uncovers context effects in which the preferred residues at one position depend on residues elsewhere. The findings establish SIMBA as a fast and incisive approach for interrogating SLiM recognition via massively parallel quantification of protein-peptide binding strength in vivo.

A quantitative intracellular peptide binding assay reveals recognition determinants and context dependence of short linear motifs

Transient protein-protein interactions play key roles in controlling dynamic cellular responses. Many examples involve globular protein domains that bind to peptide sequences known as Short Linear Motifs (SLiMs), which are enriched in intrinsically disordered regions of proteins. Here we describe a novel functional assay for measuring SLiM binding, called Systematic Intracellular Motif Binding Analysis (SIMBA). In this method, binding of a foreign globular domain to its cognate SLiM peptide allows yeast cells to proliferate by blocking a growth arrest signal. A high-throughput application of the SIMBA method involving competitive growth and deep sequencing provides rapid quantification of the relative binding strength for thousands of SLiM sequence variants, and a comprehensive interrogation of SLiM sequence features that control their recognition and potency. We show that multiple distinct classes of SLiM-binding domains can be analyzed by this method, and that the relative binding strength of peptides in vivo correlates with their biochemical affinities measured in vitro. Deep mutational scanning provides high-resolution definitions of motif recognition determinants and reveals how sequence variations at non-core positions can modulate binding strength. Furthermore, mutational scanning of multiple parent peptides that bind human tankyrase ARC or YAP WW domains identifies distinct binding modes and uncovers context effects in which the preferred residues at one position depend on residues elsewhere. The findings establish SIMBA as a fast and incisive approach for interrogating SLiM recognition via massively parallel quantification of protein-peptide binding strength in vivo.

MHCpLogics: an interactive machine learning-based tool for unsupervised data visualization and cluster analysis of immunopeptidomes

The major histocompatibility complex (MHC) encodes a range of immune response genes, including the human leukocyte antigens (HLAs) in humans. These molecules bind peptide antigens and present them on the cell surface for T cell recognition. The repertoires of peptides presented by HLA molecules are termed immunopeptidomes. The highly polymorphic nature of the genres that encode the HLA molecules confers allotype-specific differences in the sequences of bound ligands. Allotype-specific ligand preferences are often defined by peptide-binding motifs. Individuals express up to six classical class I HLA allotypes, which likely present peptides displaying different binding motifs. Such complex datasets make the deconvolution of immunopeptidomic data into allotype-specific contributions and further dissection of binding-specificities challenging. Herein, we developed MHCpLogics as an interactive machine learning-based tool for mining peptide-binding sequence motifs and visualization of immunopeptidome data across complex datasets. We showcase the functionalities of MHCpLogics by analyzing both in-house and published mono- and multi-allelic immunopeptidomics data. The visualization modalities of MHCpLogics allow users to inspect clustered sequences down to individual peptide components and to examine broader sequence patterns within multiple immunopeptidome datasets. MHCpLogics can deconvolute large immunopeptidome datasets enabling the interrogation of clusters for the segregation of allotype-specific peptide sequence motifs, identification of sub-peptidome motifs, and the exportation of clustered peptide sequence lists. The tool facilitates rapid inspection of immunopeptidomes as a resource for the immunology and vaccine communities. MHCpLogics is a standalone application available via an executable installation at: https://github.com/PurcellLab/MHCpLogics.

Computational analysis of substrate recognition of Sars-Cov-2 Mpro main protease

Mpro main protease takes an essential role in the Sars-Cov-2 viral life cycle by releasing the individual protein from the single poly-peptide chain via proteolytic cleavage in the beginning of the viral infection. Interfering with this step by inhibiting the protease with small compound-based inhibitors has been proven to be an effective strategy to treat the infection. Thus, understanding the substrate recognition mechanism of the Mpro main protease has gained great interest from the beginning of the pandemic. Here, we have studied the substrate recognition mechanism of the protease by means of the molecular dynamic methods. We have found that the glutamine residue at P1 has paramount effect in the interaction with the substrates as expected. In addition, we also have shown that for the first time, the arginine amino acid at the P3-P5 along with P4' can strengthen the interaction.

The MHC Motif Atlas: a database of MHC binding specificities and ligands

The highly polymorphic Major Histocompatibility Complex (MHC) genes are responsible for the binding and cell surface presentation of pathogen or cancer specific T-cell epitopes. This process is fundamental for eliciting T-cell recognition of infected or malignant cells. Epitopes displayed on MHC molecules further provide therapeutic targets for personalized cancer vaccines or adoptive T-cell therapy. To help visualizing, analyzing and comparing the different binding specificities of MHC molecules, we developed the MHC Motif Atlas (http://mhcmotifatlas.org/). This database contains information about thousands of class I and class II MHC molecules, including binding motifs, peptide length distributions, motifs of phosphorylated ligands, multiple specificities or links to X-ray crystallography structures. The database further enables users to download curated datasets of MHC ligands. By combining intuitive visualization of the main binding properties of MHC molecules together with access to more than a million ligands, the MHC Motif Atlas provides a central resource to analyze and interpret the binding specificities of MHC molecules.

NetTIME: a Multitask and Base-pair Resolution Framework for Improved Transcription Factor Binding Site Prediction

Motivation

Machine learning models for predicting cell-type-specific transcription factor (TF) binding sites have become increasingly more accurate thanks to the increased availability of next-generation sequencing data and more standardized model evaluation criteria. However, knowledge transfer from data-rich to data-limited TFs and cell types remains crucial for improving TF binding prediction models because available binding labels are highly skewed towards a small collection of TFs and cell types. Transfer prediction of TF binding sites can potentially benefit from a multitask learning approach; however, existing methods typically use shallow single-task models to generate low-resolution predictions. Here, we propose NetTIME, a multitask learning framework for predicting cell-type-specific TF binding sites with base-pair resolution.

Results

We show that the multitask learning strategy for TF binding prediction is more efficient than the single-task approach due to the increased data availability. NetTIME trains high-dimensional embedding vectors to distinguish TF and cell-type identities. We show that this approach is critical for the success of the multitask learning strategy and allows our model to make accurate transfer predictions within and beyond the training panels of TFs and cell types. We additionally train a linear-chain conditional random field (CRF) to classify binding predictions and show that this CRF eliminates the need for setting a probability threshold and reduces classification noise. We compare our method’s predictive performance with two state-of-the-art methods, Catchitt and Leopard, and show that our method outperforms previous methods under both supervised and transfer learning settings.

Availability and implementation

NetTIME is freely available at https://github.com/ryi06/NetTIME and the code is also archived at https://doi.org/10.5281/zenodo.6994897.

Supplementary information

Supplementary data are available at Bioinformatics online.

Protein context shapes the specificity of SH3 domain-mediated interactions in vivo

Protein–protein interactions (PPIs) between modular binding domains and their target peptide motifs are thought to largely depend on the intrinsic binding specificities of the domains. The large family of SRC Homology 3 (SH3) domains contribute to cellular processes via their ability to support such PPIs. While the intrinsic binding specificities of SH3 domains have been studied in vitro, whether each domain is necessary and sufficient to define PPI specificity in vivo is largely unknown. Here, by combining deletion, mutation, swapping and shuffling of SH3 domains and measurements of their impact on protein interactions in yeast, we find that most SH3s do not dictate PPI specificity independently from their host protein in vivo. We show that the identity of the host protein and the position of the SH3 domains within their host are critical for PPI specificity, for cellular functions and for key biophysical processes such as phase separation. Our work demonstrates the importance of the interplay between a modular PPI domain such as SH3 and its host protein in establishing specificity to wire PPI networks. These findings will aid understanding how protein networks are rewired during evolution and in the context of mutation-driven diseases such as cancer.

The SRC Homology 3 (SH3) domains mediate protein–protein interactions (PPIs). Here, the authors assess the SH3-mediated PPIs in yeast, and show that the identity of the protein itself and the position of the SH3 both affect the interaction specificity and thus the PPI-dependent cellular functions.

Deciphering the Unexpected Binding Capacity of the Third PDZ Domain of Whirlin to Various Cochlear Hair Cell Partners

Hearing is a mechanical and neurochemical process, which occurs in the hair cells of inner ear that converts the sound vibrations into electrical signals transmitted to the brain. The multi-PDZ scaffolding protein whirlin plays a critical role in the formation and function of stereocilia exposed at the surface of hair cells. In this article, we reported seven stereociliary proteins that encode PDZ binding motifs (PBM) and interact with whirlin PDZ3, where four of them are first reported. We solved the atomic resolution structures of complexes between whirlin PDZ3 and the PBMs of myosin 15a, CASK, harmonin a1 and taperin. Interestingly, the PBM of CASK and taperin are rare non-canonical PBM, which are not localized at the extreme C terminus. This large capacity to accommodate various partners could be related to the distinct functions of whirlin at different stages of the hair cell development.

Learning peptide recognition rules for a low-specificity protein

Many proteins interact with short linear regions of target proteins. For some proteins, however, it is difficult to identify a well-defined sequence motif that defines its target peptides. To overcome this difficulty, we used supervised machine learning to train a model that treats each peptide as a collection of easily-calculated biochemical features rather than as an amino acid sequence. As a test case, we dissected the peptide-recognition rules for human S100A5 (hA5), a low-specificity calcium binding protein. We trained a Random Forest model against a recently released, high-throughput phage display dataset collected for hA5. The model identifies hydrophobicity and shape complementarity, rather than polar contacts, as the primary determinants of peptide binding specificity in hA5. We tested this hypothesis by solving a crystal structure of hA5 and through computational docking studies of diverse peptides onto hA5. These structural studies revealed that peptides exhibit multiple binding modes at the hA5 peptide interface-all of which have few polar contacts with hA5. Finally, we used our trained model to predict new, plausible binding targets in the human proteome. This revealed a fragment of the protein α-1-syntrophin that binds to hA5. Our work helps better understand the biochemistry and biology of hA5, as well as demonstrating how high-throughput experiments coupled with machine learning of biochemical features can reveal the determinants of binding specificity in low-specificity proteins.

Protocol for Peptide Synthesis on Spectrally Encoded Beads for MRBLE-pep Assays

Every living cell relies on signal transduction pathways comprised of protein-protein interactions (PPIs). In many cases, these PPIs are between a folded protein domain and a short linear motif (SLiM) within an unstructured region of a protein. As a result of this small interaction interface (3-10 amino acids), the affinities of SLiM-mediated interactions are typically weak (K _ds of ~1-10 µM), allowing physiologically relevant changes in cellular concentrations of either protein partner to dictate changes in occupancy and thereby transmit cellular signals. However, these weak affinities also render detection and quantitative measurement of these interactions challenging and labor intensive. To address this, we recently developed MRBLE-pep, a technology that employs peptide libraries synthesized on spectrally encoded hydrogel beads to allow multiplexed affinity measurements between a protein and many different peptides in parallel. This approach dramatically reduces both the amount of protein and peptide as well as the time required to measure protein-peptide affinities compared to traditional methods. Here, we provide a detailed protocol describing how to: (1) functionalize polyethylene glycol diacrylate (PEG-DA) MRBLE beads with free amine groups, (2) synthesize peptide libraries on functionalized MRBLEs, (3) validate synthesized peptide sequences via MALDI mass spectrometry and quantify evenness of peptide coverage on MRBLEs, (4) use MRBLE-bound peptide libraries in multiplexed protein binding assays, and (5) analyze binding data to determine binding affinities. We anticipate that this protocol should prove useful for other researchers seeking to use MRBLE-pep in their own laboratories as well as for researchers broadly interested in solid-phase peptide synthesis and protein-protein binding assay development.

Specificity in PDZ-peptide interaction networks: Computational analysis and review

Globular PDZ domains typically serve as protein-protein interaction modules that regulate a wide variety of cellular functions via recognition of short linear motifs (SLiMs). Often, PDZ mediated-interactions are essential components of macromolecular complexes, and disruption affects the entire scaffold. Due to their roles as linchpins in trafficking and signaling pathways, PDZ domains are attractive targets: both for controlling viral pathogens, which bind PDZ domains and hijack cellular machinery, as well as for developing therapies to combat human disease. However, successful therapeutic interventions that avoid off-target effects are a challenge, because each PDZ domain interacts with a number of cellular targets, and specific binding preferences can be difficult to decipher. Over twenty-five years of research has produced a wealth of data on the stereochemical preferences of individual PDZ proteins and their binding partners. Currently the field lacks a central repository for this information. Here, we provide this important resource and provide a manually curated, comprehensive list of the 271 human PDZ domains. We use individual domain, as well as recent genomic and proteomic, data in order to gain a holistic view of PDZ domains and interaction networks, arguing this knowledge is critical to optimize targeting selectivity and to benefit human health.

Alteration of Amino Acid Profiling Influenced by the Active Ingredients of DanHong Injection After Prescription Optimization

Introduction

The aim of this work was to optimize the formulation composition of DanHong injection and to study the disturbance of microscopic components of cerebral ischemia in amino acid metabolites and metabolic pathways. The subtle relationship among these three substances and the influence of metabolic pathways were also studied.

Methods

In this study, the central composite design (CCD) matrix and response surface methodology (RSM) were used to design the experiments and to evaluate the interactive effects of three substances. Targeted metabolomics was used to detect the amino acid variation in CCD sets.

Results

Response surfaces were generated, and the formulation was optimized by superimposing the contour plots. It was found that the optimum values of the responses could be obtained at an SAB concentration (x1) of 8–9 mg/kg, a TSN concentration (x2) of 14–16 mg/kg, and an HSYA yellow A concentration (x3) of 6 mg/kg. Statistical analysis showed that the three independent variables had significant effects (p < 0.05) on the responses. A total of 22 experimental runs were performed, and the kinetic data were analyzed using a second-order polynomial. Model algorithm calculation indicated that glutamic acid, serine, leucine, glycine, and valine had a very close correlation with the active ingredients. Methionine, aspartic acid, asparagine, glutamic acid, and valine were important for distinguishing different groups, and they were identified as potential biomarkers. Cluster analysis and pathway analysis indicated that the valine, leucine, and isoleucine degradation (VLI degradation) pathway was the major metabolic pathway. Arginine and proline metabolites were most frequently detected, and they were closely associated with other networks according to the network analysis results. VLI degradation pathway and arginine and proline metabolism pathway had a significant influence on cerebral ischemia.

Discussion

The integration of CCD and metabolomics may be an effective strategy for optimizing the formulation composition and identifying the mechanism of action of traditional chinese medicine.

Robust prediction of HLA class II epitopes by deep motif deconvolution of immunopeptidomes

Predictions of epitopes presented by class II human leukocyte antigen molecules (HLA-II) have limited accuracy, restricting vaccine and therapy design. Here we combined unbiased mass spectrometry with a motif deconvolution algorithm to profile and analyze a total of 99,265 unique peptides eluted from HLA-II molecules. We then trained an epitope prediction algorithm with these data and improved prediction of pathogen and tumor-associated class II neoepitopes.

Quantitative mapping of protein-peptide affinity landscapes using spectrally encoded beads

Transient, regulated binding of globular protein domains to Short Linear Motifs (SLiMs) in disordered regions of other proteins drives cellular signaling. Mapping the energy landscapes of these interactions is essential for deciphering and perturbing signaling networks but is challenging due to their weak affinities. We present a powerful technology (MRBLE-pep) that simultaneously quantifies protein binding to a library of peptides directly synthesized on beads containing unique spectral codes. Using MRBLE-pep, we systematically probe binding of calcineurin (CN), a conserved protein phosphatase essential for the immune response and target of immunosuppressants, to the PxIxIT SLiM. We discover that flanking residues and post-translational modifications critically contribute to PxIxIT-CN affinity and identify CN-binding peptides based on multiple scaffolds with a wide range of affinities. The quantitative biophysical data provided by this approach will improve computational modeling efforts, elucidate a broad range of weak protein-SLiM interactions, and revolutionize our understanding of signaling networks.

GibbsCluster: unsupervised clustering and alignment of peptide sequences

Receptor interactions with short linear peptide fragments (ligands) are at the base of many biological signaling processes. Conserved and information-rich amino acid patterns, commonly called sequence motifs, shape and regulate these interactions. Because of the properties of a receptor-ligand system or of the assay used to interrogate it, experimental data often contain multiple sequence motifs. GibbsCluster is a powerful tool for unsupervised motif discovery because it can simultaneously cluster and align peptide data. The GibbsCluster 2.0 presented here is an improved version incorporating insertion and deletions accounting for variations in motif length in the peptide input. In basic terms, the program takes as input a set of peptide sequences and clusters them into meaningful groups. It returns the optimal number of clusters it identified, together with the sequence alignment and sequence motif characterizing each cluster. Several parameters are available to customize cluster analysis, including adjustable penalties for small clusters and overlapping groups and a trash cluster to remove outliers. As an example application, we used the server to deconvolute multiple specificities in large-scale peptidome data generated by mass spectrometry. The server is available at http://www.cbs.dtu.dk/services/GibbsCluster-2.0.

Accurate PDZ/Peptide Binding Specificity with Additive and Polarizable Free Energy Simulations

PDZ domains contain 80-100 amino acids and bind short C-terminal sequences of target proteins. Their specificity is essential for cellular signaling pathways. We studied the binding of the Tiam1 PDZ domain to peptides derived from the C-termini of its Syndecan-1 and Caspr4 targets. We used free energy perturbation (FEP) to characterize the binding energetics of one wild-type and 17 mutant complexes by simulating 21 alchemical transformations between pairs of complexes. Thirteen complexes had known experimental affinities. FEP is a powerful tool to understand protein/ligand binding. It depends, however, on the accuracy of molecular dynamics force fields and conformational sampling. Both aspects require continued testing, especially for ionic mutations. For six mutations that did not modify the net charge, we obtained excellent agreement with experiment using the additive, AMBER ff99SB force field, with a root mean square deviation (RMSD) of 0.37 kcal/mol. For six ionic mutations that modified the net charge, agreement was also good, with one large error (3 kcal/mol) and an RMSD of 0.9 kcal/mol for the other five. The large error arose from the overstabilization of a protein/peptide salt bridge by the additive force field. Four of the ionic mutations were also simulated with the polarizable Drude force field, which represents the first test of this force field for protein/ligand binding free energy changes. The large error was eliminated and the RMS error for the four mutations was reduced from 1.8 to 1.2 kcal/mol. The overall accuracy of FEP indicates it can be used to understand PDZ/peptide binding. Importantly, our results show that for ionic mutations in buried regions, electronic polarization plays a significant role.

Scribble: A master scaffold in polarity, adhesion, synaptogenesis, and proliferation

Key events ranging from cell polarity to proliferation regulation to neuronal signaling rely on the assembly of multiprotein adhesion or signaling complexes at particular subcellular sites. Multidomain scaffolding proteins nucleate assembly and direct localization of these complexes, and the protein Scribble and its relatives in the LAP protein family provide a paradigm for this. Scribble was originally identified because of its role in apical-basal polarity and epithelial integrity in Drosophila melanogaster It is now clear that Scribble acts to assemble and position diverse multiprotein complexes in processes ranging from planar polarity to adhesion to oriented cell division to synaptogenesis. Here, we explore what we have learned about the mechanisms of action of Scribble in the context of its multiple known interacting partners and discuss how this knowledge opens new questions about the full range of Scribble protein partners and their structural and signaling roles.

Allosteric Modulation of Binding Specificity by Alternative Packing of Protein Cores

Hydrophobic cores are often viewed as tightly packed and rigid, but they do show some plasticity and could thus be attractive targets for protein design. Here we explored the role of different functional pressures on the core packing and ligand recognition of the SH3 domain from human Fyn tyrosine kinase. We randomized the hydrophobic core and used phage display to select variants that bound to each of three distinct ligands. The three evolved groups showed remarkable differences in core composition, illustrating the effect of different selective pressures on the core. Changes in the core did not significantly alter protein stability, but were linked closely to changes in binding affinity and specificity. Structural analysis and molecular dynamics simulations revealed the structural basis for altered specificity. The evolved domains had significantly reduced core volumes, which in turn induced increased backbone flexibility. These motions were propagated from the core to the binding surface and induced significant conformational changes. These results show that alternative core packing and consequent allosteric modulation of binding interfaces could be used to engineer proteins with novel functions.

Mapping protein selectivity landscapes using multi-target selective screening and next-generation sequencing of combinatorial libraries

Characterizing the binding selectivity landscape of interacting proteins is crucial both for elucidating the underlying mechanisms of their interaction and for developing selective inhibitors. However, current mapping methods are laborious and cannot provide a sufficiently comprehensive description of the landscape. Here, we introduce a novel and efficient strategy for comprehensively mapping the binding landscape of proteins using a combination of experimental multi-target selective library screening and in silico next-generation sequencing analysis. We map the binding landscape of a non-selective trypsin inhibitor, the amyloid protein precursor inhibitor (APPI), to each of the four human serine proteases (kallikrein-6, mesotrypsin, and anionic and cationic trypsins). We then use this map to dissect and improve the affinity and selectivity of APPI variants toward each of the four proteases. Our strategy can be used as a platform for the development of a new generation of target-selective probes and therapeutic agents based on selective protein-protein interactions.

Combinatorial engineering of N-TIMP2 variants that selectively inhibit MMP9 and MMP14 function in the cell

Developing selective inhibitors for proteolytic enzymes that share high sequence homology and structural similarity is important for achieving high target affinity and functional specificity. Here, we used a combination of yeast surface display and dual-color selective library screening to obtain selective inhibitors for each of the matrix metalloproteinases (MMPs) MMP14 and MMP9 by modifying the non-specific N-terminal domain of the tissue inhibitor of metalloproteinase-2 (N-TIMP2). We generated inhibitor variants with 30- to 1175-fold improved specificity to each of the proteases, respectively, relative to wild type N-TIMP2. These biochemical results accurately predicted the selectivity and specificity obtained in cell-based assays. In U87MG cells, the activation of MMP2 by MMP14 was inhibited by MMP14-selective blockers but not MMP9-specific inhibitors. Target specificity was also demonstrated in MCF-7 cells stably expressing either MMP14 or MMP9, with only the MMP14-specific inhibitors preventing the mobility of MMP14-expressing cells. Similarly, the mobility of MMP9-expressing cells was inhibited by the MMP9-specific inhibitors, yet was not altered by the MMP14-specific inhibitors. The strategy developed in this study for improving the specificity of an otherwise broad-spectrum inhibitor will likely enhance our understanding of the basis for target specificity of inhibitors to proteolytic enzymes, in general, and to MMPs, in particular. We, moreover, envision that this study could serve as a platform for the development of next-generation, target-specific therapeutic agents. Finally, our methodology can be extended to other classes of proteolytic enzymes and other important target proteins.

Predicting Antigen Presentation-What Could We Learn From a Million Peptides?

Antigen presentation lies at the heart of immune recognition of infected or malignant cells. For this reason, important efforts have been made to predict which peptides are more likely to bind and be presented by the human leukocyte antigen (HLA) complex at the surface of cells. These predictions have become even more important with the advent of next-generation sequencing technologies that enable researchers and clinicians to rapidly determine the sequences of pathogens (and their multiple variants) or identify non-synonymous genetic alterations in cancer cells. Here, we review recent advances in predicting HLA binding and antigen presentation in human cells. We argue that the very large amount of high-quality mass spectrometry data of eluted (mainly self) HLA ligands generated in the last few years provides unprecedented opportunities to improve our ability to predict antigen presentation and learn new properties of HLA molecules, as demonstrated in many recent studies of naturally presented HLA-I ligands. Although major challenges still lie on the road toward the ultimate goal of predicting immunogenicity, these experimental and computational developments will facilitate screening of putative epitopes, which may eventually help decipher the rules governing T cell recognition.

The C-terminal extension landscape of naturally presented HLA-I ligands

HLA-I molecules play a central role in antigen presentation. They typically bind 9- to 12-mer peptides, and their canonical binding mode involves anchor residues at the second and last positions of their ligands. To investigate potential noncanonical binding modes, we collected in-depth and accurate HLA peptidomics datasets covering 54 HLA-I alleles and developed algorithms to analyze these data. Our results reveal frequent (442 unique peptides) and statistically significant C-terminal extensions for at least eight alleles, including the common HLA-A03:01, HLA-A31:01, and HLA-A68:01. High resolution crystal structure of HLA-A68:01 with such a ligand uncovers structural changes taking place to accommodate C-terminal extensions and helps unraveling sequence and structural properties predictive of the presence of these extensions. Scanning viral proteomes with the C-terminal extension motifs identifies many putative epitopes and we demonstrate direct recognition by human CD8⁺ T cells of a 10-mer epitope from cytomegalovirus predicted to follow the C-terminal extension binding mode.

Deciphering HLA-I motifs across HLA peptidomes improves neo-antigen predictions and identifies allostery regulating HLA specificity

The precise identification of Human Leukocyte Antigen class I (HLA-I) binding motifs plays a central role in our ability to understand and predict (neo-)antigen presentation in infectious diseases and cancer. Here, by exploiting co-occurrence of HLA-I alleles across ten newly generated as well as forty public HLA peptidomics datasets comprising more than 115,000 unique peptides, we show that we can rapidly and accurately identify many HLA-I binding motifs and map them to their corresponding alleles without any a priori knowledge of HLA-I binding specificity. Our approach recapitulates and refines known motifs for 43 of the most frequent alleles, uncovers new motifs for 9 alleles that up to now had less than five known ligands and provides a scalable framework to incorporate additional HLA peptidomics studies in the future. The refined motifs improve neo-antigen and cancer testis antigen predictions, indicating that unbiased HLA peptidomics data are ideal for in silico predictions of neo-antigens from tumor exome sequencing data. The new motifs further reveal distant modulation of the binding specificity at P2 for some HLA-I alleles by residues in the HLA-I binding site but outside of the B-pocket and we unravel the underlying mechanisms by protein structure analysis, mutagenesis and in vitro binding assays.

Predicting the differences between cancer and normal cells that are visible to the immune system is of central importance for cancer immunotherapy. Here we introduce a novel computational framework to harness the wealth of data from in-depth HLA peptidomics studies, including ten novel high quality (<1% FDR) datasets generated for this work, to improve predictions of peptides displayed on HLA-I molecules. These high-throughput and unbiased data enable us to refine models of HLA-I binding specificity for many alleles (including some that had no ligand until this study) and improve predictions of neo-antigens from exome sequencing data in melanoma and lung cancer samples. Moreover, the refined description of HLA-I binding specificity reveals cases of allosteric modulation of HLA-I binding specificity at the second amino acid position (P2) of their ligands by residues that are part of the HLA-I binding site but outside of the B pocket.

Comparing pairwise-additive and many-body generalized Born models for acid/base calculations and protein design

Generalized Born (GB) solvent models are common in acid/base calculations and protein design. With GB, the interaction between a pair of solute atoms depends on the shape of the protein/solvent boundary and, therefore, the positions of all solute atoms, so that GB is a many-body potential. For compute-intensive applications, the model is often simplified further, by introducing a mean, native-like protein/solvent boundary, which removes the many-body property. We investigate a method for both acid/base calculations and protein design that uses Monte Carlo simulations in which side chains can explore rotamers, bind/release protons, or mutate. The fluctuating protein/solvent dielectric boundary is treated in a way that is numerically exact (within the GB framework), in contrast to a mean boundary. Its originality is that it captures the many-body character while retaining the residue-pairwise complexity given by a fixed boundary. The method is implemented in the Proteus protein design software. It yields a slight but systematic improvement for acid/base constants in nine proteins and a significant improvement for the computational design of three PDZ domains. It eliminates a source of model uncertainty, which will facilitate the analysis of other model limitations. © 2017 Wiley Periodicals, Inc.

Computational Design of the Tiam1 PDZ Domain and Its Ligand Binding

PDZ domains direct protein-protein interactions and serve as models for protein design. Here, we optimized a protein design energy function for the Tiam1 and Cask PDZ domains that combines a molecular mechanics energy, Generalized Born solvent, and an empirical unfolded state model. Designed sequences were recognized as PDZ domains by the Superfamily fold recognition tool and had similarity scores comparable to natural PDZ sequences. The optimized model was used to redesign the two PDZ domains, by gradually varying the chemical potential of hydrophobic amino acids; the tendency of each position to lose or gain a hydrophobic character represents a novel hydrophobicity index. We also redesigned four positions in the Tiam1 PDZ domain involved in peptide binding specificity. The calculated affinity differences between designed variants reproduced experimental data and suggest substitutions with altered specificities.

Exhaustive search of linear information encoding protein-peptide recognition

High-throughput in vitro methods have been extensively applied to identify linear information that encodes peptide recognition. However, these methods are limited in number of peptides, sequence variation, and length of peptides that can be explored, and often produce solutions that are not found in the cell. Despite the large number of methods developed to attempt addressing these issues, the exhaustive search of linear information encoding protein-peptide recognition has been so far physically unfeasible. Here, we describe a strategy, called DALEL, for the exhaustive search of linear sequence information encoded in proteins that bind to a common partner. We applied DALEL to explore binding specificity of SH3 domains in the budding yeast Saccharomyces cerevisiae. Using only the polypeptide sequences of SH3 domain binding proteins, we succeeded in identifying the majority of known SH3 binding sites previously discovered either in vitro or in vivo. Moreover, we discovered a number of sites with both non-canonical sequences and distinct properties that may serve ancillary roles in peptide recognition. We compared DALEL to a variety of state-of-the-art algorithms in the blind identification of known binding sites of the human Grb2 SH3 domain. We also benchmarked DALEL on curated biological motifs derived from the ELM database to evaluate the effect of increasing/decreasing the enrichment of the motifs. Our strategy can be applied in conjunction with experimental data of proteins interacting with a common partner to identify binding sites among them. Yet, our strategy can also be applied to any group of proteins of interest to identify enriched linear motifs or to exhaustively explore the space of linear information encoded in a polypeptide sequence. Finally, we have developed a webserver located at http://michnick.bcm.umontreal.ca/dalel, offering user-friendly interface and providing different scenarios utilizing DALEL.

Here we describe the first strategy for the exhaustive search of the linear information encoding protein-peptide recognition; an approach that has previously been physically unfeasible because the combinatorial space of polypeptide sequences is too vast. The search covers the entire space of sequences with no restriction on motif length or composition, and includes all possible combinations of amino acids at distinct positions of each sequence, as well as positions with correlated preferences for amino acids.

Improved affinity at the cost of decreased specificity: a recurring theme in PDZ-peptide interactions

The E6 protein from human papillomavirus (HPV) plays an important role during productive infection and is a potential drug target. We have previously designed a high affinity bivalent protein binder for the E6 protein, a fusion between a helix from the E6 associated protein and PDZØ9, an engineered variant (L391F/K392M) of the second PDZ domain from synapse associated protein 97 (SAP97 PDZ2). How the substitutions improve the affinity of SAP97 PDZ2 for HPV E6 is not clear and it is not known to what extent they affect the specificity for cellular targets. Here, we explore the specificity of wild type SAP97 PDZ2 and PDZØ9 through proteomic peptide phage display. In addition, we employ a double mutant cycle of SAP97 PDZ2 in which the binding kinetics for nine identified potential cellular peptide ligands are measured and compared with those for the C-terminal E6 peptide. The results demonstrate that PDZØ9 has an increased affinity for all peptides, but at the cost of specificity. Furthermore, there is a peptide dependent coupling free energy between the side chains at positions 391 and 392. This corroborates our previous allosteric model for PDZ domains, involving sampling of intramolecular energetic pathways.

Unsupervised HLA Peptidome Deconvolution Improves Ligand Prediction Accuracy and Predicts Cooperative Effects in Peptide-HLA Interactions

Ag presentation on HLA molecules plays a central role in infectious diseases and tumor immunology. To date, large-scale identification of (neo-)Ags from DNA sequencing data has mainly relied on predictions. In parallel, mass spectrometry analysis of HLA peptidome is increasingly performed to directly detect peptides presented on HLA molecules. In this study, we use a novel unsupervised approach to assign mass spectrometry-based HLA peptidomics data to their cognate HLA molecules. We show that incorporation of deconvoluted HLA peptidomics data in ligand prediction algorithms can improve their accuracy for HLA alleles with few ligands in existing databases. The results of our computational analysis of large datasets of naturally processed HLA peptides, together with experimental validation and protein structure analysis, further reveal how HLA-binding motifs change with peptide length and predict new cooperative effects between distant residues in HLA-B07:02 ligands.

Saturation scanning of ubiquitin variants reveals a common hot spot for binding to USP2 and USP21

A detailed understanding of the molecular mechanisms whereby ubiquitin (Ub) recognizes enzymes in the Ub proteasome system is crucial for understanding the biological function of Ub. Many structures of Ub complexes have been solved and, in most cases, reveal a large structural epitope on a common face of the Ub molecule. However, owing to the generally weak nature of these interactions, it has been difficult to map in detail the functional contributions of individual Ub side chains to affinity and specificity. Here we took advantage of Ub variants (Ubvs) that bind tightly to particular Ub-specific proteases (USPs) and used phage display and saturation scanning mutagenesis to comprehensively map functional epitopes within the structural epitopes. We found that Ubvs that bind to USP2 or USP21 contain a remarkably similar core functional epitope, or "hot spot," consisting mainly of positions that are conserved as the wild type sequence, but also some positions that prefer mutant sequences. The Ubv core functional epitope contacts residues that are conserved in the human USP family, and thus it is likely important for the interactions of Ub across many family members.

Evolution of domain-peptide interactions to coadapt specificity and affinity to functional diversity

Evolution of complexity in eukaryotic proteomes has arisen, in part, through emergence of modular independently folded domains mediating protein interactions via binding to short linear peptides in proteins. Over 30 years, structural properties and sequence preferences of these peptides have been extensively characterized. Less successful, however, were efforts to establish relationships between physicochemical properties and functions of domain-peptide interactions. To our knowledge, we have devised the first strategy to exhaustively explore the binding specificity of protein domain-peptide interactions. We applied the strategy to SH3 domains to determine the properties of their binding peptides starting from various experimental data. The strategy identified the majority (∼70%) of experimentally determined SH3 binding sites. We discovered mutual relationships among binding specificity, binding affinity, and structural properties and evolution of linear peptides. Remarkably, we found that these properties are also related to functional diversity, defined by depth of proteins within hierarchies of gene ontologies. Our results revealed that linear peptides evolved to coadapt specificity and affinity to functional diversity of domain-peptide interactions. Thus, domain-peptide interactions follow human-constructed gene ontologies, which suggest that our understanding of biological process hierarchies reflect the way chemical and thermodynamic properties of linear peptides and their interaction networks, in general, have evolved.

PHASTpep: Analysis Software for Discovery of Cell-Selective Peptides via Phage Display and Next-Generation Sequencing

Next-generation sequencing has enhanced the phage display process, allowing for the quantification of millions of sequences resulting from the biopanning process. In response, many valuable analysis programs focused on specificity and finding targeted motifs or consensus sequences were developed. For targeted drug delivery and molecular imaging, it is also necessary to find peptides that are selective—targeting only the cell type or tissue of interest. We present a new analysis strategy and accompanying software, PHage Analysis for Selective Targeted PEPtides (PHASTpep), which identifies highly specific and selective peptides. Using this process, we discovered and validated, both in vitro and in vivo in mice, two sequences (HTTIPKV and APPIMSV) targeted to pancreatic cancer-associated fibroblasts that escaped identification using previously existing software. Our selectivity analysis makes it possible to discover peptides that target a specific cell type and avoid other cell types, enhancing clinical translatability by circumventing complications with systemic use.

Hammock: a hidden Markov model-based peptide clustering algorithm to identify protein-interaction consensus motifs in large datasets

Motivation: Proteins often recognize their interaction partners on the basis of short linear motifs located in disordered regions on proteins’ surface. Experimental techniques that study such motifs use short peptides to mimic the structural properties of interacting proteins. Continued development of these methods allows for large-scale screening, resulting in vast amounts of peptide sequences, potentially containing information on multiple protein-protein interactions. Processing of such datasets is a complex but essential task for large-scale studies investigating protein-protein interactions.

Results: The software tool presented in this article is able to rapidly identify multiple clusters of sequences carrying shared specificity motifs in massive datasets from various sources and generate multiple sequence alignments of identified clusters. The method was applied on a previously published smaller dataset containing distinct classes of ligands for SH3 domains, as well as on a new, an order of magnitude larger dataset containing epitopes for several monoclonal antibodies. The software successfully identified clusters of sequences mimicking epitopes of antibody targets, as well as secondary clusters revealing that the antibodies accept some deviations from original epitope sequences. Another test indicates that processing of even much larger datasets is computationally feasible.

Availability and implementation: Hammock is published under GNU GPL v. 3 license and is freely available as a standalone program (from http://www.recamo.cz/en/software/hammock-cluster-peptides/) or as a tool for the Galaxy toolbox (from https://toolshed.g2.bx.psu.edu/view/hammock/hammock). The source code can be downloaded from https://github.com/hammock-dev/hammock/releases.

Contact:muller@mou.cz

Supplementaryinformation:Supplementary data are available at Bioinformatics online.

High-throughput methods for identification of protein-protein interactions involving short linear motifs

Interactions between modular domains and short linear motifs (3–10 amino acids peptide stretches) are crucial for cell signaling. The motifs typically reside in the disordered regions of the proteome and the interactions are often transient, allowing for rapid changes in response to changing stimuli. The properties that make domain-motif interactions suitable for cell signaling also make them difficult to capture experimentally and they are therefore largely underrepresented in the known protein-protein interaction networks. Most of the knowledge on domain-motif interactions is derived from low-throughput studies, although there exist dedicated high-throughput methods for the identification of domain-motif interactions. The methods include arrays of peptides or proteins, display of peptides on phage or yeast, and yeast-two-hybrid experiments. We here provide a survey of scalable methods for domain-motif interaction profiling. These methods have frequently been applied to a limited number of ubiquitous domain families. It is now time to apply them to a broader set of peptide binding proteins, to provide a comprehensive picture of the linear motifs in the human proteome and to link them to their potential binding partners. Despite the plethora of methods, it is still a challenge for most approaches to identify interactions that rely on post-translational modification or context dependent or conditional interactions, suggesting directions for further method development.

Evolution of the SH3 Domain Specificity Landscape in Yeasts

To explore the conservation of Src homology 3 (SH3) domain-mediated networks in evolution, we compared the specificity landscape of these domains among four yeast species, Saccharomyces cerevisiae, Ashbya gossypii, Candida albicans, and Schizosaccharomyces pombe, encompassing 400 million years of evolution. We first aligned and catalogued the families of SH3-containing proteins in these four species to determine the relationships between homologous domains. Then, we tagged and purified all soluble SH3 domains (82 in total) to perform a quantitative peptide assay (SPOT) for each SH3 domain. All SPOT readouts were hierarchically clustered and we observed that the organization of the SH3 specificity landscape in three distinct profile classes remains conserved across these four yeast species. We also produced a specificity profile for each SH3 domain from manually aligned top SPOT hits and compared the within-family binding motif consensus. This analysis revealed a striking example of binding motif divergence in a C. albicans Rvs167 paralog, which cannot be explained by overall SH3 sequence or interface residue divergence, and we validated this specificity change with a yeast two-hybrid (Y2H) assay. In addition, we show that position-weighted matrices (PWM) compiled from SPOT assays can be used for binding motif screening in potential binding partners and present cases where motifs are either conserved or lost among homologous SH3 interacting proteins. Finally, by comparing pairwise SH3 sequence identity to binding profile correlation we show that for ~75% of all analyzed families the SH3 specificity profile was remarkably conserved over a large evolutionary distance. Thus, a high sequence identity within an SH3 domain family predicts conserved binding specificity, whereas divergence in sequence identity often coincided with a change in binding specificity within this family. As such, our results are important for future studies aimed at unraveling complex specificity networks of peptide recognition domains in higher eukaryotes, including mammals.

Computational design of selective peptides to discriminate between similar PDZ domains in an oncogenic pathway

Reagents that target protein-protein interactions to rewire signaling are of great relevance in biological research. Computational protein design may offer a means of creating such reagents on demand, but methods for encoding targeting selectivity are sorely needed. This is especially challenging when targeting interactions with ubiquitous recognition modules--for example, PDZ domains, which bind C-terminal sequences of partner proteins. Here we consider the problem of designing selective PDZ inhibitor peptides in the context of an oncogenic signaling pathway, in which two PDZ domains (NHERF-2 PDZ2-N2P2 and MAGI-3 PDZ6-M3P6) compete for a receptor C-terminus to differentially modulate oncogenic activities. Because N2P2 has been shown to increase tumorigenicity and M3P6 to decreases it, we sought to design peptides that inhibit N2P2 without affecting M3P6. We developed a structure-based computational design framework that models peptide flexibility in binding yet is efficient enough to rapidly analyze tradeoffs between affinity and selectivity. Designed peptides showed low-micromolar inhibition constants for N2P2 and no detectable M3P6 binding. Peptides designed for reverse discrimination bound M3P6 tighter than N2P2, further testing our technology. Experimental and computational analysis of selectivity determinants revealed significant indirect energetic coupling in the binding site. Successful discrimination between N2P2 and M3P6, despite their overlapping binding preferences, is highly encouraging for computational approaches to selective PDZ targeting, especially because design relied on a homology model of M3P6. Still, we demonstrate specific deficiencies of structural modeling that must be addressed to enable truly robust design. The presented framework is general and can be applied in many scenarios to engineer selective targeting.

The development and application of a quantitative peptide microarray based approach to protein interaction domain specificity space

Protein interaction domain (PID) linear peptide motif interactions direct diverse cellular processes in a specific and coordinated fashion. PID specificity, or the interaction selectivity derived from affinity preferences between possible PID-peptide pairs is the basis of this ability. Here, we develop an integrated experimental and computational cellulose peptide conjugate microarray (CPCMA) based approach for the high throughput analysis of PID specificity that provides unprecedented quantitative resolution and reproducibility. As a test system, we quantify the specificity preferences of four Src Homology 2 domains and 124 physiological phosphopeptides to produce a novel quantitative interactome. The quantitative data set covers a broad affinity range, is highly precise, and agrees well with orthogonal biophysical validation, in vivo interactions, and peptide library trained algorithm predictions. In contrast to preceding approaches, the CPCMAs proved capable of confidently assigning interactions into affinity categories, resolving the subtle affinity contributions of residue correlations, and yielded predictive peptide motif affinity matrices. Unique CPCMA enabled modes of systems level analysis reveal a physiological interactome with expected node degree value decreasing as a function of affinity, resulting in minimal high affinity binding overlap between domains; uncover that Src Homology 2 domains bind ligands with a similar average affinity yet strikingly different levels of promiscuity and binding dynamic range; and parse with unprecedented quantitative resolution contextual factors directing specificity. The CPCMA platform promises broad application within the fields of PID specificity, synthetic biology, specificity focused drug design, and network biology.

Prediction and experimental characterization of nsSNPs altering human PDZ-binding motifs

Single nucleotide polymorphisms (SNPs) are a major contributor to genetic and phenotypic variation within populations. Non-synonymous SNPs (nsSNPs) modify the sequence of proteins and can affect their folding or binding properties. Experimental analysis of all nsSNPs is currently unfeasible and therefore computational predictions of the molecular effect of nsSNPs are helpful to guide experimental investigations. While some nsSNPs can be accurately characterized, for instance if they fall into strongly conserved or well annotated regions, the molecular consequences of many others are more challenging to predict. In particular, nsSNPs affecting less structured, and often less conserved regions, are difficult to characterize. Binding sites that mediate protein-protein or other protein interactions are an important class of functional sites on proteins and can be used to help interpret nsSNPs. Binding sites targeted by the PDZ modular peptide recognition domain have recently been characterized. Here we use this data to show that it is possible to computationally identify nsSNPs in PDZ binding motifs that modify or prevent binding to the proteins containing the motifs. We confirm these predictions by experimentally validating a selected subset with ELISA. Our work also highlights the importance of better characterizing linear motifs in proteins as many of these can be affected by genetic variations.

Structure-based multiscale approach for identification of interaction partners of PDZ domains

PDZ domains are peptide recognition modules which mediate specific protein-protein interactions and are known to have a complex specificity landscape. We have developed a novel structure-based multiscale approach which identifies crucial specificity determining residues (SDRs) of PDZ domains from explicit solvent molecular dynamics (MD) simulations on PDZ-peptide complexes and uses these SDRs in combination with knowledge-based scoring functions for proteomewide identification of their interaction partners. Multiple explicit solvent simulations ranging from 5 to 50 ns duration have been carried out on 28 PDZ-peptide complexes with known binding affinities. MM/PBSA binding energy values calculated from these simulations show a correlation coefficient of 0.755 with the experimental binding affinities. On the basis of the SDRs of PDZ domains identified by MD simulations, we have developed a simple scoring scheme for evaluating binding energies for PDZ-peptide complexes using residue based statistical pair potentials. This multiscale approach has been benchmarked on a mouse PDZ proteome array data set by calculating the binding energies for 217 different substrate peptides in binding pockets of 64 different mouse PDZ domains. Receiver operating characteristic (ROC) curve analysis indicates that, the area under curve (AUC) values for binder vs nonbinder classification by our structure based method is 0.780. Our structure based method does not require experimental PDZ-peptide binding data for training.

Characterization of diverse internal binding specificities of PDZ domains by yeast two-hybrid screening of a special peptide library

Protein-protein interactions (PPIs) are essential events to play important roles in a series of biological processes. There are probably more ways of PPIs than we currently realized. Structural and functional investigations of weak PPIs have lagged behind those of strong PPIs due to technical difficulties. Weak PPIs are often short-lived, which may result in more dynamic signals with important biological roles within and/or between cells. For example, the characteristics of PSD-95/Dlg/ZO-1 (PDZ) domain binding to internal sequences, which are primarily weak interactions, have not yet been systematically explored. In the present study, we constructed a nearly random octapeptide yeast two-hybrid library. A total of 24 PDZ domains were used as baits for screening the library. Fourteen of these domains were able to bind internal PDZ-domain binding motifs (PBMs), and PBMs screened for nine PDZ domains exhibited strong preferences. Among 11 PDZ domains that have not been reported their internal PBM binding ability, six were confirmed to bind internal PBMs. The first PDZ domain of LNX2, which has not been reported to bind C-terminal PBMs, was found to bind internal PBMs. These results suggest that the internal PBMs binding ability of PDZ domains may have been underestimated. The data provided diverse internal binding properties for several PDZ domains that may help identify their novel binding partners.

Large-scale interaction profiling of PDZ domains through proteomic peptide-phage display using human and viral phage peptidomes

The human proteome contains a plethora of short linear motifs (SLiMs) that serve as binding interfaces for modular protein domains. Such interactions are crucial for signaling and other cellular processes, but are difficult to detect because of their low to moderate affinities. Here we developed a dedicated approach, proteomic peptide-phage display (ProP-PD), to identify domain-SLiM interactions. Specifically, we generated phage libraries containing all human and viral C-terminal peptides using custom oligonucleotide microarrays. With these libraries we screened the nine PSD-95/Dlg/ZO-1 (PDZ) domains of human Densin-180, Erbin, Scribble, and Disks large homolog 1 for peptide ligands. We identified several known and putative interactions potentially relevant to cellular signaling pathways and confirmed interactions between full-length Scribble and the target proteins β-PIX, plakophilin-4, and guanylate cyclase soluble subunit α-2 using colocalization and coimmunoprecipitation experiments. The affinities of recombinant Scribble PDZ domains and the synthetic peptides representing the C termini of these proteins were in the 1- to 40-μM range. Furthermore, we identified several well-established host-virus protein-protein interactions, and confirmed that PDZ domains of Scribble interact with the C terminus of Tax-1 of human T-cell leukemia virus with micromolar affinity. Previously unknown putative viral protein ligands for the PDZ domains of Scribble and Erbin were also identified. Thus, we demonstrate that our ProP-PD libraries are useful tools for probing PDZ domain interactions. The method can be extended to interrogate all potential eukaryotic, bacterial, and viral SLiMs and we suggest it will be a highly valuable approach for studying cellular and pathogen-host protein-protein interactions.

A graph kernel approach for alignment-free domain-peptide interaction prediction with an application to human SH3 domains

MOTIVATION

State-of-the-art experimental data for determining binding specificities of peptide recognition modules (PRMs) is obtained by high-throughput approaches like peptide arrays. Most prediction tools applicable to this kind of data are based on an initial multiple alignment of the peptide ligands. Building an initial alignment can be error-prone, especially in the case of the proline-rich peptides bound by the SH3 domains.

RESULTS

Here, we present a machine-learning approach based on an efficient graph-kernel technique to predict the specificity of a large set of 70 human SH3 domains, which are an important class of PRMs. The graph-kernel strategy allows us to (i) integrate several types of physico-chemical information for each amino acid, (ii) consider high-order correlations between these features and (iii) eliminate the need for an initial peptide alignment. We build specialized models for each human SH3 domain and achieve competitive predictive performance of 0.73 area under precision-recall curve, compared with 0.27 area under precision-recall curve for state-of-the-art methods based on position weight matrices. We show that better models can be obtained when we use information on the noninteracting peptides (negative examples), which is currently not used by the state-of-the art approaches based on position weight matrices. To this end, we analyze two strategies to identify subsets of high confidence negative data. The techniques introduced here are more general and hence can also be used for any other protein domains, which interact with short peptides (i.e. other PRMs).

AVAILABILITY

The program with the predictive models can be found at http://www.bioinf.uni-freiburg.de/Software/SH3PepInt/SH3PepInt.tar.gz. We also provide a genome-wide prediction for all 70 human SH3 domains, which can be found under http://www.bioinf.uni-freiburg.de/Software/SH3PepInt/Genome-Wide-Predictions.tar.gz.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Peptide docking and structure-based characterization of peptide binding: from knowledge to know-how

Peptide-mediated interactions are gaining increased attention due to their predominant roles in the many regulatory processes that involve dynamic interactions between proteins. The structures of such interactions provide an excellent starting point for their characterization and manipulation, and can provide leads for targeted inhibitor design. The relatively few experimentally determined structures of peptide-protein complexes can be complemented with an outburst of modeling approaches that have been introduced in recent years, with increasing accuracy and applicability to ever more systems. We review different methods to address the considerable challenges in modeling the binding of a short yet highly flexible peptide to its partner. These methods apply an array of sampling strategies and draw from a recent amassing of knowledge about the biophysical nature of peptide-protein interactions. We elaborate on applications of these structure-based approaches and in particular on the characterization of peptide binding specificity to different peptide-binding domains and enzymes. Such applications can identify new biological targets and thus complement our current view of protein-protein interactions in living organisms. Accurate peptide-protein docking is of particular importance in the light of increased appreciation of the crucial functional roles of disordered regions and the many linear binding motifs embedded within.

The SH2 domain interaction landscape

Members of the SH2 domain family modulate signal transduction by binding to short peptides containing phosphorylated tyrosines. Each domain displays a distinct preference for the sequence context of the phosphorylated residue. We have developed a high-density peptide chip technology that allows for probing of the affinity of most SH2 domains for a large fraction of the entire complement of tyrosine phosphopeptides in the human proteome. Using this technique, we have experimentally identified thousands of putative SH2-peptide interactions for more than 70 different SH2 domains. By integrating this rich data set with orthogonal context-specific information, we have assembled an SH2-mediated probabilistic interaction network, which we make available as a community resource in the PepspotDB database. A predicted dynamic interaction between the SH2 domains of the tyrosine phosphatase SHP2 and the phosphorylated tyrosine in the extracellular signal-regulated kinase activation loop was validated by experiments in living cells.

Exploring the binding diversity of intrinsically disordered proteins involved in one-to-many binding

Molecular recognition features (MoRFs) are intrinsically disordered protein regions that bind to partners via disorder-to-order transitions. In one-to-many binding, a single MoRF binds to two or more different partners individually. MoRF-based one-to-many protein-protein interaction (PPI) examples were collected from the Protein Data Bank, yielding 23 MoRFs bound to 2-9 partners, with all pairs of same-MoRF partners having less than 25% sequence identity. Of these, 8 MoRFs were bound to 2-9 partners having completely different folds, whereas 15 MoRFs were bound to 2-5 partners having the same folds but with low sequence identities. For both types of partner variation, backbone and side chain torsion angle rotations were used to bring about the conformational changes needed to enable close fits between a single MoRF and distinct partners. Alternative splicing events (ASEs) and posttranslational modifications (PTMs) were also found to contribute to distinct partner binding. Because ASEs and PTMs both commonly occur in disordered regions, and because both ASEs and PTMs are often tissue-specific, these data suggest that MoRFs, ASEs, and PTMs may collaborate to alter PPI networks in different cell types. These data enlarge the set of carefully studied MoRFs that use inherent flexibility and that also use ASE-based and/or PTM-based surface modifications to enable the same disordered segment to selectively associate with two or more partners. The small number of residues involved in MoRFs and in their modifications by ASEs or PTMs may simplify the evolvability of signaling network diversity.