Proteomic screening method for phosphopeptide motif binding proteins using peptide libraries

Phosphopeptide binding domains mediate the directed and localized assembly of protein complexes essential to intracellular kinase signaling. To identify phosphopeptide binding proteins, we developed a proteomic screening method using immobilized partially degenerate phosphopeptide mixtures combined with SILAC and microcapillary LC-MS/MS. The method was used to identify proteins that specifically bound to phosphorylated peptide library affinity matrices, including pTyr, and the motifs pSer/pThr-Pro, pSer/pThr-X-X-X-pSer/pThr, pSer/pThr-Glu/Asp, or pSer/pThr-pSer/pThr in degenerate sequence contexts. Heavy and light SILAC lysates were applied to columns containing these phosphorylated and nonphosphorylated (control) peptide libraries respectively, and bound proteins were eluted, combined, digested, and analyzed by LC-MS/MS using a hybrid quadrupole-TOF mass spectrometer. Heavy/light peptide ion ratios were calculated, and peptides that yielded ratios greater than ～3:1 were considered as being from potential phosphopeptide binding proteins since this ratio represents the lowest ratio from a known positive control. Many of those identified were known phosphopeptide-binding proteins, including the SH2 domain containing p85 subunit of PI3K bound to pTyr, 14-3-3 bound to pSer/pThr-Asp/Glu, polo-box domain containing PLK1 and Pin1 bound to pSer/pThr-Pro, and pyruvate kinase M2 binding to pTyr. Approximately half of the hits identified by the peptide library screens were novel. Protein domain enrichment analysis revealed that most pTyr hits contain SH2 domains, as expected, and to a lesser extent SH3, C1, STAT, Tyr phosphatase, Pkinase, C2, and PH domains; however, pSer/pThr motifs did not reveal enriched domains across hits.     

Phosphopeptide binding specificities of BRCA1 COOH-terminal (BRCT) domains

Protein phosphorylation by protein kinases may generate docking sites for other proteins. It thus allows the assembly of signaling complexes in response to kinase activation. Several protein domains that bind phosphoserine or phosphothreonine residues have been identified, including the 14-3-3, PIN1, FHA, KIX, WD-40 domain, and polo box (Yaffe, M. B., and Elia, A. E. (2001) Curr. Opin. Cell Biol. 13, 131-138; Elia, A. E., Cantley, L. C., and Yaffe, M. B. (2003) Science 299, 1228-1231). The BRCA1 COOH-terminal (BRCT) domains are protein modules found in many proteins that regulate DNA damage responses (Koonin, E. V., Altschul, S. F., and Bork, P. (1996) Nat. Genet. 13, 266-268). Whether BRCT domains can mediate phosphorylation-dependent interactions has not been systematically investigated. We report here that the BRCT domains also recognize phosphopeptides. Oriented peptide library analysis indicated that the BRCT domains from BRCA1, MDC1, BARD1, and DNA Ligase IV preferred distinct phosphoserine-containing peptides. In addition, the interaction between BRCA1 and the BRCT binding motif of BACH1 was required for BACH1 checkpoint activity. Furthermore, BRCT domains of the yeast DNA repair protein Rad9 could bind phosphopeptides, suggesting that the BRCT domains represent a class of ancient phosphopeptide-binding modules. Potential targets of BRCT domains were identified through data base search. Structural analysis of BRCA1 BRCT repeats also predicted conserved residues that may form the phosphopeptide-binding pocket. Thus, the BRCT repeats are a new family of phosphopeptide-binding domains in DNA damage responses.     

Crystal Structure of the Pml1p Subunit of the Yeast Precursor mRNA Retention and Splicing Complex

The precursor mRNA retention and splicing (RES) complex mediates nuclear retention and enhances splicing of precursor mRNAs. The RES complex from yeast comprises three proteins, Snu17p, Bud13p and Pml1p. Snu17p acts as a central platform that concomitantly binds the Bud13p and Pml1p subunits via short peptide epitopes. As a step to decipher the molecular architecture of the RES complex, we have determined crystal structures of full-length Pml1p and N-terminally truncated Pml1p. The first 50 residues of full-length Pml1p, encompassing the Snu17p-binding region, are disordered, showing that Pml1p binds to Snu17p via an intrinsically unstructured region. The remainder of Pml1p folds as a forkhead-associated (FHA) domain, which is expanded by a number of noncanonical elements compared with known FHA domains from other proteins. An atypical N-terminal appendix runs across one β-sheet and thereby stabilizes the domain as shown by deletion experiments. FHA domains are thought to constitute phosphopeptide-binding elements. Consistently, a sulfate ion was found at the putative phosphopeptide-binding loops of full-length Pml1p. The N-terminally truncated version of the protein lacked a similar phosphopeptide mimic but retained an almost identical structure. A long loop neighboring the putative phosphopeptide-binding site was disordered in both structures. Comparison with other FHA domain proteins suggests that this loop adopts a defined conformation upon ligand binding and thereby confers ligand specificity. Our results show that in the RES complex, an FHA domain of Pml1p is flexibly tethered via an unstructured N-terminal region to Snu17p.     

