Dimethyl-Labeling-Based Quantification of the Lysine Acetylome and Proteome of Plants

Coordinated metabolic adaptation of Arabidopsis thaliana to high light

In plants, exposure to high light irradiation induces various stress responses, which entail complex metabolic rearrangements. To explore these dynamics, we conducted time‐course experiments spanning 2 min to 72 h with Arabidopsis thaliana under high and control light. Comparative metabolomics, transcriptomics, redox proteomics, and stable isotope labeling on leaf rosettes identified a series of synchronous and successive responses that provide a deeper insight into well‐orchestrated mechanisms contributing to high‐light acclimation. We observed transient transcriptome downregulation related to light harvesting and electron flow before the profound remodeling of the photosynthetic apparatus. Throughout the entire time course, redox homeostasis is tightly balanced between downregulation of production and enhanced transformation of NADPH accompanied by redistribution of reducing equivalents across several subcellular compartments. In both light conditions, C4 acids such as malate and fumarate are produced via anaplerosis. In carbon units, their accumulation in vacuoles surpasses plastidic levels of starch and intensifies notably under high light. In parallel, citrate synthesis from pyruvate is significantly hindered diurnally. Isotopic labeling in 2‐oxoglutarate and glutamate suggests a moderate de novo synthesis of C5 acids from a vacuolar citrate reservoir during the light phase while they are largely renewed during the night. In the absence of a diurnal clockwise flow through the tricarboxylic acid (TCA) cycle, increased oxidation of photorespiratory glycine takes over as a source of reductants to fuel mitochondrial ATP production. These findings, along with previous research, contribute to a model integrating redox balance and linking increased carbon assimilation and nitrogen metabolism, especially in the context of an incomplete TCA cycle.

Nα-acetyltransferase NAA50 mediates plant immunity independent of the Nα-acetyltransferase A complex

In humans and plants, 40% of the proteome is cotranslationally acetylated at the N-terminus by a single Nα-acetyltransferase (Nat) termed NatA. The core NatA complex is comprised of the catalytic subunit Nα-acetyltransferase 10 (NAA10) and the ribosome-anchoring subunit NAA15. The regulatory subunit Huntingtin Yeast Partner K (HYPK) and the acetyltransferase NAA50 join this complex in humans. Even though both are conserved in Arabidopsis (Arabidopsis thaliana), only AtHYPK is known to interact with AtNatA. Here we uncover the AtNAA50 interactome and provide evidence for the association of AtNAA50 with NatA at ribosomes. In agreement with the latter, a split-luciferase approach demonstrated close proximity of AtNAA50 and AtNatA in planta. Despite their interaction, AtNatA/HYPK and AtNAA50 exerted different functions in vivo. Unlike NatA/HYPK, AtNAA50 did not modulate drought tolerance or promote protein stability. Instead, transcriptome and proteome analyses of a novel AtNAA50-depleted mutant (amiNAA50) implied that AtNAA50 negatively regulates plant immunity. Indeed, amiNAA50 plants exhibited enhanced resistance to oomycetes and bacterial pathogens. In contrast to what was observed in NatA-depleted mutants, this resistance was independent of an accumulation of salicylic acid prior to pathogen exposure. Our study dissects the in vivo function of the NatA interactors HYPK and NAA50 and uncovers NatA-independent roles for NAA50 in plants.

The interplay of post-translational protein modifications in Arabidopsis leaves during photosynthesis induction

Diurnal dark to light transition causes profound physiological changes in plant metabolism. These changes require distinct modes of regulation as a unique feature of photosynthetic lifestyle. The activities of several key metabolic enzymes are regulated by light-dependent post-translational modifications (PTM) and have been studied at depth at the level of individual proteins. In contrast, a global picture of the light-dependent PTMome dynamics is lacking, leaving the response of a large proportion of cellular function undefined. Here, we investigated the light-dependent metabolome and proteome changes in Arabidopsis rosettes in a time resolved manner to dissect their kinetic interplay, focusing on phosphorylation, lysine acetylation, and cysteine-based redox switches. Of over 24 000 PTM sites that were detected, more than 1700 were changed during the transition from dark to light. While the first changes, as measured 5 min after onset of illumination, occurred mainly in the chloroplasts, PTM changes at proteins in other compartments coincided with the full activation of the Calvin-Benson cycle and the synthesis of sugars at later timepoints. Our data reveal connections between metabolism and PTM-based regulation throughout the cell. The comprehensive multiome profiling analysis provides unique insight into the extent by which photosynthesis reprograms global cell function and adds a powerful resource for the dissection of diverse cellular processes in the context of photosynthetic function.

