The plant kinome

Plant kinases are one of the largest protein families in Arabidopsis. There are almost 600 membrane-located receptor kinases and almost 400 soluble kinases with distinct functions in signal transduction. In this minireview we discuss phylogeny and functional context of prominent members from major protein kinase subfamilies in plants.

Receptor kinase-mediated control of primary active proton pumping at the plasma membrane

Acidification of the cell wall space outside the plasma membrane is required for plant growth and is the result of proton extrusion by the plasma membrane-localized H+-ATPases. Here we show that the major plasma membrane proton pumps in Arabidopsis, AHA1 and AHA2, interact directly in vitro and in planta with PSY1R, a receptor kinase of the plasma membrane that serves as a receptor for the peptide growth hormone PSY1. The intracellular protein kinase domain of PSY1R phosphorylates AHA2/AHA1 at Thr-881, situated in the autoinhibitory region I of the C-terminal domain. When expressed in a yeast heterologous expression system, the introduction of a negative charge at this position caused pump activation. Application of PSY1 to plant seedlings induced rapid in planta phosphorylation at Thr-881, concomitant with an instantaneous increase in proton efflux from roots. The direct interaction between AHA2 and PSY1R observed might provide a general paradigm for regulation of plasma membrane proton transport by receptor kinases.

The Arabidopsis Kinome: phylogeny and evolutionary insights into functional diversification

Background

Protein kinases constitute a particularly large protein family in Arabidopsis with important functions in cellular signal transduction networks. At the same time Arabidopsis is a model plant with high frequencies of gene duplications. Here, we have conducted a systematic analysis of the Arabidopsis kinase complement, the kinome, with particular focus on gene duplication events. We matched Arabidopsis proteins to a Hidden-Markov Model of eukaryotic kinases and computed a phylogeny of 942 Arabidopsis protein kinase domains and mapped their origin by gene duplication.

Results

The phylogeny showed two major clades of receptor kinases and soluble kinases, each of which was divided into functional subclades. Based on this phylogeny, association of yet uncharacterized kinases to families was possible which extended functional annotation of unknowns. Classification of gene duplications within these protein kinases revealed that representatives of cytosolic subfamilies showed a tendency to maintain segmentally duplicated genes, while some subfamilies of the receptor kinases were enriched for tandem duplicates. Although functional diversification is observed throughout most subfamilies, some instances of functional conservation among genes transposed from the same ancestor were observed. In general, a significant enrichment of essential genes was found among genes encoding for protein kinases.

Conclusions

The inferred phylogeny allowed classification and annotation of yet uncharacterized kinases. The prediction and analysis of syntenic blocks and duplication events within gene families of interest can be used to link functional biology to insights from an evolutionary viewpoint. The approach undertaken here can be applied to any gene family in any organism with an annotated genome.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-548) contains supplementary material, which is available to authorized users.

A kinase-phosphatase signaling module with BSK8 and BSL2 involved in regulation of sucrose-phosphate synthase

External supply of sucrose to carbon-starved Arabidopsis seedlings induced changes in phosphorylation of Brassinosteroid Signaling Kinase 8 (BSK8) at two different sites. Serine S(20) lies within a phosphorylation hotspot at the N-terminal region of the protein, while S(213) is located within the kinase domain of BSK8. Upon sucrose supply phosphorylation of BSK8(S20) and BSK8(S213) showed opposite behavior with increasing phosphorylation of S(213) and decreased phosphorylation of S(20) at 5 min after sucrose supply. Here we aim to systematically analyze the effects of BSK8 mutations on downstream cellular regulatory events and characterize molecular functions of BSK8 and its phosphorylation. Comparative phosphoproteomic profiling of a bsk8 knockout mutant and wild type revealed potential targets in sucrose metabolism. Activity of sucrose-phosphate synthase (SPS) was decreased by phosphorylation at S(152), and SPS phosphorylation inversely correlated with sucrose-induced BSK8 activity. Furthermore, BSK8 was found to interact with BSL2, a Kelch-type phosphatase. On the basis of a combination of kinase activity measurements, SPS activity assays, and phosphorylation site mutations in BSK8 at S(20) and S(213), we conclude that regulation of SPS by BSK8 occurs through activation of a phosphatase that in turn may dephosphorylate SPS and thus activates the enzyme.

