Phosphoproteomics of the Arabidopsis Plasma Membrane and a New Phosphorylation Site Database[W]

Abstract

Functional genomic technologies are generating vast amounts of data describing the presence of transcripts or proteins in plant cells. Together with classical genetics, these approaches broaden our understanding of the gene products required for specific responses. Looking to the future, the focus of research must shift to the dynamic aspects of biology: molecular mechanisms of function and regulation. Phosphorylation is a key regulatory factor in all aspects of plant biology; but it is difficult, if not impossible, for most researchers to identify in vivo phosphorylation sites within their proteins of interest. We have developed a large-scale strategy for the isolation of phosphopeptides and identification by mass spectrometry (Nühse et al., 2003b). Here, we describe the identification of more than 300 phosphorylation sites from Arabidopsis thaliana plasma membrane proteins. These data will be a valuable resource for many fields of plant biology and overcome a major impediment to the elucidation of signal transduction pathways. We present an analysis of the characteristics of phosphorylation sites, their conservation among orthologs and paralogs, and the existence of putative motifs surrounding the sites. These analyses yield general principles for predicting other phosphorylation sites in plants and provide indications of specificity determinants for responsible kinases. In addition, more than 50 sites were mapped on receptor-like kinases and revealed an unexpected complexity of regulation. Finally, the data also provide empirical evidence on the topology of transmembrane proteins. This information indicates that prediction programs incorrectly identified the cytosolic portion of the protein in 25% of the transmembrane proteins found in this study. All data are deposited in a new searchable database for plant phosphorylation sites maintained by PlantsP (http://plantsp.sdsc.edu) that will be updated as the project expands to encompass additional tissues and organelles.

INTRODUCTION

After the completion of several eukaryotic genome sequences and with the automation of gene expression studies, biologists face the prospect that the next step toward understanding cellular processes is far more complex. The functions of proteins are regulated by localization, binding partners, posttranslational modifications, and stability. Most of these parameters are poorly predictable from their primary sequence—a postgenomic challenge. It is thus necessary to examine the protein itself, which can be a formidable task for regulatory proteins of low abundance. In the past few years, rapid progress in mass spectrometric technology has allowed analysis and identification of proteins with unprecedented levels of sensitivity, and numerous proteomic studies have tackled the aforementioned challenges. The protein complements of subcellular structures have been explored by organellar proteomics (Taylor et al., 2003). The interaction partners of hundreds of yeast proteins identified with affinity purification techniques (Gavin et al., 2002; Ho et al., 2002) and the sites of posttranslational modifications were determined using specific enrichment techniques (Ficarro et al., 2002; Elortza et al., 2003; Peng et al., 2003).

Protein phosphorylation is one of the most widespread and arguably best understood posttranslational modifications. Virtually all cellular processes are regulated in one or multiple ways by phosphorylation and dephosphorylation, and the identification of kinases, their substrates, and the specific site of phosphorylation is the key to a molecular understanding of signaling. In vivo labeling (Peck et al., 2001), phosphorylation-specific antibodies, and/or band shifts in immunoblots are all methods to show that a particular protein becomes phosphorylated in response to a stimulus. The identification of the actual site(s) of phosphorylation, however, is necessary for a functional characterization and represents a major bottleneck for analysis of individual proteins. With the described advances in mass spectrometric technology, we can now approach the problem from a different angle. Phosphopeptides can be isolated from complex protein digests and identified in a large scale by mass spectrometry in a shotgun approach (Ficarro et al., 2002). In a sufficiently large study, the phosphorylation sites of interest may be found directly. Alternatively, the identified sites from related proteins provide enough knowledge to guide further experiments.

We know from the Arabidopsis thaliana genome that plants have evolved signal transduction components that are substantially different from those of animals, in particular Ser/Thr rather than Tyr receptor kinases (Arabidopsis Genome Initiative, 2000). Uncritical transfer of mammalian signaling paradigms may therefore be misleading, and it will be necessary to establish mechanisms that are unique in plants. Two major obstacles, however, impede the analysis of integral membrane signaling proteins. First, the proteins are difficult to work with biochemically because of their hydrophobic nature and typically low abundance. Second, identification of phosphorylation sites requires an enrichment method for phosphopeptides because suppression effects during ionization prevent the detection of these peptides by mass spectrometry when in complex mixtures.

The problems of membrane protein insolubility are circumvented by proteolytic digestion of the intact membranes and analysis of peptides released from extramembrane domains (described in Nühse et al., 2003b). Because the proteins are digested in situ and only the fragments from intracellular loops and N and C termini recovered, no selection against large proteins or those with multiple transmembrane domains is apparent. This strategy, then, overcomes the first limitation.

For enrichment of phosphopeptides, immobilized metal-ion affinity chromatography (IMAC) for phosphopeptides showed great promise for large-scale studies but had a reputation for poor specificity. We investigated the potential of IMAC in combination with capillary liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) for the identification of plasma membrane phosphoproteins of Arabidopsis (Nühse et al., 2003b). Using either simple batch binding of peptides or prefractionation via strong anion exchange (SAX) chromatography before IMAC, a distinct set of phosphopeptides was detected. The degree of enrichment over nonphosphorylated peptides (75 to 90% purity) allowed efficient analysis and solved the second challenge.

Here, we present an analysis of several hundred phosphorylation sites on ∼200 putative plasma membrane proteins of Arabidopsis collected from multiple LC-MS/MS analyses of independent plasma membrane preparations. This phosphoproteome study is among the largest performed in any organism. Besides providing new insights into the functional regulation of numerous proteins, the size of the data set enabled us to find characteristics and patterns that form the basis for developing predictive models of phosphorylation in plants.

RESULTS