Three-Dimensional Models of Immune Cell Surface Proteins and Identification of Binding Sites

The extracellular regions of many cell surface proteins of the immune system contain distinct domains that may be linked in many different ways and are often only loosely tethered to the transmembrane segment. In efforts to identify regions critical for binding, molecular models of these domains are used to select residues for mutagenesis and to map binding sites. Many immune cell surface proteins belong to protein superfamilies and display only limited sequence identity compared to proteins of known three-dimensional (3D) structure, often 30% or less. Therefore, detailed 3D structures are difficult to predict, and structure-based sequence analysis and model assessment are particularly important components of the model building process. In some cases, experimentally determined structures have made it possible to assess the accuracy of predictions, which illustrates the opportunities and shortcomings of the approach. Herein the model-based identification of binding sites in cell surface proteins is described and representative examples are discussed.     

14-3-3 proteins, FHA domains and BRCT domains in the DNA damage response

The DNA damage response depends on the concerted activity of protein serine/threonine kinases and modular phosphoserine/threonine-binding domains to relay the damage signal and recruit repair proteins. The PIKK family of protein kinases, which includes ATM/ATR/DNA-PK, preferentially phosphorylate Ser-Gln sites, while their basophilic downstream effecter kinases, Chk1/Chk2/MK2 preferentially phosphorylate hydrophobic-X-Arg-X-X-Ser/Thr-hydrophobic sites. A subset of tandem BRCT domains act as phosphopeptide binding modules that bind to ATM/ATR/DNA-PK substrates after DNA damage. Conversely, 14-3-3 proteins interact with substrates of Chk1/Chk2/MK2. FHA domains have been shown to interact with substrates of ATM/ATR/DNA-PK and CK2. In this review we consider how substrate phsophorylation together with BRCT domains, FHA domains and 14-3-3 proteins function to regulate ionizing radiation-induced nuclear foci and help to establish the G₂/M checkpoint. We discuss the role of MDC1 a molecular scaffold that recruits early proteins to foci, such as NBS1 and RNF8, through distinct phosphodependent interactions. In addition, we consider the role of 14-3-3 proteins and the Chk2 FHA domain in initiating and maintaining cell cycle arrest.     

Diversity in the SH2 domain family phosphotyrosyl peptide binding site

Src homology 2 (SH2) domains are approximately 100 residue phosphotyrosyl peptide binding modules found in signalling proteins and are important targets for therapeutic intervention. The peptide binding site is evolutionarily well conserved, particularly at the two major binding pockets, pTyr and pTyr + 3. We present a computational analysis of diversity within the peptide binding region and discuss molecular recognition beyond the conventional binding motif, drawing attention to novel conserved ligand interaction sites which may be exploitable in ligand binding studies. The peptide binding site is defined by selecting crystal contacts and domains are clustered according to binding site residue similarity. Comparison with a classification based on experimental peptide screening reveals a high level of qualitative agreement, indicating that the method is able independently to generate functional information. A conservation scoring method reveals extensive patches of conservation in some groups not present across the whole family, challenging the notion that the domains recognise only a linear phosphopeptide sequence. Conservation difference maps determine group-dependent clusters of conserved residues that are not seen when considering a larger experimentally determined group. Many of these residues contact the peptide outside the pTyr to pTyr + 3 motif, challenging the conventional view that this motif is largely responsible for ligand recognition and discrimination.     

Microfluidic,bead-based assay: Theory and experiments

