Identification of preferred protein interactions by phage-display of the human Src homology-3 proteome

We have determined the human genome to contain 296 different Src homology-3 (SH3) domains and cloned them into a phage-display vector. This provided a powerful and unbiased system for simultaneous assaying of the complete human SH3 proteome for the strongest binding to target proteins of interest, without the limitations posed by short linear peptide ligands or confounding variables of more indirect methods for protein interaction screening. Studies involving three ligand proteins, human immunodeficiency virus-1 Nef, p21-activated kinase (PAK)2 and ADAM15, showed previously reported as well as novel SH3 partners with nanomolar affinities specific for them. This argues that SH3 domains may have a more dominant role in directing cellular protein interactions than has been assumed. Besides showing potentially important new SH3-directed interactions, these studies also led to the discovery of novel signalling proteins, such as the PAK2-binding adaptor protein POSH2 and the ADAM15-binding sorting nexin family member SNX30.     

Accurate Prediction of Peptide Binding Sites on Protein Surfaces

Many important protein–protein interactions are mediated by the binding of a short peptide stretch in one protein to a large globular segment in another. Recent efforts have provided hundreds of examples of new peptides binding to proteins for which a three-dimensional structure is available (either known experimentally or readily modeled) but where no structure of the protein–peptide complex is known. To address this gap, we present an approach that can accurately predict peptide binding sites on protein surfaces. For peptides known to bind a particular protein, the method predicts binding sites with great accuracy, and the specificity of the approach means that it can also be used to predict whether or not a putative or predicted peptide partner will bind. We used known protein–peptide complexes to derive preferences, in the form of spatial position specific scoring matrices, which describe the binding-site environment in globular proteins for each type of amino acid in bound peptides. We then scan the surface of a putative binding protein for sites for each of the amino acids present in a peptide partner and search for combinations of high-scoring amino acid sites that satisfy constraints deduced from the peptide sequence. The method performed well in a benchmark and largely agreed with experimental data mapping binding sites for several recently discovered interactions mediated by peptides, including RG-rich proteins with SMN domains, Epstein-Barr virus LMP1 with TRADD domains, DBC1 with Sir2, and the Ago hook with Argonaute PIWI domain. The method, and associated statistics, is an excellent tool for predicting and studying binding sites for newly discovered peptides mediating critical events in biology.     

Ubiquitin binds to and regulates a subset of SH3 domains

SH3 domains are modules of 50-70 amino acids that promote interactions among proteins, often participating in the assembly of large dynamic complexes. These domains bind to peptide ligands, which usually contain a core Pro-X-X-Pro (PXXP) sequence. Here we identify a class of SH3 domains that bind to ubiquitin. The yeast endocytic protein Sla1, as well as the mammalian proteins CIN85 and amphiphysin, carry ubiquitin-binding SH3 domains. Ubiquitin and peptide ligands bind to the same hydrophobic groove on the SH3 domain surface, and ubiquitin and a PXXP-containing protein fragment compete for binding to SH3 domains. We conclude that a subset of SH3 domains constitutes a distinct type of ubiquitin-binding domain and that ubiquitin binding can negatively regulate interaction of SH3 domains with canonical proline-rich ligands.     

Biochemical and computational analysis of LNX1 interacting proteins

PDZ (Post-synaptic density, 95 kDa, Discs large, Zona Occludens-1) domains are protein interaction domains that bind to the carboxy-terminal amino acids of binding partners, heterodimerize with other PDZ domains, and also bind phosphoinositides. PDZ domain containing proteins are frequently involved in the assembly of multi-protein complexes and clustering of transmembrane proteins. LNX1 (Ligand of Numb, protein X 1) is a RING (Really Interesting New Gene) domain-containing E3 ubiquitin ligase that also includes four PDZ domains suggesting it functions as a scaffold for a multi-protein complex. Here we use a human protein array to identify direct LNX1 PDZ domain binding partners. Screening of 8,000 human proteins with isolated PDZ domains identified 53 potential LNX1 binding partners. We combined this set with LNX1 interacting proteins identified by other methods to assemble a list of 220 LNX1 interacting proteins. Bioinformatic analysis of this protein list was used to select interactions of interest for future studies. Using this approach we identify and confirm six novel LNX1 binding partners: KCNA4, PAK6, PLEKHG5, PKC-alpha1, TYK2 and PBK, and suggest that LNX1 functions as a signalling scaffold.     

