Cysteine oxidation as a regulatory mechanism of Arabidopsis plastidial NAD-dependent malate dehydrogenase.

Malate dehydrogenases (MDHs) catalyze a reversible NAD(P)-dependent-oxidoreductase reaction that plays an important role in central metabolism and redox homeostasis of plant cells. Recent studies suggest a moonlighting function of plastidial NAD-dependent MDH (plNAD-MDH; EC 1.1.1.37) in plastid biogenesis, independent of its enzyme activity. In this study, redox effects on activity and conformation of recombinant plNAD-MDH from Arabidopsis thaliana were investigated. We show that reduced plNAD-MDH is active while it is inhibited upon oxidation. Interestingly, the presence of its cofactors NAD+ and NADH could prevent oxidative inhibition of plNAD-MDH. In addition, a conformational change upon oxidation could be observed via non-reducing SDS-PAGE. Both effects, its inhibition and conformational change, were reversible by re-reduction. Further investigation of single cysteine substitutions and mass spectrometry revealed that oxidation of plNAD-MDH leads to oxidation of all four cysteine residues. However, cysteine oxidation of C129 leads to inhibition of plNAD-MDH activity and oxidation of C147 induces its conformational change. In contrast, oxidation of C190 and C333 does not affect plNAD-MDH activity or structure. Our results demonstrate that plNAD-MDH activity can be reversibly inhibited, but not inactivated, by cysteine oxidation and might be co-regulated by the availability of its cofactors in vivo.

Specificity and dynamics of H2O2 detoxification by the cytosolic redox regulatory network as revealed by in vitro reconstitution.

The thiol redox state is a decisive functional characteristic of proteins in cell biology. Plasmatic cell compartments maintain a thiol-based redox regulatory network linked to the glutathione/glutathione disulfide couple (GSH/GSSG) and the NAD(P)H system. The basic network constituents are known and in vivo cell imaging with gene-encoded probes have revealed insight into the dynamics of the [GSH]2/[GSSG] redox potential, cellular H2O2 and NAD(P)H+H+ amounts in dependence on metabolic and environmental cues. Less understood is the contribution and interaction of the network components, also because of compensatory reactions in genetic approaches. Reconstituting the cytosolic network of Arabidopsis thaliana in vitro from fifteen recombinant proteins at in vivo concentrations, namely glutathione peroxidase-like (GPXL), peroxiredoxins (PRX), glutaredoxins (GRX), thioredoxins, NADPH-dependent thioredoxin reductase A and glutathione reductase and applying Grx1-roGFP2 or roGFP2-Orp1 as dynamic sensors, allowed for monitoring the response to a single H2O2 pulse. The major change in thiol oxidation as quantified by mass spectrometry-based proteomics occurred in relevant peptides of GPXL, and to a lesser extent of PRX, while other Cys-containing peptides only showed small changes in their redox state and protection. Titration of ascorbate peroxidase (APX) into the system together with dehydroascorbate reductase lowered the oxidation of the fluorescent sensors in the network but was unable to suppress it. The results demonstrate the power of the network to detoxify H2O2, the partially independent branches of electron flow with significance for specific cell signaling and the importance of APX to modulate the signaling without suppressing it and shifting the burden to glutathione oxidation.

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

In humans and plants, 40% of the proteome is co-translationally 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.

Chemoselective umpolung of thiols to episulfoniums for cysteine bioconjugation.

Cysteine conjugation is an important tool in protein research and relies on fast, mild and chemoselective reactions. Cysteinyl thiols can either be modified with prefunctionalized electrophiles, or converted into electrophiles themselves for functionalization with selected nucleophiles in an independent step. Here we report a bioconjugation strategy that uses a vinyl thianthrenium salt to transform cysteine into a highly reactive electrophilic episulfonium intermediate in situ, to enable conjugation with a diverse set of bioorthogonal nucleophiles in a single step. The reactivity profile can connect several nucleophiles to biomolecules through a short and stable ethylene linker, ideal for introduction of infrared labels, post-translational modifications or NMR probes. In the absence of reactive exogenous nucleophiles, nucleophilic amino acids can react with the episulfonium intermediate for native peptide stapling and protein-protein ligation. Ready synthetic access to isotopologues of vinyl thianthrenium salts enables applications in quantitative proteomics. Such diverse applications demonstrate the utility of vinyl-thianthrenium-based bioconjugation as a fast, selective and broadly applicable tool for chemical biology.