Microbeads are frequently used as a solid support for biomolecules such as proteins and nucleic acids in heterogeneous microfluidic assays. However, relatively few studies investigate the binding kinetics on modified bead surfaces in a microfluidics context. In this study, a customized hot embossing technique is used to stamp microwells in a thin plastic substrate where streptavidin-coated agarose beads are selectively placed and subsequently immobilized within a conduit. Biotinylated quantum dots are used as a label to monitor target analyte binding to the bead's surface. Three-dimensional finite element simulations are carried out to model the binding kinetics on the bead's surface. The model accounts for surface exclusion effects resulting from a single quantum dot occluding multiple receptor sites. The theoretical predictions are compared and favorably agree with experimental observations. The theoretical simulations provide a useful tool to predict how varying parameters affect microbead reaction kinetics and sensor performance. This study enhances our understanding of bead-based microfluidic assays and provides a design tool for developers of point-of-care, lab-on-chip devices for medical diagnosis, food and water quality inspection, and environmental monitoring.     

SH2 domains from suppressor of cytokine signaling-3 and protein tyrosine phosphatase SHP-2 have similar binding specificities

Suppressor of cytokine signaling-3 (SOCS-3) and the protein tyrosine phosphatase SHP-2 both regulate signaling by cytokines of the interleukin-6 family, and this is dependent upon recruitment to tyrosine 757 in the shared cytokine receptor subunit gp130. To better explore the overlap in ligand binding specificities exhibited by these two signaling regulators, we have mapped the phosphopeptide binding preferences of the SH2 domains from SOCS-3 and SHP-2. Degenerate phosphopeptide libraries were screened against recombinantly produced SH2 domains to determine the sequences of optimal phosphopeptide ligands. We found that the consensus ligand binding motif for SOCS-3 was pY-(S/A/V/Y/F)-hydrophobic-(V/I/L)-hydrophobic-(H/V/I/Y), while the consensus motif for SHP-2 was pY-(S/T/A/V/I)-X-(V/I/L)-X-(W/F). We validated these data through the design of phosphopeptide ligands based on the consensus motifs and found that these bound to SOCS-3 and SHP-2 with high affinity. Finally, we have compared the affinity of SOCS-3 for binding to phosphopeptides representing putative docking sites in the gp130, leptin and erythropoietin receptors. While SOCS-3 binds with much higher affinity to a gp130 phosphopeptide than to phosphopeptides derived from the other receptors, multiple SOCS-3 binding sites are predicted to exist in the leptin and erythropoietin receptors which may compensate for weaker binding to individual sites.     

Differentiation of peptide molecular recognition by phospholipase C gamma-1 Src homology-2 domain and a mutant Tyr phosphatase PTP1bC215S.

Activated epidermal growth factor receptor (EGFR) undergoes autophosphorylation on several cytoplasmic tyrosine residues, which may then associate with the src homology-2 (SH2) domains of effector proteins such as phospholipase C gamma-1 (PLC gamma-1). Specific phosphotyrosine (pTyr)-modified EGFR fragment peptides can inhibit this intermolecular binding between activated EGFR and a tandem amino- and carboxy-terminal (N/C) SH2 protein construct derived from PLC gamma-1. In this study, we further explored the molecular recognition of phosphorylated EGFR988-998 (Asp-Ala-Asp-Glu-pTyr-Leu-Ile-Pro-Gln-Gln-Gly, I) by PLC gamma-1 N/C SH2 in terms of singular Ala substitutions for amino acid residues N- and C-terminal to the pTyr (P site) of phosphopeptide I. Comparison of the extent to which these phosphopeptides inhibited binding of PLC gamma-1 N/C SH2 to activated EGFR showed the critical importance of amino acid side chains at positions P+2 (Ile994), P+3 (Pro995), and P+4 (Gln996). Relative to phosphopeptide I, multiple Ala substitution throughout the N-terminal sequence, N-terminal sequence, N-terminal truncation, or dephosphorylation of pTyr each resulted in significantly decreased binding to PLC gamma-1 N/C SH2. These structure-activity results were analyzed by molecular modeling studies of the predicted binding of phosphopeptide I to each the N- and C-terminal SH2 domains of PLC gamma-1.(ABSTRACT TRUNCATED AT 250 WORDS)     

Structural insights into the innate immune recognition specificities of L- and H-ficolins

Innate immunity relies critically upon the ability of a few pattern recognition molecules to sense molecular markers on pathogens, but little is known about these interactions at the atomic level. Human L- and H-ficolins are soluble oligomeric defence proteins with lectin-like activity, assembled from collagen fibers prolonged by fibrinogen-like recognition domains. The X-ray structures of their trimeric recognition domains, alone and in complex with various ligands, have been solved to resolutions up to 1.95 and 1.7 A, respectively. Both domains have three-lobed structures with clefts separating the distal parts of the protomers. Ca(2+) ions are found at sites homologous to those described for tachylectin 5A (TL5A), an invertebrate lectin. Outer binding sites (S1) homologous to the GlcNAc-binding pocket of TL5A are present in the ficolins but show different structures and specificities. In L-ficolin, three additional binding sites (S2-S4) surround the cleft. Together, they define an unpredicted continuous recognition surface able to sense various acetylated and neutral carbohydrate markers in the context of extended polysaccharides such as 1,3-beta-D-glucan, as found on microbial or apoptotic surfaces.     

Common EF-hand motifs in cholinesterases and neuroligins suggest a role for Ca2+ binding in cell surface associations

Comparisons of protein sequence via cyclic training of Hidden Markov Models (HMMs) in conjunction with alignments of three-dimensional structure, using the Combinatorial Extension (CE) algorithm, reveal two putative EF-hand metal binding domains in acetylcholinesterase. Based on sequence similarity, putative EF-hands are also predicted for the neuroligin family of cell surface proteins. These predictions are supported by experimental evidence. In the acetylcholinesterase crystal structure from Torpedo californica, the first putative EF-hand region binds the Zn2+ found in the heavy metal replacement structure. Further, the interaction of neuroligin 1 with its cognate receptor neurexin depends on Ca2+. Thus, members of the alpha,beta hydrolase fold family of proteins contain potential Ca2+ binding sites, which in some family members may be critical for heterologous cell associations.     

Targeted Disruption of an EH‐domain Protein Endocytic Complex,Pan1‐End3

Pan1 is a multi‐domain scaffold that enables dynamic interactions with both structural and regulatory components of the endocytic pathway. Pan1 is composed of Eps15 Homology (EH) domains which interact with adaptor proteins, a central region that is responsible for its oligomerization and C‐terminal binding sites for Arp2/3, F‐actin, and type‐I myosin motors. In this study, we have characterized the binding sites between Pan1 and its constitutive binding partner End3, another EH domain containing endocytic protein. The C‐terminal End3 Repeats of End3 associate with the N‐terminal part of Pan1's central coiled‐coil region. These repeats appear to act independently of one another as tandem, redundant binding sites for Pan1. The end3‐1 allele was sequenced, and corresponds to a C‐terminal truncation lacking the End3 Repeats. Mutations of the End3 Repeats highlight that those residues which are identical between these repeats serve as contact sites for the interaction with Pan1.      

Caged phosphopeptides reveal a temporal role for 14-3-3 in G1 arrest and S-phase checkpoint function

Using classical genetics to study modular phosphopeptide-binding domains within a family of proteins that are functionally redundant is difficult when other members of the domain family compensate for the product of the knocked-out gene. Here we describe a chemical genetics approach that overcomes this limitation by using UV-activatable caged phosphopeptides. By incorporating a caged phosphoserine residue within a consensus motif, these reagents simultaneously and synchronously inactivate all phosphoserine/phosphothreonine-binding domain family members in a rapid and temporally regulated manner. We applied this approach to study the global function of 14-3-3 proteins in cell cycle control. Activation of the caged phosphopeptides by UV irradiation displaced endogenous proteins from 14-3-3-binding, causing premature cell cycle entry, release of G1 cells from interphase arrest and loss of the S-phase checkpoint after DNA damage, accompanied by high levels of cell death. This class of reagents will greatly facilitate molecular dissection of kinase-dependent signaling pathways when applied to other phosphopeptide-binding domains including SH2, Polo-box and tandem BRCT domains.     

In vitro evolution of recognition specificity mediated by SH3 domains reveals target recognition rules

We have designed a repertoire of 10(7) different SH3 domains by grafting the residues that are represented in the binding surfaces of natural SH3 domains onto the scaffold of the human Abl-SH3 domain. This phage-displayed library was screened by affinity selection for SH3 domains that bind to the synthetic peptides, APTYPPPLPP and LSSRPLPTLPSP, which are peptide ligands for the human Abl or Src SH3 domains, respectively. By characterizing the isolates, we have observed that as few as two or three amino acid substitutions lead to dramatic changes in recognition specificity. We propose that the ability to shift recognition specificity with a small number of amino acid replacements is an important evolutionary characteristic of protein binding modules. Furthermore, we have used the information obtained by these in vitro evolution experiments to generate a scoring matrix that evaluates the probability that any SH3 domain binds to the peptide ligands for the Abl and Src SH3 domains. A table of predictions for the 28 SH3 domains of baker's yeast is presented.     

Prediction of carbohydrate binding sites on protein surfaces with 3-dimensional probability density distributions of interacting atoms

Non-covalent protein-carbohydrate interactions mediate molecular targeting in many biological processes. Prediction of non-covalent carbohydrate binding sites on protein surfaces not only provides insights into the functions of the query proteins; information on key carbohydrate-binding residues could suggest site-directed mutagenesis experiments, design therapeutics targeting carbohydrate-binding proteins, and provide guidance in engineering protein-carbohydrate interactions. In this work, we show that non-covalent carbohydrate binding sites on protein surfaces can be predicted with relatively high accuracy when the query protein structures are known. The prediction capabilities were based on a novel encoding scheme of the three-dimensional probability density maps describing the distributions of 36 non-covalent interacting atom types around protein surfaces. One machine learning model was trained for each of the 30 protein atom types. The machine learning algorithms predicted tentative carbohydrate binding sites on query proteins by recognizing the characteristic interacting atom distribution patterns specific for carbohydrate binding sites from known protein structures. The prediction results for all protein atom types were integrated into surface patches as tentative carbohydrate binding sites based on normalized prediction confidence level. The prediction capabilities of the predictors were benchmarked by a 10-fold cross validation on 497 non-redundant proteins with known carbohydrate binding sites. The predictors were further tested on an independent test set with 108 proteins. The residue-based Matthews correlation coefficient (MCC) for the independent test was 0.45, with prediction precision and sensitivity (or recall) of 0.45 and 0.49 respectively. In addition, 111 unbound carbohydrate-binding protein structures for which the structures were determined in the absence of the carbohydrate ligands were predicted with the trained predictors. The overall prediction MCC was 0.49. Independent tests on anti-carbohydrate antibodies showed that the carbohydrate antigen binding sites were predicted with comparable accuracy. These results demonstrate that the predictors are among the best in carbohydrate binding site predictions to date.     

Mutational analysis reveals a single binding interface between RhoA and its effector, PRK1

Protein kinase C-related kinases (PRKs) are serine/threonine kinases that are members of the protein kinase C superfamily and can be activated by binding to members of the Rho family of small G proteins via a Rho binding motif known as an HR1 domain. The PRKs contain three tandem HR1 domains at their N-termini. The structure of the HR1a domain from PRK1 in complex with RhoA [Maesaki, R., et al. (1999) Mol. Cell 4, 793-803] identified two potential contact interfaces between the G protein and the HR1a domain. In this work, we have used an alanine scanning mutagenesis approach to identify whether both contact sites are used when the two proteins interact in solution and also whether HR1b, the second HR1 domain from PRK1, plays a role in binding to RhoA. The mutagenesis identified just one contact site as being relevant for binding of RhoA and HR1a in solution, and the HR1b domain was found not to contribute to RhoA binding. The folded state and thermal stability of the HR1a and HR1b domains were also investigated. HR1b was found to be more thermally stable than HR1a, and it is hypothesized that the differences in the biophysical properties of these two domains govern their interaction with small G proteins.     

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.     

The BARD1 C-terminal domain structure and interactions with polyadenylation factor CstF-50

The BARD1 N-terminal RING domain binds BRCA1 while the BARD1 C-terminal ankyrin and tandem BRCT repeat domains bind CstF-50 to modulate mRNA processing and RNAP II stability in response to DNA damage. Here we characterize the BARD1 structural biochemistry responsible for CstF-50 binding. The crystal structure of the BARD1 BRCT domain uncovers a degenerate phosphopeptide binding pocket lacking the key arginine required for phosphopeptide interactions in other BRCT proteins. Small angle X-ray scattering together with limited proteolysis results indicates that ankyrin and BRCT domains are linked by a flexible tether and do not adopt a fixed orientation relative to one another. Protein pull-down experiments utilizing a series of purified BARD1 deletion mutants indicate that interactions between the CstF-50 WD-40 domain and BARD1 involve the ankyrin-BRCT linker but do not require ankyrin or BRCT domains. The structural plasticity imparted by the ANK-BRCT linker helps to explain the regulated assembly of different protein BARD1 complexes with distinct functions in DNA damage signaling including BARD1-dependent induction of apoptosis plus p53 stabilization and interactions. BARD1 architecture and plasticity imparted by the ANK-BRCT linker are suitable to allow the BARD1 C-terminus to act as a hub with multiple binding sites to integrate diverse DNA damage signals directly to RNA polymerase.