Proteogenomics and in silico structural and functional annotation of the barley powdery mildew Blumeria graminis f. sp. hordei

Blumeria graminis is an economically important obligate plant-pathogenic fungus, whose entire genome was recently sequenced and manually annotated using ab initio in silico predictions (Spanu et al. 2010, Science 330, 1543-1546). Employing large scale proteogenomic analysis we are now able to verify independently the existence of proteins predicted by ∼24% of open reading frame models. We compared the haustoria and sporulating hyphae proteomes and identified 71 proteins exclusively in haustoria, the feeding and effector-delivery organs of the pathogen. These proteins are significantly smaller than the rest of the protein pool and predicted to be secreted. Most do not share any similarities with Swiss–Prot or Trembl entries nor possess any identifiable Pfam domains. We used a novel automated prediction pipeline to model the 3D structures of the proteins, identify putative ligand binding sites and predict regions of intrinsic disorder. This revealed that the protein set found exclusively in haustoria is significantly less disordered than the rest of the identified Blumeria proteins or random (and representative) protein sets generated from the yeast proteome. For most of the haustorial proteins with unknown functions no good templates could be found, from which to generate high quality models. Thus, these unknown proteins present potentially new protein folds that can be specific to the interaction of the pathogen with its host.     

Genome-Wide Prediction and Validation of Peptides That Bind Human Prosurvival Bcl-2 Proteins

Programmed cell death is regulated by interactions between pro-apoptotic and prosurvival members of the Bcl-2 family. Pro-apoptotic family members contain a weakly conserved BH3 motif that can adopt an alpha-helical structure and bind to a groove on prosurvival partners Bcl-x_L, Bcl-w, Bcl-2, Mcl-1 and Bfl-1. Peptides corresponding to roughly 13 reported BH3 motifs have been verified to bind in this manner. Due to their short lengths and low sequence conservation, BH3 motifs are not detected using standard sequence-based bioinformatics approaches. Thus, it is possible that many additional proteins harbor BH3-like sequences that can mediate interactions with the Bcl-2 family. In this work, we used structure-based and data-based Bcl-2 interaction models to find new BH3-like peptides in the human proteome. We used peptide SPOT arrays to test candidate peptides for interaction with one or more of the prosurvival proteins Bcl-x_L, Bcl-w, Bcl-2, Mcl-1 and Bfl-1. For the 36 most promising array candidates, we quantified binding to all five human receptors using direct and competition binding assays in solution. All 36 peptides showed evidence of interaction with at least one prosurvival protein, and 22 peptides bound at least one prosurvival protein with a dissociation constant between 1 and 500 nM; many peptides had specificity profiles not previously observed. We also screened the full-length parent proteins of a subset of array-tested peptides for binding to Bcl-x_L and Mcl-1. Finally, we used the peptide binding data, in conjunction with previously reported interactions, to assess the affinity and specificity prediction performance of different models.     

Proteome scanning to predict PDZ domain interactions using support vector machines

Background  

PDZ domains mediate protein-protein interactions involved in important biological processes through the recognition of short linear motifs in their target proteins. Two recent independent studies have used protein microarray or phage display technology to detect PDZ domain interactions with peptide ligands on a large scale. Several computational predictors of PDZ domain interactions have been developed, however they are trained using only protein microarray data and focus on limited subsets of PDZ domains. An accurate predictor of genomic PDZ domain interactions would allow the proteomes of organisms to be scanned for potential binders. Such an application would require an accurate and precise predictor to avoid generating too many false positive hits given the large amount of possible interactors in a given proteome. Once validated these predictions will help to increase the coverage of current PDZ domain interaction networks and further our understanding of the roles that PDZ domains play in a variety of biological processes.     

