Email updates

Keep up to date with the latest news and content from BMC Bioinformatics and BioMed Central.

This article is part of the supplement: Proceedings of the Tenth Annual MCBIOS Conference

Open Access Proceedings

Circles within circles: crosstalk between protein Ser/Thr/Tyr-phosphorylation and Met oxidation

R Shyama Prasad Rao12, Dong Xu23, Jay J Thelen12 and Ján A Miernyk124*

Author Affiliations

1 Department of Biochemistry, University of Missouri, Columbia, MO 65211, USA

2 Interdisciplinary Plant Group, University of Missouri, Columbia, MO 65211, USA

3 Computer Science, University of Missouri, Columbia, MO 65211, USA

4 Plant Genetics Research Unit, USDA, Agricultural Research Service, University of Missouri, Columbia, MO 65211, USA

For all author emails, please log on.

BMC Bioinformatics 2013, 14(Suppl 14):S14  doi:10.1186/1471-2105-14-S14-S14


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2105/14/S14/S14


Published:9 October 2013

© 2013 Prasad Rao et al.; licensee 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.

Abstract

Background

Reversible posttranslational protein modifications such as phosphorylation of Ser/Thr/Tyr and Met oxidation are critical for both metabolic regulation and cellular signalling. Although these modifications are typically studied individually, herein we describe the potential for cross-talk and hierarchical regulation.

Results

The proximity of Met to Ser/Thr/Tyr within the proteome has not previously been addressed. In order to consider the possibility of a generalized interaction, we performed a trans-kingdom sequence analysis of known phosphorylation sites in proteins from bacteria, fungi, plants, and animals. The proportion of phosphorylation sites that include a Met within a 13-residue window centered upon Ser/Thr/Tyr is significantly less than the occurrence of Met in proximity to all Ser/Thr/Tyr residues. Met residues are present at all positions (-6 to +6, inclusive) within the 13-residue window that we have considered. Detailed analysis of sequences from eight disparate plant taxa revealed that many conserved phosphorylation sites have a Met residue in the proximity. Results from GO enrichment analysis indicated that the potential for phosphorylation and Met oxidation crosstalk is most prevalent in kinases and proteins involved in signalling.

Conclusion

The large proportion of known phosphorylation sites with Met in the proximity fulfils the necessary condition for cross-talk. Kinases/signalling proteins are enriched for Met around phosphorylation sites. These proteins/sites are likely candidates for cross-talk between oxidative signalling and reversible phosphorylation.

Keywords:
Methionine; methionine sulfoxide; oxidation; protein phosphorylation; regulation; reversible posttranslational modifications; signalling

Background

The proteome includes the combination of at least three components; genome-encoded proteins, the products of alternative initiation and splicing, and of posttranslational modifications (PTM). This combination is orders of magnitude larger than could be solely genome-encoded [1]. The diversity of PTM is extraordinarily large; to date as many as 435 different PTMs are known [2-5]. They can be reversible (e.g., acetylation) or irreversible (e.g., proteolytic cleavage), enzyme catalysed (e.g., kinase-mediated phosphorylation) or the result direct chemical reaction (e.g., oxidation), and individual or combinatorial [6,7]. In toto, these factors comprise a well-adapted basis for signalling, regulation, targeting, and interaction, all at least potentially in the absence of de novo protein synthesis [3,8].

Among PTM, reversible phosphorylation has been the most extensively studied [9-11]. This is in part because reversible phosphorylation is a component of both cellular signalling [8,12] and direct regulatory control of protein function [13-15]. Protein phosphorylation can be stoichiometric at an individual site, or combinatorial and sequential/hierarchical at multiple sites [16,17]. During the 'genomic era,' an enormous number of different protein kinases and phospho-protein (P-protein) phosphatases were identified, and reversible protein phosphorylation has been shown to be involved in regulating many fundamental cellular processes. It has been reported that the combination of genes encoding protein kinases and P-protein phosphatases constitute between 2 and 4% of the total number of genes in a typical eukaryotic genome [18]. Saccharomyces cerevisiae has genes for 113 protein kinases while the Homo sapiens has 518 such genes [19]. The genome of the reference dicot plant Arabidopsis thaliana includes genes encoding 1055 protein kinases [19] and 112 P-protein phosphatase catalytic subunits [20]. Clearly, this PTM is ubiquitous across all branches of the tree of life [21].

Reactive oxygen species (ROS) are an inevitable consequence of aerobic metabolism [22-24]. High levels of ROS cause oxidative damage to proteins, which is lethal unless repaired or reversed [25,26]. In contrast, low levels of ROS have evolved functional roles in many aspects of cellular signalling [27-29]. The reversible chemical oxidation of Met to methionine sulfoxide (MetSO) spans the two extremes in response to ROS. Oxidation of Met to MetSO can damage/inactivate proteins, but because this PTM is enzymatically reversible [30,31] it can participate in the types of regulation [32] and signalling [33,34] typified by reversible phosphorylation. Met oxidation is a particularly important component of cellular responses to oxidative stress [35,36].

In addition to Ser/Thr/Tyr phosphorylation and Met oxidation, other common reversible PTMs include Ser/Thr O-glycosylation [37], Lys/Arg methylation [38], and Lys acetylation [39]. There are examples of each directly regulating protein activities as well as playing roles in cellular signalling [40-43]. While a single PTM can modulate protein function, combinatorial and sequential/hierarchical interplay among two or more types of PTM can integrate signals from multiple pathways [6,44,45]. Dynamic interactions among PTM can take place at the same target residue [5,46], or there can be cross-talk between or among multiple residues [47,48]. In some instances the PTM of a specific residue can require a prior (primed) PTM at another site [6,17,49].