Auxin efflux by PIN-FORMED proteins is activated by two different protein kinases, D6 PROTEIN KINASE and PINOID

The development and morphology of vascular plants is critically determined by synthesis and proper distribution of the phytohormone auxin. The directed cell-to-cell distribution of auxin is achieved through a system of auxin influx and efflux transporters. PIN-FORMED (PIN) proteins are proposed auxin efflux transporters, and auxin fluxes can seemingly be predicted based on the--in many cells--asymmetric plasma membrane distribution of PINs. Here, we show in a heterologous Xenopus oocyte system as well as in Arabidopsis thaliana inflorescence stems that PIN-mediated auxin transport is directly activated by D6 PROTEIN KINASE (D6PK) and PINOID (PID)/WAG kinases of the Arabidopsis AGCVIII kinase family. At the same time, we reveal that D6PKs and PID have differential phosphosite preferences. Our study suggests that PIN activation by protein kinases is a crucial component of auxin transport control that must be taken into account to understand auxin distribution within the plant.

Auxin efflux by PIN-FORMED proteins is activated by two different protein kinases, D6 PROTEIN KINASE and PINOID

The development and morphology of vascular plants is critically determined by synthesis and proper distribution of the phytohormone auxin. The directed cell-to-cell distribution of auxin is achieved through a system of auxin influx and efflux transporters. PIN-FORMED (PIN) proteins are proposed auxin efflux transporters, and auxin fluxes can seemingly be predicted based on the—in many cells—asymmetric plasma membrane distribution of PINs. Here, we show in a heterologous Xenopus oocyte system as well as in Arabidopsis thaliana inflorescence stems that PIN-mediated auxin transport is directly activated by D6 PROTEIN KINASE (D6PK) and PINOID (PID)/WAG kinases of the Arabidopsis AGCVIII kinase family. At the same time, we reveal that D6PKs and PID have differential phosphosite preferences. Our study suggests that PIN activation by protein kinases is a crucial component of auxin transport control that must be taken into account to understand auxin distribution within the plant.

DOI:http://dx.doi.org/10.7554/eLife.02860.001

In plants, a hormone called auxin controls the growth of the stems and roots. This chemical is transported from cell to cell, and its flow though the plant is redirected continuously as the plant is developing. Auxin is pumped out of cells by proteins in the cell membrane called ‘auxin efflux carriers’. These proteins are usually found on one side of each cell and this is what gives the direction to auxin transport.

Zourelidou, Absmanner et al. now report that being positioned on the correct side of a plant cell is not enough to enable an efflux carrier to do its job—it must also be turned on by kinases before it can pump auxin out of cells. Kinases are enzymes that add phosphate groups to specific sites on other proteins, and plants without certain kinases are unable to transport auxin.

When Zourelidou, Absmanner et al. produced the efflux carrier and a plant kinase—which turns the efflux carrier on—in immature egg cells from frogs, auxin was rapidly pumped out of the cells. However, cells that contained the efflux carrier but not the kinase could not transport the hormone. Importantly egg cells from frogs do not normally transport auxin, but these cells are commonly used in experiments because they are large, which makes them easier to work with in the lab.

One of at least two kinases must tag a number of sites on the efflux carrier to ensure that it is switched on. It was already known that some of these sites are involved in making sure that the efflux carrier is located on the correct side of the cell. Zourelidou, Absmanner et al. also found that auxin itself encourages the addition of phosphate groups onto the efflux carrier.

Though it was thought that knowing where the auxin transporters are was enough to explain the direction of auxin transport in plants, it is now clear that activation by the kinases needs to be taken into account too. And since these kinases may activate the transporters to different extents, identifying how these proteins are controlled, for example by auxin itself, will be the next challenge in the field.

DOI:http://dx.doi.org/10.7554/eLife.02860.002

Meta-Analysis of Arabidopsis thaliana Phospho-Proteomics Data Reveals Compartmentalization of Phosphorylation Motifs