Distinct binding specificity of the multiple PDZ domains of INADL, a human protein with homology to INAD from Drosophila melanogaster

PDZ domains are protein-protein interaction modules that typically bind to short peptide sequences at the carboxyl terminus of target proteins. Proteins containing multiple PDZ domains often bind to different trans-membrane and intracellular proteins, playing a central role as organizers of multimeric complexes. To characterize the rules underlying the binding specificity of different PDZ domains, we have assembled a novel repertoire of random peptides that are displayed at high density at the carboxyl terminus of the capsid D protein of bacteriophage lambda. We have exploited this combinatorial library to determine the peptide binding preference of the seven PDZ domains of human INADL, a multi-PDZ protein that is homologous to the INAD protein of Drosophila melanogaster. This approach has permitted the determination of the consensus ligand for each PDZ domain and the assignment to class I, class II, and to a new specificity class, class IV, characterized by the presence of an acidic residue at the carboxyl-terminal position. Homology modeling and site-directed mutagenesis experiments confirmed the involvement of specific residues at contact positions in determining the domain binding preference. However, these experiments failed to reveal simple rules that would permit the association of the chemical characteristics of any given residue in the peptide binding pocket to the preference for specific amino acid sequences in the ligand peptide. Rather, they suggested that to infer the binding preference of any PDZ domain, it is necessary to simultaneously take into account all contact positions by using computational procedures. For this purpose we extended the SPOT algorithm, originally developed for SH3 domains, to evaluate the probability that any peptide would bind to any given PDZ domain.     

EVH1 domains: structure, function and interactions

Drosophila enabled/vasodilator-stimulated phosphoprotein homology 1 (EVH1) domains are 115 residue protein-protein interaction modules which provide essential links for their host proteins to various signal transduction pathways. Many EVH1-containing proteins are associated closely with actin-based structures and are involved in re-organization of the actin cytoskeleton. EVH1 domains are also present in proteins enriched in neuronal tissue, thus implicating them as potential mediators of synaptic plasticity, linking them to memory formation and learning. Like Src homology 3, WW and GYF domains and profilin, EVH1 domains recognize and bind specific proline-rich sequences (PRSs). The binding is of low affinity, but tightly regulated by the high specificity encoded into residues in the protein:peptide interface. In general, a small (3-6 residue) 'core' PRS in the target protein binds a 'recognition pocket' on the domain surface. Further affinity- and specificity-increasing interactions are then formed between additional domain epitopes and peptide 'core-flanking' residues. The three-dimensional structures of EVH1:peptide complexes now reveal, in great detail, some of the most important features of these interactions and allow us to better understand the origins of specificity, ligand orientation and sequence degeneracy of target peptides, in low affinity signalling complexes.     

Systematic discovery of new recognition peptides mediating protein interaction networks

Many aspects of cell signalling, trafficking, and targeting are governed by interactions between globular protein domains and short peptide segments. These domains often bind multiple peptides that share a common sequence pattern, or “linear motif” (e.g., SH3 binding to PxxP). Many domains are known, though comparatively few linear motifs have been discovered. Their short length (three to eight residues), and the fact that they often reside in disordered regions in proteins makes them difficult to detect through sequence comparison or experiment. Nevertheless, each new motif provides critical molecular details of how interaction networks are constructed, and can explain how one protein is able to bind to very different partners. Here we show that binding motifs can be detected using data from genome-scale interaction studies, and thus avoid the normally slow discovery process. Our approach based on motif over-representation in non-homologous sequences, rediscovers known motifs and predicts dozens of others. Direct binding experiments reveal that two predicted motifs are indeed protein-binding modules: a DxxDxxxD protein phosphatase 1 binding motif with a K_D of 22 μM and a VxxxRxYS motif that binds Translin with a K_D of 43 μM. We estimate that there are dozens or even hundreds of linear motifs yet to be discovered that will give molecular insight into protein networks and greatly illuminate cellular processes.     

