PPSP: prediction of PK-specific phosphorylation site with Bayesian decision theory
© Xue et al; licensee BioMed Central Ltd. 2006
Received: 21 August 2005
Accepted: 20 March 2006
Published: 20 March 2006
As a reversible and dynamic post-translational modification (PTM) of proteins, phosphorylation plays essential regulatory roles in a broad spectrum of the biological processes. Although many studies have been contributed on the molecular mechanism of phosphorylation dynamics, the intrinsic feature of substrates specificity is still elusive and remains to be delineated.
In this work, we present a novel, versatile and comprehensive program, PPSP (Prediction of PK-specific Phosphorylation site), deployed with approach of Bayesian decision theory (BDT). PPSP could predict the potential phosphorylation sites accurately for ~70 PK (Protein Kinase) groups. Compared with four existing tools Scansite, NetPhosK, KinasePhos and GPS, PPSP is more accurate and powerful than these tools. Moreover, PPSP also provides the prediction for many novel PKs, say, TRK, mTOR, SyK and MET/RON, etc. The accuracy of these novel PKs are also satisfying.
Taken together, we propose that PPSP could be a potentially powerful tool for the experimentalists who are focusing on phosphorylation substrates with their PK-specific sites identification. Moreover, the BDT strategy could also be a ubiquitous approach for PTMs, such as sumoylation and ubiquitination, etc.
Protein phosphorylation, as one of the most common post-translational modifications (PTM), is reversibly and transiently performed by protein kinases (PKs). It plays crucial regulatory roles in a variety of biological cellular processes, including transcription , translation , mitosis/cell cycle , neurite outgrowth [4, 5] and signal transductions , etc. Many previous researches have contributed to increase our knowledge on phosphorylation. However, the intrinsic features of phosphorylation dynamics are still cryptic and remain to be dissected. Biochemically, the catalytic site of a PK hydrolyzes adenosine triphosphate (ATP) and transfers a phosphate moiety to the acceptor residue (S/T, Y in eukaryotes) in the substrate. Each PK only modifies a defined subset of substrates specifically to ensure signaling fidelity, and defects of PK function often induce a variety of diseases and cancers .
There is an extensively-adopted hypothesis that PKs phosphorylate their substrates at the specific sites (consensus sequence) flanking with canonical motif [8–10]. To date, the consensus motifs of ~30 PKs have been reported . However, there is still a large number of PKs with their specific target motifs remained to be identified. Therefore, elucidating PK-specific phosphorylation sites on the substrates is the foundation of understanding the molecular mechanism of substrates specificity and important for the biomedical drug design. However, it has been described that only consensus motif is not enough for providing the specificity of PK recognition in vivo . There are numerous mechanisms have been proposed to contribute specificity for PKs, such as co-complex of PKs with their substrates, subcellular co-localization, interacting through modular docking sites, phosphopeptide-binding mechanisms, etc [12–17]. In a cell, protein kinase usually forms a tight complex with its target either through a third scaffold protein, or by recognizing and binding short sequence of the substrate, known as a docking site [12, 18]. Moreover, phosphopeptide-binding domains (PBDs) are also important to achieve substrate specificity. Numerous PBDs (PTB, WW, SH2, SH3, FHA, MH2, WD40, Polo-box, and 14-3-3, etc) bind the phosphorylated forms of specific proteins, with recognizing distinct peptides surrounding the phosphorylated sites (pS/T, or pY) [14–17, 19]. However, how these mechanisms achieve the additional specificity for PKs beyond phosphorylated motifs is still elusive, and there are very few computational studies published on this area [13, 16, 19]. In addition, many docking sites and PBDs still remain to be dissected. Thus, in this work, we focus on the prediction of PK-specific phosphorylation sites based on profiles/features of the surrounding primary sequences, as previously described [8–10].
Conventional experimental identifications of PK-specific phosphorylation sites on substrates in vivo and in vitro have provided the foundation of understanding the mechanisms of phosphorylation dynamics. However, these experiments are often time-consuming and expensive. And the enzymatic activity of the PKs are usually diminished or impeded in vitro, hampering on the studies of phosphorylation greatly. Recently, phospho-proteomic studies with mass spectrometry (MS) approaches have generated numerous data in yeast , mouse , and human , etc. But in these cases, it's still difficult to distinguish the PK-specific sites on the substrates. With regard of this, it is of note that the in silico prediction of PK-specific phosphorylation sites is in urgent need for the further experimental manipulation. To address this question, several excellent predictors have been implemented and reported [13, 22–25]. For example, NetPhos has employed the consensus-motif-based approaches implemented in the artificial neural networks (ANNs) algorithm . The enhanced version, NetPhosK can predict PK-specific phosphorylation sites for ~17 PKs . Another online tool Scansite  has constructed the motif profiles of phosphorylation sites for ~20 PKs, and could predict their target sites, respectively. Previously, we have reported a web server GPS, which has been implemented in GPS (Group-based phosphorylation Predicting and Scoring) algorithm [26, 27]. GPS could predict ~70 kinds of PK-specific phosphorylation sites, and gain excellent performance on several PK groups, especially for kinase Aurora-B. Recently, a novel and excellent web tool of KinasePhos has been incorporated with HMM (Hidden Markov Models) algorithm and constructed for phosphorylation sites predicting of 18 PK-specific groups [24, 25].
In this study, we present a novel, convenient and comprehensive program, PPSP (Prediction of PK-specific Phosphorylation site), implemented in an algorithm of Bayesian decision theory (BDT). An online PPSP web service has been also constructed, accurately predicting PK-specific phosphorylation sites for 68 PK groups. The prediction performances of PPSP are satisfactory on several well-studied PKs and comparable with the other existing tools NetPhosK, Scansite, KinasePhos and GPS. Moreover, PPSP also provides the accurate prediction for many novel PKs, such as TRK, mTOR, SyK, and MET/RON, etc. Obviously, PPSP is more accurate and powerful. Therefore, we propose that PPSP could be useful and insightful for further experimental design. In addition, the prediction results of PPSP combined with delicate experiments verifications will propel our understanding of the mechanisms of phosphorylation into a new phase.
Preparation of training data set
Firstly, we obtained the data set of phosphorylation sites from Phospho.ELM (Ver 2.0, Sep. 2004)  and filtered the phosphorylation sites without information of PKs. There were ~1,400 sites preserved. We also manually curated the recent literature and acquired ~660 items (Before Nov. 2004). These newly curated data has been submitted to Phospho.ELM for further integration. The two data sets were integrated, and the redundant items were removed if two items exactly pinned point to the same phosphorylation site from one protein sequence. Then the total training data set contained >2,000 non-redundant positive data with very few homologous sites (see additional file 1).
Since there were several PKs with too few known phosphorylation sites, we clustered them into distinct sub-groups based on sequence homology. For example, eight ribosomal protein S6 kinases (RSK1, Q15418; RSK3, Q15349; RSK2, P51812; MSK1, O75676; MSK1, O75582; RSK4, Q9UK32; S6K1, P23443; STK14B, Q9UBS0) are homologous with high similarity, so we clustered these PKs into a unique PK group of S6K (Ribosomal protein S6 kinase, or RSK). In total, we have enabled 68 PK grouped.
Although Swiss-Prot also curates a huge amount of phosphorylation sites, we have found ~69% of the annotation to be ambiguous (7,924 of 11,520) (see additional file 2). There are only 842 items to be kinase-specific sites, and only 18 PKs with not less than ten sites (see additional file 3). Phospho.ELM has been constructed based on the rationale of allowing both experimentalists and bioinformatists to easily access extensive information of phosphoproteins with their sites, i.e., tracking the primary reference to find whether the site is really phosphorylated, identified in vivo or in vitro, and the relationship between the phosphorylation with physiological response . And these data has been collected from literature manually with high quality. Taken together, although other resources also have collected some phosphorylation sites, we chose Phospho.ELM for its comprehensiveness.
Positive & negative control for evaluation
The sequence information of these phosphorylation substrates was retrieved from ExPASy. As previously described , we adopted the experimental phosphorylation sites as the positive control, while all other residues (S/T or Y) in the phosphorylation substrates were regarded as the negative control. The detailed statistics of the positive and negative data sets categorized by PK groups is available (see additional file 4).
Bayesian Decision Theory (BDT)
Supposed that we have an unclassified data x that belongs to one of two certain categories: C1 (defined as phosphorylated sites in this work) and C2 (defined as non-phosphorylated sites). In addition, suppose the posterior probability of x for these two categories can be denoted as: p(C1|x) and p(C2|x). Then the probability of wrong prediction is:
To minimize the expectation of error probability that is defined as :
P(error) = ∫P(error|x)p(x)dx (2)
It is obvious that one should choose the more probable category as the prediction result, which can be formulated by the Bayesian Decision Rule :
Furthermore, by definition we can introduce the loss function λ(α i |C j ), where α i ,i = 1,2 is the finite set of possible solution. Thus the expected loss (risk) of taking action α i is:
In this condition, the goal of optimization becomes to minimize the overall risk for every x. Similar to the rationale of Bayesian Decision Rule, we can obtain the best performance by computing R(α i |x) for each solution α i and choose that for which has the minimal overall risk (also named as Bayes Risk) .
Training and prediction procedure
In this study, a local ennea-peptide (9aa) is deployed to define a candidate phosphorylation site, which has 4 upstream and 4 downstream residues of the potential phosphorylation site and can be denoted as = (x1,x2,...,x9)'. Given some positive (training) data, there are many ways to estimate R(α i |x) (where α1 and α2 denote different prediction results: true and false phosphorylation sites, respectively). One simple way is to assume that all flanking residues are mutual independent, and then the Bayes Risk can be formulated as:
Here p(C l |x j ) is the posterior probability of x j belonging to category C l and can be further described by the Bayesian formula:
Here p(C l ) is the prior probability of category C l and p(x j |C l ) can be estimated by observing the occurrence of each residue in training data given the hypothesis of equation (5). Although there are much more false phosphorylation site in data set, we give equal prior probability for each category (no prior information), which can avoid bias prediction results. The loss function we construct is based on BLOSUM62 matrix  by considering the biochemical difference of residues, which can be formulated as:
Although other matrices could be also adopted, the BLOSUM62 matrix is chosen in this work. Moreover, we introduce a trade-off threshold b as the only parameter in this method to control the performance for different categories. Thus the final Discriminant function for prediction is:
Results and discussions
Prediction performance of PPSP
Three measurements, i.e., Sensitivity (Sn), Specificity (Sp), and Mathew correlation coefficient (MCC) are widely employed to evaluate the performance of prediction with definitions as below:
Among the data with positive predictions by PPSP, the real positives are regarded as true positives (TP), while the others are defined as false positives (FP). Among the data with negative predictions by PPSP, the real positives are regarded as false negatives (FN), while the others are defined as true negatives (TN). The Sensitivity (Sn) and Specificity (Sp) illustrate the correct prediction ratios of positive and negative data sets respectively. But when the number of positive data and negative data differ too much from each other, the Mathew correlation coefficient (MCC) should be calculated to assess the prediction performance. The value of MCC ranges from -1 to 1, and bigger MCC stands for better prediction performance.
To assess whether PPSP is unbiased and robust for prediction, we adopt the standard method "Jack-Knife" validation. We perform the Jack-Knife validation for these PKs by removing one real site from the training data set at a time and re-calculating the Sn &Sp, respectively. The final results are the average of the all Sn &Sp of the Jack-Knife validation. Although "Jack-Knife" validation does make sense when the size of the data set is small (i.e., N < 30), we have also taken an additional test with n-fold (4-, 6-, 8- and 10-fold in this work) cross-validation for 22 PK groups with larger positive data set (N ≥ 30). As previously proposed , the tests are repeated for 20 times and the Sn &Sp is computed each time. Then the average Sn & Sp are calculated as the final results (see additional file 5).
The performances of self-consistency, Jack-knife validation and n -fold (4-, 6-, 8-, 10-fold in this work) cross-validation for four well-studied PKs of PKA, CK2, ATM and S6K. The n- fold cross-validation has been performed for the large data sets (N ≥ 30).
n -fold cross-validation
Data set (No.)
The self-consistency performance and Jack-knife validation for four novel PKs of TRK, mTOR, SyK and MET/RON.
With the default cut-off of PPSP, the percentile of the sites predicted to be potential true positive hits is listed. Both random ennea-peptides and data sets from human proteome have been computed, separately.
Random ennea- peptides
Ennea-peptides from human proteome
Comparison of PPSP with Scansite, NetPhosK, KinasePhos and GPS
The prediction performance of Scansite, NetPhosK, KinasePhos and GPS for four well-studied PKs of PKA, CK2, ATM and S6K.
The experimental verified vs. predicted CK2-specific phosphorylation sites of Bluetongue virus (BTV) nonstructural protein 2 (NS2), Drosophila transcription factor GAGA and human Calmodulin protein.
T79, S81, S101, T117
T247, S249, S259
T123, S335, S337, S380, S388
T5, T79, S81, T117
T87, T88, S204, T247, S249, S259, T266
T5, T28, T44, T62, T79,
S81, S101, T117
T88, S182, T247, S249, S259
T79, S81, T117
T247, S249, S259
S240, S241, S339, T378, S386, S389,
S391, S393, S397, S518, S521, S523
T5, T44, T79, S81, S101, T117
S240, S339, T378, S386,
S391, S397, S518, S521
T5, T44, T79, S81,
T5, T79, S81, S101, T117
T123, S337, S339, S385, S386,
S388, S518, S521
T5, T44, T79, S81, T117
Application of PPSP to the novel PKs
For application of PPSP to the novel PKs, here we employ PPSP to predict the phosphorylation sites of TRK. TRK is a sub-family of receptor tyrosine kinases (RTK), consisting three highly similar homologs, TRK-A, -B, and -C . TRK-A, -B, and -C could be activated specifically by nerve growth factor (NGF), brain-derived neurotrophic factor (BDNF) and NT-4/-5, and NT-3, respectively. Under activated state, TRK could regulate a variety of biological processes including cell survival, embryo, differentiation, proliferation, axon and dendrite growth and patterning, and apoptosis, etc .
In this work, we present a novel computational program–PPSP (prediction of PK-specific phosphorylation sites) based on Bayesian decision theory (BDT). We model a candidate phosphorylation motif as an ennea/nona-peptide (9aa) flanking with 4 upstream and 4 downstream residues of a potential phosphorylation site (S/T, or Y). With the BDT algorithm, we estimate the probability distributions of true and false phosphorylation sites and make prediction based on a loss function constructed with BLOSUM62 matrix . We have evaluated the sensitivity and specificity of PPSP by "Jack-knife" validation. An online PPSP web service has been also constructed, accurately predicting PK-specific phosphorylation sites for 68 PK groups. For comparison with four reported systems Scansite, NetPhosK, KinasePhos and GPS, we take four well-studied PKs of PKA, CK2, ATM and S6K as model kinases. The prediction performances of PPSP are satisfactory judged using these well-studied PKs and comparable with the other existing tools. Moreover, PPSP also provides the accurate prediction for many novel PKs, such as TRK, mTOR, SyK, and MET/RON, etc. Thus, comparison with the previous work, PPSP provides more accurate and powerful ability. Moreover, the BDT approach could also be an extensive method for PTMs prediction, such as sumoylation and ubiquitination, etc. In addition, although many phospho-proteomic researches have generated numerous data [8, 20, 21], however, the up-regulated PKs still remain to be dissected. Despite the demonstration of phosphor-regulation of protein kinases and their respective substrates, the exact phosphorylation sites are unclear [4, 5]. Taken together, the prediction results of PPSP should be insightful and important for further experiments. The combination of computational and experimental identifications will propel our understanding of phosphorylation dynamics into a new phase.
Availability and requirements
PPSP has been implemented in Linux + Apache + PHP, and is freely available at: http://bioinformatics.lcd-ustc.org/PPSP. A latest web browser (eg. Internet Explorer, Netscape, or Mozilla, etc) is required.
We thank Dr. Fengfeng Zhou (UGA) for helpful discussions and critical reading of this manuscript. The authors thank Drs T.J. Gibson and F. Diella for providing the data set of Phospho.ELM for this work. We thank the anonymous referees for their many helpful comments. The work is supported by grants from Chinese Natural Science Foundation (39925018 and 30121001), Chinese Academy of Science (KSCX2-2-01), Chinese 973 project (2002CB713700), Beijing Office for Science (H020220020220) and American Cancer Society (RPG-99-173-01) to X. Yao. X. Yao is a GCC Distinguished Cancer Research Scholar.
- Schafmeier T, Haase A, Kaldi K, Scholz J, Fuchs M, Brunner M: Transcriptional feedback of neurospora circadian clock gene by phosphorylation-dependent inactivation of its transcription factor. Cell 2005, 122(2):235–246. 10.1016/j.cell.2005.05.032View ArticlePubMedGoogle Scholar
- Singh CR, Curtis C, Yamamoto Y, Hall NS, Kruse DS, He H, Hannig EM, Asano K: Eukaryotic translation initiation factor 5 is critical for integrity of the scanning preinitiation complex and accurate control of GCN4 translation. Mol Cell Biol 2005, 25(13):5480–5491. 10.1128/MCB.25.13.5480-5491.2005PubMed CentralView ArticlePubMedGoogle Scholar
- Lou Y, Yao J, Zereshki A, Dou Z, Ahmed K, Wang H, Hu J, Wang Y, Yao X: NEK2A interacts with MAD1 and possibly functions as a novel integrator of the spindle checkpoint signaling. J Biol Chem 2004, 279(19):20049–20057. 10.1074/jbc.M314205200View ArticlePubMedGoogle Scholar
- Liu HY, MacDonald JI, Hryciw T, Li C, Meakin SO: Human tumorous imaginal disc 1 (TID1) associates with Trk receptor tyrosine kinases and regulates neurite outgrowth in nnr5-TrkA cells. J Biol Chem 2005, 280(20):19461–19471. 10.1074/jbc.M500313200View ArticlePubMedGoogle Scholar
- Robinson KN, Manto K, Buchsbaum RJ, MacDonald JI, Meakin SO: Neurotrophin-dependent tyrosine phosphorylation of Ras guanine-releasing factor 1 and associated neurite outgrowth is dependent on the HIKE domain of TrkA. J Biol Chem 2005, 280(1):225–235. 10.1074/jbc.M505720200View ArticlePubMedGoogle Scholar
- Pawson T: Specificity in signal transduction: from phosphotyrosine-SH2 domain interactions to complex cellular systems. Cell 2004, 116(2):191–203. 10.1016/S0092-8674(03)01077-8View ArticlePubMedGoogle Scholar
- Ma L, Chen Z, Erdjument-Bromage H, Tempst P, Pandolfi PP: Phosphorylation and functional inactivation of TSC2 by Erk implications for tuberous sclerosis and cancer pathogenesis. Cell 2005, 121(2):179–193. 10.1016/j.cell.2005.02.031View ArticlePubMedGoogle Scholar
- Beausoleil SA, Jedrychowski M, Schwartz D, Elias JE, Villen J, Li J, Cohn MA, Cantley LC, Gygi SP: Large-scale characterization of HeLa cell nuclear phosphoproteins. Proc Natl Acad Sci U S A 2004, 101(33):12130–12135. 10.1073/pnas.0404720101PubMed CentralView ArticlePubMedGoogle Scholar
- Kreegipuu A, Blom N, Brunak S: PhosphoBase, a database of phosphorylation sites: release 2.0. Nucleic Acids Res 1999, 27(1):237–239. 10.1093/nar/27.1.237PubMed CentralView ArticlePubMedGoogle Scholar
- Manning BD, Cantley LC: Hitting the target: emerging technologies in the search for kinase substrates. Sci STKE 2002, 2002(162):PE49.PubMedGoogle Scholar
- Kim JH, Lee J, Oh B, Kimm K, Koh I: Prediction of phosphorylation sites using SVMs. Bioinformatics 2004, 20(17):3179–3184. 10.1093/bioinformatics/bth382View ArticlePubMedGoogle Scholar
- Biondi RM, Nebreda AR: Signalling specificity of Ser/Thr protein kinases through docking-site-mediated interactions. Biochem J 2003, 372(Pt 1):1–13. 10.1042/BJ20021641PubMed CentralView ArticlePubMedGoogle Scholar
- Obenauer JC, Cantley LC, Yaffe MB: Scansite 2.0: Proteome-wide prediction of cell signaling interactions using short sequence motifs. Nucleic Acids Res 2003, 31(13):3635–3641. 10.1093/nar/gkg584PubMed CentralView ArticlePubMedGoogle Scholar
- Uhlik MT, Temple B, Bencharit S, Kimple AJ, Siderovski DP, Johnson GL: Structural and evolutionary division of phosphotyrosine binding (PTB) domains. J Mol Biol 2005, 345(1):1–20. 10.1016/j.jmb.2004.10.038View ArticlePubMedGoogle Scholar
- Yaffe MB, Elia AE: Phosphoserine/threonine-binding domains. Curr Opin Cell Biol 2001, 13(2):131–138. 10.1016/S0955-0674(00)00189-7View ArticlePubMedGoogle Scholar
- Yaffe MB, Leparc GG, Lai J, Obata T, Volinia S, Cantley LC: A motif-based profile scanning approach for genome-wide prediction of signaling pathways. Nat Biotechnol 2001, 19(4):348–353. 10.1038/86737View ArticlePubMedGoogle Scholar
- Yaffe MB, Smerdon SJ: The use of in vitro peptide-library screens in the analysis of phosphoserine/threonine-binding domain structure and function. Annu Rev Biophys Biomol Struct 2004, 33: 225–244. 10.1146/annurev.biophys.33.110502.133346View ArticlePubMedGoogle Scholar
- Holland PM, Cooper JA: Protein modification: docking sites for kinases. Curr Biol 1999, 9(9):R329–31. 10.1016/S0960-9822(99)80205-XView ArticlePubMedGoogle Scholar
- Joughin BA, Tidor B, Yaffe MB: A computational method for the analysis and prediction of protein:phosphopeptide-binding sites. Protein Sci 2005, 14(1):131–139. 10.1110/ps.04964705PubMed CentralView ArticlePubMedGoogle Scholar
- Ficarro SB, McCleland ML, Stukenberg PT, Burke DJ, Ross MM, Shabanowitz J, Hunt DF, White FM: Phosphoproteome analysis by mass spectrometry and its application to Saccharomyces cerevisiae. Nat Biotechnol 2002, 20(3):301–305. 10.1038/nbt0302-301View ArticlePubMedGoogle Scholar
- Ballif BA, Villen J, Beausoleil SA, Schwartz D, Gygi SP: Phosphoproteomic analysis of the developing mouse brain. Mol Cell Proteomics 2004, 3(11):1093–1101. 10.1074/mcp.M400085-MCP200View ArticlePubMedGoogle Scholar
- Blom N, Gammeltoft S, Brunak S: Sequence and structure-based prediction of eukaryotic protein phosphorylation sites. J Mol Biol 1999, 294(5):1351–1362. 10.1006/jmbi.1999.3310View ArticlePubMedGoogle Scholar
- Blom N, Sicheritz-Ponten T, Gupta R, Gammeltoft S, Brunak S: Prediction of post-translational glycosylation and phosphorylation of proteins from the amino acid sequence. Proteomics 2004, 4(6):1633–1649. 10.1002/pmic.200300771View ArticlePubMedGoogle Scholar
- Huang HD, Lee TY, Tzeng SW, Horng JT: KinasePhos: a web tool for identifying protein kinase-specific phosphorylation sites. Nucleic Acids Res 2005, 33(Web Server issue):W226–9. 10.1093/nar/gki471PubMed CentralView ArticlePubMedGoogle Scholar
- Huang HD, Lee TY, Tzeng SW, Wu LC, Horng JT, Tsou AP, Huang KT: Incorporating hidden Markov models for identifying protein kinase-specific phosphorylation sites. J Comput Chem 2005, 26(10):1032–1041. 10.1002/jcc.20235View ArticlePubMedGoogle Scholar
- Xue Y, Zhou F, Zhu M, Ahmed K, Chen G, Yao X: GPS: a comprehensive www server for phosphorylation sites prediction. Nucleic Acids Res 2005, 33(Web Server issue):W184–7. 10.1093/nar/gki393PubMed CentralView ArticlePubMedGoogle Scholar
- Zhou FF, Xue Y, Chen GL, Yao X: GPS: a novel group-based phosphorylation predicting and scoring method. Biochem Biophys Res Commun 2004, 325(4):1443–1448. 10.1016/j.bbrc.2004.11.001View ArticlePubMedGoogle Scholar
- Diella F, Cameron S, Gemund C, Linding R, Via A, Kuster B, Sicheritz-Ponten T, Blom N, Gibson TJ: Phospho.ELM: a database of experimentally verified phosphorylation sites in eukaryotic proteins. BMC Bioinformatics 2004, 5(1):79. 10.1186/1471-2105-5-79PubMed CentralView ArticlePubMedGoogle Scholar
- Duda RO, Hart PE, Stork DG: Pattern classification. 2nd edition. Beijing, China Machine Press; 2004:680.Google Scholar
- Henikoff S, Henikoff JG: Amino acid substitution matrices from protein blocks. Proc Natl Acad Sci U S A 1992, 89(22):10915–10919.PubMed CentralView ArticlePubMedGoogle Scholar
- Modrof J, Lymperopoulos K, Roy P: Phosphorylation of bluetongue virus nonstructural protein 2 is essential for formation of viral inclusion bodies. J Virol 2005, 79(15):10023–10031. 10.1128/JVI.79.15.10023-10031.2005PubMed CentralView ArticlePubMedGoogle Scholar
- Horscroft NJ, Roy P: NTP binding and phosphohydrolase activity associated with purified bluetongue virus non-structural protein NS2. J Gen Virol 2000, 81(Pt 8):1961–1965.View ArticlePubMedGoogle Scholar
- Lymperopoulos K, Wirblich C, Brierley I, Roy P: Sequence specificity in the interaction of Bluetongue virus non-structural protein 2 (NS2) with viral RNA. J Biol Chem 2003, 278(34):31722–31730. 10.1074/jbc.M301072200View ArticlePubMedGoogle Scholar
- Taraporewala ZF, Chen D, Patton JT: Multimers of the bluetongue virus nonstructural protein, NS2, possess nucleotidyl phosphatase activity: similarities between NS2 and rotavirus NSP2. Virology 2001, 280(2):221–231. 10.1006/viro.2000.0764View ArticlePubMedGoogle Scholar
- Bonet C, Fernandez I, Aran X, Bernues J, Giralt E, Azorin F: The GAGA Protein of Drosophila is Phosphorylated by CK2. J Mol Biol 2005, 351(3):562–572. 10.1016/j.jmb.2005.06.039View ArticlePubMedGoogle Scholar
- Arrigoni G, Marin O, Pagano MA, Settimo L, Paolin B, Meggio F, Pinna LA: Phosphorylation of calmodulin fragments by protein kinase CK2. Mechanistic aspects and structural consequences. Biochemistry 2004, 43(40):12788–12798. 10.1021/bi049365cView ArticlePubMedGoogle Scholar
- Huang EJ, Reichardt LF: Trk receptors: roles in neuronal signal transduction. Annu Rev Biochem 2003, 72: 609–642. 10.1146/annurev.biochem.72.121801.161629View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.