Protein (de)phosphorylation plays an important role in plants. To provide a robust foundation for subcellular phosphorylation signaling network analysis and kinase-substrate relationships, we performed a meta-analysis of 27 published and unpublished in-house mass spectrometry-based phospho-proteome data sets for Arabidopsis thaliana covering a range of processes, (non)photosynthetic tissue types, and cell cultures. This resulted in an assembly of 60,366 phospho-peptides matching to 8141 nonredundant proteins. Filtering the data for quality and consistency generated a set of medium and a set of high confidence phospho-proteins and their assigned phospho-sites. The relation between single and multiphosphorylated peptides is discussed. The distribution of p-proteins across cellular functions and subcellular compartments was determined and showed overrepresentation of protein kinases. Extensive differences in frequency of pY were found between individual studies due to proteomics and mass spectrometry workflows. Interestingly, pY was underrepresented in peroxisomes but overrepresented in mitochondria. Using motif-finding algorithms motif-x and MMFPh at high stringency, we identified compartmentalization of phosphorylation motifs likely reflecting localized kinase activity. The filtering of the data assembly improved signal/noise ratio for such motifs. Identified motifs were linked to kinases through (bioinformatic) enrichment analysis. This study also provides insight into the challenges/pitfalls of using large-scale phospho-proteomic data sets to nonexperts.

Dissecting the subcellular compartmentation of proteins and metabolites in arabidopsis leaves using non-aqueous fractionation

Non-aqueous fractionation is a technique for the enrichment of different subcellular compartments derived from lyophilized material. It was developed to study the subcellular distribution of metabolites. Here we analyzed the distribution of about 1,000 proteins and 70 metabolites, including 22 phosphorylated intermediates in wild-type Arabidopsis rosette leaves, using non-aqueous gradients divided into 12 fractions. Good separation of plastidial, cytosolic, and vacuolar metabolites and proteins was achieved, but cytosolic, mitochondrial, and peroxisomal proteins clustered together. There was considerable heterogeneity in the fractional distribution of transcription factors, ribosomal proteins, and subunits of the vacuolar-ATPase, indicating diverse compartmental location. Within the plastid, sub-organellar separation of thylakoids and stromal proteins was observed. Metabolites from the Calvin-Benson cycle, photorespiration, starch and sucrose synthesis, glycolysis, and the tricarboxylic acid cycle grouped with their associated proteins of the respective compartment. Non-aqueous fractionation thus proved to be a powerful method for the study of the organellar, and in some cases sub-organellar, distribution of proteins and their association with metabolites. It remains the technique of choice for the assignment of subcellular location to metabolites in intact plant tissues, and thus the technique of choice for doing combined metabolite-protein analysis on a single tissue sample.

Plasma membrane lipid-protein interactions affect signaling processes in sterol-biosynthesis mutants in Arabidopsis thaliana

The plasma membrane is an important organelle providing structure, signaling and transport as major biological functions. Being composed of lipids and proteins with different physicochemical properties, the biological functions of membranes depend on specific protein-protein and protein-lipid interactions. Interactions of proteins with their specific sterol and lipid environment were shown to be important factors for protein recruitment into sub-compartmental structures of the plasma membrane. System-wide implications of altered endogenous sterol levels for membrane functions in living cells were not studied in higher plant cells. In particular, little is known how alterations in membrane sterol composition affect protein and lipid organization and interaction within membranes. Here, we conducted a comparative analysis of the plasma membrane protein and lipid composition in Arabidopsis sterol-biosynthesis mutants smt1 and ugt80A2;B1. smt1 shows general alterations in sterol composition while ugt80A2;B1 is significantly impaired in sterol glycosylation. By systematically analyzing different cellular fractions and combining proteomic with lipidomic data we were able to reveal contrasting alterations in lipid-protein interactions in both mutants, with resulting differential changes in plasma membrane signaling status.

SILAC and alternatives in studying cellular proteomes of plants

Quantitative proteomics by metabolic labeling has a high impact on the growing field of plant systems biology. SILAC has been pioneered and optimized for plant cell culture systems allowing for SILAC-based quantitative experiments in specialized experimental setups. In comparison to other model organisms, the application of SILAC to whole plants is challenging. As autotrophic organisms, plants under their natural growth conditions can hardly be fully labeled with stable isotope-coded amino acids. The metabolic labeling with inorganic nitrogen is therefore the method of choice for most whole-plant physiological questions. Plants can easily metabolize different inorganic nitrogen isotopes. The incorporation of the labeled inorganic nitrogen then results in proteins and metabolites with distinct molecular mass, which can be detected on a mass spectrometer. In comparative quantitative experiments, similarly as in SILAC experiments, treated and untreated samples are differentially labeled by nitrogen isotopes and jointly processed, thereby minimizing sample-to-sample variation. In recent years, heavy nitrogen labeling has become a widely used strategy in quantitative proteomics and novel approaches were developed for metabolite identification. Here we present a typical hydroponics setup, the workflow for processing of samples, mass spectrometry and data analysis for large-scale metabolic labeling experiments of whole plants.