Interactions between peptides containing nucleobase amino acids and T7 phages displaying S. cerevisiae proteins

The importance of high-throughput analyses of protein abundances and functions is interestingly increasing in genomic/proteomic studies. In such postgenome sequencing era, a protein-detecting chip, in which a large number of molecules specifically capturing target proteins (capturing agents) such as antibodies, recombinant proteins, and small molecules are arrayed onto solid, wet, or semi-wet substrates, enables comprehensive analysis of proteomes by a single experiment. However, whole proteomes are generally complicated for comprehensive analyses so that alternative approaches to subproteome analysis categorized by protein functions and binding properties (focused proteome) would be effective. Approaching the goal of development of designed peptide chip for protein analysis, diversity increases in peptide structures and validation of target proteins are needed. We herein describe design and synthesis of nucleobase amino acid (NBA)-containing peptides, selection of nucleic acid-related proteins derived from S. cerevisiae, and detection of interactions between NBA-containing peptides and T7 phages displaying proteins by both enzyme-linked immunosorbent assays (ELISA) and label-free anomalous reflection of gold (AR) measurements. Twenty-eight phage clones were obtained by the phage-display method and sequenced. Ten of 28 clones were expected to be nucleic acid-related proteins including initiation factor, TYB protein, ribosomal proteins, elongation factor, ATP synthase subunit, GTP-binding protein, and ribonuclease. Other phage clones encoded several classes of enzymes such as reductase, oxidase, aldolase, metalloprotease, and hexokinase. Both ELISA and AR measurements suggested that the methodology of in vitro selection for recognition of the NBA-containing peptide presented in this study was successfully established. Such a combination of NBA and phage display technologies would be potential to efficiently confirm valuable target proteins binding specifically to capturing agents, to be arrayed onto solid surfaces to develop the designed peptide chip.     

The extracellular and cytoplasmic proteomes of the non-virulent Bacillus anthracis strain UM23C1-2

The recently published genome sequence of Bacillus anthracis Ames has facilitated the prediction of proteins associated with the virulence of this bacterium. The aim of this study was to define reference maps for the extracellular and cytoplasmic proteomes of the avirulent B. anthracis strain UM23C1-2 that are useful for physiological studies and the development of improved vaccines. Using 2-DE and subsequent MALDI-TOF-TOF MS, 64 proteins were identified in the extracellular proteome, only 29 of which were predicted to be exported into the culture medium. The latter included chitinases, proteases, nucleotidases, sulfatases, phosphatases and proteins of unknown function. Of the remaining proteins in the culture medium, 18 were predicted to be associated with the cell wall or anchored on the trans side of the cytoplasmic membrane while 17 other proteins lacked identifiable export signals and were predicted to be cytoplasmic proteins. Among the S-layer proteins, Sap and Eag account for 10% of the total extracellular proteome. Many of the proteins are predicted to contribute to the virulence and antigenic signature of B. anthracis. We have also studied the composition of the cytoplasmic proteome, identifying 300 distinct proteins. The most abundant cytoplasmic proteins are primarily those involved in glycolysis, amino acid metabolism, protein translation, protein folding and stress adaptation. The presence of a variety of proteases, peptidases, peptide binding proteins, as well as enzymes required for the metabolism of amino acids, suggests that B. anthracis is adapted to life in a protein-rich environment rather than the soil. We therefore speculate that proteases and peptidases could be useful targets for the development of improved vaccines. In addition, both of these B. anthracis compartment-specific proteomes can be used as reference maps to monitor changes in the production of secreted and cytosolic proteins that occur, for example, during growth in macrophages.     