Overcoming Limitations in Plasma Membrane Phosphoproteome Analysis

Two benchmark parameters were used to test whether the methodology delivers comprehensive proteome coverage: the representation of proteins with multiple transmembrane-spanning domains and the sequence coverage of cytoplasmic domains that can be phosphorylated. Analyzing the data from this large-scale phosphoproteomic study, we found that approximately two-thirds of the identified phosphoproteins had at least one transmembrane (TM) domain (Figure 1A  

Figure 1.

Coverage and Functional Classification of PM Phosphoproteins.

(A) Proportional representation of transmembrane helices within proteins found in this study as predicted by TMHMM (http://www.cbs.dtu.dk/services/TMHMM/).

(B) Graphical representation of sequence coverage of AHA2 in a shotgun proteomic experiment with total PM peptides. Structural model adapted from Palmgren (2001).

(C) Representation of functional protein classes among PM phosphoproteins.

). One of them, a calpain-like protein, contained a predicted 23 TM domains. It has been suggested that precleavage of membrane proteins with cyanogen bromide is necessary for complete tryptic digestion (Washburn et al., 2001). In an analysis of a total plasma membrane (PM) digest, however, we found very good coverage of even short cytoplasmic loops for the proton ATPase AHA2 (Figure 1B). We postulate that membrane shaving of tryptic peptides is facilitated by the fact that TM domains are often anchored in the membrane by polybasic residues (i.e., trypsin cleavage sites) on the cytosolic face. Thus, the accessibility of PM proteins to trypsin does not appear to be a limiting factor. Moreover, all fragments were from regions predicted to reside on the cytoplasmic face of the PM (Figure 1B), indicating that this type of analysis provides information on the topography of proteins within the membrane (Wu et al., 2003).

The functional categories of the 200 identified phosphoproteins (Nühse et al., 2003b; this study) show a high representation of proteins with a role in signal transduction. Receptor-like kinases represent a quarter of all phosphoproteins and all signaling-related proteins, as far as could be annotated by similarity, more than a third. Proteins with a role in transport (channels and ABC transporters) and ion homeostasis (H⁺ and Ca²⁺ ATPases) constitute another large family. Although the plasma membranes were of high purity (Nühse et al., 2003a), a small number of proteins were identified as clear contaminants, among them several nuclear proteins, such as DNA polymerases and putative transcription factors.

Confirmed and Novel Phosphorylation Sites

Unsurprisingly, all but a few phosphorylation sites described here are novel. Among the few well-characterized plant plasma membrane phosphoproteins are the H⁺-ATPases and the aquaporins. We reported the identification of previously known and novel phosphorylation sites on several isoforms of the H⁺-ATPase (AHA1, AHA2, AHA3, and AHA4) in an earlier publication (Nühse et al., 2003b). Vacuolar and plasma membrane aquaporins are known to be regulated by posttranslational modification to adapt to rapid changes in the apoplastic water potential (Johansson et al., 1996; Maurel, 1997). The phosphorylation sites we identified on two PIP2 isoforms (see Supplemental Table 1 online) have also been detected in another proteomic study (Santoni et al., 2003). They are equivalent to Ser262 of soybean (Glycine max) nodulin26 and Ser274 of spinach (Spinacia oleracea) PM28A (Johansson et al., 1998; Guenther et al., 2003), both of which are phosphorylated in response to increasing water potential, probably by a calcium-dependent kinase.

In several cases, phosphorylation of a particular protein was known or predicted, but no data on the precise site of modification were published. The Na⁺/H⁺-antiporter SOS1 is phosphorylated by the SOS2 kinase in complex with the calcium sensor protein SOS3 on (an) unknown site(s) (Quintero et al., 2002). The site we identified is consistent with both the anticipated SNF1/AMPK1 motif and the localization on the hydrophilic C terminus (Halfter et al., 2000), and thus may be (one of) the authentic site(s) targeted by SOS2/SOS3.

For the remainder of the phosphorylation sites, however, this report represents the first insights into potential posttranslational regulation of these proteins. In addition to the importance in studying individual proteins, the data provide the basis of drawing some general predictions and conclusions about protein phosphorylation in plants.

Phosphorylation Sites Tend to Cluster Outside Known Functional Domains

Proteins are generally composed of catalytic/functional domains surrounded or interspersed by less defined (or not yet understood) variable regions. We analyzed our data to investigate if there is a trend in the location of phosphorylation sites in relation to the architecture of the protein using protein domain predictions from the SMART Web site (http://smart.embl-heidelberg.de/).

Figure 2  

Figure 2.

Graphical Representation of Protein Architecture, Sequence Conservation within the Protein Family, and Localization of Phosphorylation Sites.

Arrowheads indicate the location of phosphorylation sites within the protein. Protein domain architecture was adapted from the SMART Web site (http://smart.embl-heidelberg.de/). Sequence conservation was plotted from ClustalW alignments using MDR family members as listed in PlantsT (http://plantst.sdsc.edu/) and CDPK members from PlantsP (http://plantsp.sdsc.edu/). Arabidopsis CesA sequences were downloaded from http://cellwall.stanford.edu/php/download.php. ABC_membr., TM domain of ABC transporters; AAA, ATPases associated with various cellular activities; RING, RING-type zinc finger; EF, EF hand/calcium binding motif.

shows the position of phosphorylation sites relative to functional domains for three multidomain protein families represented in our study. In the majority of cases, the sites were outside catalytic and known regulatory domains, and sequences surrounding the sites were poorly conserved within a gene family. Nevertheless, such isoform-specific regulatory sites can still be clustered in equivalent regions of the protein. In the case of the multiple drug resistance (MDR) family of ABC transporters, the poorly conserved linker region between the duplicated ATPase domains is a well-known target for phosphorylation in yeast and human transporters (Decottignies et al., 1999; Idriss et al., 2000). Three sites in plant ABC transporters share an [R,K]-X-pS motif (see Supplemental Table 1 online), suggesting that these sites are phosphorylated by the same or a related kinase. In the case of the cellulose synthase family, all but one of the sites are in the unconserved region upstream of the first TM domain, whereas one phosphorylation motif in the catalytic domain of CesA1 is loosely conserved among several other CesA proteins. Among the calcium-dependent protein kinases, we found two sites within the C-terminal calmodulin-like domains (where in one case the phosphoserine replaces a conserved Asp, see below), two in the kinase domain, and only two in the highly diverse N termini. Taking into account the data from receptor-like kinases (see below) and other families (data not shown), the majority of phosphorylation sites are located outside of defined catalytic/functional domains.

An additional outcome of this analysis was that phosphorylation sites in most of the multispanning TM proteins were close to either the N or C terminus. Assuming that phosphorylation is exclusively intracellular, the identification of a phosphopeptide establishes the corresponding terminus as cytoplasmic. This knowledge is useful for the determination of the inside/outside topology of unknown membrane proteins. We found a Web-based prediction tool (http://www.cbs.dtu.dk/services/TMHMM/) in conflict with our empirical data for more than one-quarter of proteins with more than one TM domain (data not shown). Thus, these data have significance on a structural as well as functional level.

Conservation of Phosphorylation Sites between Related Proteins

Related to the question addressed above, we also investigate if the amino acid sequence surrounding a phosphorylation site is conserved between members of functional families of proteins. It could be envisaged that sites involved in the direct regulation of enzyme activity might be more likely to occur within conserved regions of a protein family. On the other hand, phosphorylation of unique sites would allow isoform-specific regulation. In the latter case, these sites may still be conserved among true orthologs in related plant species.

Figure 3  

Figure 3.

Comparison of Phosphorylation Sites Conserved in Paralogs and Orthologs of Cellulose Synthases.

Three regions surrounding CesA3 phosphorylation sites (marked with arrows) have been aligned with all paralogous CesA sequences from Arabidopsis (top panel) and with putative orthologs in other species. The orthologs were found by BLAST search in EST databases using 40 to 60 amino acids from CesA3. The location of the short sequences relative to the functional domains is indicated above the alignments. Identical residues are in black boxes, and similar residues are in gray.

shows an alignment of the N-terminal phosphorylation sites of the cellulose synthase CesA3 with the corresponding sequences of paralogous CesA proteins (top panel) and putative orthologs (bottom panel), the latter found by BLAST search of EST databases. It is striking that whereas the CesA3 sites are unique among the paralogous CesA proteins, highly similar sequences were found in putative orthologs from the EST database, indicating that isoform-specific regulation of this cellulose synthase protein is conserved even in monocots.

To get a more representative overview of sequence conservation around phosphorylation sites, we did a similar comparison for all putative transporter proteins. Putative paralogs were taken from the Aramemnon database (http://aramemnon.botanik.uni-koeln.de/). Putative orthologs were found by BLAST search of ∼40 to 60 amino acids surrounding the phosphorylation site against an EST database, including a few amino acids of sequence that are conserved in the protein family when possible to reduce random low-scoring hits. As shown in Table 1  

Table 1.

Conservation of Phosphorylation Sites in Membrane Transporters

PlantsT	Family	Subfamily/Name	Protein	Phosphorylation Site	Conservation within Protein Family	Conservation among Putative Orthologs (EST Database)
2.A.1	STP	Monosaccharide transporter	STP1	RGVDDVpSQEFDDL	–	(+)
			STP13	FVNGEKpSNGKSNG	–	+
	PHT	Phosphate transporter	PHT1-4	NEDNENpSNNDSRT	–	n.f.
2.A.2	SUC	Sucrose-proton symporter	SUC5	AANNApTALEpTQpSpSPEDLGQ	3 (9)	+
2.A.7	PUP	Purine permease	PUP18	EPDQILpSPRRSLE	2 (13)	+
2.A.16	TDT	C4-dicarboxylate/malic acid transporter (telluride resistance)	At4g13800	PVFINpSGpSSRSSNS	2 (7)	n.f.
			At5g24030	DQLRNVpSpSENIENYLK	–	n.f.
2.A.17	POT	Proton-dependent oligopeptide transporter	At3g47960	NTDVVDpSFEEEQR	–	+
			At1g69870	GSFSKSpSPpSELDV..YKRIpSpSPGSILD	–	(+)
2.A.18	AAAT	amino acid/auxin permease	AtLTH4	PATPRVpSpTPEILTP	–	+
			At5g41800	PVTRLDpSDAGALF	–	(+)
			At5g01240	GRKVEDpSAAEEDI	–	–
2.A.29		Mitochondrial carrier	At3g55640	KFMYMVpTGMENHK	–	+
2.A.31	AE	Anion exchangers/boron transporter	At3g06450	EAPAILpSFNLKPE	–	+
2.A.36	NHX	Sodium-proton exchanger	SOS1/NHX7	VNVYRRpSVSFGGI	–	+
			NHX1,2	VPFVPGpSPpTERNPPD	–	+
2.A.40	NCS2	Nucleobase:cation symporter 2	At2g27810	GFRPKFpSGEpTTApTDpSpSSGQLSL	–	+
2.A.49	AMT	Ammonium transporter	AMT1	MAGMDMpTRHGGFA	3 (6)	+
2.A.55	NRAMP	Metal ion (Fe/Mn) transporter	NRAMP1	QLPCRVpSTSDVD	–	+
2.A.57	ENT	Equilibrative nucleoside transporter	ENT6	AGIQNQpSDLpSDDDSKN	4 (7)	+
2.A.72	KUP	K+ uptake permease	KUP2	EDDNARpSVQpSNEpSSSESR	–	(+)
			KUP11	DQKLDQpSMDEEAG	–	–

The conservation of a phosphorylatable residue (S or T) in paralogs and putative orthologs of phosphorylated transporters was analyzed. Phosphoproteins were grouped into families according to PlantsT (http://plantst.sdsc.edu/). Proteins listed in the Aramemnon database (http://aramemnon.botanik.uni-koeln.de/) as isospecific homologs of the phosphoprotein with at least 30% similarity were called protein family. Under the heading Conservation within Protein Family, the number of family members that share the phosphorylation site is listed. Under the heading Conservation among Putative Orthologs, (+) denotes one BLAST hit in the EST database with conserved phosphorylation site; + denotes two or more hits; n.f., no ambiguous ortholog was found in the EST databases.

, the majority of phosphorylation sites are conserved among putative orthologs and, to a lesser extent, among some members of the same protein family. Because of potential errors in EST sequence data, a negative result cannot be taken as evidence for absence of related phosphorylation sites. Conversely, however, the conservation of a putative phosphorylation motif may be taken as additional supporting evidence for the exact phosphorylation site if the mass spectrometric data are ambiguous. The SOS1 phosphorylation site could not be unambiguously assigned to S1 or S3 in the peptide SVSFGGIYNNK, but S3 is more highly conserved in other species. Thus, it is the more likely candidate (Figure 4A  

Figure 4.

Phosphorylation of Ser Residues May Act as a Reversible Mimic of Functionally Important Acidic Residues.

(A) ClustalW alignment of the experimentally determined phosphorylation site on the Na⁺/H⁺ transporter SOS1 with putative orthologous sequences. Potential phosphorylation sites in SOS1 are indicated by arrows.

(B) and (C) Alignment of sequences corresponding to phosphorylation sites (white on black background) in the glucose transporter STP1 and the calcium-dependent protein kinase CPK32, respectively. Acidic residues in the position equivalent to experimentally determined phosphorylation sites have a gray background.

). In addition, careful analysis of aligned sequences may also help to define specificity determinants of the kinase that targets the sequence. In SOS1, the Phe residue in the +1 position is also conserved and may be important in addition to the known basic-X-X-Ser motif for the SNF1/AMP-related kinase SOS2 (Weekes et al., 1993; Qiu et al., 2002).

In the glucose transporter STP1 (Figure 4B) and, most strikingly, the calcium-dependent protein kinase CPK32 (Figure 4C), we found the interesting case of a conserved acidic residue being replaced by a phosphoamino acid in one or two sequences. In CPK32, the phosphoserine takes the place of an absolutely conserved Asp that coordinates the calcium ion in the EF hand. It is tempting to speculate that phosphorylation of this Ser residue represents the switching on of a charge-dependent function in one isoform where the others are constitutively on.

Phosphorylation of Receptor-Like Kinases Implies Complex Regulation

Receptor-like Ser/Thr kinases (RLKs) form an extremely diverse protein family in Arabidopsis. Because the genome sequence has revealed their number to be in the hundreds, the question has arisen of how signaling specificity is ensured. Even if only a fraction of all RLK genes is expressed at any given time in one cell, it is puzzling how unique output signals are generated from receptors with highly similar kinase domains.

Our experimental data on RLK phosphorylation sites are perhaps the most striking of this study. With ∼50 identified phosphoproteins, this protein family is by far the best represented and accounts for a quarter of all proteins (Figure 1C). Figure 5A  

Figure 5.

Phosphorylation Sites of Receptor-Like Kinases.

A ClustalW alignment of phosphorylated receptor-like kinases was rendered with Boxshade. Experimentally determined phosphorylation sites are represented by white squares with black frames. A schematic drawing of the secondary structure (below) shows conserved α helices and β sheets of the kinase domain as cylinders and arrows, respectively. An enlargement of this sequence alignment is provided in Supplemental Figure 1 online.

shows a graphic representation of the position of the phosphorylation sites relative to the conserved kinase subdomains (for a full-size alignment, see Supplemental Figure 1 online). Unexpectedly, three-quarters of the phosphopeptides came from either the juxtamembrane region ([JM]; the linker between the TM domain and the kinase core) or the domain C-terminal of the kinase core. The JM and C-terminal sequences are highly divergent in length and sequence, and the identified phosphorylation sites are generally unique for a single RLK (data not shown). The implications for signaling mechanisms are discussed below.

Only eight phosphopeptides were found in the activation segment or A loop (i.e., the sequence between the conserved DFG and APE motifs of the subdomains VII and VIII, respectively), the best-characterized domain for kinase activation. Several of them cover the same conserved residues, and one of them (GpSFGYLD…) corresponds to one of the very few previously known RLK phosphorylation sites, Thr468 in SERKa (Shah et al., 2001). Most, but not all, RLKs with phosphorylated A loops belong to the group of RD kinases that require this modification for charge neutralization in the active site (Johnson et al., 1996).

Outside of the activation segment, some more unusual phosphorylation sites were found. Several peptides cover an [RK]-A-S-A-E motif just upstream of domain I, which is highly conserved among the LRRK III and VII families and probably represents a common regulatory site. Modification of residues in the unconserved linker regions between the conserved subdomains I/II, II/III, and IX/X, respectively, or within domains X/XI (Figure 5) have to our knowledge not been reported for any protein kinase. This certainly reflects a bias in hypothesis-driven research toward previously known sites (particularly in the A loop), but these unusual sites may indeed prove to be unique for plant RLKs or even the particular protein itself.

Phosphorylation Site Motifs

Although kinases may have additional protein–protein docking or interaction domains, the enzymatic active site shows specificity toward its substrates based on the linear sequence surrounding the residue to be phosphorylated. We sorted all unambiguously assigned phosphorylation sites according to single residues in the −6 to +6 positions to detect the frequency of residues within specific positions and to look for common motifs (Figure 6A  

Figure 6.

Frequency Distribution of Amino Acid Residues Surrounding Phosphorylation Sites.

(A) Frequency of amino acids in single positions relative to the total frequency in positions −6 to +6.

(B) Alignment of phosphorylation site sequences (residues −6 to +6) with the experimentally found G-pS-[FYW] motif. The color scheme is simplified from the scheme in (Taylor, 1997).

). Although most residues were evenly distributed across all positions, Pro was seldom in the −1 position and overrepresented in the +1. Aromatic residues were also found with disproportionate frequency in the +1 position but were virtually excluded from position +2. Taking these analyses further, we found restraint on the residues in the −1 position when an aromatic residue was in the +1 position. Overall, ∼12% of total Gly residues are present in the −1 position. However, when an aromatic residue is present in the +1 position, Gly is the −1 residue 50% of the time (13/26; Figure 6B). Conversely, Pro or aromatic residues are never present at the −1 position with an aromatic residue in the +1 position. These indications of bias probably reflect spatial restraints of the substrate binding site of the relevant kinases. Information from these types of analyses will be useful both for designing peptide libraries to test kinase specificity determinants and for refining programs for predicting phosphorylation sites.

Finally, we investigated the possibility of using patterns of residues surrounding the phosphorylation site to predict kinase specificity determinants or to cluster proteins that may be coordinately regulated. We found an unusually large number of motifs with acidic residues in the −1 to +4 positions (see Supplemental Figure 2D online), many even with a [DE]-pS-[DE] motif (see Supplemental Figure 2E online). The second largest group were protein kinase A-like motifs [RK]-x₁₋₃-pS (see Supplemental Figures 2A and 2B online). The high frequency of these motifs reflects the fact that acidic- or basic-directed specificity is common among the well-characterized kinases. The G – pS – aromatic motif discussed above is conserved among most PIP2 (plasma membrane intrinsic protein) and all Syp1 (syntaxin family 1) proteins (Figure 6A). Interestingly, phosphorylation of the homologous site of the aquaporin PM28A has been shown to be calcium dependent (Weaver and Roberts, 1992), and uncharacterized calcium-dependent activities also phosphorylated a H⁺-ATPase (Schaller and Sussman, 1988) and the syntaxin Syp122 (Nühse et al., 2003a). It is conceivable that this unusual phosphorylation motif is subject to regulation by the same or closely related kinases. If these kinases are indeed CDPKs, the substrate motifs are more variable than anticipated (Cheng et al., 2002). Several motifs have been described as preferred substrates for the related Snf1-like and calcium-dependent protein kinases (Weekes et al., 1993; McMichael et al., 1995; Bachmann et al., 1996; Huber and Huber, 1996; Huang and Huber, 2001; Huang et al., 2001). Many of the experimentally determined sites correspond to the simple pattern basic-X₂-pSer, but only very few match the ideal motif determined in peptide-based assays (see Supplemental Figure 2 online). It is possible that the plasma membrane samples contain only very few substrates of Snf1-like kinases, such as metabolic enzymes, but a more likely explanation is that full-length proteins in vivo reveal a more complex recognition motif than a short stretch of primary structure. Thus, the amino acid positions described here for sorting may not actually be the sole specificity determinant for the respective kinase. It is likely, however, that the inclusion of structural and other parameters as well as empirical data on kinase specificity will allow us to use these data sets for predictions as was successfully done for mammalian proteins (Blom et al., 1999).

DISCUSSION

With the set of Arabidopsis plasma membrane phosphorylation sites presented in this work, we have made two major achievements. First, as a database resource, the data will facilitate signaling research on the proteins identified. Technically, it remains difficult, if not impossible, for most laboratories to identify phosphorylation sites within their proteins of interest. These data, and those that follow, should accelerate studies into the mechanistic regulation of biological processes in plants in a manner comparable to the way genome and EST sequencing projects have facilitated the cloning of genes of interest. All data will be released as a free access, searchable database through PlantsP (http://plantsp.sdsc.edu). Second, in a wider scope, the data provide insight into general emerging patterns of regulation through phosphorylation as well as an appreciation of the complexity of these processes.

It is important to stress that the identification of phosphorylation sites as described in this work relies on MS/MS spectra from a single peptide. Thus, the analysis requires a sequenced genome to ensure a high-confidence assignment and currently restricts such large-scale plant studies to Arabidopsis and rice (Oryza sativa). We found, however, that phosphorylation sites identified in Arabidopsis generally are conserved in a broad range of species (Table 1, Figures 3 and 4; data not shown). These results indicate that the evolutionary conservation of regulatory mechanisms allows for the direct transfer of phosphorylation site data from model species to crop species.

Concerning individual phosphorylation sites, our key findings are posttranslational regulation of many novel proteins and protein families and a comprehensive insight into the complexity of receptor-like kinase phosphorylation. For most of the identified protein families, such as receptor-like kinases, ABC transporters, and ion channels, posttranslational regulation in yeast and mammals is well characterized. Therefore, finding phosphorylation sites on analogous Arabidopsis proteins was not surprising. For the family of cellulose synthases and the functionally associated cellulase KORRIGAN, however, regulation by phosphorylation was unanticipated. Understanding of the regulation of cellulose biosynthesis is still vague, but only roles for transcriptional regulation, protein stability, and redox potential have been shown (Doblin et al., 2002). Cellulose synthases, especially those involved in building the primary cell wall, have proven difficult to work with on a biochemical level. Mutational analysis of the phosphorylation sites described in this work may lead to new understanding of the mechanisms regulating these important proteins.

Following the paradigms of receptor Tyr kinases of animals (Schlessinger, 2000), the regulation of plant RLKs by phosphorylation comes as no surprise. Investigation of the catalytic properties of many receptor-like kinases have relied on bacterial expression of fusion proteins and most display strong autophosphorylation activity with varying specificity for Ser or Thr residue. Where these sites were analyzed, 10 or more phosphorylated tryptic peptides were found (Stone et al., 1998; Coello et al., 1999; Lease et al., 2001; Liu et al., 2002). When the mechanism of this autophosphorylation was analyzed, some kinases showed an intermolecular mechanism (Coello et al., 1999; van der Knaap et al., 1999), whereas BRI1, Xa21, and CrRLK phosphorylated mostly or exclusively via an intramolecular mechanism (Schulze-Muth et al., 1996; Oh et al., 2000; Liu et al., 2002). It is worth noting that experiments with kinase domain fusions and full-length proteins can come to remarkably different results. Li et al. (2002)found that Escherichia coli–expressed MBP and GST fusions of the BRI1 and BAK1 kinase domains, respectively, have autophosphorylation activity, whereas the full-length proteins expressed in yeast do not (auto)phosphorylate unless coexpressed (Nam and Li, 2002). Similarly, soluble kinase-dead BRI1 binds soluble BAK1 much less strongly than does wild-type BRI1 (Li et al., 2002), whereas no such effect was seen with the full-length proteins (Nam and Li, 2002) in accordance with the earlier observation of brassinolide binding sites in kinase-dead BRI1 mutant plants (Wang et al., 2001). These differences may reflect the absence of important regulatory factors, either from truncations of the protein or from expression in heterologous systems. Thus, it is extremely important that conclusions made about RLK function are based upon in vivo phosphorylation data.

What does the occurrence of the identified phosphorylation sites tell us about the activation status of the corresponding RLK? For the most common type of regulation, the phosphorylation of residues in the activation loop, it is safe to assume that the modification induces the well-characterized structural rearrangements that render the kinase active (Johnson et al., 1996; Huse and Kuriyan, 2002). Whether this is a regulated step or part of the protein maturation process as for many AGC kinases (Newton, 2003) remains to be established. Several mutations in the activation loops of plant RLKs that lead to aberrant phenotypes have been described (Dievart and Clark, 2003). It has also been shown that the phosphatase KAPP binds to the activation loop in a phosphorylation-dependent manner (Braun et al., 1997). These results underscore the key function of this domain at least in some proteins, but activation mechanisms independent of A loop modification, such as relief of autoinhibition by a domain that folds back onto the kinase core, are very common (Huse and Kuriyan, 2002).

Possible roles for phosphorylation in the highly divergent JM and C-terminal domains are intriguing. For both receptor Tyr and Ser/Thr kinases, a dual function of the JM domain has emerged in the last few years (Hubbard, 2001). The unphosphorylated JM domain contacts the N-terminal lobe of the kinase and keeps it in an inactive conformation, either on its own as for the ephrin receptor EphB2 (Wybenga-Groot et al., 2001) or with the help of an additional protein, such as immunophilin FKBP12 in the case of the TGFβRI (Huse et al., 2001). After kinase activation and phosphorylation of the JM, inhibition of the kinase is released. In addition, these phosphorylated regions then become docking sites for substrate proteins with domains that bind phosphotyrosine (SH2 or PTB domain proteins) or phosphoserine/threonine residues (MH2 or FHA domain proteins). The recent elucidation of the autoinhibited FLT3 kinase structure (Griffith et al., 2004) has added another mechanistically different version of autoinhibition by the JM. It will be interesting to test if any of these mechanisms are realized in plant RLKs. No details are known about intramolecular (autoinhibitory) interactions in RLKs, but two proteins binding to RLKs in a phosphorylation-dependent manner are known: the phosphatase KAPP (Li et al., 1999) and the ubiquitin ligase ARC (Stone et al., 2003). Additional and alternative mechanisms should not be overlooked. For example, the C-terminal hydrophobic motif of AGC kinases folds back onto the kinase domain upon phosphorylation and acts as an allosteric activator (Yang et al., 2002).

We have not been able to find common motifs among the JM and C-terminal phosphorylation sites, and they may indeed be highly individual regulatory sites. Although this is not encouraging in the quest for general paradigms about how plant RLKs work, it may be a requirement for signaling specificity among the hundreds of RLKs in Arabidopsis. Regardless of the current level of understanding, the discovery of phosphorylation in the JM and C-terminal domains of plant RLKs clearly indicates that these regions should be included in experiments designed to elucidate the mechanistic regulation of these proteins.

Few plant proteins in the databases are annotated with bona fide phosphorylation sites. Like in other organisms, the site determination from in vivo phosphorylated protein is very challenging. Many, or perhaps even most, sites are predicted by similarity to related proteins or through computer algorithms and later verified by site-directed mutagenesis. Our experimental data have shown that both strategies are potentially misleading. Although phosphorylation sites are conserved among most putative orthologs and some Arabidopsis paralogs, they do not necessarily stand out to facilitate the choice of sites to mutagenize. Concerning in silico predictions, we have found that approximately one-third of experimentally found sites would have been missed by Netphos (data not shown), one of the best available prediction tools (Blom et al., 1999) (http://www.cbs.dtu.dk/services/NetPhos/). This success rate is significantly lower than that estimated for the prediction of sites in mammalian proteins. In addition, the algorithms used in Netphos have been trained on well-characterized mammalian kinase substrates, and several of the best studied mammalian kinases (PKA and PKC) have no clear homologs in Arabidopsis. Thus, the algorithm probably yields a large number of false positives. To make in silico predictions for plants, we clearly need a much larger number of experimentally verified in vivo phosphorylation sites and (in vitro) determined kinase specificities.

Unlike the open reading frame of a protein, its phosphorylation status is dynamic, and any stocktaking of phosphorylated sites can only be a snapshot during the experimental conditions. Our experimental design, with samples pooled from control and elicitor-treated suspension-cultured cells, should encompass a wide variety of basal, stress-, and elicitor-induced phosphorylation sites. Still, kinases that are only activated by very specific stimuli or only in differentiated tissue would have been inactive and their substrates not found in this data set. The consequences are as follows: (1) none of the cataloged sites can be regarded as the regulatory site of the respective protein, but they may serve as a good first guess for mutagenesis and further biological characterization; (2) this study is only the start of a much larger-scale phosphoproteomic project that takes developmental stages, variously treated plants, differentiated tissues, and other subcellular structures into account. We currently are not able to perform comparisons of phosphoproteomes between two treatment conditions to identify differential components, but this achievement is the goal of our ongoing research.

Criticism has occasionally been raised about the descriptive nature of such large-scale projects, misjudging the complementary nature of hypothesis-driven and data-driven research (Kell and Oliver, 2004). Omic data may not reveal biological mechanisms per se but allow the formulation of new testable hypotheses that are not restrained by established ideas. For example, hypotheses about RLK regulation focused mostly on the activation loop. With the unanticipated complexity of phosphorylation sites presented here, new models will have to take other intramolecular interactions into account. The daunting task ahead of us is to understand the function of all these phosphorylation sites; but on the way to this goal, a careful analysis of the existing data can spark many new ideas.

METHODS

Materials

POROS chromatography materials (Self Pack OligoR3, MC 20, and HQ 20) were purchased from Applied Biosystems (Foster City, CA) and modified porcine trypsin from Promega (Southampton, UK). Microcolumns were packed in Gelloader tips (Eppendorf, Cambridge, UK) as described by Gobom et al. (1999).

Phosphopeptide Preparation

The isolation of plasma membranes from suspension-cultured Arabidopsis thaliana cells has been described by Nühse et al. (2003a). For this large-scale study, plasma membranes from control cells and cells that have been treated with 100 nM flg22 peptide for 8 min were pooled. Before harvesting of plasma membranes from the last upper phase of the dextran/polyethylene glycol partitioning, 0.01% Brij58 was added to invert the membranes to inside-out orientation. After centrifugation (60 min, 150,000g), the membranes were resuspended in 100 mM sodium carbonate, left on ice for 10 min, and then collected by centrifugation (30 min, 100,000g). The wash was repeated with 500 mM and then 50 mM ammonium hydrogen carbonate.

The membranes were digested at a 1:50 ratio of trypsin:protein overnight at 37°C. Phosphopeptides were isolated as described by Nühse et al. (2003b). Briefly, the supernatant containing the released peptides was either used directly for IMAC after addition of acetic acid (final concentration 0.2 M) or formic acid was added (final concentration 5%) to the digest and the supernatant purified over an OligoR3 column (i.e., washed with 5% formic acid and eluted with 100 mM acetic acid and 50% acetonitrile).

POROS MC material was prepared according to the manufacturer's instructions and loaded with 100 mM FeCl₃ in 100 mM acetic acid. Chromatography was performed as described by Stensballe et al. (2001)with minor modifications. Peptides were batch-bound to the IMAC material by shaking at room temperature in a typical volume of 20 to 50 μL, containing ∼2 to 5 μL (settled volume) of POROS MC material. After the incubation, the slurry was packed into Gelloader pipette tips with constricted tips and washed once each with 15 μL of 0.1 M acetic acid and 0.1 M acetic acid/30% acetonitrile, respectively, before elution with dilute ammonia, pH 10.5, or 50 mM ammonium phosphate, pH 9. The majority of phosphopeptides were identified from multiple runs of a two-dimensional liquid chromatographic separation scheme as described by Nühse et al. (2003b). Briefly, 500 μg of plasma membranes were trypsin digested and the released peptides purified over an OligoR3 microcolumn as described above but washed with water before elution with 50% acetonitrile. The eluate was diluted with buffer and pH adjusted to final concentrations of 30% acetonitrile/20 mM NH₄HCO₃, pH ∼7. A microcolumn was packed with POROS HQ20 (SAX) and preequilibrated with 30% acetonitrile/25 mM NH₄HCO₃, pH 7.0 (SAX buffer). The sample was slowly loaded onto the column and the flow-through collected. Twelve fractions were collected by step eluting with 20 μL each of 40 to 500 mM NaCl in SAX buffer. Flow-through and eluate fractions were briefly concentrated in a speed vac to reduce the acetonitrile concentration, brought to 5% formic acid, and desalted on R3 microcolumns. IMAC purification of phosphopeptides was done as described above.

Mass Spectrometry

Automated nanoflow LC-MS/MS analysis was performed using a QTOF Ultima mass spectrometer (Waters/Micromass, Manchester, UK) employing automated data-dependent acquisition. A nanoflow-HPLC system (Ultimate; Switchos2; Famos; LC Packings, Amstersdam, The Netherlands) was used to deliver a flow rate of 175 nL min⁻¹ to the mass spectrometer. Chromatographic separation was accomplished using a 2-cm fused silica precolumn (75-μm i.d. and 360-μm o.d.; Zorbax SB-C18 5 μm; Agilent, Wilmington, DE) connected to an 8-cm analytical column (50-μm i.d. and 360-μm o.d.; Agilent Zorbax SB-C18 3.5 μm). Peptides were eluted by a gradient of 5 to 32% acetonitrile for 35 min.

The mass spectrometer was operated in positive ion mode with a source temperature of 80°C and a countercurrent gas flow rate of 150 L h⁻¹. Data-dependent analysis was employed (three most abundant ions in each cycle): 1 s MS mass-to-charge ratio (m/z) 350 to 1500 and max 4 s MS/MS m/z 50 to 2000 (continuum mode), 30 s dynamic exclusion.

Raw data were processed using MassLynx 3.5 ProteinLynx (smooth 3/2 Savitzky Golay and center 4 channels/80% centroid) and the resulting MS/MS data set exported in the Micromass pkl format. We performed the peptide identification and assignment of partial posttranslational modifications using an in-house version of Mascot version 1.9. All data sets were searched twice, first with relatively large peptide mass tolerances, followed by internal mass recalibration by an in-house software algorithm using theoretical masses from unambiguously identified peptides obtained from the first search. The recalibrated data sets were searched against NCBInr (all species) using the following constraints: only tryptic peptides with up to three missed cleavage sites were allowed; 0.1 D mass tolerances for MS and MS/MS fragment ions. Phosphorylation (STY), deamidation (NQ), and oxidation (M) were specified as variable modifications. The results were filtered for non-Arabidopsis peptide assignments. External mass calibration using NaI resulted generally in mass errors of <50 ppm, typically 5 to 15 ppm in the m/z range 50 to 2000.

Approximately 50 of the MS/MS spectra were inspected for neutral loss, sequence-specific fragmentation patterns, and database match. Most of the spectra showed prominent loss of phosphate. The Mascot search results calculated all hits with a score above 40 to be significant; that is, a nonrandom match at P < 0.05 (the Mascot score is a probabilistic implementation [http://www.matrixscience.com/help/scoring_help.html] of the MOWSE score; Pappin et al., 1993). However, for all manually inspected spectra with a score above 30 and many spectra with a score above 20, the top search hit was in agreement with the spectrum on the basis of characteristic fragmentation. Correct hits may have a low score because of poorer quality spectra. Phosphorylation sites are only marked as unambiguous (bold red type in Supplemental Table 1 online) if sequences with alternative phosphate positions scored at least 10 units lower or were not reported as significant hits. Tryptic peptides that are not unique to one protein are marked in the table. All protein hits were verified in The Arabidopsis Information Resource database (http://www.arabidopsis.org/); no exclusions were made on the basis of predicted function, localization, etc.

Bioinformatics

The following Web-based bioinformatic tools were used: TMHMM (http://www.cbs.dtu.dk/services/TMHMM/) for the prediction of transmembrane helices; SMART (http://smart.embl-heidelberg.de/) for functional domains; ClustalW (http://www.ebi.ac.uk/clustalw/) for protein sequence alignments; Boxshade (http://www.ch.embnet.org/software/BOX_form.html) for rendering ClustalW results (setting 50% of sequences must agree for shading).

Putative orthologs of phosphoproteins were found by BLAST search of (http://www.ncbi.nlm.nih.gov/BLAST/; translated search tblastn) 40 to 60 amino acid sequences surrounding the empirical phosphorylation sites against an EST database (setting est_others; expect 10,000). Where possible, a few amino acids were included in the search that are conserved among several homologs of the phosphoprotein to distinguish unrelated sequences more easily.

Sequence data from this article have been deposited with the EMBL/GenBank data libraries under the following accession numbers: AHA2 (At4g30190), MDR1 (At2g36910), MDR4 (At2g47000), MDR8 (At1g02520), MDR17 (At3g62150), AtCesA1 (At4g32410), AtCesA2 (At4g39350), AtCesA3 (At5g05170), AtCesA4 (At5g44030), AtCesA5 (At5g09870), AtCesA6 (At5g64740), AtCesA7 (At5g17420), AtCesA8 (At4g18780), AtCesA9 (At2g21770), and AtCesA10 (At2g25540). For CDPKs, see PlantsP database (http://plantsp.sdsc.edu); for the STP family, see PlantsT (http://plantst.sdsc.edu). CesA3 putative orthologs are as follows: B. napus (CD837398); I. nil (BJ553825); M. truncatula (BQ139762 and BF634138); G. hirsutum (AF150630); P. balsamifera (BU879114); T. aestivum (BQ578769 and BT009438); Z. mays (ZmCesA1, AF200525; ZmCesA2, AF200526; ZmCesA4, AF200528; ZmCesA5, AF200529; ZmCesA9, AF200533); L. japonicus (AV408250); O. sativa (AK072356.1, AK098978, and D48636); G. max (BE660209); L. sativa (BQ869850); Rosa sp (BQ104792); V. vinifera (CD005820); S. propinquum (BG465505); S. bicolor (BG273410); H. vulgare (BQ471177); T. monococcum (BQ802778). SOS1 (At2g01980); putative orthologs: P. trem (CA924192); G. max (BE609663); B. vulgaris (BQ591881); S. tuberosum (BQ510147); M. chrysanthemum (BE577626); O. sativa (CB676883).

ACKNOWLEDGMENTS

This work was supported by Biotechnology and Biological Science Research Council Grant 83/C17990 (T.S.N. and S.C.P.), by the Gatsby Charitable Foundation (T.S.N. and S.C.P.), by a grant from the Danish Natural Sciences Research Council (O.N.J.), by an EMBO short-term fellowship (T.S.N.), and by a Danish Industrial PhD fellowship (A.S.).

Author notes

[W]

Online version contains Web-only data.

Article, publication date, and citation information can be found at www.plantcell.org/cgi/doi/10.1105/tpc.104.023150.

1

To whom correspondence should be addressed. E-mail ; fax [phone].

The author responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (www.plantcell.org) is: Scott Peck ().