Plasma membrane lipid-protein interactions affect signaling processes in sterol-biosynthesis mutants in Arabidopsis thaliana

The plasma membrane is an important organelle providing structure, signaling and transport as major biological functions. Being composed of lipids and proteins with different physicochemical properties, the biological functions of membranes depend on specific protein-protein and protein-lipid interactions. Interactions of proteins with their specific sterol and lipid environment were shown to be important factors for protein recruitment into sub-compartmental structures of the plasma membrane. System-wide implications of altered endogenous sterol levels for membrane functions in living cells were not studied in higher plant cells. In particular, little is known how alterations in membrane sterol composition affect protein and lipid organization and interaction within membranes. Here, we conducted a comparative analysis of the plasma membrane protein and lipid composition in Arabidopsis sterol-biosynthesis mutants smt1 and ugt80A2;B1. smt1 shows general alterations in sterol composition while ugt80A2;B1 is significantly impaired in sterol glycosylation. By systematically analyzing different cellular fractions and combining proteomic with lipidomic data we were able to reveal contrasting alterations in lipid-protein interactions in both mutants, with resulting differential changes in plasma membrane signaling status.

Reconstruction and analysis of nutrient-induced phosphorylation networks in Arabidopsis thaliana

Elucidating the dynamics of molecular processes in living organisms in response to external perturbations is a central goal in modern systems biology. We investigated the dynamics of protein phosphorylation events in Arabidopsis thaliana exposed to changing nutrient conditions. Phosphopeptide expression levels were detected at five consecutive time points over a time interval of 30 min after nutrient resupply following prior starvation. The three tested inorganic, ionic nutrients NH(+) 4, NO(-) 3, PO(3-) 4 elicited similar phosphosignaling responses that were distinguishable from those invoked by the sugars mannitol, sucrose. When embedded in the protein-protein interaction network of Arabidopsis thaliana, phosphoproteins were found to exhibit a higher degree compared to average proteins. Based on the time-series data, we reconstructed a network of regulatory interactions mediated by phosphorylation. The performance of different network inference methods was evaluated by the observed likelihood of physical interactions within and across different subcellular compartments and based on gene ontology semantic similarity. The dynamic phosphorylation network was then reconstructed using a Pearson correlation method with added directionality based on partial variance differences. The topology of the inferred integrated network corresponds to an information dissemination architecture, in which the phosphorylation signal is passed on to an increasing number of phosphoproteins stratified into an initiation, processing, and effector layer. Specific phosphorylation peptide motifs associated with the distinct layers were identified indicating the action of layer-specific kinases. Despite the limited temporal resolution, combined with information on subcellular location, the available time-series data proved useful for reconstructing the dynamics of the molecular signaling cascade in response to nutrient stress conditions in the plant Arabidopsis thaliana.

Sucrose-induced receptor kinase SIRK1 regulates a plasma membrane aquaporin in Arabidopsis

The transmembrane receptor kinase family is the largest protein kinase family in Arabidopsis, and it contains the highest fraction of proteins with yet uncharacterized functions. Here, we present functions of SIRK1, a receptor kinase that was previously identified with rapid transient phosphorylation after sucrose resupply to sucrose-starved seedlings. SIRK1 was found to be an active kinase with increasing activity in the presence of an external sucrose supply. In sirk1 T-DNA insertional mutants, the sucrose-induced phosphorylation patterns of several membrane proteins were strongly reduced; in particular, pore-gating phosphorylation sites in aquaporins were affected. SIRK1-GFP fusions were found to directly interact with aquaporins in affinity pull-down experiments on microsomal membrane vesicles. Furthermore, protoplast swelling assays of sirk1 mutants and SIRK1-GFP expressing lines confirmed a direct functional interaction of receptor kinase SIRK1 and aquaporins as substrates for phosphorylation. A lack of SIRK1 expression resulted in the failure of mutant protoplasts to control water channel activity upon changes in external sucrose concentrations. We propose that SIRK1 is involved in the regulation of sucrose-specific osmotic responses through direct interaction with and activation of an aquaporin via phosphorylation and that the duration of this response is controlled by phosphorylation-dependent receptor internalization.