Mitochondrial alternative NADH dehydrogenases NDA1 and NDA2 promote survival of reoxygenation stress in Arabidopsis by safeguarding photosynthesis and limiting ROS generation

Plant submergence stress is a growing problem for global agriculture. During desubmergence, rising O₂ concentrations meet a highly reduced mitochondrial electron transport chain (mETC) in the cells. This combination favors the generation of reactive oxygen species (ROS) by the mitochondria, which at excess can cause damage. The cellular mechanisms underpinning the management of reoxygenation stress are not fully understood. We investigated the role of alternative NADH dehydrogenases (NDs), as components of the alternative mETC in Arabidopsis, in anoxia-reoxygenation stress management. Simultaneous loss of the matrix-facing NDs, NDA1 and NDA2, decreased seedling survival after reoxygenation, while overexpression increased survival. The absence of NDAs led to reduced maximum potential quantum efficiency of photosystem II linking the alternative mETC to photosynthetic function in the chloroplast. NDA1 and NDA2 were induced upon reoxygenation, and transcriptional activation of NDA1 was controlled by the transcription factors ANAC016 and ANAC017 that bind to the mitochondrial dysfunction motif (MDM) in the NDA1 promoter. The absence of NDA1 and NDA2 did not alter recovery of cytosolic ATP levels and NADH : NAD⁺ ratio at reoxygenation. Rather, the absence of NDAs led to elevated ROS production, while their overexpression limited ROS. Our observations indicate that the control of ROS formation by the alternative mETC is important for photosynthetic recovery and for seedling survival of anoxia-reoxygenation stress.

Proteome-wide lysine acetylation profiling to investigate the involvement of histone deacetylase HDA5 in the salt stress response of Arabidopsis leaves

Post-translational modifications (PTMs) of proteins play important roles in the acclimation of plants to environmental stress. Lysine acetylation is a dynamic and reversible PTM, which can be removed by histone deacetylases. Here we investigated the role of lysine acetylation in the response of Arabidopsis leaves to 1 week of salt stress. A quantitative mass spectrometry analysis revealed an increase in lysine acetylation of several proteins from cytosol and plastids, which was accompanied by altered histone deacetylase activities in the salt-treated leaves. While activities of HDA14 and HDA15 were decreased upon salt stress, HDA5 showed a mild and HDA19 a strong increase in activity. Since HDA5 is a cytosolic-nuclear enzyme from the class II histone deacetylase family with yet unknown protein substrates, we performed a lysine acetylome analysis on hda5 mutants and characterized its substrate proteins. Next to histone H2B, the salt stress-responsive transcription factor GT2L and the dehydration-related protein ERD7 were identified as HDA5 substrates. In addition, in protein-protein interaction studies, HDA18 was discovered, among other interacting proteins, to work in a complex together with HDA5. Altogether, this study revealed the substrate proteins of HDA5 and identified new lysine acetylation sites which are hyperacetylated upon salt stress. The identification of specific histone deacetylase substrate proteins, apart from histones, will be important to unravel the acclimation response of Arabidopsis to salt stress and their role in plant physiology.

Light acclimation interacts with thylakoid ion transport to govern the dynamics of photosynthesis in Arabidopsis

Understanding photosynthesis in natural, dynamic light environments requires knowledge of long-term acclimation, short-term responses, and their mechanistic interactions. To approach the latter, we systematically determined and characterized light-environmental effects on thylakoid ion transport-mediated short-term responses during light fluctuations. For this, Arabidopsis thaliana wild-type and mutants of the Cl^- channel VCCN1 and the K⁺ exchange antiporter KEA3 were grown under eight different light environments and characterized for photosynthesis-associated parameters and factors in steady state and during light fluctuations. For a detailed characterization of selected light conditions, we monitored ion flux dynamics at unprecedented high temporal resolution by a modified spectroscopy approach. Our analyses reveal that daily light intensity sculpts photosynthetic capacity as a main acclimatory driver with positive and negative effects on the function of KEA3 and VCCN1 during high-light phases, respectively. Fluctuations in light intensity boost the accumulation of the photoprotective pigment zeaxanthin (Zx). We show that KEA3 suppresses Zx accumulation during the day, which together with its direct proton transport activity accelerates photosynthetic transition to lower light intensities. In summary, both light-environment factors, intensity and variability, modulate the function of thylakoid ion transport in dynamic photosynthesis with distinct effects on lumen pH, Zx accumulation, photoprotection, and photosynthetic efficiency.

A Proteomic Approach for the Quantification of Posttranslational Protein Lysine Acetylation in Candida albicans

Candida albicans is a normal component of the human microflora that colonizes mucosal/epithelial surfaces and the gastrointestinal tract as a commensal organism. However, in an immunocompromised host, it can cause life-threatening infections of high mortality and morbidity. Virulence as well as antifungal drug resistance of C. albicans is often regulated by posttranslational modifications (PTM) of proteins via lysine acetylation by lysine acetyltransferases. Here, we report an experimental approach using tandem mass tag (TMT) labeling for the detection and quantification of lysine acetylation of histone and nonhistone proteins in C. albicans.

Dynamic light- and acetate-dependent regulation of the proteome and lysine acetylome of Chlamydomonas

The green alga Chlamydomonas reinhardtii is one of the most studied microorganisms in photosynthesis research and for biofuel production. A detailed understanding of the dynamic regulation of its carbon metabolism is therefore crucial for metabolic engineering. Post-translational modifications can act as molecular switches for the control of protein function. Acetylation of the ɛ-amino group of lysine residues is a dynamic modification on proteins across organisms from all kingdoms. Here, we performed mass spectrometry-based profiling of proteome and lysine acetylome dynamics in Chlamydomonas under varying growth conditions. Chlamydomonas liquid cultures were transferred from mixotrophic (light and acetate as carbon source) to heterotrophic (dark and acetate) or photoautotrophic (light only) growth conditions for 30 h before harvest. In total, 5863 protein groups and 1376 lysine acetylation sites were identified with a false discovery rate of <1%. As a major result of this study, our data show that dynamic changes in the abundance of lysine acetylation on various enzymes involved in photosynthesis, fatty acid metabolism, and the glyoxylate cycle are dependent on acetate and light. Exemplary determination of acetylation site stoichiometries revealed particularly high occupancy levels on K175 of the large subunit of RuBisCO and K99 and K340 of peroxisomal citrate synthase under heterotrophic conditions. The lysine acetylation stoichiometries correlated with increased activities of cellular citrate synthase and the known inactivation of the Calvin-Benson cycle under heterotrophic conditions. In conclusion, the newly identified dynamic lysine acetylation sites may be of great value for genetic engineering of metabolic pathways in Chlamydomonas.

Acetylation of conserved lysines fine-tunes mitochondrial malate dehydrogenase activity in land plants

Plants need to rapidly and flexibly adjust their metabolism to changes of their immediate environment. Since this necessity results from the sessile lifestyle of land plants, key mechanisms for orchestrating central metabolic acclimation are likely to have evolved early. Here, we explore the role of lysine acetylation as a post-translational modification to directly modulate metabolic function. We generated a lysine acetylome of the moss Physcomitrium patens and identified 638 lysine acetylation sites, mostly found in mitochondrial and plastidial proteins. A comparison with available angiosperm data pinpointed lysine acetylation as a conserved regulatory strategy in land plants. Focusing on mitochondrial central metabolism, we functionally analyzed acetylation of mitochondrial malate dehydrogenase (mMDH), which acts as a hub of plant metabolic flexibility. In P. patens mMDH1, we detected a single acetylated lysine located next to one of the four acetylation sites detected in Arabidopsis thaliana mMDH1. We assessed the kinetic behavior of recombinant A. thaliana and P. patens mMDH1 with site-specifically incorporated acetyl-lysines. Acetylation of A. thaliana mMDH1 at K169, K170, and K334 decreases its oxaloacetate reduction activity, while acetylation of P. patens mMDH1 at K172 increases this activity. We found modulation of the malate oxidation activity only in A. thaliana mMDH1, where acetylation of K334 strongly activated it. Comparative homology modeling of MDH proteins revealed that evolutionarily conserved lysines serve as hotspots of acetylation. Our combined analyses indicate lysine acetylation as a common strategy to fine-tune the activity of central metabolic enzymes with likely impact on plant acclimation capacity.

Functional characterization of proton antiport regulation in the thylakoid membrane

During photosynthesis, energy is transiently stored as an electrochemical proton gradient across the thylakoid membrane. The resulting proton motive force (pmf) is composed of a membrane potential (ΔΨ) and a proton concentration gradient (ΔpH) and powers the synthesis of ATP. Light energy availability for photosynthesis can change very rapidly and frequently in nature. Thylakoid ion transport proteins buffer the effects that light fluctuations have on photosynthesis by adjusting pmf and its composition. Ion channel activities dissipate ΔΨ, thereby reducing charge recombinations within photosystem II. The dissipation of ΔΨ allows for increased accumulation of protons in the thylakoid lumen, generating the signal that activates feedback downregulation of photosynthesis. Proton export from the lumen via the thylakoid K+ exchange antiporter 3 (KEA3), instead, decreases the ΔpH fraction of the pmf and thereby reduces the regulatory feedback signal. Here, we reveal that the Arabidopsis (Arabidopsis thaliana) KEA3 protein homo-dimerizes via its C-terminal domain. This C-terminus has a regulatory function, which responds to light intensity transients. Plants carrying a C-terminus-less KEA3 variant show reduced feed-back downregulation of photosynthesis and suffer from increased photosystem damage under long-term high light stress. However, during photosynthetic induction in high light, KEA3 deregulation leads to an increase in carbon fixation rates. Together, the data reveal a trade-off between long-term photoprotection and a short-term boost in carbon fixation rates, which is under the control of the KEA3 C-terminus.

Isotopically Dimethyl Labeling-Based Quantitative Proteomic Analysis of Phosphoproteomes of Soybean Cultivars

Isotopically dimethyl labeling was applied in a quantitative post-translational modification (PTM) proteomic study of phosphoproteomic changes in the drought responses of two contrasting soybean cultivars. A total of 9457 phosphopeptides were identified subsequently, corresponding to 4571 phosphoprotein groups and 3889 leading phosphoproteins, which contained nine kinase families consisting of 279 kinases. These phosphoproteins contained a total of 8087 phosphosites, 6106 of which were newly identified and constituted 54% of the current soybean phosphosite repository. These phosphosites were converted into the highly conserved kinase docking sites by bioinformatics analysis, which predicted six kinase families that matched with those newly found nine kinase families. The overly post-translationally modified proteins (OPP) occupies 2.1% of these leading phosphoproteins. Most of these OPPs are photoreceptors, mRNA-, histone-, and phospholipid-binding proteins, as well as protein kinase/phosphatases. The subgroup population distribution of phosphoproteins over the number of phosphosites of phosphoproteins follows the exponential decay law, Y = 4.13e-0.098X - 0.04. Out of 218 significantly regulated unique phosphopeptide groups, 188 phosphoproteins were regulated by the drought-tolerant cultivar under the water loss condition. These significantly regulated phosphoproteins (SRP) are mainly enriched in the biological functions of water transport and deprivation, methionine metabolic processes, photosynthesis/light reaction, and response to cadmium ion, osmotic stress, and ABA response. Seventeen and 15 SRPs are protein kinases/phosphatases and transcription factors, respectively. Bioinformatics analysis again revealed that three members of the calcium dependent protein kinase family (CAMK family), GmSRK2I, GmCIPK25, and GmAKINβ1 kinases, constitute a phosphor-relay-mediated signal transduction network, regulating ion channel activities and many nuclear events in this drought-tolerant cultivar, which presumably contributes to the development of the soybean drought tolerance under water deprivation process.

Two mitochondrial phosphatases, PP2c63 and Sal2, are required for posttranslational regulation of the TCA cycle in Arabidopsis

Protein phosphorylation is a well-established post-translational mechanism that regulates protein functions and metabolic pathways. It is known that several plant mitochondrial proteins are phosphorylated in a reversible manner. However, the identities of the protein kinases/phosphatases involved in this mechanism and their roles in the regulation of the tricarboxylic acid (TCA) cycle remain unclear. In this study, we isolated and characterized plants lacking two mitochondrially targeted phosphatases (Sal2 and PP2c63) along with pyruvate dehydrogenase kinase (PDK). Protein-protein interaction analysis, quantitative phosphoproteomics, and enzymatic analyses revealed that PDK specifically regulates pyruvate dehydrogenase complex (PDC), while PP2c63 nonspecifically regulates PDC. When recombinant PP2c63 and Sal2 proteins were added to mitochondria isolated from mutant plants, protein-protein interaction and enzymatic analyses showed that PP2c63 directly phosphorylates and modulates the activity of PDC, while Sal2 only indirectly affects TCA cycle enzymes. Characterization of steady-state metabolite levels and fluxes in the mutant lines further revealed that these phosphatases regulate flux through the TCA cycle, and that altered metabolism in the sal2 pp2c63 double mutant compromises plant growth. These results are discussed in the context of current models of the control of respiration in plants.

Dual lysine and N-terminal acetyltransferases reveal the complexity underpinning protein acetylation

Protein acetylation is a highly frequent protein modification. However, comparatively little is known about its enzymatic machinery. N-α-acetylation (NTA) and ε-lysine acetylation (KA) are known to be catalyzed by distinct families of enzymes (NATs and KATs, respectively), although the possibility that the same GCN5-related N-acetyltransferase (GNAT) can perform both functions has been debated. Here, we discovered a new family of plastid-localized GNATs, which possess a dual specificity. All characterized GNAT family members display a number of unique features. Quantitative mass spectrometry analyses revealed that these enzymes exhibit both distinct KA and relaxed NTA specificities. Furthermore, inactivation of GNAT2 leads to significant NTA or KA decreases of several plastid proteins, while proteins of other compartments were unaffected. The data indicate that these enzymes have specific protein targets and likely display partly redundant selectivity, increasing the robustness of the acetylation process in vivo. In summary, this study revealed a new layer of complexity in the machinery controlling this prevalent modification and suggests that other eukaryotic GNATs may also possess these previously underappreciated broader enzymatic activities.

High-Resolution Lysine Acetylome Profiling by Offline Fractionation and Immunoprecipitation

Acetylation of lysine side chains at their ε-amino group is a reversible posttranslational modification (PTM), which can affect diverse protein functions. Lysine acetylation was first described on histones, and nowadays gains more and more attention due to its more general occurrence in proteomes, and its possible crosstalk with other protein modifications. Here we describe a workflow to investigate the acetylation of lysine-containing peptides on a large scale. For this high-resolution lysine acetylome analysis, dimethyl-labeled peptide samples are pooled and offline-fractionated using hydrophilic interaction liquid chromatography (HILIC). The offline fractionation is followed by an immunoprecipitation and liquid chromatography-tandem mass spectrometry (LC-MS/MS) for data acquisition and subsequent data analysis.

New beginnings and new ends: methods for large-scale characterization of protein termini and their use in plant biology

Dynamic regulation of protein function and abundance plays an important role in virtually every aspect of plant life. Diversifying mechanisms at the RNA and protein level result in many protein molecules with distinct sequence and modification, termed proteoforms, arising from a single gene. Distinct protein termini define proteoforms arising from translation of alternative transcripts, use of alternative translation initiation sites, and different co- and post-translational modifications of the protein termini. Also site-specific proteolytic processing by endo- and exoproteases generates truncated proteoforms, defined by distinct protease-generated neo-N- and neo-C-termini, that may exhibit altered activity, function, and localization compared with their precursor proteins. In eukaryotes, the N-degron pathway targets cytosolic proteins, exposing destabilizing N-terminal amino acids and/or destabilizing N-terminal modifications for proteasomal degradation. This enables rapid and selective removal not only of unfolded proteins, but also of substrate proteoforms generated by proteolytic processing or changes in N-terminal modifications. Here we summarize current protocols enabling proteome-wide analysis of protein termini, which have provided important new insights into N-terminal modifications and protein stability determinants, protein maturation pathways, and protease-substrate relationships in plants.

Establishment of Dimethyl Labeling-based Quantitative Acetylproteomics in Arabidopsis

Protein acetylation, one of many types of post-translational modifications (PTMs), is involved in a variety of biological and cellular processes. In the present study, we applied both CsCl density gradient (CDG) centrifugation-based protein fractionation and a dimethyl-labeling-based 4C quantitative PTM proteomics workflow in the study of dynamic acetylproteomic changes in Arabidopsis. This workflow integrates the dimethyl chemical labeling with chromatography-based acetylpeptide separation and enrichment followed by mass spectrometry (MS) analysis, the extracted ion chromatogram (XIC) quantitation-based computational analysis of mass spectrometry data to measure dynamic changes of acetylpeptide level using an in-house software program, named Stable isotope-based Quantitation-Dimethyl labeling (SQUA-D), and finally the confirmation of ethylene hormone-regulated acetylation using immunoblot analysis. Eventually, using this proteomic approach, 7456 unambiguous acetylation sites were found from 2638 different acetylproteins, and 5250 acetylation sites, including 5233 sites on lysine side chain and 17 sites on protein N termini, were identified repetitively. Out of these repetitively discovered acetylation sites, 4228 sites on lysine side chain (i.e. 80.5%) are novel. These acetylproteins are exemplified by the histone superfamily, ribosomal and heat shock proteins, and proteins related to stress/stimulus responses and energy metabolism. The novel acetylproteins enriched by the CDG centrifugation fractionation contain many cellular trafficking proteins, membrane-bound receptors, and receptor-like kinases, which are mostly involved in brassinosteroid, light, gravity, and development signaling. In addition, we identified 12 highly conserved acetylation site motifs within histones, P-glycoproteins, actin depolymerizing factors, ATPases, transcription factors, and receptor-like kinases. Using SQUA-D software, we have quantified 33 ethylene hormone-enhanced and 31 hormone-suppressed acetylpeptide groups or called unique PTM peptide arrays (UPAs) that share the identical unique PTM site pattern (UPSP). This CDG centrifugation protein fractionation in combination with dimethyl labeling-based quantitative PTM proteomics, and SQUA-D may be applied in the quantitation of any PTM proteins in any model eukaryotes and agricultural crops as well as tissue samples of animals and human beings.

Lysine acetylome profiling uncovers novel histone deacetylase substrate proteins in Arabidopsis

Histone deacetylases have central functions in regulating stress defenses and development in plants. However, the knowledge about the deacetylase functions is largely limited to histones, although these enzymes were found in diverse subcellular compartments. In this study, we determined the proteome‐wide signatures of the RPD3/HDA1 class of histone deacetylases in Arabidopsis. Relative quantification of the changes in the lysine acetylation levels was determined on a proteome‐wide scale after treatment of Arabidopsis leaves with deacetylase inhibitors apicidin and trichostatin A. We identified 91 new acetylated candidate proteins other than histones, which are potential substrates of the RPD3/HDA1‐like histone deacetylases in Arabidopsis, of which at least 30 of these proteins function in nucleic acid binding. Furthermore, our analysis revealed that histone deacetylase 14 (HDA14) is the first organellar‐localized RPD3/HDA1 class protein found to reside in the chloroplasts and that the majority of its protein targets have functions in photosynthesis. Finally, the analysis of HDA14 loss‐of‐function mutants revealed that the activation state of RuBisCO is controlled by lysine acetylation of RuBisCO activase under low‐light conditions.