In silico analysis of protein Lys-N𝜀-acetylation in plants
- 1Division of Biochemistry, University of Missouri, Columbia, MO, USA
- 2Interdisciplinary Plant Group, University of Missouri, Columbia, MO, USA
- 3Plant Genetics Research Unit, United States Department of Agriculture – Agricultural Research Service, Columbia, MO, USA
Among post-translational modifications, there are some conceptual similarities between Lys-N𝜀-acetylation and Ser/Thr/Tyr O-phosphorylation. Herein we present a bioinformatics-based overview of reversible protein Lys-acetylation, including some comparisons with reversible protein phosphorylation. The study of Lys-acetylation of plant proteins has lagged behind studies of mammalian and microbial cells; 1000s of acetylation sites have been identified in mammalian proteins compared with only hundreds of sites in plant proteins. While most previous emphasis was focused on post-translational modifications of histones, more recent studies have addressed metabolic regulation. Being directly coupled with cellular CoA/acetyl-CoA and NAD/NADH, reversible Lys-N𝜀-acetylation has the potential to control, or contribute to control, of primary metabolism, signaling, and growth and development.
Protein regulation is highly choreographed, and encompasses multilayered and interconnected transcriptional and translational processes. Subsequent PTM allow the fine-tuning of function. Many PTM are reversible, and modulate the activity, localization, turnover, and interactions of proteins (Mann and Jensen, 2003). The dynamics of PTM allow a much more rapid response to cellular signals than transcriptional or translational regulation.
More than 461 distinct PTM have been described1 (Khoury et al., 2011), and it is increasingly clear that many, perhaps most, proteins are decorated with multiple PTM (Hunter, 2007). This yields a proteome vastly larger and more diverse than the translated genome (Khoury et al., 2011; Minguez et al., 2012). The combinatorial diversity of PTM results in enormous flexibility for control of structure, function, and interaction. For example, Ser/Thr/Tyr-phosphorylation introduces a negative charge while Lys-N𝜀-acetylation neutralizes a positive charge.
Reversible protein phosphorylation was first reported in 1954, and has remained the hallmark of PTM-based regulation and signaling (Pawson and Scott, 2005). Approximately 130 protein kinases are encoded by the Saccharomyces cerevisiae genome, 518 in Homo sapiens, 1085 in Arabidopsis thaliana, and as many as 2500 in some land plants (Lehti-Shiu and Shiu, 2012). In yeast there are >2500 known phosphoproteins (P-proteins; Rao and Møller, 2012), more than 31,480 phosphosites (P-sites) have been described in 7280 Rattus norvegicus proteins (Lundby et al., 2012b), and ∼14000 P-sites in >4000 proteins have been described from the reference eudicot plant A. thaliana2. Considering all plants, these numbers swell to ∼33,300 sites on ∼11,700 P-proteins2, and more are discovered virtually daily. The recent reports of thousands of Lys-N𝜀-acetylation sites in bacterial, yeast, and animal proteins have led to the suggestion that this PTM could be as common as phosphorylation (Kouzarides, 2000; Choudhary et al., 2009; Hu et al., 2010; Filippakopoulos and Knapp, 2012; Lundby et al., 2012a,b; Xiong and Guan, 2012). An overview comparison between these two major systems for PTM is presented in Table 1.
TABLE 1. The components of the Ser/Thr/Tyr O-phosphorylation and Lys-N𝜀-acetylation systems from Homo sapiens, Saccharomyces cerevisiae, and Arabidopsis thaliana.
There are multiple forms of protein acetylation; O-acetylation of Ser/Thr residues, non-reversible Nα-acetylation of protein N-termini (Martinez et al., 2008), and reversible protein N𝜀-acetylation of internal Lys residues (Soppa, 2010; Xing and Poirier, 2012). While important in their own right, the first two types of Lys-acetylation will not be further addressed herein. Instead, “acetylation” should be understood to mean protein Lys-N𝜀-acetylation (PKA). Histones were first reported to be acetylated nearly 50 years ago (Phillips, 1963). It was more than two decades later that the first cytoplasmic protein, α-tubulin, was identified as being Lys-acetylated (Piperno and Fuller, 1985). Recently, a plethora of acetylated proteins have been described encompassing all subcellular compartments (Choudhary et al., 2009; Table 2). During these few years, PKA has changed from an obscure histone PTM to a mechanism for controlling (or contributing to control of) many aspects of primary metabolism, gene expression, signaling, and development (Zhao et al., 2010; Rao et al., 2014).
Use of high-throughput, high-resolution tandem MS has led to detection of PKA in all three domains of life; archaea, bacteria, and eukaryotes (Yang, 2004a; Hu et al., 2010; Soppa, 2010; Jones and O’Connor, 2011). Recent MS-based studies revealed 1070 acetylation sites in 349 proteins in Escherichia coli (Zhang et al., 2013), and 332 acetylation sites on 185 proteins in Bacillus subtilis (Kim et al., 2013). In S. cerevisiae, more than 4000 sites of Lys-N𝜀-acetylation of 1059 proteins have been described (Henriksen et al., 2012), and in mammalian systems >21,000 sites on >7000 proteins have been described (Choudhary et al., 2009; Zhao et al., 2010; Lundby et al., 2012a). A large majority of the yeast and mammalian Lys-acetylation sites have been described since 2009 (Choudhary et al., 2009; Smith and Workman, 2009), suggesting that our understanding and appreciation for this PTM is in a relatively early stage especially in comparison with phosphorylation (Tables 1 and 2).
The results from large-scale secondary structure analyses have led to the conclusion that O-phosphorylation is substantially enriched in regions of ISD (Iakoucheva et al., 2004). It has been reported that PKA sites are significantly enriched in ordered regions of mammalian proteins and depleted in regions of ISD (Choudhary et al., 2009), that PKA sites are equally distributed in ordered and disordered regions (Gao and Xu, 2011), and that PKA sites preferentially occur in regions of ISD in Toxoplasma gondii tachyzoite proteins (Xue et al., 2013). We observed that versus all Lys residues in our soybean database, PKA was approximately twice as likely to occur in long ISD-regions (Smith-Hammond et al., 2014b). If, as has recently been proposed (Cumberworth et al., 2013), regions of ISD are important in mediating protein interactions, then PTM of residues within regions of ISD might explain the basis for multiple layers of regulation (Nishi et al., 2013). Clearly this aspect of PKA deserves a greater focus.
Lys-N𝜀-Acetylation in Plants
While thousands of acetylated proteins are known in animals and yeast, when preparing this manuscript fewer than 400 plant proteins (∼500 sites) have been reported. These include 131 A. thaliana proteins (155 sites; Finkemeier et al., 2011; Wu et al., 2011), 97 Vitis vinifera berry proteins (138 sites; Melo-Braga et al., 2012), 121 proteins (190 sites) from developing Glycine max cotyledons (Smith-Hammond et al., 2014b), and 31 proteins (35 sites) from Solanum tuberosum mitochondrial proteins (Salvato et al., 2014). Even after results from work in progress are added; 90 A. thaliana mitochondrial proteins (174 sites; König et al., 2014a), and 358 Pisum sativum proteins (664 sites; Smith-Hammond et al., 2014a), the extant number of PKA proteins is small in comparison to P-proteins and only five species have been examined.
“Writers, Erasers, and Readers” of Lys-N𝜀-Acetylation
The occurrence of a mammalian enzyme specific for histone acetylation was first reported by Racey and Byvoet (1971), and the first histone acetyltransferase (HAT1) gene was cloned from yeast (Kleff et al., 1995). For over 40 years, protein Lys-acetylation has meant “histone Lys-acetylation.” The relatively recent discovery that transcription factors, structural proteins, metabolic enzymes, and a host of other non-histone proteins are Lys-N𝜀-acetylated has led to some confusion about specificity and terminology (Josling et al., 2012; Xing and Poirier, 2012). There are even instances where a “HAT” has been reported to acetylate histones, non-histone proteins, and even small molecules (Gu and Roeder, 1997). At this point, we favor erring on the side of generalization rather than claiming unsupported specificity. With this caveat, herein we refer to the “writers” of the PTM code (Muntean and Hess, 2009), protein Lys-N𝜀-acetyltransferases, as KATs, and the “erasers,” protein Lys-N𝜀-deacetylases as KDACs.
Based upon comparative sequence analyses, it has been concluded that there are four distinct KAT families: GNAT (GCN5-related N-terminal acetyltransferases); MYST [MOZ, Ybf2 (Sas3), Sas2, and Tip60; Sapountzi and Côté, 2011]; p300/CREB (CBP); and the nuclear receptor coactivator family (Roth et al., 2001). The first three families are widespread in eukaryotes [including A. thaliana (Pandey et al., 2002)], eubacteria, and archaea. The nuclear receptor coactivator family is thought to be absent from plants, fungi, and lower animals (Pandey et al., 2002). The different KAT families have distinct kinetic and catalytic mechanisms (Yan et al., 2002; Röttig and Steinbüchel, 2013), but are characteristically large, complex, multi-domain proteins (Roth et al., 2001; Figures 1A,C).
FIGURE 1. Phylogeny and domain architecture of the Arabidopsis thaliana KATs and KDACs. The phylogenetic trees of KATs (A) and KDACs (B) are shown with distinct groups colored differently. The KAT and KDAC sequences were aligned using MUSCLE, and trees were constructed in MEGA5 (http://www.megasoftware.net/) using maximum-likelihood method with 1000 bootstrap replicates (Values of more than 50% are shown). Scale bars indicate the number of substitutions per site. For the domain architecture of selected KATs (C) and KDACs (D), information was collected from PFAM (http://pfam.sanger.ac.uk/). Horizontal bars are proportional to sequence length (129–1919 residues for KATs, 203–682 residues for KDACs).
Perhaps not all protein acetylation is KAT-mediated? Similar to non-kinase based protein phosphorylation (autophosphorylation; Miernyk et al., 1992; Bae and Schlessinger, 2010), there have also been reported examples of autoacetylation (e.g., Yang et al., 2012). An intriguing recent publication described “widespread and enzyme-independent N𝜀-acetylation” of a number of proteins in the mitochondrial matrix (Wagner and Payne, 2013). This mechanism might also be extended to include peroxisomal and plastidial proteins? It is noteworthy that König et al. (2014b) found N𝜀-acetylated proteins, and a sirtuin, within A. thaliana mitochondria, but found no evidence of a KAT. Nor was there any evidence for a KAT within highly purified potato tuber mitochondria (Salvato et al., 2014).
If the KATs write the PTM code, then it is the deacetylases that are responsible for “erasing” it. The genome of the reference eudicot plant A. thaliana includes 18 genes encoding KDAC proteins; at least two members each of the RPD3-like (reduced potassium deficiency 3), HD-tuin, and sirtuin families (Hollender and Liu, 2008). Much like the KATs, the KDAC proteins display complex domain organization (Figures 1B,D), tissue-specific expression, and physiological functions. Members of the RPD3-like family are apparently present in all eukaryotes and have been the most widely studied KDACs (Murfett et al., 2001; Rossi et al., 2003). The HD-tuins appear to be present only in plants (Dangl et al., 2001; Luo et al., 2012) and have been the least studied.
The sirtuins (Silent Information Regulator 2 proteins) are a ubiquitous family of NAD++-dependent KDACs. It has been reported that mammalian cells contain seven sirtuin homologs (SIRT1–7) with diverse cellular localization [for example, some proteins of the Srt3, 4, and 5 families are mitochondrial (Huang et al., 2010; Sol et al., 2012)] and physiological functions. The sirtuins have been demonstrated to deacetylate a wide spectrum of clients (Finkemeier et al., 2011; Duan, 2013). König et al. (2014b) recently described a Srt2-orthologous protein localized within the matrix of A. thaliana mitochondria that deacetylates a specific cohort of mitochondrial client proteins. In contrast to the seven sirtuin genes present in mammalian genomes, A. thaliana and G. max (Glyma04g38730.1 and Glyma06g16260.1) have only two sirtuin-encoding genes. König et al. (2014b), however, detected more than seven alternative splicing variants of Atsrt2.
Like a molecular barcode, the information present in acetylated-Lys must be recognized and decoded by some sort of “reader.” Originally discovered as a component of histone-binding proteins, bromodomains (BRD) are conserved structural motifs (Figure 2B) that recognize and bind PKA (Dhalluin et al., 1999; Zeng and Zhou, 2002; Sanchez and Zhou, 2009). The term “BRD” comes from brahma, a regulatory protein in Drosophila melanogaster. The human genome encodes at least 46 BRD-proteins (each of which has between one and six BRDs) which have been sorted into eight classes (Filippakopoulos et al., 2012). The yeast genome encodes at least nine BRD-proteins (Sanchez and Zhou, 2009; Table 1).
FIGURE 2. Plant bromodomain proteins. (A) Phylogenetic analysis of the A. thaliana BRD-proteins. Distinct clusters (more than 50% bootstrap value) are colored differently in the tree. Sequences were aligned with MUSCLE and the tree was constructed using maximum likelihood method with 1000 bootstraps in MEGA5.The scale bar indicates the number of substitutions per site. (B) A representative BRD structure (AT3G54610.1) has four helices (Z, A, B, and C, each ∼110 residues in length). The ZA plus BC loops form the acetyl-Lys binding pocket. Unlike many mammalian and yeast proteins which include two or more BRD, most plant proteins have a single BRD. (C) An exception is the product of the Chlamydomonas reinhardtii (gi| 159485810) locus which includes three BRD (plus a PHD domain). The C. reinhardtii protein structure was predicted using I-TASSER (http://zhanglab.ccmb.med.umich.edu/I-TASSER).
In A. thaliana there are 29 BRD-proteins3, which can be separated into multiple groups (Figure 2A). The number of BRD-proteins varies considerably among plants, from as many as 57 in G. max to as few as nine in the red nanoalga Cyanidioschyzon merolae. There are only a few instances of plant proteins that include more than a single BRD (Figure 2C). The relationship between BRDs and Lys-acetylated client proteins (Figure 3) has been compared with the recognition and binding of O-phosphorylated client proteins with the SH2 domain or with 14-3-3 proteins (Yang, 2004a; de Boer et al., 2013). It is not yet clear if recognition and binding involve only acetylated-Lys residues or if these residues must be in a particular context/domain/environment.
FIGURE 3. A plant Lys-N𝜀-acetyl-protein interactome. (A) The BRD-proteins are presented as squares, KATs as arrow-heads, KDACs as diamonds, and PKA proteins as circles. Some of the PKA proteins [histone H4 (AT1G07660) and H3.2 (AT1G09200), GTP-binding EF-Tu (AT1G07930), and nuclear PGK1 (AT3G12780)] interact with BRD-proteins. However, either no interactions or only self-interactions are known for the majority of plant PKA-proteins (∼69%; B). The interactome information was collected from TAIR (http://www.arabidopsis.org/), AtPID (http://www.megabionet.org/atpid/), AtPIN (http://bioinfo.esalq.usp.br/), CCSB (http://interactome.dfci.harvard.edu/), and PAIR (http://www.cls.zju.edu.cn/pair/), and visualized using Cytoscape. Edge color (red is negative and green is positive) and width are proportional to gene co-expression (correlation coefficient) of genes based on GSE3011 from NCBI-GEO (http://www.ncbi.nlm.nih.gov/geo/). Nodes with no known interactions have not been included. Node size is proportional to the number of interactions. Node colors indicate subcellular protein localization, based on information from the TAIR database. (C) The percentage of currently known plant PKA-proteins in different subcellular localization categories is shown. The A. thaliana homologs of known PKA-proteins from other plants species were obtained using BLAST.
An additional possible PKA-reader is the YEATS domain, which is an evolutionarily conserved structural feature found in a variety of proteins in chromatin-modifying and transcription complexes ranging from yeast to humans (Schulze et al., 2010). Recently, the 3D structure of the YEATS domain from Yaf9 protein has been resolved, which is shown to have a region containing a shallow groove that might constitute aPKA-binding pocket (Zhang et al., 2011). The YEATS-domain containing protein family comprises more than 100 members in over 70 eukaryotic species (Schulze et al., 2009), including A. thaliana (Zacharaki et al., 2012), Oryza sativa, and Medicago truncatula.
It is worthwhile to note that not all effects of acetylation related downstream functions require BRD proteins, or any kind of mediating protein binding. Just as with O-phosphorylation, acetylation can exert direct effects on protein function or enzyme activity. For example, in Escherichia coli 85 PKA proteins (125 sites) have been identified (Yu et al., 2008), which seemingly function in the absence of any known BRD-proteins3.
Although the total number of P-sites and PKA sites are converging (Table 2), there remain large differences in the numbers of components of each system. While there are hundreds of Ser/Thr/Tyr-kinases (518/130/1085 in H. sapiens/S. cerevisiae/A. thaliana) and phosphatases (147/38/150), there are relatively few KATs (22/11/16) and KDACs (18/4/18). How can so few enzymes specifically regulate thousands of acetylation events in a cell (Kouzarides, 2000; Xing and Poirier, 2012; Xiong and Guan, 2012)? How does PKA accomplish specific and dynamic control with a relatively small complement of writers, erasers, and readers? Choi and Bakkenist (2013) recently suggested an intriguing possibility, that a PTM-binding protein might function in part by blocking access of a kinase and/or acetyltransferase to a different, potentially regulatory-site. Additionally, there are several reports that both KAT and KDAC are present in vivo as components of multi-protein complexes. Associated proteins can both stimulate (Berndsen et al., 2008) and inhibit (Kim et al., 2014) acetyl-transferase activity, and might also specify or modify client protein selectivity.
The known writers, erasers, and readers are all relatively large, multi-domain proteins (Figure 1). Multiple domains often function as protein interaction modules or scaffolds. This implies a profusion of protein interactions, multi-protein complexes, and networks (Figure 3). With the exception of BRD-proteins binding to PKA-histones (Dhalluin et al., 1999), most protein interactions have been addressed from a computational perspective (Lu et al., 2011). Protein interactions, signaling networks, etc. have now been well-studied as related to reversible O-phosphorylation (Trost et al., 2010; Newman et al., 2013; Uhrig et al., 2013), and it will be important to extend our understanding of PKA in this direction.
Lys-N𝜀-Acetylation and Metabolic Control
Lys-acetylation has rapidly become established as an important PTM involved with metabolic regulation in mammalian and microbial systems. In yeast and mammalian systems, virtually every enzyme of glycolysis, gluconeogenesis, the Krebs cycle, and urea, fatty acid, and glycogen metabolism has been reported to be acetylated (Wang et al., 2010; Zhao et al., 2010; Oliveira and Sauer, 2012). However, distinguishing bona fide regulatory sites among the thousands of PKA sites detected by contemporary high-resolution MS, and elucidating the mechanisms by which the modifications alter protein function remain a primary challenge.
The multiple mechanisms by which reversible-O-phosphorylation can directly affect protein activity includes effects on catalysis, the binding or release of substrates/products/regulators, protein complex formation, localization, turnover, etc. (Zhang et al., 2007; Trost et al., 2010; Oliveira et al., 2012; Fischer, 2013). As a rule of thumb, a PTM directly involved with control of protein function will be dynamic (e.g., have a shorter lifetime than the protein itself). A successful experimental strategy has been to use phosphatase or deacetylase inhibitors to treat tissues, cells, or organelles. Alternatively, it is also useful to employ knockout or knock-down phosphatase or deacetylase mutants. Finally, it is possible to use recombinant phosphatases or deacetylases to treat modified proteins coupled with direct measurements of activity.
While our current understanding of reversible PKA in plants is both preliminary and fragmentary, there is support for regulation in a few instances. Deacetylation of 3PGA-kinase by incubation with a heterologous sirtuin led to a significant increase in catalytic activity (Finkemeier et al., 2011). Likewise, RuBisCO LSU is Lys-acetylated in vivo, which reduces activity (Finkemeier et al., 2011). There are several reports that PKA inhibits/prevents/reverses protein interactions. Plastidial glycolytic/Calvin cycle enzymes form a complex/metabolon in the stroma (Graciet et al., 2004) mediated by PKA of GAPDH (Winkel, 2004). The position and acetylation of the specific GAPDH Lys residue are conserved in animal and bacterial sequences (Zhao et al., 2010; Zhang et al., 2013), as is, presumably, the role in mediating protein interactions. Another example of PKA affecting protein interactions is LHC subunit trimer formation and association with the thylakoid membranes (Wu et al., 2011). Unfortunately, thus far there have not yet been many instances where protein interactions have been directly related to protein activity.
Lys-N𝜀-Acetylation and Signaling
Changes in the complex pattern of histone PKA as a mechanism for controlling gene expression is being studied extensively (e.g., Cigliano et al., 2013; Gu et al., 2013), and is not addressed herein. Instead, we will focus on signaling between the nucleus and the cytoplasmic organelles. Plastids and mitochondria are specialized for both production and utilization of ATP and reduced pyridine nucleotides. While both plastids and mitochondria are “semi-autonomous,” the vast majority of proteins resident in these organelles is nuclear-encoded, translated in the cytoplasm, and imported post-translationally. The protein complements of these organelles are dynamic, and must be regulated to match cellular energy demands. Mechanisms for this regulation include sensing metabolic states and signaling the nucleus of changes.
Nuclear regulation of organellar protein composition and concentration is extensive and complex, and signaling is considered anterograde or forward. However, the idea of organellar signaling leading to changes in nuclear gene expression is a newer concept and is referred to as retrograde signaling (Ng et al., 2014). While our understanding of retrograde signaling is at a very early nascent stage, it has nevertheless been surprising that the best understood signaling molecules are simple metabolic intermediates (e.g., Czarnecki et al., 2012).
In a remarkable example of both flexibility and economy, acetyl-CoA is both a central metabolite and the substrate for PKA (Hartl and Finkemeier, 2012; Wellen and Thompson, 2012). Acetyl-CoA is a key component of major metabolic pathways in the cytoplasm, peroxisomes, plastids, and mitochondria, and examples of PKA of enzymes have been identified in all of these subcellular compartments (Finkemeier et al., 2011; Wu et al., 2011). Thus the metabolic status of these organelles, and of the cell in general, might easily be signaled to the nucleus. In order to maintain compartmental specificity, it would be necessary that the signaling molecules be either an up- or downstream component of the specific pathways. It has been recently suggested that citrate and possibly malate could be the signals for mitochondrial retrograde signaling (Finkemeier et al., 2013).
Another potentially important aspect of PKA and signaling involves subcellular dynamics. Both organelles and cytoplasmic protein complexes employ molecular motors to move along the cytoskeleton. Actin filaments, intermediate filaments, and microtubules are all subject to PKA, which can affect both intracellular trafficking and protein interactions (Zencheck et al., 2012). Cytoplasmic GAPDH can function either as part of a soluble, non-associated glycolytic pathway or as a component of a glycolytic metabolon at the mitochondrial outer membrane (Graham et al., 2007). All aspects of this dynamic microcompartmentation are potentially controlled by PKA, including association of the GAPDH subunits, association with the metabolon (several if not all components of which are subject to PKA) and positioning in the cell via association with actin (Wojtera-Kwiczor et al., 2013). It seems reasonable to assume that the metabolic signals sent to the nucleus would differ under each of these conditions. Finally, it is important to consider interactive and hierarchical interactions among PKA, other PTM, and other aspects of signaling (e.g., oxidative signaling).
Crosstalk Between KAC and Other PTM
As yet we have only a nascent understanding and appreciation of the complexity of various interacting PTM (van Noort et al., 2012; Rao et al., 2013). Decoding the various levels of crosstalk patterns is critical to appreciating the role of PTM in protein regulation, signaling and plant development, and controlling gene expression. The possibilities are manifold, and include multiple instances of the same PTM at different sites [e.g., a priming modification at site A is necessary for subsequent phosphorylation or acetylation of site B (Lu et al., 2011; Woodsmith et al., 2013)], to hierarchical responses to multiple PTM of the same site (Minguez et al., 2012; Zauner et al., 2012), differential responses to multiple different PTM at different sites within the same protein (cis-crosstalk), and ultimately to crosstalk between PTM of different proteins (trans-crosstalk).
An instance of PTM cis-crosstalk is the phosphorylation of Ser10 residue of histone H3, which subsequently leads to acetylation of Lys14 residue (Roth et al., 2001). The amino- and carboxy-terminal tails of the core histones are decorated with multisite-modifications including methylation, acetylation, phosphorylation, ADP-ribosylation, ubiquitination, and sumoylation (Lau and Cheung, 2013). The “histone code” was introduced as an explanation of how combinatorial systems of histone PTM regulate transcription (Jenuwein and Allis, 2001). The histone-code hypothesis was subsequently modified, extended to include transcription factors, and referred to as the PTM code (Benayoun and Veitia, 2009). With the widespread occurrence of a plethora of PTM, we propose further extension of the “PTM code” to include regulation, signaling, and development, as well as control of gene expression.
While reversible O-phosphorylation has received the most attention of any PTM, there are an increasing number of reports of Lys-N𝜀-acetylation. Based upon results from analyses of mammalian systems, it is reasonable to expect discovery of many additional sites of PKA of plant and microbial proteins. It will be important to shift research emphasis from descriptive to quantitative and to determine the stoichiometry and dynamics of PKA rather than only sites of acetylation. In O-phosphorylation there are large differences between the number of kinases/client proteins, and phosphatases. In contrast, in PKA there are large differences in the numbers of acetyltransferases/deacetylases, and client proteins. The bases for these differences are obscure. Network analyses of the writers, erasers, and readers of O-phosphorylation are relatively well-developed in contrast to corresponding network analyses of N-acetylation. In the absence of any apparent candidates for KAT in mitochondria, plastids or peroxisomes, how and where are proteins resident within these organelles Lys-acetylated? A very complex and multifaceted question addresses the nature and extent of PKA crosstalk with other PTM. Finally, it will be important to achieve an improved understanding of the roles of PKA in the long-range signaling pathways.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Research in the Thelen Lab is supported by a grant from the National Science Foundation Plant Genome Research Program Grant No. DBI-0604439.
BRD, bromodomains; HAT, histone acetyltransferase; ISD, intrinsic sequence disorder; KAC, N𝜀-acetyl-Lys; KAT, protein Lys-N𝜀-acetyltransferases; KDAC, protein Lys-N𝜀-deacetylases as KDACs; MS, mass spectrometry; PKA, protein Lys-N𝜀-acetylation; PTM, post-translational modifications.
Bae, J. H., and Schlessinger, J. (2010). Asymmetric tyrosine kinase arrangements in activation or autophosphorylation of receptor tyrosine kinases. Mol. Cells 29, 443–448. doi: 10.1007/s10059-010-0080-5
Benayoun, B. A., and Veitia, R. A. (2009). A post-translational modification code for transcription factors: sorting through a sea of signals. Trends Cell Biol. 19, 189–197. doi: 10.1016/j.tcb.2009.02.003
Berndsen, C. E., Tsubota, T., Lindner, S. E., Lee, S., Holton, J. M., Kaufman, P. D.,et al. (2008). Molecular functions of the histone acetyltransferase chaperone complex Rtt109-Vps75. Nat. Struct. Mol. Biol. 15, 948–956. doi: 10.1038/nsmb.1459
Breitkreutz, A., Choi, H., Sharom, J. R., Boucher, L., Neduva, V., Larsen, B.,et al. (2010). A global protein kinase and phosphatase interaction network in yeast. Science 328, 1043–1046. doi: 10.1126/science.1176495
Choudhary, C., Kumar, C., Gnad, F., Nielsen, M. L., Rehman, M., Walther, T. C.,et al. (2009). Lysine acetylation targets protein complexes and co-regulates major cellular functions. Science 325, 834–840. doi: 10.1126/science.1175371
Cigliano, R. A., Sanseverino, W., Cremona, G., Ercolano, M. R., Conicella, C., and Consiglio, F. M. (2013). Genome-wide analysis of histone modifiers in tomato: gaining an insight into their developmental roles. BMC Genomics 14:57. doi: 10.1186/1471-2164-14-57
Cumberworth, A., Lamour, G., Babu, M., and Gsponer, J. (2013). Promiscuity as a functional trait: intrinsically disordered regions as central players of interactomes. Biochem. J. 454, 361–369. doi: 10.1042/BJ20130545
Czarnecki, O., Gläßer, C., Chen, J.-G., Mayer, K. F. X., and Grimm, B. (2012). Evidence for a contribution of ALA synthesis to plastid-to-nucleus signaling. Front. Plant Sci. 3:236. doi: 10.3389/fpls.2012.00236
Earley, K. W., Shook, M. S., Brower-Toland, B., Hicks, L., and Pikaard, C. S. (2007). In vitro specificities of Arabidopsis co-activator histone acetyltransferases: implications for histone hyperacetylation in gene activation. Plant J. 52, 615–626. doi: 10.1111/j.1365-313X.2007.03264.x
Filippakopoulos, P., Picaud, S., Mangos, M., Keates, T., Lambert, J. P., Barsyte-Lovejoy, D.,et al. (2012). Histone recognition and large-scale structural analysis of the human bromodomain family. Cell 149, 214–231. doi: 10.1016/j.cell.2012.02.013
Finkemeier, I., König, A. C., Heard, W., Nunes-Nesi, A., Pham, P. A., Leister, D.,et al. (2013). Transcriptomic analysis of the role of carboxylic acids in metabolite signaling in Arabidopsis leaves. Plant Physiol. 162, 239–253. doi: 10.1104/pp.113.214114
Finkemeier, I., Laxa, M., Miguet, L., Howden, A. J., and Sweetlove, L. J. (2011). Proteins of diverse function and subcellular location are lysine acetylated in Arabidopsis. Plant Physiol. 155, 1779–1790. doi: 10.1104/pp.110.171595
Graciet, E., Lebreton, S., and Gontero, B. (2004). Emergence of new regulatory mechanisms in the Benson-Calvin pathway via protein-protein interactions: a glyceraldehyde-3-phosphate dehydrogenase/CP12/phosphoribulokinase complex. J. Exp. Bot. 55, 1245–1254. doi: 10.1093/jxb/erh107
Graham, J. W. A., Williams, T. C. R., Morgan, M., Fernie, A. R., Ratcliffe, R. G., and Sweetlove, L. J. (2007). Glycolytic enzymes associate dynamically with mitochondria in response to respiratory demand and support substrate channeling. Plant Cell 19, 3723–3738. doi: 10.1105/tpc.107.053371
Gu, X., Wang, Y., and He, Y. (2013). Photoperiodic regulation of flowering time through periodic histone deacetylation of the florigen gene FT. PLoS Biol. 11:e1001649. doi: 10.1371/journal.pbio.1001649
Henriksen, P., Wagner, S. A., Weinert, B. T., Sharma, S., Bacinskaja, G., Rehman, M.,et al. (2012). Proteome-wide analysis of lysine acetylation suggests its broad regulatory scope in Saccharomyces cerevisiae. Mol. Cell. Proteomics 11, 1510–1522. doi: 10.1074/mcp.M112.017251
Hrabak, E. M., Chan, C. W., Gribskov, M., Harper, J. F., Choi, J. H., Halford, J. F.,et al. (2003). The Arabidopsis CDPK-SnRK superfamily of protein kinases. Plant Physiol. 132, 666–680. doi: 10.1104/pp.102.011999
Iakoucheva, L. M., Radivojac, P., Brown, C. J., O’Connor, T. R., Sikes, J. G., Obradovic, Z.,et al. (2004). The importance of intrinsic disorder for protein phosphorylation. Nucleic Acids Res. 32, 1037–1049. doi: 10.1093/nar/gkh253
Jeffers, V., and Sullivan, W. J. Jr. (2012). Lysine acetylation is widespread on proteins of diverse function and localization in the protozoan parasite Toxoplasma gondii. Eukaryot. Cell 11, 735–742. doi: 10.1128/EC.00088-12
Kerk, D., Templeton, G., and Moorhead, G. B. (2008). Evolutionary radiation pattern of novel protein phosphatases revealed by analysis of protein data from the completely sequenced genomes of humans, green algae, and higher plants. Plant Physiol. 146, 351–367. doi: 10.1104/pp.107.111393
Khoury, G. A., Baliban, R. C., and Floudasa, C. A. (2011). Proteome-wide post-translational modification statistics: frequency analysis and curation of the swiss-prot database. Sci. Rep. 1:90. doi: 10.1038/srep00090
Kim, D., Yu, B. J., Kim, J. A., Lee, Y. J., Choi, S. G., Kang, S.,et al. (2013). The acetylproteome of gram-positive model bacterium Bacillus subtilis. Proteomics 13, 1726–1736. doi: 10.1002/pmic.201200001
Kim, K. B., Kim, D. W., Park, J. W., Jeon, Y. J., Kim, D., Rhee, S.,et al. (2014). Inhibition of Ku70 acetylation by INHAT subunit SET/TAF-Iβ regulates Ku70-mediated DNA damage response. Cell. Mol. Life Sci. 71, 2731–2745. doi: 10.1007/s00018-013-1525-8
Kim, S. C., Sprung, R., Chen, Y., Xu, Y., Ball, H., Pei, J.,et al. (2006). Substrate and functional diversity of lysine acetylation revealed by a proteomics survey. Mol. Cell 23, 607–618. doi: 10.1016/j.molcel.2006.06.026
Kleff, S., Andrulis, E. D., Anderson, C. W., and Sternglanz, R. (1995). Identification of a gene encoding a yeast histone H4 acetyltransferase. J. Biol. Chem. 270, 24674–24677. doi: 10.1074/jbc.270.42.24674
König, A.-C., Hartl, M., Boersemab, P. J., Manna, M., and Finkemeiera, I. (2014a). The mitochondrial lysine acetylome of Arabidopsis. Mitochondrion doi: 10.1016/j.mito.2014.03.004 [Epub ahead of print].
König, A.-C., Hartl, M., Pham, P. A., Laxa, M., Boersema, P. J., Orwat, A.,et al. (2014b). The Arabidopsis class II sirtuin is a lysine deacetylase and interacts with mitochondrial energy metabolism. Plant Physiol. 164, 1401–1414. doi: 10.1104/pp.113.232496
Lau, P. N., and Cheung, P. (2013). Elucidating combinatorial histone modifications and crosstalk by coupling histone-modifying enzyme with biotin ligase activity. Nucleic Acids Res. 41:e49. doi: 10.1093/nar/gks1247
Liu, B. A., Shah, E., Jablonowski, K., Stergachis, A., Engelmann, B., and Nash, P. D. (2011). The SH2 domain-containing proteins in 21 species establish the provenance and scope of phosphotyrosine signaling in eukaryotes. Sci. Signal. 4:ra83. doi: 10.1126/scisignal.2002105
Lu, Z., Cheng, Z., Zhao, Y., and Volchenboum, S. L. (2011). Bioinformatic analysis and post-translational modification crosstalk prediction of lysine acetylation. PLoS ONE 6:e28228. doi: 10.1371/journal.pone.0028228
Lundby, A., Lage, K., Weinert, B. T., Bekker-Jensen, D. B., Secher, A., Skovgaard, T.,et al. (2012a). Proteomic analysis of lysine acetylation sites in rat tissues reveals organ specificity and subcellular patterns. Cell Rep. 2, 419–431. doi: 10.1016/j.celrep.2012.07.006
Lundby, A., Secher, A., Lage, K., Nordsborg, N. B., Dmytriyev, A., Lundby, C.,et al. (2012b). Quantitative maps of protein phosphorylation sites across 14 different rat organs and tissues. Nat. Commun. 3:876. doi: 10.1038/ncomms1871
Luo, M., Wang, Y. Y., Liu, X., Yang, S., Lu, Q., Cui, Y.,et al. (2012). HD2C interacts with HDA6 and is involved in ABA and salt stress response in Arabidopsis. J. Exp. Bot. 63, 3297–3306. doi: 10.1093/jxb/ers059
Martinez, A., Traverso, J. A., Valot, B., Ferro, M., Espagne, C., Ephritikhine, G.,et al. (2008). Extent of N-terminal modifications in cytosolic proteins from eukaryotes. Proteomics 8, 2809–2831. doi: 10.1002/pmic.200701191
Melo-Braga, M. N., Verano-Braga, T., León, I. R., Antonacci, D., Nogueira, F. C., Thelen, J. J.,et al. (2012). Modulation of protein phosphorylation, N-glycosylation and Lys-acetylation in grape (Vitis vinifera) mesocarp and exocarp owing to Lobesia botrana infection. Mol. Cell. Proteomics 11, 945–956. doi: 10.1074/mcp.M112.020214
Miao, J., Lawrence, M., Jeffers, V., Zhao, F., Parker, D., Ge, Y.,et al. (2013). Extensive lysine acetylation occurs in evolutionarily conserved metabolic pathways and parasite-specific functions during Plasmodium falciparum intraerythrocytic development. Mol. Microbiol. 89, 660–675. doi: 10.1111/mmi.12303
Minguez, P., Parca, L., Diella, F., Mende, D. R., Kumar, R., Helmer-Citterich, M.,et al. (2012). Deciphering a global network of functionally associated post-translational modifications. Mol. Syst. Biol. 8:599. doi: 10.1038/msb.2012.31
Murfett, J., Wang, X. J., Hagen, G., and Guilfoyle, T. J. (2001). Identification of Arabidopsis histone deacetylase HDA6 mutants that affect transgene expression. Plant Cell 13, 1047–1061. doi: 10.1105/tpc.13.5.1047
Ng, S., De Clercq, I., Van Aken, O., Law, S. R., Ivanova, A., Willems, P.,et al. (2014). Anterograde and retrograde regulation of nuclear genes encoding mitochondrial proteins during growth, development, and stress. Mol. Plant. 7, 1075–1093. doi: 10.1093/mp/ssu037
Nishi, H., Fong, J. H., Chang, C., Teichmann, S. A., and Panchenko, A. R. (2013). Regulation of protein-protein binding by coupling between phosphorylation and intrinsic disorder: analysis of human protein complexes. Mol. Biosyst. 9, 1620–1626. doi: 10.1039/c3mb25514j
Oliveira, A. P., Ludwig, C., Picotti, P., Kogadeeva, M., Aebersold, R., and Sauer, U. (2012). Regulation of yeast central metabolism by enzyme phosphorylation. Mol. Syst. Biol. 8:623. doi: 10.1038/msb.2012.55
Oliveira, A. P., and Sauer, U. (2012). The importance of post-translational modifications in regulating Saccharomyces cerevisiae metabolism. FEMS Yeast Res. 12, 104–117. doi: 10.1111/j.1567-1364.2011.00765.x
Pandey, R., Müller, A., Napoli, C. A., Selinger, D. A., Pikaard, C. S., Richards, E. J.,et al. (2002). Analysis of histone acetyltransferase and histone deacetylase families of Arabidopsis thaliana suggests functional diversification of chromatin modification among multicellular eukaryotes. Nucleic Acids Res. 30, 5036–5055. doi: 10.1093/nar/gkf660
Piperno, G., and Fuller, M. T. (1985). Monoclonal antibodies specific for an acetylated form of α-tubulin recognize the antigen in cilia and flagella from a variety of organisms. J. Cell Biol. 101, 2085–2094. doi: 10.1083/jcb.101.6.2085
Racey, L. A., and Byvoet, P. (1971). Histone acetyltransferase in chromatin: evidence for in vitro enzymatic transfer of acetate from acetyl-coenzyme a to histones. Exp. Cell Res. 64, 366–370. doi: 10.1016/0014-4827(71)90089-9
Rao, R. S. P., Thelen, J. J., and Miernyk, J. A. (2014). Is Lys-N𝜀-acetylation the next big thing in post-translational modifications? Trends Plant Sci. doi: 10.1016/j.tplants.2014.05.001 [Epub ahead of print].
Rao, R. S. P., Xu, D., Thelen, J. J., and Miernyk, J. A. (2013). Circles within circles: crosstalk between protein Ser/Thr/Tyr-phosphorylation and Met oxidation. BMC Bioinformatics 14(Suppl. 14):S14. doi: 10.1186/1471-2105-14-S14-S14
Rossi, V., Locatelli, S., Lanzanova, C., Boniotti, M. B., Varotto, S., Pipal, A.,et al. (2003). A maize histone deacetylase and retinoblastoma-related protein physically interact and cooperate in repressing gene transcription. Plant Mol. Biol. 51, 401–413. doi: 10.1023/A:1022090916446
Salvato, F., Havelund, J. F., Chen, M., Rao, R. S. P., Rogowska-Wrzesinska, A., Jensen, O. N.,et al. (2014). The potato tuber mitochondrial proteome. Plant Physiol. 164, 637–653. doi: 10.1104/pp.113.229054
Schulze, J. M., Wang, A. Y., and Kobor, M. S. (2009). YEATS domain proteins: a diverse family with many links to chromatin modification and transcription. Biochem. Cell Biol. 87, 65–75. doi: 10.1139/O08-111
Smith-Hammond, C. L., Swatek, K., Johnston, M. L., Thelen, J. J., and Miernyk, J. A. (2014b). Initial description of the developing soybean seed protein-Lys-N𝜀-acetylome. J. Proteomics 96, 56–66. doi: 10.1016/j.jprot.2013.10.038
Sol, E. M., Wagner, S. A., Weinert, B. T., Kumar, A., Kim, H. S., Deng, C. X.,et al. (2012). Proteomic investigations of lysine acetylation identify diverse substrates of mitochondrial deacetylase sirt3. PLoS ONE 7:e50545. doi: 10.1371/journal.pone.0050545
Tinti, M., Johnson, C., Toth, R., Ferrier, D. E., and Mackintosh, C. (2012). Evolution of signal multiplexing by 14-3-3-binding 2R-ohnologue protein families in the vertebrates. Open Biol. 2:120103. doi: 10.1098/rsob.120103
Tyler, R. C., Bitto, E., Berndsen, C. E., Bingman, C. A., Singh, S., Lee, M. S.,et al. (2006). Structure of Arabidopsis thaliana At1g77540 protein, a minimal acetyltransferase from the COG2388 family. Biochemistry 45, 14325–14336. doi: 10.1021/bi0612059
Uhrig, R. G., Labandera, A.-M., and Moorhead, G. B. (2013). Arabidopsis PPP family of serine/threonine protein phosphatases: many targets but few engines. Trends Plant Sci. 18, 505–513. doi: 10.1016/j.tplants.2013.05.004
van Noort, V., Seebacher, J., Bader, S., Mohammed, S., Vonkova, I., Betts, M. J.,et al. (2012). Cross-talk between phosphorylation and lysine acetylation in a genome-reduced bacterium. Mol. Syst. Biol. 8:571. doi: 10.1038/msb.2012.4
Wagner, G. R., and Payne, R. M. (2013). Widespread and enzyme-independent N𝜀-acetylation and N𝜀-succinylation in the chemical conditions of the mitochondrial matrix. J. Biol. Chem. 288, 29036–29045. doi: 10.1074/jbc.M113.486753
Wang, Q., Zhang, Y., Yang, C., Xiong, H., Lin, Y., Yao, J.,et al. (2010). Acetylation of metabolic enzymes coordinates carbon source utilization and metabolic flux. Science 327, 1004–1007. doi: 10.1126/science.1179687
Wojtera-Kwiczor, J., Groß, F., Leffers, H. M., Kang, M., Schneider, M., and Scheibe, R. (2013). Transfer of a redox-signal through the cytosol by redox-dependent microcompartmentation of glycolytic enzymes at mitochondria, and actin cytoskeleton. Front. Plant. Sci. 9:284. doi: 10.3389/fpls.2012.00284
Wu, X., Oh, M.-H., Schwarz, E. M., Larue, C. T., Sivaguru, M., Imai, B. S.,et al. (2011). Lysine acetylation is a widespread protein modification for diverse proteins in Arabidopsis. Plant Physiol. 155, 1769–1778. doi: 10.1104/pp.110.165852
Xue, B., Jeffers, V., Sullivan, W. J., and Uversky, V. N. (2013). Protein intrinsic disorder in the acetylome of intracellular and extracellular Toxoplasma gondii. Mol. Biosyst. 9, 645–657. doi: 10.1039/c3mb25517d
Yan, Y., Harper, S., Speicher, D. W., and Marmorstein, R. (2002). The catalytic mechanism of the ESA1 histone acetyltransferase involves a self-acetylated intermediate. Nat. Struct. Biol. 9, 862–869. doi: 10.1038/nsb849
Zacharaki, V., Benhamed, M., Poulios, S., Latrasse, D., Papoutsoglou, P., Delarue, M.,et al. (2012). The Arabidopsis ortholog of the YEATS domain containing protein YAF9a regulates flowering by controlling H4 acetylation levels at the FLC locus. Plant Sci. 196, 44–52. doi: 10.1016/j.plantsci.2012.07.010
Zhang, J., Sprung, R., Pei, J., Tan, X., Kim, S., Zhu, H.,et al. (2009). Lysine acetylation is a highly abundant and evolutionarily conserved modification in Escherichia coli. Mol. Cell. Proteomics 8, 215–225. doi: 10.1074/mcp.M800187-MCP200
Zhang, Q., Chakravarty, S., Ghersi, D., Zeng, L., Plotnikov, A. N., Sanchez, R.,et al. (2010). Biochemical profiling of histone binding selectivity of the yeast bromodomain family. PLoS ONE 5:e8903. doi: 10.1371/journal.pone.0008903
Zhang, W., Zhang, J., Zhang, X., Xu, C., and Tu, X. (2011). Solution structure of the Taf14 YEATS domain and its roles in cell growth of Saccharomyces cerevisiae. Biochem. J. 436, 83–90. doi: 10.1042/BJ20110004
Keywords: acetyltransferase, bromodomain, deacetylase, interactome, post-translational modification, regulation, signaling
Citation: Rao RSP, Thelen JJ and Miernyk JA (2014) In silico analysis of protein Lys-N𝜀-acetylation in plants. Front. Plant Sci. 5:381. doi: 10.3389/fpls.2014.00381
Received: 18 March 2014; Accepted: 17 July 2014;
Published online: 04 August 2014.
Edited by:Steven Carl Huber, United States Department of Agriculture – Agricultural Research Service, USA
Reviewed by:Ingo Dreyer, Universidad Politecnica de Madrid, Spain
Iris Finkemeier, Max Planck Institute for Plant Breeding Research, Germany
Copyright © 2014 Rao, Thelen and Miernyk. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jan A. Miernyk, Division of Biochemistry, University of Missouri, 102 Curtis Hall, Columbia, MO 65211, USA e-mail: email@example.com