It was recently suggested that reversible Met oxidation might serve as a rheostat to control proximal Ser/Thr/Tyr phosphorylation [50,51], and that this could be a mechanism for integrating information from both metabolic and ROS-based signalling pathways [52,53]. However, it is not known how often phosphorylation sites contain Met residues in their proximity. It is important to identify candidate proteins with coupled Met oxidation and phosphorylation sites. The results described herein are from a study undertaken to test the possibility that phosphorylation and Met oxidation crosstalk has a broad occurrence.

Methods

Proteome sequences and phosphorylation sites

In order to consider our analysis of potential crosstalk in as broad a context as possible we have employed a data-mining strategy. Selected reference proteome or complete proteome sequences were downloaded from Uniprot [54], NCBI [55], or TAIR [56]. Experimentally confirmed phosphorylation (Ser/Thr/Tyr) site information was obtained from multiple databases [21,57]. Plant-specific phosphorylation site data was obtained from P3DB [58]. Our aim was to explore the proportion of phosphorylated (and potentially phosphorylatable) sites that contain a Met residue within 13 residue window in different taxa. We used sequences and sites from several species of plants, yeast, animals, bacteria, and archaea (Additional file 1). Because Met oxidation is an aerobic process, we also included species that have an anaerobic metabolism/lifestyle such as Bacteroides fragilis, Clostridium botulinum, and Ascaris suum.

Additional file 1. S1.pdf. Percentage of Ser/Thr/Tyr sites (Ser, Thr, and Tyr from top to bottom) with Met within a ± 6 window in different species.

Format: PDF Size: 225KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Phosphorylation site predictions

Phosphorylation sites were predicted for complete proteomes using Musite [59], NetPhosBac [60], or Disphos [61]. Predictions were based on taxa-specific Ser/Thr/Tyr models at 95% specificity for the respective proteomes. Predictions were made to compare the differences, if any, between the proportion of Met residues in the vicinity of experimentally known phosphorylation sites and predicted phosphorylation sites.

Sequence analyses

The proportion of Ser/Thr/Tyr sites that include at least one Met-residue in a window ranging up to 21-residues was scanned. As most known phosphorylation motifs are less than 13-residues in length [62], a 13-residue window (six residues either side of the phosphorylation site) was selected for our analyses. Amino acid frequencies for individual proteomes were computed. Amino-terminal Met-residues were discounted in both calculations of Met frequency and sequence-distribution analyses. The hydropathy score for each site (Ser/Thr/Tyr centered in a 13-residue window; the phosphorylation site itself was not included) was calculated using the Kyte-Doolittle scale [57,63].

In order to evaluate the extent of conservation of phosphorylation sites that contain a Met, and to identify candidate proteins for Phos-Met oxidation crosstalk in plants, a PSI-BLAST search was conducted for Arabidopsis thaliana, Glycine max, and Oryza sativa proteins that include a Met-residue within the 13-residue phosphorylation-site window. Proteins from these three taxa account for ~90% of known phosphorylation sites in plants. The results were used to identify related sequences in five additional proteomes (Pinus taeda, Selaginella moellendorffii, Physcomitrella patens, Chlamydomonas reinhardtii, and Cyanidioschyzon merolae), using a significance E-value cutoff of 0.00001. The total number of phosphorylation sites with Met-conservation was then calculated among the eight disparate taxa. Additionally, the proportion of hydrophobic residues (Met/Leu/Ile/Val/Phe/) was enumerated for each phosphorylation site.

Gene ontology enrichment analysis

A gene ontology (GO) term-enrichment analysis was performed for proteins that contain at least one known phosphorylation site including a Met within our specified ± 6-residue window, using the agriGO webserver [64] for A. thaliana or the g:Profiler webserver [65] for S. cerevisiae and H. sapiens. A hypergeometric test (with Bonferroni correction for multiple comparisons) was used to select significantly enriched (p < 0.05) terms [64,65]. A log2 ratio of numbers of each significantly enriched GO term in the test set (phosphoproteins that have Met residues in the vicinity of phosphorylation site) to background set (all phosphoproteins) was calculated. All data analyses were performed using Python ver. 2.7. A Z-test was used to calculate the difference between two group proportions; p < 0.05 was considered significant.

Results

Evaluating the translated proteome from diverse taxa, it is clear that Met belongs to the group of relatively low-abundance amino acids along with Cys, His, Trp, and Tyr (Figure 1). In plants, yeast, and animals Ser has a relatively high abundance, and Thr is of intermediate abundance in all taxa. The low Met abundance is in contrast to the other members of the hydrophobic group (Leu, Ile, Val, and Phe) all of which are moderately to highly abundant.

thumbnailFigure 1. Proportions (%) of amino acids present in the translated proteomes from different taxa. Met is the least abundant hydrophobic residue. Proportions of Ser and Thr are nearly equal in bacteria and archaea.

The proportion of Ser/Thr/Tyr sites which include a Met within the window of ± 6 residues is shown in Figure 2. On average 14% of phosphorylated sites contained a Met. Overall, the proportions of phosphorylated Ser/Thr/Tyr sites (experimentally determined and predicted) including a Met was significantly less (p ~ 0, Z-test for two proportions) compared to proteome-wide Ser/Thr/Tyr sites. Proportionally, Tyr sites contain more Met than Ser and Thr sites in eukaryotes. The apparent deviation from this pattern seen with archaeal and bacterial proteins can likely be attributed to the small total number of phosphorylated sites in these organisms. Proteome specific proportions within individual taxa are shown in Additional file 1.

thumbnailFigure 2. The percentage of Ser/Thr/Tyr sites (individually from left to right in each group) including a Met residue within ± 6 window in different taxa. On average, phosphorylated sites contain less Met compared to all Ser/Thr/Tyr sites in proteomes. In eukaryotes, Tyr sites contain proportionally more Met compared to Ser and Thr sites. The large deviation from the otherwise 'normal' pattern seen with bacteria and archaea is likely due to the small number of phosphorylated sites.