Metabolic labeling and membrane fractionation for comparative proteomic analysis of Arabidopsis thaliana suspension cell cultures

Plasma membrane microdomains are features based on the physical properties of the lipid and sterol environment and have particular roles in signaling processes. Extracting sterol-enriched membrane microdomains from plant cells for proteomic analysis is a difficult task mainly due to multiple preparation steps and sources for contaminations from other cellular compartments. The plasma membrane constitutes only about 5-20% of all the membranes in a plant cell, and therefore isolation of highly purified plasma membrane fraction is challenging. A frequently used method involves aqueous two-phase partitioning in polyethylene glycol and dextran, which yields plasma membrane vesicles with a purity of 95% (1). Sterol-rich membrane microdomains within the plasma membrane are insoluble upon treatment with cold nonionic detergents at alkaline pH. This detergent-resistant membrane fraction can be separated from the bulk plasma membrane by ultracentrifugation in a sucrose gradient (2). Subsequently, proteins can be extracted from the low density band of the sucrose gradient by methanol/chloroform precipitation. Extracted protein will then be trypsin digested, desalted and finally analyzed by LC-MS/MS. Our extraction protocol for sterol-rich microdomains is optimized for the preparation of clean detergent-resistant membrane fractions from Arabidopsis thaliana cell cultures. We use full metabolic labeling of Arabidopsis thaliana suspension cell cultures with K(15)NO3 as the only nitrogen source for quantitative comparative proteomic studies following biological treatment of interest (3). By mixing equal ratios of labeled and unlabeled cell cultures for joint protein extraction the influence of preparation steps on final quantitative result is kept at a minimum. Also loss of material during extraction will affect both control and treatment samples in the same way, and therefore the ratio of light and heave peptide will remain constant. In the proposed method either labeled or unlabeled cell culture undergoes a biological treatment, while the other serves as control (4).

Unraveling sterol-dependent membrane phenotypes by analysis of protein abundance-ratio distributions in different membrane fractions under biochemical and endogenous sterol depletion

During the last decade, research on plasma membrane focused increasingly on the analysis of so-called microdomains. It has been shown that function of many membrane-associated proteins involved in signaling and transport depends on their conditional segregation within sterol-enriched membrane domains. High throughput proteomic analysis of sterol-protein interactions are often based on analyzing detergent resistant membrane fraction enriched in sterols and associated proteins, which also contain proteins from these microdomain structures. Most studies so far focused exclusively on the characterization of detergent resistant membrane protein composition and abundances. This approach has received some criticism because of its unspecificity and many co-purifying proteins. In this study, by a label-free quantitation approach, we extended the characterization of membrane microdomains by particularly studying distributions of each protein between detergent resistant membrane and detergent-soluble fractions (DSF). This approach allows a more stringent definition of dynamic processes between different membrane phases and provides a means of identification of co-purifying proteins. We developed a random sampling algorithm, called Unicorn, allowing for robust statistical testing of alterations in the protein distribution ratios of the two different fractions. Unicorn was validated on proteomic data from methyl-β-cyclodextrin treated plasma membranes and the sterol biosynthesis mutant smt1. Both, chemical treatment and sterol-biosynthesis mutation affected similar protein classes in their membrane phase distribution and particularly proteins with signaling and transport functions.

Calcium-dependent protein kinase/NADPH oxidase activation circuit is required for rapid defense signal propagation

In animals and plants, pathogen recognition triggers the local activation of intracellular signaling that is prerequisite for mounting systemic defenses in the whole organism. We identified that Arabidopsis thaliana isoform CPK5 of the plant calcium-dependent protein kinase family becomes rapidly biochemically activated in response to pathogen-associated molecular pattern (PAMP) stimulation. CPK5 signaling resulted in enhanced salicylic acid-mediated resistance to the bacterial pathogen Pst DC3000, differential plant defense gene expression, and synthesis of reactive oxygen species (ROS). Using selected reaction monitoring MS, we identified the plant NADPH oxidase, respiratory burst oxidase homolog D (RBOHD), as an in vivo phosphorylation target of CPK5. Remarkably, CPK5-dependent in vivo phosphorylation of RBOHD occurs on both PAMP- and ROS stimulation. Furthermore, rapid CPK5-dependent biochemical and transcriptional activation of defense reactions at distal sites is compromised in cpk5 and rbohd mutants. Our data not only identify CPK5 as a key regulator of innate immune responses in plants but also support a model of ROS-mediated cell-to-cell communication, where a self-propagating mutual activation circuit consisting of the protein kinase, CPK5, and the NADPH oxidase RBOHD facilitates rapid signal propagation as a prerequisite for defense response activation at distal sites within the plant.

Calcium-dependent protein kinase/NADPH oxidase activation circuit is required for rapid defense signal propagation

In animals and plants, pathogen recognition triggers the local activation of intracellular signaling that is prerequisite for mounting systemic defenses in the whole organism. We identified that Arabidopsis thaliana isoform CPK5 of the plant calcium-dependent protein kinase family becomes rapidly biochemically activated in response to pathogen-associated molecular pattern (PAMP) stimulation. CPK5 signaling resulted in enhanced salicylic acid-mediated resistance to the bacterial pathogen Pst DC3000, differential plant defense gene expression, and synthesis of reactive oxygen species (ROS). Using selected reaction monitoring MS, we identified the plant NADPH oxidase, respiratory burst oxidase homolog D (RBOHD), as an in vivo phosphorylation target of CPK5. Remarkably, CPK5-dependent in vivo phosphorylation of RBOHD occurs on both PAMP- and ROS stimulation. Furthermore, rapid CPK5-dependent biochemical and transcriptional activation of defense reactions at distal sites is compromised in cpk5 and rbohd mutants. Our data not only identify CPK5 as a key regulator of innate immune responses in plants but also support a model of ROS-mediated cell-to-cell communication, where a self-propagating mutual activation circuit consisting of the protein kinase, CPK5, and the NADPH oxidase RBOHD facilitates rapid signal propagation as a prerequisite for defense response activation at distal sites within the plant.

Function of calcium-dependent protein kinase CPK28 of Arabidopsis thaliana in plant stem elongation and vascular development

After a period of vegetative growth, plants undergo a developmental switch to the reproductive phase, inducing the transition to bolting, elongation of the inflorescence and flowering. We have identified calcium-dependent protein kinase CPK28 from Arabidopsis thaliana as a regulatory component that controls stem elongation and vascular development. In two independent mutant alleles of cpk28, a reduction of stem elongation, accompanied by shorter leaf petioles and enhanced anthocyanin levels, is observed upon the transition to the generative phase. Anatomical analysis revealed an altered vascular pattern characterised by fewer xylem tracheary elements but at the same time increased lignification and secondary growth. Coincident with these morphological changes, cpk28 mutants showed altered expression of NAC transcriptional regulators NST1 and NST3 as well as of GA3ox1, a key regulator of gibberellic acid homeostasis. In vitro protein kinase activity of CPK28 is strictly calcium-dependent. Furthermore, CPK28 is phosphorylated in vivo at several sites. Site-specific amino acid substitutions at these phosphorylation sites resulted in reduced in vitro activity. However, when introduced into a cpk28 mutant background, wild-type and phosphorylation site variants, but not kinase-inactive variants of CPK28 complemented the morphological and developmental defects. Our data identify CPK28 as a developmentally controlled regulator for coordinated stem elongation and secondary growth.

PhosPhAt goes kinases--searchable protein kinase target information in the plant phosphorylation site database PhosPhAt