Differential proteome profiling of bacterial culture supernatants reveals candidates for the induction of oral immune priming in the red flour beetle.

Most organisms are host to symbionts and pathogens, which led to the evolution of immune strategies to prevent harm. Whilst the immune defences of vertebrates are classically divided into innate and adaptive, insects lack specialized cells involved in adaptive immunity, but have been shown to exhibit immune priming: the enhanced survival upon infection after a first exposure to the same pathogen or pathogen-derived components. An important piece of the puzzle are the pathogen-associated molecules that induce these immune priming responses. Here, we make use of the model system consisting of the red flour beetle (Tribolium castaneum) and its bacterial pathogen Bacillus thuringiensis, to compare the proteomes of culture supernatants of two closely related B. thuringiensis strains that either induce priming via the oral route, or not. Among the proteins that might be immunostimulatory to T. castaneum, we identify the Cry3Aa toxin, an important plasmid-encoded virulence factor of B. thuringiensis. In further priming-infection assays we test the relevance of Cry-carrying plasmids for immune priming. Our findings provide valuable insights for future studies to perform experiments on the mechanisms and evolution of immune priming.

Eukaryote-specific assembly factor DEAP2 mediates an early step of photosystem II assembly in Arabidopsis.

The initial step of oxygenic photosynthesis is the thermodynamically challenging extraction of electrons from water and the release of molecular oxygen. This light-driven process, which is the basis for most life on Earth, is catalyzed by photosystem II (PSII) within the thylakoid membrane of photosynthetic organisms. The biogenesis of PSII requires a controlled step-wise assembly process of which the early steps are considered to be highly conserved between plants and their cyanobacterial progenitors. This assembly process involves auxiliary proteins, which are likewise conserved. In the present work, we used Arabidopsis (Arabidopsis thaliana) as a model to show that in plants, a eukaryote-exclusive assembly factor facilitates the early assembly step, during which the intrinsic antenna protein CP47 becomes associated with the PSII reaction center (RC) to form the RC47 intermediate. This factor, which we named DECREASED ELECTRON TRANSPORT AT PSII (DEAP2), works in concert with the conserved PHOTOSYNTHESIS AFFECTED MUTANT 68 (PAM68) assembly factor. The deap2 and pam68 mutants showed similar defects in PSII accumulation and assembly of the RC47 intermediate. The combined lack of both proteins resulted in a loss of functional PSII and the inability of plants to grow photoautotrophically on the soil. While overexpression of DEAP2 partially rescued the pam68 PSII accumulation phenotype, this effect was not reciprocal. DEAP2 accumulated at 20-fold higher levels than PAM68, together suggesting that both proteins have distinct functions. In summary, our results uncover eukaryotic adjustments to the PSII assembly process, which involve the addition of DEAP2 for the rapid progression from RC to RC47.

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.

Peptide CoA conjugates for in situ proteomics profiling of acetyltransferase activities.

The acetylation of protein N-termini is a co- or posttranslational modification that plays important roles in protein homeostasis and stability. N-terminal acetyltransferases (NATs) catalyze the introduction of this modification using acetyl-coenzyme A (acetyl-CoA) as source of the acetyl-group. NATs operate in complex with auxiliary proteins that impact activity and specificity of these enzymes. Proper function of NATs is essential for development in plants and mammals alike. High resolution mass spectrometry (MS) is a powerful tool for investigating NATs and protein complexes in general. However, efficient methods for enriching NAT complexes ex vivo from cellular extracts are needed for the subsequent analysis. Based on bisubstrate analog inhibitors of lysine acetyltransferases, peptide-CoA conjugates have been developed as capture compounds of NATs. The N-terminal residue of these probes, serving as attachment site of the CoA moiety, was shown to impact NAT binding according to the respective amino acid specificity of these enzymes. This chapter reports the detailed protocols for the synthesis of peptide-CoA conjugates, the experimental procedures for NAT enrichment as well as the MS and data analysis. Collectively, these protocols provide a set of tools for profiling NAT complexes in cell lysates of healthy or diseases backgrounds.

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.

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 O2 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.

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.

Alternative splicing of Arabidopsis G6PD5 recruits NADPH-producing OPPP reactions to the endoplasmic reticulum.

Glucose-6-phosphate dehydrogenase is the rate-limiting enzyme of the oxidative pentose-phosphate pathway (OPPP). The OPPP mainly provides NADPH and sugar-phosphate building blocks for anabolic pathways and is present in all eukaryotes. In plant cells, the irreversible part of the OPPP is found in several compartments. Among the isoforms catalyzing the first OPPP step in Arabidopsis, G6PD1 to G6PD4 target plastids (with G6PD1 being also directed to peroxisomes), whereas G6PD5 and G6PD6 operate in the cytosol. We noticed that alternative splice forms G6PD5.4 and G6PD5.5 encode N-terminally extended proteoforms. Compared to G6PD5.1, RT-PCR signals differed and fluorescent reporter fusions expressed in Arabidopsis protoplasts accumulated in distinct intracellular sites. Co-expression with organelle-specific markers revealed that the G6PD5.4 and G6PD5.5 proteoforms label different subdomains of the endoplasmic reticulum (ER), and analysis of C-terminal roGFP fusions showed that their catalytic domains face the cytosol. In g6pd5-1 g6pd6-2 mutant protoplasts lacking cytosolic G6PDH activity, the ER-bound proteoforms were both active and thus able to form homomers. Among the Arabidopsis 6-phosphogluconolactonases (catalyzing the second OPPP step), we noticed that isoform PGL2 carries a C-terminal CaaX motif that may be prenylated for membrane attachment. Reporter-PGL2 fusions co-localized with G6PD5.4 in ER subdomains, which was abolished by Cys-to-Ser exchange in the 256CSIL motif. Among the Arabidopsis 6-phosphogluconate dehydrogenases (catalyzing the third OPPP step), S-acylated peptides were detected for all three isoforms in a recent palmitoylome, with dual cytosolic/peroxisomal PGD2 displaying three sites. Co-expression of GFP-PGD2 diminished crowding of OFP-G6PD5.4 at the ER, independent of PGL2's presence. Upon pull-down of GFP-G6PD5.4, not only unlabeled PGD2 and PGL2 were enriched, but also enzymes that depend on NADPH provision at the ER, indicative of physical interaction with the OPPP enzymes. When membrane-bound G6PD5.5 and 5.4 variants were co-expressed with KCR1 (ketoacyl-CoA reductase, involved in fatty acid elongation), ATR1 (NADPH:cytochrome-P450 oxidoreductase), or pulled C4H/CYP73A5 (cinnamate 4-hydroxylase) as indirectly (via ATR) NADPH-dependent cytochrome P450 enzyme, co-localization in ER subdomains was observed. Thus, alternative splicing of G6PD5 can direct the NADPH-producing OPPP reactions to the cytosolic face of the ER, where they may operate as membrane-bound metabolon to support several important biosynthetic pathways of plant cells.

Glutamate 1-semialdehyde aminotransferase is connected to GluTR by GluTR-binding protein and contributes to the rate-limiting step of 5-aminolevulinic acid synthesis.

Tetrapyrroles play fundamental roles in crucial processes including photosynthesis, respiration, and catalysis. In plants, 5-aminolevulinic acid (ALA) is the common precursor of tetrapyrroles. ALA is synthesized from activated glutamate by the enzymes glutamyl-tRNA reductase (GluTR) and glutamate-1-semialdehyde aminotransferase (GSAAT). ALA synthesis is recognized as the rate-limiting step in this pathway. We aimed to explore the contribution of GSAAT to the control of ALA synthesis and the formation of a protein complex with GluTR. In Arabidopsis thaliana, two genes encode GSAAT isoforms: GSA1 and GSA2. A comparison of two GSA knockout mutants with the wild-type revealed the correlation of reduced GSAAT activity and ALA-synthesizing capacity in leaves with lower chlorophyll content. Growth and green pigmentation were more severely impaired in gsa2 than in gsa1, indicating the predominant role of GSAAT2 in ALA synthesis. Interestingly, GluTR accumulated to higher levels in gsa2 than in the wild-type and was mainly associated with the plastid membrane. We propose that the GSAAT content modulates the amount of soluble GluTR available for ALA synthesis. Several different biochemical approaches revealed the GSAAT-GluTR interaction through the assistance of GluTR-binding protein (GBP). A modeled structure of the tripartite protein complex indicated that GBP mediates the stable association of GluTR and GSAAT for adequate ALA synthesis.

Investigating Peptide-Coenzyme A Conjugates as Chemical Probes for Proteomic Profiling of N-Terminal and Lysine Acetyltransferases.

Acetyl groups are transferred from acetyl-coenzyme A (Ac-CoA) to protein N-termini and lysine side chains by N-terminal acetyltransferases (NATs) and lysine acetyltransferases (KATs), respectively. Building on lysine-CoA conjugates as KAT probes, we have synthesized peptide probes with CoA conjugated to N-terminal alanine (α-Ala-CoA), proline (α-Pro-CoA) or tri-glutamic acid (α-3Glu-CoA) units for interactome profiling of NAT complexes. The α-Ala-CoA probe enriched the majority of NAT catalytic and auxiliary subunits, while a lysine CoA-conjugate bound only a subset of endogenous KATs. Interactome profiling with the α-Pro-CoA probe showed reduced NAT recruitment in favor of metabolic CoA binding proteins and α-3Glu-CoA steered the interactome towards NAA80 and NatB. These findings agreed with the inherent substrate specificities of the target proteins and showed that N-terminal CoA-conjugated peptides are versatile probes for NAT complex profiling in lysates of physiological and pathological backgrounds.

Lysine acetylation regulates moonlighting activity of the E2 subunit of the chloroplast pyruvate dehydrogenase complex in Chlamydomonas.

The dihydrolipoamide acetyltransferase subunit DLA2 of the chloroplast pyruvate dehydrogenase complex (cpPDC) in the green alga Chlamydomonas reinhardtii has previously been shown to possess moonlighting activity in chloroplast gene expression. Under mixotrophic growth conditions, DLA2 forms part of a ribonucleoprotein particle (RNP) with the psbA mRNA that encodes the D1 protein of the photosystem II (PSII) reaction center. Here, we report on the characterization of the molecular switch that regulates shuttling of DLA2 between its functions in carbon metabolism and D1 synthesis. Determination of RNA-binding affinities by microscale thermophoresis demonstrated that the E3-binding domain (E3BD) of DLA2 mediates psbA-specific RNA recognition. Analyses of cpPDC formation and activity, as well as RNP complex formation, showed that acetylation of a single lysine residue (K197) in E3BD induces the release of DLA2 from the cpPDC, and its functional shift towards RNA binding. Moreover, Förster resonance energy transfer microscopy revealed that psbA mRNA/DLA2 complexes localize around the chloroplast's pyrenoid. Pulse labeling and D1 re-accumulation after induced PSII degradation strongly suggest that DLA2 is important for D1 synthesis during de novo PSII biogenesis.

Mass Spectrometry-Based Quantitative Cysteine Redox Proteome Profiling of Isolated Mitochondria Using Differential iodoTMT Labeling.

Mitochondria are central hubs of redox biochemistry in the cell. An important role of mitochondrial carbon metabolism is to oxidize respiratory substrates and to pass the electrons down the mitochondrial electron transport chain to reduce oxygen and to drive oxidative phosphorylation. During respiration, reactive oxygen species are produced as a side reaction, some of which in turn oxidize cysteine thiols in proteins. Hence, the redox status of cysteine-containing mitochondrial proteins has to be controlled by the mitochondrial glutathione and thioredoxin systems, which draw electrons from metabolically derived NADPH. The redox status of mitochondrial cysteines can undergo fast transitions depending on the metabolic status of the cell, as for instance at early seed germination. Here, we describe a state-of-the-art method to quantify redox state of protein cysteines in isolated Arabidopsis seedling mitochondria of controlled metabolic and respiratory state by MS2-based redox proteomics using the isobaric thiol labeling reagent Iodoacetyl Tandem Mass Tag™ (iodoTMT). The procedure is also applicable to isolated mitochondria of other plant and nonplant systems.

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.