The kinetics of PDZ domain-ligand interactions and implications for the binding mechanism

PDZ domains are protein adapter modules present in a few hundred human proteins. They play important roles in scaffolding and signal transduction. PDZ domains usually bind to the C termini of their target proteins. To assess the binding mechanism of this interaction we have performed the first in-solution kinetic study for PDZ domains and peptides corresponding to target ligands. Both PDZ3 from postsynaptic density protein 95 and PDZ2 from protein tyrosine phosphatase L1 bind their respective target peptides through an apparent A + B --> A.B mechanism without rate-limiting conformational changes. But a mutant with a fluorescent probe (Trp) outside of the binding pocket suggests that slight changes in the structure take place upon binding in protein tyrosine phosphatase-L1 PDZ2. For PDZ3 from postsynaptic density protein 95 the pH dependence of the binding reaction is consistent with a one-step mechanism with one titratable group. The salt dependence of the interaction shows that the formation of electrostatic interactions is rate-limiting for the association reaction but not for dissociation of the complex.     

HDAC6 and Ubp-M BUZ domains recognize specific C-terminal sequences of proteins

The BUZ/Znf-UBP domain is a protein module found in the cytoplasmic deacetylase HDAC6, E3 ubiquitin ligase BRAP2/IMP, and a subfamily of ubiquitin-specific proteases. Although several BUZ domains have been shown to bind ubiquitin with high affinity by recognizing its C-terminal sequence (RLRGG-COOH), it is currently unknown whether the interaction is sequence-specific or whether the BUZ domains are capable of binding to proteins other than ubiquitin. In this work, the BUZ domains of HDAC6 and Ubp-M were subjected to screening against a one-bead-one-compound (OBOC) peptide library that exhibited random peptide sequences with free C-termini. Sequence analysis of the selected binding peptides as well as alanine scanning studies revealed that the BUZ domains require a C-terminal Gly-Gly motif for binding. At the more N-terminal positions, the two BUZ domains have distinct sequence specificities, allowing them to bind to different peptides and/or proteins. A database search of the human proteome on the basis of the BUZ domain specificities identified 11 and 24 potential partner proteins for Ubp-M and HDAC6 BUZ domains, respectively. Peptides corresponding to the C-terminal sequences of four of the predicted binding partners (FBXO11, histone H4, PTOV1, and FAT10) were synthesized and tested for binding to the BUZ domains by fluorescence polarization. All four peptides bound to the HDAC6 BUZ domain with low micromolar K(D) values and less tightly to the Ubp-M BUZ domain. Finally, in vitro pull-down assays showed that the Ubp-M BUZ domain was capable of binding to the histone H3-histone H4 tetramer protein complex. Our results suggest that BUZ domains are sequence-specific protein-binding modules, with each BUZ domain potentially binding to a different subset of proteins.     

A motif-based profile scanning approach for genome-wide prediction of signaling pathways

The rapid increase in genomic information requires new techniques to infer protein function and predict protein-protein interactions. Bioinformatics identifies modular signaling domains within protein sequences with a high degree of accuracy. In contrast, little success has been achieved in predicting short linear sequence motifs within proteins targeted by these domains to form complex signaling networks. Here we describe a peptide library-based searching algorithm, accessible over the World Wide Web, that identifies sequence motifs likely to bind to specific protein domains such as 14-3-3, SH2, and SH3 domains, or likely to be phosphorylated by specific protein kinases such as Src and AKT. Predictions from database searches for proteins containing motifs matching two different domains in a common signaling pathway provides a much higher success rate. This technology facilitates prediction of cell signaling networks within proteomes, and could aid in the identification of drug targets for the treatment of human diseases.     

Structural properties and peptide ligand binding of the capsid homology domains of human Arc

The activity-regulated cytoskeleton-associated protein (Arc) is important for synaptic plasticity and the normal function of the brain. Arc interacts with neuronal postsynaptic proteins, but the mechanistic details of its function have not been fully established. The C-terminal domain of Arc consists of tandem domains, termed the N- and C-lobe. The N-lobe harbours a peptide binding site, able to bind multiple targets. By measuring the affinity of human Arc towards various peptides from stargazin and guanylate kinase-associated protein (GKAP), we have refined its specificity determinants. We found two sites in the GKAP repeat region that bind to Arc and confirmed these interactions by X-ray crystallography. Phosphorylation of the stargazin peptide did not affect binding affinity but caused changes in thermodynamic parameters. Comparison of the crystal structures of three high-resolution human Arc-peptide complexes identifies three conserved C–H…π interactions at the binding cavity, explaining the sequence specificity of short linear motif binding by Arc. We further characterise central residues of the Arc lobe fold, show the effects of peptide binding on protein dynamics, and identify acyl carrier proteins as structures similar to the Arc lobes. We hypothesise that Arc may affect protein-protein interactions and phase separation at the postsynaptic density, affecting protein turnover and re-modelling of the synapse. The present data on Arc structure and ligand binding will help in further deciphering these processes.     

Structure of the binding site for inositol phosphates in a PH domain.

Phosphatidylinositol bisphosphate has been found to bind specifically to pleckstrin homology (PH) domains that are commonly present in signalling proteins but also found in cytoskeleton. We have studied the complexes of the beta-spectrin PH domain and soluble inositol phosphates using both circular dichroism and nuclear magnetic resonance spectroscopy, and X-ray crystallography. The specific binding site is located in the centre of a positively charged surface patch of the domain. The presence of 4,5-bisphosphate group on the inositol ring is critical for binding. In the crystal structure that has been determined at 2.0 A resolution, inositol-1,4,5-trisphosphate is bound with salt bridges and hydrogen bonds through these phosphate groups whereas the 1-phosphate group is mostly solvent-exposed and the inositol ring has virtually no interactions with the protein. We propose a model in which PH domains are involved in reversible anchoring of proteins to membranes via their specific binding to phosphoinositides. They could also participate in a response to a second messenger such as inositol trisphosphate, organizing cross-roads in cellular signalling.     

Probability-based evaluation of peptide and protein identifications from tandem mass spectrometry and SEQUEST analysis: the human proteome

Large-scale protein identifications from highly complex protein mixtures have recently been achieved using multidimensional liquid chromatography coupled with tandem mass spectrometry (LC/LC-MS/MS) and subsequent database searching with algorithms such as SEQUEST. Here, we describe a probability-based evaluation of false positive rates associated with peptide identifications from three different human proteome samples. Peptides from human plasma, human mammary epithelial cell (HMEC) lysate, and human hepatocyte (Huh)-7.5 cell lysate were separated by strong cation exchange (SCX) chromatography coupled offline with reversed-phase capillary LC-MS/MS analyses. The MS/MS spectra were first analyzed by SEQUEST, searching independently against both normal and sequence-reversed human protein databases, and the false positive rates of peptide identifications for the three proteome samples were then analyzed and compared. The observed false positive rates of peptide identifications for human plasma were significantly higher than those for the human cell lines when identical filtering criteria were used, suggesting that the false positive rates are significantly dependent on sample characteristics, particularly the number of proteins found within the detectable dynamic range. Two new sets of filtering criteria are proposed for human plasma and human cell lines, respectively, to provide an overall confidence of >95% for peptide identifications. The new criteria were compared, using a normalized elution time (NET) criterion (Petritis et al. Anal. Chem. 2003, 75, 1039-1048), with previously published criteria (Washburn et al. Nat. Biotechnol. 2001, 19, 242-247). The results demonstrate that the present criteria provide significantly higher levels of confidence for peptide identifications from mammalian proteomes without greatly decreasing the number of identifications.