The overall positional distribution of Met in the proteome of eukaryotic organisms is nearly uniform. In Figure 3 we have plotted the distribution of Met within the ± 6 residue window centered on; any amino acid ("non-STY"), Ser, Thr, or Tyr (STY), or P-Ser, P-Thr, or P-Tyr (Phospho-STY). The distribution of Met flanking P-Tyr residues is enriched at positions -4, -3, and +5 in plant proteins, and to a lesser degree at position +3 in all but yeast proteins. Considering the large sample size, the position-specific proportional differences are significant (p ~ 0, χ2 test for independence). In contrast, the apparent enrichment of Met at several positions in bacterial or archaeal proteins is based upon much smaller sample numbers and is not significant.

thumbnailFigure 3. Distribution of Met residues (%) around non-STY and STY (individually from left to right, all and phosphorylated) sites in plants (A), yeast (B), animals (C), bacteria (D) and archaea (E) (n = 10353, 3645, 21562, 126 and 63 phosphoproteins respectively).

The patterns are somewhat more complex when comparing the distribution of other members of the hydrophobic amino acid homology group (Leu, Ile, Val, and Phe) within the ± 6 residue window used to evaluate Met distribution (Additional file 2). With few exceptions, plant, yeast, and animal proteins have a lower occurrence of a hydrophobic residue within the ± 6 residue window centered on a P-Ser/Thr/Tyr residue. The exceptions include yeast position +1, and mammalian positions -1 and +3 from a P-Tyr residue (Additional file 2). Overall distribution of hydrophobic residues in the ± 6 window of prokaryotic and archaeal proteins is quite uneven, but it must be considered that the total number of phosphorylation sites in these proteins is much lower than in eukaryotic proteins (e.g., E. coli has been reported to have only 32 P-Tyr-containing proteins and two tyrosine-kinases [66].

Additional file 2. S2.pdf. Distribution of other hydrophobic residues (Val/Ile/Leu/Phe, %) around non-Ser/Thr/Tyr and Ser/Thr/Tyr (individually from left to right, all and phosphorylated) sites in plant (A), yeast (B), animal (C), bacterial (D) and archaeal (E) P-proteins.

Format: PDF Size: 39KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

We conducted GO analysis of the sequences of all the plant, yeast, and animal proteins that have a Met-containing phosphorylated-site, in order to test for enrichment in functional domains/motifs. The analyses were conducted versus a total P-proteins background (Figure 4). Bacteria and archaeal proteins were not included in this analysis because of the small sample sizes. Only the GO terms that are significantly enriched (p < 0.05, hypergeometric test) in our dataset are displayed; most are related to protein kinases, phosphorylation, signalling, and regulation. There is also enrichment in proteins annotated as having an involvement in cell division in yeast or animals, but not in plants (Figure 4). For example, a log2 value of 0.5 indicates 41% increase in a specific GO term in the test set (proteins that contain a Met near phosphorylation site) as against the background set (all phosphoproteins).

thumbnailFigure 4. Results from GO analyses. The proteins (test set) that contain a Met residue within ± 6 window around known phosphorylation sites are enriched (p < 0.05, hypergeometric test and z-score > 1.96; X-axis shows the log2 ratio of number of each significantly enriched GO term in test set against background set) for several GO terms against the total phosphoproteins (background set) in plant (A. thaliana, A), yeast (S. cerevisiae, B) and animal (H. sapiens, C). In general Met oxidation and phosphorylation cross-talk are more likely in protein kinases, as several GO terms related to kinase activity as enriched in all three taxa. (BP: Biological Process, MF: Molecular Function, and CC: Cellular Component).

We assume that a phosphorylation-site that includes a Met at the same position within the ± 6 residue-window among diverse taxa would have the greatest probability for regulatory crosstalk. We specifically searched proteins from 8 plant taxa in order to evaluate potential candidates. A list of the consensus 13-residue sequences for potentially conserved P-Ser/Thr/Tyr sites which include a Met residue is presented in Additional file 3. The extent of Met positional-conservation is presented as a heat map in Figure 5A. Several sites with a "high degree" of positional conservation for Met (≥ 6/8 homologs) were detected. At higher Met conservation, most of the non-Met containing remaining sequences have a hydrophobic (Val/Ile/Leu/Phe) residue at that position (Figure 5B).

Additional file 3. S3.pdf. Consensus sequences (13-mers) for phosphorylation sites with highly conserved Met (≥ 7 Met among eight taxa).

Format: PDF Size: 724KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

thumbnailFigure 5. The occurrence of Met versus other hydrophobic amino acids vicinal to phosphorylation sites. (A) The heat map indicates the proportion of phosphorylation sites at each level of Met conservation. Many phosphorylation sites with conserved (e.g., ≥ 6/8) Met are present across different (eight) taxa. (B) Proportion (%) of sequences (in disparate taxa) including a hydrophobic (Val/Ile/Leu/Phe) residue. At higher Met conservation, a higher proportion (%) of remaining sequences includes a hydrophobic residue in place of Met.

The proteins that contain an experimentally verified phosphorylation site with a vicinal Met show significant enrichment (p < 0.05, hypergeometric test) for many GO terms when the Met is conserved (panel A ≥7, B ≥6, C ≥5 and D ≥4) in eight disparate plant taxa (Figure 6). Many significant GO terms are enriched several-fold; kinase activity (GO:0004712), for example. However, as the Met conservation level decreases the GO term enrichment also decreases until it is no longer significant when compared to the background of total phosphoproteins. The results from GO term analysis indicate that Met and phosphorylation site cross-talk is likely to be more prevalent in kinases and various stress-related proteins (Figure 6).

thumbnailFigure 6. Results from GO term enrichment analyses. The A. thaliana proteins (test set) that contain a known phosphorylation site (based upon published data from A. thaliana, G. max, or O. sativa) including a Met residue within the ± 6 residue window shows significant enrichment (p < 0.005, hypergeometric test, X-axis shows the log2 ratio of number of each significantly enriched GO term in test set against background set) for several GO terms when the Met is conserved (panel A ≥7, B ≥6, C ≥5 and D ≥4) in eight disparate taxa. As the Met conservation level decreases the significance of the GO term enrichment decreases compared to the background of A. thaliana P-proteins. The GO term enrichment results agree with the proposal that Met oxidation:phosphorylation-site cross-talk is most likely in kinases (green bars) and various stress-related proteins (red bars).

Figure 7 shows some of the possible scenarios in the phosphorylation and Met-oxidation crosstalk. Depending on the extent of evolutionary conservation of Met (Figure 5, 7), it could play a regulatory role if it is well conserved across species (for example, senescence-associate protein: AT1G66580) or as a structural role if it serves merely as a hydrophobic residue along with VILF residues (for example, as in ATPase F1º: AT4G04640). The relevance of the crosstalk may also depend on the distance between the phosphorylation site and Met residue. The farther away the Met residue is, the less likely it is involved in crosstalk.

thumbnailFigure 7. Crosstalk between Met and Ser/Thr/Tyr might be related to the nature/extent of conservation of Met in proteins. (I) Met could play a regulatory role (crosstalk) (for example, senescence associate protein: AT1G66580) when it is conserved across species. (II) It might play only a structural role in a protein (for example, ATPase F1º AT4G04640), present as a hydrophobic residue. (III-V) If either the phosphorylation site or the Met residue itself is not conserved, it is unlikely that crosstalk will evolve in such proteins (for example, Poly(A) binding protein: AT3G16380.1; peptidylprolyl cis/trans isomerase: AT2G18040.1; P450 reductase: AT4G24520.1). Crosstalk can also depend on the distance between phosphorylation site and Met residue in the 1° sequence, and the 3° structural context. (ara: Arabidopsis thaliana, gly: Glycine max, ory: Oryza sativa, pin: Pinus Taeda, Sel: Selaginella moellendorffii, phy: Physcomitrella patens, chl: Chlamydomonas reinhardtii, and cya: Cyanidioschyzon merolae).

Discussion

Two unrelated plant enzymes, cytoplasmic nitrate reductase [50], and mitochondrial pyruvate dehydrogenase [51], have a Met residue proximal to a regulatory phosphorylation site, and in both cases oxidation of this Met to MetSO inhibits phosphorylation. It has been suggested that the redox status of these Met residues allows the enzymes to monitor oxidative stress, and crosstalk between these two PTMs can serve to fine-tune activity/metabolism [52,53]. The study described herein was undertaken to evaluate how widespread the potential is for this sort of regulation by crosstalk between PTMs.

A 13-amino acid motif centered on the phosphorylatable Ser/Thr/Tyr-residue is apparently adequate for recognition by protein kinases [62,67]. Multiple distinct factors likely contributed to evolution of this motif, including the occurrence of each of the amino acids [68,69] and the constraints inherent to protein structure [70].

Given that thousands of phosphorylation sites are known in proteomes spanning taxa [57], the presence of Met in ~14% of these sites (Figure 2, Additional file 1) translates into a number of potential sites and proteins for phosphorylation and Met-oxidation crosstalk. Many phosphorylation motifs contain hydrophobic residues (Met/Val/Ile/Leu/Phe). This presumably reflects some essential aspect of protein structure or folding. However, Met is less abundant than the other hydrophobic residues [71], and is less likely to accommodate a structural requirement [72]. Furthermore, biosynthesis of Met is more energetically expensive than other members of the hydrophobic group [69,71,73]. Considering all these constraints, the presence of Met in the phosphorylation motif would need to be of some evolutionary benefit otherwise it would have been lost during subsequent optimizing of the regulatory network [74]. We propose that the benefit might have been acquisition of the ability to monitor ROS status by reversible Met oxidation, a capability not available with other hydrophobic residues [75]. Even non-surface exposed Met residues are susceptible to oxidation because of local polypeptide flexibility [76].

Met residue in the neighbourhood of phosphorylation sites can have two distinct roles (Figure 7). It could have a solely structural role in which Met acts as a hydrophobic residue in a kinase recognition motif. If Met is replaced by other hydrophobic (VILF) residues in related species, then it is likely to have a structural role. On the other hand, Met could play a regulatory role if it is conserved near phosphorylation site in a protein from diverse species [75]. We identified a large number of plant proteins and phosphorylation sites including Met that are highly conserved (Figure 5, Figure 7 and Additional file 3). These proteins are potential candidates wherein Met oxidation may be involved in coupling oxidative signals with phosphorylation and regulation of protein function.

This proposed evolutionary acquisition of Met/MetSO as a mechanism for fine-control of phosphorylation status might be addressed by comparison of phosphorylation sites in proteins from aerobic microorganisms with those from closely-related but obligate anaerobes. If Met oxidation exerts an influence (selection) on phosphorylation then it would be predicted that one group would have a significantly different proportion of (potential) phosphorylation sites including a Met residue. At least superficially this is not the case. In the very limited data available to us, aerobes (e.g., Escherichia coli and Staphylococcus aurius) and anaerobes (e.g., Bacillus fragilis and Clostridium botulinum) have a similar proportion of sites with Met (Additional file 1). It would be particularly interesting to evaluate the P-proteome of an organism such as the nematode Ascaris suum which has both aerobic and anaerobic components of its' life-cycle. Unfortunately, no P-proteomic data are currently available.

While the regulatory importance of protein PTM has been long known, it has only been more recently that roles for multiple and interacting PTM's have been appreciated [77]. Examples range from interaction between multiple instances of the same PTM at different sites (e.g., prior or priming phosphorylation of site A is necessary for subsequent phosphorylation of site B [14-17]), through hierarchical responses to multiple PTM of the same site [5,37,46], or differential responses to multiple different PTM at different sites [78], and ultimately to crosstalk among different PTM at different sites [79]. A protein with 10 instances of a single PTM could give rise to 1024 (2n) different molecular species, but if there are two types of PTM there is the potential to generate 59049 (3n) species [46]. The enormity of this molecular diversity would generate a gradient in response that will allow an exquisitely fine control over metabolic activity, signalling, and cellular function [46].

Our initial interest in crosstalk between Ser/Thr/Tyr-phosphorylation and Met oxidation was stimulated by observations that two important metabolic enzymes could be responsive to cellular redox signalling through this mechanism [50-52]. While this study was underway, it came to our attention that there have been additional reports of what now appears to be crosstalk between phosphorylation and Met oxidation. Decreased phosphorylation is correlated with increased Met oxidation of α-synuclein in MetSO knockdown mutants of both yeast and mice [80,81]. Site-specific oxidation of Met in calmodulin was shown to affect its structure, and thus its interaction with Ca2+/calmodulin-dependent protein kinase II and subsequent downstream signalling [82]. There is also an analogous example [83] where activity of a protein P-Tyr-phosphatase is regulated by reversible Cys-oxidation. Finally, there is the intriguing recent report that Met oxidation in A. thaliana is responsive to cGMP [84].

Results from GO analysis indicated that the occurrence of Met near phosphorylation sites could be more prevalent in proteins related to signalling such as kinases and stress related proteins, and could be common to all three taxa - plants, yeast, and animals (Figure 4 and 6). Many phosphorylation motifs - kinase recognition signatures - contain hydrophobic Met residue [57,62]. Oxidation of Met to MetSO causes a shift from hydrophobic to hydrophilic in nature, and would thus likely disrupt kinase recognition [50]. So it is potentially useful for kinases and stress signalling proteins to have evolutionary selection for Met near phosphorylation sites for potential crosstalk. This would allow direct communication of oxidative signals to other mainstream processes.

Conclusion

Our study identified a large number of phosphorylation sites which include a vicinal Met residue. The proteins containing these sites can potentially function as redox-sensors capable of transducing input from ROS signalling to regulation of phosphorylation. These observations should stimulate further research on PTM crosstalk and control of protein function in response to oxidative signalling.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

RSPR participated in design of the study and conducted the computational analyses. JJT and DX participated in design and coordination of the study. JAM conceived the study, participated in its design and coordination, and drafted the manuscript. All authors read and approved the final manuscript.

Acknowledgements

This work was supported in part by National Science Foundation-Plant Genome Research Program Grant DBI-0604439 to J.J.T.

Declarations

Publication costs for this article were funded by the United Stated Department of Agriculture, Agricultural Research Service, Project number 3622-21000-032-00D.

This article has been published as part of BMC Bioinformatics Volume 14 Supplement 14, 2013: Proceedings of the Tenth Annual MCBIOS Conference. Discovery in a sea of data. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcbioinformatics/supplements/14/S14.

References

  1. Prabakaran S, Lippens G, Steen H, Gunawardena J: Post-translational modification: nature's escape from genetic imprisonment and the basis for dynamic information encoding.

    Wiley Interdiscip Rev Syst Biol Med 2012, 4:565-583. PubMed Abstract | Publisher Full Text OpenURL

  2. Mann M, Jensen ON: Proteomic analysis of post-translational modifications.

    Nature Biotechnol 2003, 21:255-261. Publisher Full Text OpenURL

  3. Walsh CT, Garneau-Tsodikova S, Gatto GJ Jr: Protein posttranslational modifications: The chemistry of proteome diversifications.

    Angew Chem Int Ed Engl 2005, 44:7342-7372. PubMed Abstract | Publisher Full Text OpenURL

  4. Khoury GA, Baliban RC, Floudasa CA: Proteome-wide post-translational modification statistics: Frequency analysis and curation of the swiss-prot database.

    Sci Rep 2011, 1:90. OpenURL

  5. Minguez P, Parca L, Diella F, Mende DR, Kumar R, Helmer-Citterich M, Gavin AC, van Noort V, Bork P: Deciphering a global network of functionally associated post-translational modifications.

    Molec Syst Biol 2012, 8:599. OpenURL

  6. Hunter T: The age of crosstalk: Phosphorylation, ubiquitination and beyond.

    Mol Cell 2007, 28:730-738. PubMed Abstract | Publisher Full Text OpenURL

  7. Beltrao P, Albanèse V, Kenner LR, Swaney DL, Burlingame A, Villén J, Lim WA, Fraser JS, Frydman J, Krogan NJ: Systematic functional prioritization of protein posttranslational modifications.

    Cell 2012, 150:413-425. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Deribe YL, Pawson T, Dikic I: Post-translational modifications in signal integration.

    Nature Struct Mol Biol 2010, 17:666-672. Publisher Full Text OpenURL

  9. Loroch S, Dickhut C, Zahedi RP, Sickmann A: Phosphoproteomics - More than meets the eye.

    Electrophoresis 2013, in press.

    doi: 10.1002/elps.201200710.

    OpenURL

  10. Havelund JF, Thelen JJ, Møller IM: Biochemistry, proteomics, and phosphoproteomics of plant mitochondria from non-photosynthetic cells.

    Front Plant Sci 2013, 4:51. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  11. Mijakovic I, Macek B: Impact of phosphoproteomics on studies of bacterial physiology.

    FEMS Microbiol Rev 2012, 36:877-892. PubMed Abstract | Publisher Full Text OpenURL

  12. Rogne M, Taskén K: Cell signalling analyses in the functional genomics era.

    Nature Biotechnol 2013, 30:333-338. OpenURL

  13. Furuya E, Yokoyama M, Uyeda K: Regulation of fructose-6-phosphate 2-kinase by phosphorylation and dephosphorylation: possible mechanism for coordinated control of glycolysis and glycogenolysis.

    Proc Natl Acad Sci USA 1982, 79:325-329. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Dexter JP, Gunawardena J: Dimerization and bifunctionality confer robustness to the isocitrate dehydrogenase regulatory system in Escherichia coli.

    J Biol Chem 2013, 288:5770-5778. PubMed Abstract | Publisher Full Text OpenURL

  15. Oliveira AP, Ludwig C, Picotti P, Kogadeeva M, Aebersold R, Sauer U: Regulation of yeast central metabolism by enzyme phosphorylation.

    Mol Syst Biol 2012, 8:623. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Just S, Illing S, Trester-Zedlitz M, Lau EK, Kotowski SJ, Miess E, Mann A, Doll C, Trinidad JC, Burlingame AL, von Zastrow M, Schulz S: Differentiation of opioid drug effects by hierarchical multi-site phosphorylation.

    Mol Pharmacol 2013, 83:633-639. PubMed Abstract | Publisher Full Text OpenURL

  17. Ahsan N, Huang Y, Tovar-Mendez A, Swatek KN, Zhang J, Miernyk JA, Xu D, Thelen JJ: A versatile mass spectrometry-based method to both identify kinase client-relationships and characterize signaling network topology.

    J Proteome Res 2013, 12:937-948. PubMed Abstract | Publisher Full Text OpenURL

  18. Moorhead GBG, De Wever V, Templeton G, Kerk D: Evolution of protein phosphatases in plants and animals.

    Biochem J 2002, 417:401-409. OpenURL

  19. Niedner RH, Buzko OV, Haste NM, Taylor A, Gribskov M, Taylor SS: Protein kinase resource: an integrated environment for phosphorylation research.

    Proteins 2006, (63):78-86. OpenURL

  20. Kerk D, Bulgrien J, Smith DW, Barsam B, Veretnik S, Gribskov M: The complement of protein phosphatase catalytic subunits encoded in the genome of Arabidopsis.

    Plant Physiol 2002, 129:908-925. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. Aivaliotis M, Macek B, Gnad F, Reichelt P, Mann M, Oesterhelt D: Ser/Thr/Tyr protein phosphorylation in the archaeon Halobacterium salinarum - A representative of the third domain of life.

    PLoS ONE 2009, 4:e4777. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Hatanaka E, Dermargos A, Hirata AE, Vinolo MA, Carpinelli AR, Newsholme P, Armelin HA, Curi R: Oleic, linoleic and linolenic acids increase ROS production by fibroblasts via NADPH oxidase activation.

    PLoS One 2013, 8:e58626. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. Grivennikova VG, Vinogradov AD: Partitioning of superoxide and hydrogen peroxide production by mitochondrial respiratory complex I.

    Biochim Biophys Acta 2013, 1827:446-454. PubMed Abstract | Publisher Full Text OpenURL

  24. Panieri E, Gogvadze V, Norberg E, Venkatesh R, Orrenius S, Zhivotovsky B: Reactive oxygen species generated in different compartments induce cell death, survival, or senescence.

    Free Radic Biol Med 2013, 57:176-187. PubMed Abstract | Publisher Full Text OpenURL

  25. Møller IM, Jensen PE, Hansson A: Oxidative modifications to cellular components in plants.

    Annu Rev Plant Biol 2007, 58:459-481. PubMed Abstract | Publisher Full Text OpenURL

  26. Höhn A, König J, Grune T: Protein oxidation in aging and the removal of oxidized proteins.

    J Proteomics 2013, in press.

    doi:pii: S1874-3919(13)00021-3.

    OpenURL

  27. Møller IM, Sweetlove LJ: ROS signalling--specificity is required.

    Trends Plant Sci 2010, 15:370-374. PubMed Abstract | Publisher Full Text OpenURL

  28. Aiken CT, Kaake RM, Wang X, Huang L: Oxidative stress-mediated regulation of proteasome complexes.

    Mol Cell Proteomics 2011, 10:R110-006924. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Wang P, Du Y, Zhao X, Miao Y, Song CP: The MPK6-ERF6-ROS-responsive cis-acting Element7/GCC box complex modulates oxidative gene transcription and the oxidative response in Arabidopsis.

    Plant Physiol 2013, 161:1392-1408. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Boschi-Muller S, Olry A, Antoine M, Branlant G: The enzymology and biochemistry of methionine sulfoxide reductases.

    Biochim Biophys Acta 2005, 1703:231-238. PubMed Abstract | Publisher Full Text OpenURL

  31. Boschi-Muller S, Gand A, Branlant G: The methionine sulfoxide reductases: Catalysis and substrate specificities.

    Arch Biochem Biophys 2008, 474:266-273. PubMed Abstract | Publisher Full Text OpenURL

  32. Stadtman ER, Moskovitz J, Levine RL: Oxidation of methionine residues of proteins: biological consequences.

    Antioxid Redox Signal 2003, 5:577-582. PubMed Abstract | Publisher Full Text OpenURL

  33. Carruthers NJ, Stemmer PM: Methionine oxidation in the calmodulin-binding domain of calcineurin disrupts calmodulin binding and calcineurin activation.

    Biochemistry 2008, 47:3085-3095. PubMed Abstract | Publisher Full Text OpenURL

  34. Bigelow DJ, Squier TC: Thioredoxin-dependent redox regulation of cellular signaling and stress response through reversible oxidation of methionines.

    Mol Biosyst 2011, 7:2101-2109. PubMed Abstract | Publisher Full Text OpenURL

  35. Erickson JR, Joiner ML, Guan X, Kutschke W, Yang J, Oddis CV, Bartlett RK, et al.: A dynamic pathway for calcium-independent activation of CaMKII by methionine oxidation.

    Cell 2008, 133:462-474. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Bigelow DJ, Squier TC: Redox modulation of cellular signaling and metabolism through reversible oxidation of methionine sensors in calcium regulatory proteins.

    Biochim Biophys Acta 2005, 1703:121-134. PubMed Abstract | Publisher Full Text OpenURL

  37. Zauner G, Kozak RP, Gardner RA, Fernandes DL, Deelder AM, Wuhrer M: Protein O-glycosylation analysis.

    Biol Chem 2012, 393:687-708. PubMed Abstract | Publisher Full Text OpenURL

  38. Erce MA, Pang CN, Hart-Smith G, Wilkins MR: The methylproteome and the intracellular methylation network.

    Proteomics 2012, 12:564-586. PubMed Abstract | Publisher Full Text OpenURL

  39. Zhang K, Tian S, Fan E: Protein lysine acetylation analysis: current MS-based proteomic technologies.

    Analyst 2013, 138:1628-1636. PubMed Abstract | Publisher Full Text OpenURL

  40. Theillet FX, Smet-Nocca C, Liokatis S, Thongwichian R, Kosten J, Yoon MK, Kriwacki RW, et al.: Cell signaling, post-translational protein modifications and NMR spectroscopy.

    J Biomol NMR 2012, 54:217-236. PubMed Abstract | Publisher Full Text OpenURL

  41. Shi SP, Qiu JD, Sun XY, Suo SB, Huang SY, Liang RP: PMeS: prediction of methylation sites based on enhanced feature encoding scheme.

    PLoS One 2012, 7:e38772. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  42. Donlin LT, Andresen C, Just S, Rudensky E, Pappas CT, Kruger M, Jacobs EY, et al.: Smyd2 controls cytoplasmic lysine methylation of Hsp90 and myofilament organization.

    Genes Dev 2012, 26:114-119. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  43. Gill DJ, Clausen H, Bard F: Location, location, location: new insights into O-GalNAc protein glycosylation.

    Trends Cell Biol 2011, 21:149-158. PubMed Abstract | Publisher Full Text OpenURL

  44. Young NL, Plazas-Mayorca MD, Garcia BA: Systems-wide proteomic characterization of combinatorial post-translational modification patterns.

    Expert Rev Proteomics 2010, 7:79-92. PubMed Abstract | Publisher Full Text OpenURL

  45. Hart GW, Slawson C, Ramirez-Correa G, Lagerlof O: Cross talk between O-GlcNAcylation and phosphorylation: Roles in signaling, transcription, and chronic disease.

    Annu Rev Biochem 2011, 80:825-858. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  46. Hart GW, Housley MP, Slawson C: Cycling of O-linked beta-N-acetylglucosamine on nucleocytoplasmic proteins.

    Nature 2007, 446:1017-1022. PubMed Abstract | Publisher Full Text OpenURL

  47. Lu Z, Cheng Z, Zhao Y, Volchenboum SL: Bioinformatic analysis and post-translational modification crosstalk prediction of lysine acetylation.

    PLoS One 2011, 6:e28228. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. van Noort V, Seebacher J, Bader S, Mohammed S, Vonkova I, Betts MJ, Kühner S, et al.: Cross-talk between phosphorylation and lysine acetylation in a genome-reduced bacterium.

    Mol Syst Biol 2012, 8:571. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  49. Sharfi H, Eldar-Finkelman H: Sequential phosphorylation of insulin receptor substrate-2 by glycogen synthase kinase-3 and c-Jun NH2-terminal kinase plays a role in hepatic insulin signaling.

    Am J Physiol Endocrinol Metab 2008, 294:E307-315. PubMed Abstract | Publisher Full Text OpenURL

  50. Hardin SC, Larue CT, Oh MH, Jain V, Huber SC: Coupling oxidative signals to protein phosphorylation via methionine oxidation in Arabidopsis.

    Biochem J 2009, 422:305-312. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  51. Miernyk JA, Johnston ML, Huber SC, Tovar-Méndez A, Hoyos E, Randall DD: Oxidation of an adjacent methionine residue inhibits regulatory seryl-phosphorylation of pyruvate dehydrogenase.

    Proteomics Insights 2009, 2:15-22. OpenURL

  52. Emes MJ: Oxidation of methionine residues: the missing link between stress and signalling responses in plants.

    Biochem J 2009, 422:e1-2. PubMed Abstract | Publisher Full Text OpenURL

  53. Zhang X-H: Regulation of protein function by residue oxidation.

    Proteomics Insights 2010, 3:17-24. OpenURL

  54. Uniprot [http://www.uniprot.org/] webcite

  55. NCBI [http://www.ncbi.nlm.nih.gov/] webcite

  56. TAIR [http://www.arabidopsis.org/] webcite

  57. Rao RSP, Møller IM: Large-scale analysis of phosphorylation site occupancy in eukaryotic proteins.

    Biochim Biophys Acta 2012, 1824:405-412. PubMed Abstract | Publisher Full Text OpenURL

  58. P3DB [http://www.p3db.org/] webcite

  59. Gao J, Thelen JJ, Dunker AK, Xu D: Musite, a tool for global prediction of general and kinase-specific phosphorylation sites.

    Mol Cell Proteomics 2010, 9:2586-2600. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  60. Miller ML, Soufi B, Jers C, Blom N, Macek B, Mijakovic I: NetPhosBac - A predictor for Ser/Thr phosphorylation sites in bacterial proteins.

    Proteomics 2009, 9:116-125. PubMed Abstract | Publisher Full Text OpenURL

  61. Iakoucheva LM, Radivojac P, Brown CJ, O'Connor TR, Sikes JG, Obradovic Z, Dunker AK: The importance of intrinsic disorder for protein phosphorylation.

    Nucleic Acids Res 2004, 32:1037-1049. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  62. Beausoleil SA, Jedrychowski M, Schwartz D, Elias JE, Villén J, Li J, Cohn MA, Cantley LC, Gygi SP: Large-scale characterization of HeLa cell nuclear phosphoproteins.

    Proc Natl Acad Sci USA 2004, 101:12130-12135. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  63. Kyte J, Doolittle RF: A simple method for displaying the hydropathic character of a protein.

    J Mol Biol 1982, 157:105-132. PubMed Abstract | Publisher Full Text OpenURL

  64. Du Z, Zhou X, Ling Y, Zhang Z, Su Z: AgriGO: A GO analysis toolkit for the agricultural community.

    Nucleic Acids Res 2010, 38:W64-70. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  65. Reimand J, Kull M, Peterson H, Hansen J, Vilo J: g:Profiler - a web-based toolset for functional profiling of gene lists from large-scale experiments.

    Nucleic Acids Res 2007, 35:W193-200. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  66. Mijakovic I: Protein phosphorylation in bacteria.

    Microbe 5:21-25. OpenURL

  67. Mok J, Kim PM, Lam HYK, Piccirillo S, Zhou X, Jeschke GR, Sheridan DL, et al.: Deciphering protein kinase specificity through large-scale analysis of yeast phosphorylation site motifs.

    Science Signaling 2010, 3:ra12. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  68. Grantham R: Amino acid difference formula to help explain protein evolution.

    Science 1974, 185:862-864. PubMed Abstract | Publisher Full Text OpenURL

  69. Akashi H, Gojobori T: Metabolic efficiency and amino acid composition in the proteomes of Escherichia coli and Bacillus subtilis.

    Proc Natl Acad Sci USA 2002, 99:3695-3700. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  70. Siltberg-Liberles J, Grahnen JA, Liberles DA: The evolution of protein structures and structural ensembles under functional constraint.

    Genes 2011, 2:748-762. Publisher Full Text OpenURL

  71. Raiford DW, Heizer EM Jr, Miller RV, Akashi H, Raymer ML, Krane DE: Do amino acid biosynthetic costs constrain protein evolution in Saccharomyces cerevisiae?

    J Mol Evol 2008, 67:621-630. PubMed Abstract | Publisher Full Text OpenURL

  72. Jha AN, Vishveshwara S, Banavar JR: Amino acid interaction preferences in proteins.

    Protein Sci 2010, 19:603-616. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  73. Barton MD, Delneri D, Oliver SG, Rattray M, Bergman CM: Evolutionary systems biology of amino acid biosynthetic cost in yeast.

    PLoS ONE 2010, 5:e11935. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  74. Tan CS, Pasculescu A, Lim WA, Pawson T, Bader GD, Linding R: Positive selection of tyrosine loss in metazoan evolution.

    Science 2009, 325:1686-1688. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  75. Møller IM, Rogowska-Wrzesinska A, Rao RS: Protein carbonylation and metal-catalyzed protein oxidation in a cellular perspective.

    J Proteomics 2011, 74:2228-2242. PubMed Abstract | Publisher Full Text OpenURL

  76. Xu K, Uversky VN, Xue B: Local flexibility facilitates oxidization of buried methionine residues.

    Protein Pept Lett 2012, 19:688-697. PubMed Abstract | Publisher Full Text OpenURL

  77. Woodsmith J, Kamburov A, Stelzl U: Dual coordination of post translational modifications in human protein networks.

    PLoS Comput Biol 2013, 9:e1002933. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  78. Wellen KE, Thompson CB: A two-way street: Reciprocal regulation of metabolism and signalling.

    Nature Rev Mol Cell Biol 2012, 13:270-276. OpenURL

  79. Yang X-J, Seto E: Lysine acetylation: Codified crosstalk with other posttranslational modifications.

    Mol Cell 2008, 31:449-461. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  80. Oien DB, Carrasco GA, Moskovitz J: Decreased phosphorylation and increased methionine oxidation of α-synuclein in the methionine sulfoxide reductase A knockout mouse.

    J Amino Acids 2011, 2011:721094. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  81. Oien DB, Shinogle HE, Moore DS, Moskovitz J: Clearance and phosphorylation of alpha-synuclein are inhibited in methionine sulfoxide reductase A null yeast cells.

    J Mol Neurosci 2009, 39:323-332. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  82. Snijder J, Rose RJ, Raijmakers R, Heck AJ: Site-specific methionine oxidation in calmodulin affects structural integrity and interaction with Ca2+/calmodulin-dependent protein kinase II.

    J Struct Biol 174:187-195. OpenURL

  83. Xu S, Ding H, Su F, Zhang A, Jiang M, Ostman A, Frijhoff J, et al.: Regulation of protein tyrosine phosphatases by reversible oxidation.

    J Biochem 2011, 150:345-356. PubMed Abstract | Publisher Full Text OpenURL

  84. Marondedze C, Turek I, Parrott B, Thomas L, Jankovic B, Lilley KS, Gehring C: Structural and functional characteristics of cGMP-dependent methionine oxidation in Arabidopsis thaliana proteins.

    Cell Commun Signal 2013, 11:1. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL