|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
,
,¶
From the || Max Planck Institute of Molecular Plant Physiology, Am Mühlenberg 1, 14476 Golm, Germany,
Carnegie Institution, Stanford, California 94305, and ¶ Department of Plant Biology, University of Copenhagen, Thorvaldsenvej 40, DK-1871 Frederiksberg, Denmark
| ABSTRACT |
|---|
|
|
|---|
Similarly in plants, glucose and sucrose are important metabolites and function as signaling molecules, potentially being involved in the regulation of growth and development at the whole plant level (6–9). Although the transcription levels of many genes respond to sugars in general and changes of transcription in response to sucrose depletion and resupply have been globally described (10) with over 100 transcriptional changes after 30 min of sucrose treatment, only a few genes that are regulated by sucrose specifically have been characterized (11). Examples for sucrose-specific responses are the transcription factor ATB2, which is repressed by sucrose at physiological concentrations (12), and the sucrose transporter SUT1 from sugar beet that is regulated by sucrose at the transcriptional level (13). The regulation of sucrose transporter transcription is thought to be mediated via a protein phosphorylation cascade (14). Because hexokinase has been shown to function as a key regulator of glucose-specific responses in plants (15, 16), the question remains how specificity in characteristic sucrose responses may be achieved. Similar to yeast, one may suspect that extracellular sucrose may be directly sensed at the plasma membrane affecting sucrose-specific transport and metabolism. A plasma membrane-localized sucrose sensing mechanism involving sucrose transporters has been suggested in the past (17, 18). Even if the sucrose sensing mechanisms were cytosolic, the responses are likely to involve modulation of sucrose transport at the plasma membrane.
Phosphorylation is the most well known post-translational modification (PTM)1 involved in signaling. This principle of activation and inactivation of proteins by phosphorylation as well as the function of phosphorylated residues as docking sites for protein scaffolds and complex assembly has been well studied in the field of mammalian signal transduction (19–22). However, most approaches in plant biology so far have been focused on phosphorylation of specific proteins and protein families (23, 24) and the study of specific signaling pathways (25, 26). There have been only a few efforts to globally analyze phosphorylation sites of membrane proteins (27). Although today a vast collection of techniques is available to specifically enrich for and detect phosphorylation sites in plant proteins (28–30), this study is the first that specifically follows phosphorylation events in a plant over time in response to an external stimulus. Such time course studies are of great value for the understanding of phosphorylation during plant signal transduction because they may reveal some important insights into specificity and relevance of individual phosphorylation sites for that given response (31, 32).
In this study, we took a proteomics approach to identify plasma membrane proteins that are phosphorylated in response to sucrose resupply after depletion to Arabidopsis seedlings. By quantitative comparison of the relative intensity of phosphoproteins at different times we present first insights into sucrose-induced phosphorylation responses of Arabidopsis membrane proteins 3, 5, 10, and 30 min after sucrose resupply. We identified 67 phosphorylation sites of which 70% have not been described before, and we showed sucrose-induced changes in phosphorylation level for 40 of these sites and established in vivo roles for some of the responses. The importance of this approach for understanding transporter regulation is demonstrated for a novel phosphorylation site in the C terminus of the AHA proton ATPase, which can override regulation by the well established 14-3-3 autoregulation of activity (33). Importantly receptor kinases also exhibited a change in their phosphorylation status within the first 3 min, suggesting a role close to the core of the signal perception.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
Plasma Membrane Preparation and Phosphopeptide Enrichment—
Plasma membranes of seedlings were isolated using a two-phase partitioning system as described previously (27). Briefly plasma membrane vesicles were inverted using Brij-58, and intracellular protein parts were digested using trypsin. Protein amounts were equalized before tryptic digestion. Protein amounts for each series ranged from 50 to 150 µg between repeated experiments.
Digested peptides were resuspended in 25% acetonitrile and 0.1% TFA and incubated with IMAC resin for phosphopeptide enrichment (PhosSelect, Sigma-Aldrich). Phosphopeptides were eluted from IMAC resin using 400 mM NH4OH in 25% acetonitrile. Samples were desalted over C18 material prior to mass spectrometric analysis (34).
Mass Spectrometric Analysis—
Tryptic peptide mixtures were analyzed by LC/MS/MS using nanoflow HPLC (Proxeon Biosystems) and a linear ion trap instrument (LTQ, Thermo Electron) as mass analyzer. Peptides were eluted from a 75-µm analytical column (Reprosil C18, Dr. Maisch GmbH) on a linear gradient running from 10 to 30% acetonitrile in 50 min and sprayed directly into the LTQ mass spectrometer. Proteins were identified by tandem mass spectrometry (MS/MS) by information-dependent acquisition of fragmentation spectra of multiply charged peptides. Zoom scans were implemented after a scan over the full mass range (300–1500 m/z) to be able to determine the charge state of the peptides. Additional data-dependent fragmentation (MS3) was used if peptides displayed a loss of phosphoric acid (neutral loss, 98 Da) upon MS/MS fragmentation (35). Fragment MS/MS spectra from raw files were extracted as DTA files and then merged to peak lists using the default settings of DTASuperCharge version 1.17 (SourceForge) with a tolerance of 500 ppm for precursor ion detection. With the default settings, the m/z of the respective zoom scans was used as a precursor m/z value in the DTA files and peak list. By using the option "Precursor mass N–2," the m/z value from the respective zoom scan was assigned to the DTA file and peak lists also for the MS3 fragment spectra.
Fragmentation spectra were searched against a non-redundant Arabidopsis protein database (TIGR6, version 2006-09; 30,711 entries; www.arabidopsis.org) using the Mascot algorithm (version 2.1; Matrix Science). The database contains the full Arabidopsis proteome and commonly observed contaminants (human keratin, trypsin, and lysyl endopeptidase) so no taxonomic restrictions were applied during automated database search. The following search parameters were applied: trypsin as cleaving enzyme; peptide mass tolerance, 300 ppm; MS/MS tolerance, 0.8 Da; one missed cleavage allowed. Carbamidomethylation of cysteine was set as a fixed modification, and methionine oxidation and phosphorylation of serine, threonine, and tyrosine were set as variable modifications. Only peptides with a length of more than five amino acids were considered. Putative phosphopeptides were manually inspected using the software package MSQuant version 1.4.1 (SourceForge). Peptides were accepted as phosphopeptides if they displayed a neutral loss for serine and threonine phosphorylation or if their fragmentation spectrum displayed b or y ions indicative of phosphorylation.
Phosphorylation site assignment was done using MSQuant version 1.4.1 as described previously (32). Briefly for each peptide different combinations of phosphorylation sites were scored (PTM score), and the highest scoring match was accepted if the sum of the Mascot score and PTM score was higher than 30. If no distinction could be made on this basis, the phosphorylation site was marked as ambiguous.
In general, peptides were accepted without manual interpretation if they displayed a Mascot score greater than 40; peptides with a score greater than 24 were manually inspected, requiring a series of three y or b ions to be accepted. Using the above criteria for protein identification, the rate of false identifications as determined by a search against a reversed TIGR6 Arabidopsis database (decoy database) was 1.12%. This was calculated as the percentage of twice the number of "hits" from the decoy database to the total number of hits in both searches (36).
Quantitative Analysis—
Within one time series, label-free quantitative analysis of the changes in phosphopeptide intensities was achieved by protein correlation profiling (37). Extracted ion chromatograms of each LC/MS/MS analysis representing different times of sucrose resupply were compared based on their m/z values and retention times. Ion intensity sums of m/z values identifying the same peptide sequence at the same retention time in different liquid chromatography runs were used for quantitative comparison between the different samples representing different time points.
Label-free relative quantitation of phosphorylation was performed as described previously (38). Briefly for ion intensities of each non-phosphopeptide the deviation from their mean ion intensity across the five time points was calculated. This was performed separately for each peptide species, i.e. for each m/z. Medians of all relative deviations for peptide species (i.e. for each m/z) from each protein were then calculated, and outlier peptides were identified based on their significant deviation from the median (
2 test). The outlier peptides were excluded from further analysis (see Supplemental Fig. 1). The mean of the remaining non-phosphopeptides was calculated for each time point and used for normalization. Ion intensities of phosphopeptide species (different m/z, i.e. different charge states treated independently) were normalized against the mean of non-phosphopeptides from the same protein. Subsequently for each phosphopeptide sequence, the mean of normalized intensities was calculated based on the different phosphopeptide m/z species, and this value was used to calculate ratios between treated (time points 3, 5, 10, and 30 min) and untreated (time point 0) samples. For each protein only unique non-phosphopeptides were used for normalization.
Peptides conserved in multiple members of a protein family were identified using the "show subsets" option in Mascot, and the respective peptides present in multiple proteins were excluded from quantitative analysis. Also if ion intensities for a given peptide showed signal to noise ratios <2, the respective values were omitted from quantitative analysis. If no unphosphorylated peptides were identified for a given protein, the mean of all identified non-phosphopeptides was used for normalization. These proteins are marked by asterisks in Table II.
|
The time series experiment of sucrose depletion and resupply was repeated in three independent biological experiments, and two of the experiments were suitable for quantitative analysis. One additional experiment was carried out with carbon depletion and mannitol resupply. If not stated otherwise, mean log2 ratios ±S.D. from analytical replicates of the experiment with the most quantified phosphopeptides are presented. Reproducibility of the time course responses between two biological experiments was confirmed by linear regression analysis (r2 = 0.725; p < 0.01) of the log2 ratios of quantifiable phosphopeptides identified in both of the two independent experiments across all time points (Supplemental Fig. 2).
Cluster analysis of the responses was done using Ward's method of increase in the sum of squares as a measure using ClustanGraphics version 8. For a given classification, the increase in the sum of squares of all cases in that classification is the minimum. The clustering method is tolerant for missing values.
Constructs for Expression of Arabidopsis AHA2 in Saccharomyces cerevisiae—
A centromeric yeast expression vector (pMP1733) containing the PMA1 promoter was used for expression of different versions of the AHA2 gene. aha2T881A and aha2T947A mutant genes (33) were transferred to this expression vector by inserting XhoI-SpeI fragments from pMP720 and pMP832, respectively. The aha2T947D mutation was generated by insertion of a linker containing the mutation into a SacI-SpeI-digested vector containing the AHA2 gene (pMP1745). The aha2T881D mutant was generated by polymerase chain reaction-assisted cassette mutagenesis in which an ApaI-SacI fragment containing the desired mutation was inserted into ApaI-SacI-digested pMP1745. All mutations were confirmed by sequencing.
Yeast Complementation Tests—
S. cerevisiae strain RS-72 (Mat a; ade1-100 his4-519 leu2-3,312 pPMA1::pGAL1) was used for complementation tests (39). The experiment was replicated independently three times with cells from independent transformation events. Isolation of plasma membrane protein from the yeast strain RS-72 was performed essentially as described previously (33).
Measurement of Proton Transport by the ATPase—
Proton transport was assayed as described previously (40) by monitoring fluorescence quenching of 9-amino-6-chloro-2-methoxyacridine, a dye that upon protonation accumulates in an impermeant form inside vesicles. The initial decrease in fluorescence is therefore a direct effect of the amount of protons transported into the plasma membrane vesicles by the H+-ATPase. The reaction medium contained 20 mM MOPS-KOH, pH 7.0, 1 mM ATP, 40 mM K2SO4, 25 mM KNO3, and 60 nM valinomycin. The assay were performed in a volume of 2.4 ml using 8 µg of plasma membrane protein. The reaction was started by addition of 2 mM MgSO4. Protein concentrations were determined by the method of Bradford with the Bio-Rad protein assay reagent and
-globulin as a standard.
Overlay Assays and Western Blots—
SDS-PAGE, Western blotting, and 14-3-3 overlay assays were performed as described previously (41). Briefly equal amounts of plasma membrane protein were loaded on a regular SDS gel, blotted to a nitrocellulose membrane, and incubated with MRGSH6-tagged GF14-
14-3-3 protein. Bound 14-3-3 protein was subsequently detected immunologically using a primary anti-RGSH6 antibody (Qiagen) followed by incubation with a secondary anti-IgG antibody conjugated with alkaline phosphatase. For detection of H+-ATPase antibodies raised against either the N terminus (number 762) or C terminus (number 759) of plasma membrane H+-ATPase were used (1:5000) (42).
| RESULTS |
|---|
|
|
|---|
19%, which is in agreement with experimentally determined contamination of plasma membrane preparations using stable isotope labeling (43). The majority of proteins (36%) identified from Arabidopsis seedling plasma membrane preparations in all experiments were transport-related. These included seven isoforms of the plasma membrane proton ATPase family, five aquaporins, four peptide transporters, four ABC transporters, five monosaccharide transporters, one sucrose transporter, and two phosphate and sulfate transporters (Table I). The second largest functional category with 12% of identified proteins consisted of membrane-anchored proteins involved in cell wall architecture, such as glycosylphosphatidylinositol-anchored proteins. Roughly 10% of all identified peptides belonged to proteins with yet uncharacterized functions.
|
|
Activation of the H+-ATPase provides the driving force for proton-coupled transport processes (i.e. sucrose transport). Therefore, it is not surprising to find phosphopeptides of the most abundant ATPases AHA1 and AHA2 as highly responsive to sucrose treatment (Table II). The relative ion intensity of the phosphopeptides GLDIDTAGHHYpTV (where pT is phosphothreonine) and GLDIETPSHYpTV representing Thr-947 in AHA1 and AHA2 increased with increasing duration of the sucrose resupply and reached highest ratios after 30 min of sucrose resupply (Table II and Fig. 2A). The C-terminal 14-3-3 protein-binding site (peptides GLDIDTAGHHYpTV and GLDIETPSHYpTV) is well conserved, and this phosphorylation site is known to be essential for the regulation of proton pumps by interaction of a 14-3-3 protein with the phosphorylated threonine (41, 44, 45).
|
Another phosphopeptide for ATPases, pTLHGLQPK, identifies a previously undescribed phosphorylation site also in the C terminus of AHA1 and AHA2 at Thr-881. Because the peptide pTLHGLQPK is identical between AHA1 and AHA2, we cannot distinguish which of the two proton ATPase isoforms contributes more to the phosphorylation of pTLHGLQPK, and the observed time course of phosphorylation change is an overlay of the response in both proteins (Table II). The phosphorylation site pTLHGLQPK was clearly identified by neutral loss and subsequent MS3 fragmentation (Fig. 3A). By expressing site-directed mutants of this phosphorylation site in yeast, we were able to show that this site, Thr-881, is crucial for plasma membrane ATPase activity in AHA2. In a yeast mutant deficient in the endogenous plasma membrane H+-ATPase, pma1, expression of AHA2 to a limited degree complemented growth. When placed on medium containing galactose, the yeast endogenous H+-ATPase PMA1 was expressed. When grown on medium containing glucose as the single carbon source, plant H+-ATPase isoforms will be expressed only, and the growth of the yeast cells is dependent on the activity of the plant H+-ATPase. Thr-881 and Thr-947 were mutated into Asp and Ala to study the effect of a negative charge at these positions. Strikingly if Thr-881 was mutated to Asp, thus mimicking phosphorylation, growth was stimulated markedly and was comparable to growth of the pma1 mutant complemented with a mutant of AHA2 (aha2
92) completely devoid of its autoinhibitory C-terminal domain (Fig. 3B). Mutation of Thr-881 to Ala inhibited the growth rate of transformed cells. The fact that mutating Thr-881 to an Asp resulted in pump activation suggests that phosphorylation of Thr-881 directly produces an activated pump without the need for additional protein partners such as 14-3-3 proteins. The phosphorylation site at Thr-881 of AHA1 and AHA2 is homologous to the peptide pTLHGLQAPDAK representing Thr-889 in the C terminus of the plasma membrane ATPase AHA11, and this phosphorylation site is conserved in all 12 members of the plasma membrane ATPase family except for AHA10 and AHA12.
|
A plasma membrane Ca2+-transporting ATPase (ACA8; At5g57110) displayed a decrease of phosphorylation at Ser-22 (Table II and Fig. 4A) after 5 min. The N-terminal region of type II B calcium-ATPases (such as ACA8) is known to have regulatory functions through a calmodulin-binding site, an autoinhibitory domain covering also Ser-22, and through phosphorylation (47). Phosphorylation of Ser-45 in the N terminus of ACA2 was shown to reduce calmodulin stimulation of ACA2 (48).
|
Aquaporins—
Three phosphorylation sites were found for aquaporins (Table II and Fig. 4A). Aquaporin PIP2e was found to be phosphorylated in the N-terminal peptide TKDELTEEESLSGK and at Ser-282 (pSQLHELHA where pS is phosphoserine). A phosphorylation site of aquaporin PIP2b was found at Ser-278 (SLGpSFR). Phosphorylated Ser-278 and Ser-281 in PIP2b are among the two conserved serine residues that control pore gating. Upon dephosphorylation, the pore is closed (52). Interestingly Ser-278 appears to locate the C-terminal domain, which interacts with the neighboring subunit of the tetramer. This finding may suggest that PIP2 aquaporins use an allosteric trans-regulatory mechanism for rapid shutoff as recently found in AMT ammonium transporters (53). The phosphorylation site at Ser-282 in PIP2e is in a conserved position homologous to the regulatory Ser-281 of PIP2b and may well have regulatory function also in PIP2e. The time course of Ser-282 at PIP2e displays a significant increase of relative phosphorylation after 5 min and a strong decrease after 30 min (Fig. 4A), suggesting that addition of external sucrose leads to a rapid transient opening of the PIP2e pore and a subsequent adaptation response with pore closing. Possibly adjustments to osmotic changes are involved in these processes. The role of the N-terminal phosphorylation site of PIP2e remains unclear especially because it did not show a strong response to sucrose. Phosphorylation of the aquaporin PIP2b showed only slight response to sucrose resupply (Fig. 4A).
Protein Kinases—
Because the receptors and signaling cascades for the sugar responses are not known in plants, the identification of receptors as candidates for involvement in sucrose responses is of major interest. In total, 34 receptor-like kinases were identified with more than two tryptic peptides. Among these, six were found to be phosphorylated at a total of seven different phosphorylation sites (Table II). Two of the phosphorylation sites of the identified receptor-like kinases are novel sites, which have not been described previously (30). In addition, seven phosphorylation sites to five non-membrane protein kinases were identified. Among the receptor kinases, At5g10020 from the LRR III subfamily and At1g53730 and At3g14350 from the LRR V subfamily (54) showed a significant increase in relative phosphorylation after only 3 min at sites in the juxtamembrane region followed by a strong decrease already at 5 min of sucrose exposure. In contrast, a phosphopeptide from the C-terminal region of LRR III receptor kinase At3g02880 decreased significantly in relative phosphorylation after 30 min of sucrose resupply (Fig. 4B). Phosphopeptide ApSAEVLGK was conserved to the kinase domain of 10 different protein kinases (At3g02880, At5g53320, At5g16590, At4g23740, At3g08680, At5g05160, At3g17840, At5g58300, At1g48480, and At2g26730) and has also been identified previously (27).
One soluble protein kinase showed an increase in phosphorylation at a phosphorylation site in the kinase domain (peptide SYSTNLAFTPPEYLR; Fig. 4B) after 3 min of sucrose resupply but showed dephosphorylation at the N-terminal Thr-25 (peptide SNPDVTGLDEEGR; Fig. 4B) after 5 min. A very similar phosphopeptide identified from the kinase domain of a receptor-like cytoplasmic kinase II showed a more moderate response to sucrose treatment (peptide SYSTNLAYTPPEYLR; Fig. 4B). The protein kinase At4g35600 from the receptor-like cytoplasmic kinase IV family was phosphorylated in the kinase domain at Ser-117 after 5 min of sucrose resupply, and phosphorylation of the kinase At3g01490 increased after 10 min at a site located in the N-terminal region (Fig. 4B).
The two closely related phosphopeptides SYSTNLAFTPPEYLR and SYSTNLAYTPPEYLR of protein kinases At5g41260 and At4g35230, respectively, were clearly identified as different sequences based on their MS/MS fragmentation pattern (Fig. 5). Both peptides displayed a neutral loss of 98, indicating serine or threonine phosphorylation. However, the information of the fragmentation pattern did not allow differentiation between phosphorylation of the first or the second serine in each of the peptides.
|
Specificity of the Response to Sucrose—
To define sucrose specificity of the response, we carried out an experiment in which mannitol was resupplied to the sucrose-depleted seedlings. Interestingly the plasma membrane ATPase AHA2 showed no change in phosphorylation at the C-terminal peptide GLDIETPSHYpTV (Table II), whereas proton ATPase AHA1 did show increased phosphorylation of the C-terminal peptide GLDIDTAGHHYpTV (Table II). This indicates that there may indeed be specificity among the responses, and Thr-947 of AHA2 may more specifically be activated by sucrose. The phosphorylation site pSQLHELHA of aquaporin PIP2e (Table II) also showed an increase in relative phosphorylation after 5 min of resupply, indicating that this protein displays a general response to changes in osmotic conditions. The phosphorylation sites of the sucrose transporter and the three fast responding receptor kinases were not identified in the mannitol treatment.
| DISCUSSION |
|---|
|
|
|---|
Recently comparative phosphoprotein analysis has yielded large datasets of elicitor responses in cell cultures (60), but temporal analysis of phosphorylation sites at the level of single phosphorylated peptides has not yet been studied in detail in plants. A temporal phosphoproteomics analysis of pathogenesis-related responses has been attempted using enrichment of phosphoproteins and iTRAQ (isobaric tags for relative and absolute quantification) labeling, but no phosphorylation sites were identified in that study (61). In mammalian systems an approach involving SILAC (stable isotope labeling by amino acids in cell culture) and subsequent immunoprecipitation of phosphorylated proteins has revealed ordered phosphorylation of proteins according to their position in a signaling cascade but was still lacking precise information of the actual phosphorylation sites involved (31). The first analysis of site-specific phosphorylation dynamics was carried out in yeast to analyze the response to mating pheromone (35) and has been optimized in mammalian cells, identifying site-specific phosphorylation dynamics (32).
With respect to phosphorylation site identification, in this study rather stringent criteria for the acceptance of phosphorylation sites were applied, and less starting material for plasma membrane preparation was used than in previous large scale studies. Although the number of identified sites in our study is 5 times lower than what has been published previously for membrane preparations from cell cultures (27), 70% of the phosphorylation sites identified in this study are novel sites that had not been identified in previous work. The differences may be due to differences in tissue (cell culture versus seedlings) as well as to stress and condition-specific responses. However, there is no bias toward certain classes of proteins identified in our study versus previous work.
Because metabolic labeling of whole plant seedlings is a rather lengthy task over at least two plant generations, we chose a label-free approach for relative comparison of ion intensities between samples that was used successfully before to characterize the proteome of the human centrosome (62). The variation of retention times between the five different LC runs was around 1%, which still allowed efficient peak assignment based on m/z values and their retention times between different samples. The average relative error of quantitation per protein as calculated from the mean ratios of analytical replicates ranged from 6 to 20%, which is slightly higher than what has been noted for metabolic stable isotope labeling approaches (63). We were able to (i) reproduce the response of the most common peptides in independent biological experiments (Supplemental Fig. 2) and (ii) validate the result of the label-free quantitation in the case of the plasma membrane ATPase in which the relationship between C-terminal phosphorylation and protein activity is well described (59). By label-free quantitation we found an increase in phosphorylation of C-terminal peptides GLDIETPSHYpTV and GLDIDTAGHHYpTV, and this clearly was accompanied by an increase in proton pumping activities in sucrose-supplied plasma membrane vesicles. Thus, the quantitation method indeed yields information relevant for the biology of the system under investigation.
In this study, we clearly identified 67 phosphorylation sites. Among the identified phosphorylation sites, we also confirmed well known sites, such as the C-terminal threonine of plasma membrane ATPases or known regulatory sites of aquaporins. In those cases, the time course of the relative intensities of the respective phosphopeptides is in good agreement with known biochemical and structural properties of these proteins or their activation or inactivation status. In addition, we clearly established an in vivo relevance for one of the new phosphorylation sites of AHA2 in controlling activity of the proton ATPase. Therefore, we are convinced that temporal analysis of phosphorylation sites using a label-free protein correlation approach is a suitable strategy to gain valuable insights into the regulatory processes involved in the specific response under investigation in tissues and systems not suited for metabolic labeling. Moreover the conserved phosphorylation site Ser-278 in the C terminus of PIP2e was found to have increased phosphorylation, suggesting that addition of sucrose leads to activation potentially via a trans-regulation within the complex as recently found for the ammonium transporter AMT1;1 (53). The novel sucrose transporter phosphorylation site may play an important role in acclimating sucrose uptake to the extracellular supply. Further work using site-directed mutagenesis is required to determine the role of this site in regulation of sucrose transport.
The Arabidopsis proteome is far from being annotated completely, and a large fraction of proteins still have no function assigned. In that respect, studies of the dynamic behavior and the analysis of phosphorylation under different conditions may provide a valuable contribution to the functional categorization of yet poorly annotated proteins. For example, two "unknown" proteins (At1g15400 and At5g57110; Table II) were identified in this study that showed temporal behavior of phosphorylation similar to that of the plasma membrane Ca2+-transporting ATPase ACA8. Similarly the uncharacterized protein kinase At5g41260 (Fig. 4B) could be involved in sucrose-related or osmotic responses.
Especially the changes in relative phosphorylation of kinases may reveal interesting information about signaling pathways. For example, in this study three membrane-bound kinases were identified that showed a rapid transient change in relative phosphorylation, whereas the identified soluble kinases showed peaks in relative phosphorylation at later time points. It certainly is too early to draw a direct connection or a signaling cascade, but further and more detailed similar phosphorylation dynamics studies under different nutrient conditions may help to identify such connections. In combination with bioinformatics analysis tools about domain structure (e.g. SMART (Simple Modular Architecture Research Tool), European Molecular Biology Laboratory) or short, characteristic sequence motifs (e.g. ELM (Eukaryotic Linear Motif), the ELM Consortium) and with thorough clustering analysis of the identified responses we may gain more insight into the links between patterns of phosphorylation and possible function of unknown proteins. Clustering of the quantified phosphorylation sites identified in this study revealed four significant response classes (Fig. 6). Interestingly the protein distribution between these classes shows that the kinases identified in our study exclusively group in clusters 1 and 2, which on average show rapid increases or decreases of relative phosphorylation. In contrast, enzymes and other proteins are more represented in "slow" responding clusters 3 and 4. Clustering of the responses may therefore be an important tool to identify proteins involved in a signaling cascade in a time-resolved manner especially if this type of analysis is being extended over a wide range of conditions in a large scale fashion. Of course, detailed functional characterizations of candidate proteins will ultimately have to go hand in hand with the analysis of the respective T-DNA insertional mutants and will have to include transcriptomics and metabolomics profiling.
|
With this work we have made a first step in including a dynamic component to plant phosphoproteomics. Although in total the fraction of phosphorylation sites for which quantitative time course information could be resolved was rather low (65% of all identified sites), we believe that with more optimization of phosphopeptide enrichment, larger quantities of starting material, and alternative MS/MS fragmentation techniques the number can well be exceeded. Nevertheless the results presented in this first study of dynamic plant phosphoproteomics already indicate exciting new candidates possibly involved in the key process of sucrose transport and regulation. These will be validated in in-depth biological experiments in the future. However, it is important to point out that this unbiased proteomics survey study provided a terrific basis for novel targeted analyses in a field that has been under intense investigation for a long time in plant biology.
In conclusion, our work on dynamics of phosphorylation under sucrose depletion and resupply validates the method of protein correlation profiling as an alternative approach for quantitation in cases where metabolic labeling is not possible (i.e. tissue samples, whole plant, and seedlings). Furthermore we clearly identified several previously undescribed phosphorylation sites even in well studied proteins such as the plasma membrane proton ATPases, and most importantly, we give first insights into the dynamic behavior of individual phosphorylation sites under specific conditions. We believe that similar experiments carried out under different stress conditions will provide a strong basis for our understanding of signaling processes in plants.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
Published, MCP Papers in Press, June 23, 2007, DOI 10.1074/mcp.M700164-MCP200
1 The abbreviations used are: PTM, post-translational modification; MS3, MS/MS/MS fragmentation; ANOVA, analysis of variance; ABC, ATP-binding cassette; LRR, leucine-rich repeat. ![]()
* This work has been supported by the Human Frontier Science Program Fellowship (to T. N.), Danish Research Council for Technology and Production Grant FTP-274-05-0269 (to M. G. P.), and Department of Energy Grant DE-FG02-04ER15542 (to W. B. F.). The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. ![]()
S The on-line version of this article (available at http://www.mcponline.org) contains supplemental material. ![]()
Both authors made equal contributions to this work. ![]()
** Supported by the Emmy-Noether Program of the Deutsche Forschungsgemeinschaft. To whom correspondence should be addressed. E-mail: wschulze{at}mpimp-golm.mpg.de
| REFERENCES |
|---|
|
|
|---|