Reversible phosphorylation is a key mechanism for regulating protein function. Thus it is of high interest to know which kinase can phosphorylate which proteins. Comprehensive information about phosphorylation sites in Arabidopsis proteins is hosted within the PhosPhAt database (http://phosphat.mpimp-golm.mpg.de). However, our knowledge of the kinases that phosphorylate those sites is dispersed throughout the literature and very difficult to access, particularly for investigators seeking to interpret large scale and high-throughput experiments. Therefore, we aimed to compile information on kinase–substrate interactions and kinase-specific regulatory information and make this available via a new functionality embedded in PhosPhAt. Our approach involved systematic surveying of the literature for regulatory information on the members of the major kinase families in Arabidopsis thaliana, such as CDPKs, MPK(KK)s, AGC kinases and SnRKs, as well as individual kinases from other families. To date, we have researched more than 4450 kinase-related publications, which collectively contain information on about 289 kinases. Users can now query the PhosPhAt database not only for experimental and predicted phosphorylation sites of individual proteins, but also for known substrates for a given kinase or kinase family. Further developments include addition of new phosphorylation sites and visualization of clustered phosphorylation events, known as phosphorylation hotspots.

The protein composition of the digestive fluid from the venus flytrap sheds light on prey digestion mechanisms

The Venus flytrap (Dionaea muscipula) is one of the most well-known carnivorous plants because of its unique ability to capture small animals, usually insects or spiders, through a unique snap-trapping mechanism. The animals are subsequently killed and digested so that the plants can assimilate nutrients, as they grow in mineral-deficient soils. We deep sequenced the cDNA from Dionaea traps to obtain transcript libraries, which were used in the mass spectrometry-based identification of the proteins secreted during digestion. The identified proteins consisted of peroxidases, nucleases, phosphatases, phospholipases, a glucanase, chitinases, and proteolytic enzymes, including four cysteine proteases, two aspartic proteases, and a serine carboxypeptidase. The majority of the most abundant proteins were categorized as pathogenesis-related proteins, suggesting that the plant's digestive system evolved from defense-related processes. This in-depth characterization of a highly specialized secreted fluid from a carnivorous plant provides new information about the plant's prey digestion mechanism and the evolutionary processes driving its defense pathways and nutrient acquisition.

Proteomics wants cRacker: automated standardized data analysis of LC-MS derived proteomic data

The large-scale analysis of thousands of proteins under various experimental conditions or in mutant lines has gained more and more importance in hypothesis-driven scientific research and systems biology in the past years. Quantitative analysis by large scale proteomics using modern mass spectrometry usually results in long lists of peptide ion intensities. The main interest for most researchers, however, is to draw conclusions on the protein level. Postprocessing and combining peptide intensities of a proteomic data set requires expert knowledge, and the often repetitive and standardized manual calculations can be time-consuming. The analysis of complex samples can result in very large data sets (lists with several 1000s to 100,000 entries of different peptides) that cannot easily be analyzed using standard spreadsheet programs. To improve speed and consistency of the data analysis of LC-MS derived proteomic data, we developed cRacker. cRacker is an R-based program for automated downstream proteomic data analysis including data normalization strategies for metabolic labeling and label free quantitation. In addition, cRacker includes basic statistical analysis, such as clustering of data, or ANOVA and t tests for comparison between treatments. Results are presented in editable graphic formats and in list files.

Characterization and Prediction of Protein Phosphorylation Hotspots in Arabidopsis thaliana

The regulation of protein function by modulating the surface charge status via sequence-locally enriched phosphorylation sites (P-sites) in so called phosphorylation "hotspots" has gained increased attention in recent years. We set out to identify P-hotspots in the model plant Arabidopsis thaliana. We analyzed the spacing of experimentally detected P-sites within peptide-covered regions along Arabidopsis protein sequences as available from the PhosPhAt database. Confirming earlier reports (Schweiger and Linial, 2010), we found that, indeed, P-sites tend to cluster and that distributions between serine and threonine P-sites to their respected closest next P-site differ significantly from those for tyrosine P-sites. The ability to predict P-hotspots by applying available computational P-site prediction programs that focus on identifying single P-sites was observed to be severely compromised by the inevitable interference of nearby P-sites. We devised a new approach, named HotSPotter, for the prediction of phosphorylation hotspots. HotSPotter is based primarily on local amino acid compositional preferences rather than sequence position-specific motifs and uses support vector machines as the underlying classification engine. HotSPotter correctly identified experimentally determined phosphorylation hotspots in A. thaliana with high accuracy. Applied to the Arabidopsis proteome, HotSPotter-predicted 13,677 candidate P-hotspots in 9,599 proteins corresponding to 7,847 unique genes. Hotspot containing proteins are involved predominantly in signaling processes confirming the surmised modulating role of hotspots in signaling and interaction events. Our study provides new bioinformatics means to identify phosphorylation hotspots and lays the basis for further investigating novel candidate P-hotspots. All phosphorylation hotspot annotations and predictions have been made available as part of the PhosPhAt database at http://phosphat.mpimp-golm.mpg.de.