Multisubstrate specificity shaped the complex evolution of the aminotransferase family across the tree of life.

Aminotransferases (ATs) are an ancient enzyme family that play central roles in core nitrogen metabolism, essential to all organisms. However, many of the AT enzyme functions remain poorly defined, limiting our fundamental understanding of the nitrogen metabolic networks that exist in different organisms. Here, we traced the deep evolutionary history of the AT family by analyzing AT enzymes from 90 species spanning the tree of life (ToL). We found that each organism has maintained a relatively small and constant number of ATs. Mapping the distribution of ATs across the ToL uncovered that many essential AT reactions are carried out by taxon-specific AT enzymes due to wide-spread nonorthologous gene displacements. This complex evolutionary history explains the difficulty of homology-based AT functional prediction. Biochemical characterization of diverse aromatic ATs further revealed their broad substrate specificity, unlike other core metabolic enzymes that evolved to catalyze specific reactions today. Interestingly, however, we found that these AT enzymes that diverged over billion years share common signatures of multisubstrate specificity by employing different nonconserved active site residues. These findings illustrate that AT family enzymes had leveraged their inherent substrate promiscuity to maintain a small yet distinct set of multifunctional AT enzymes in different taxa. This evolutionary history of versatile ATs likely contributed to the establishment of robust and diverse nitrogen metabolic networks that exist throughout the ToL. The study provides a critical foundation to systematically determine diverse AT functions and underlying nitrogen metabolic networks across the ToL.

Biochemical characterization of plant aromatic aminotransferases.

Aromatic aminotransferases (Aro ATs) are pyridoxal-5-phosphate (PLP)-dependent enzymes that catalyze the transamination reactions of an aromatic amino acid (AAA) or a keto acid. Aro ATs are involved in biosynthesis or degradation of AAAs and play important functions in controlling the production of plant hormones and secondary metabolites, such as auxin, tocopherols, flavonoids, and lignin. Most Aro ATs show substrate promiscuity and can accept multiple aromatic and non-aromatic amino and keto acid substrates, which complicates and limits our understanding of their in planta functions. Considering the critical roles Aro ATs play in plant primary and secondary metabolism, it is important to accurately determine substrate specificity and kinetic properties of Aro ATs. This chapter describes various methodologies of protein expression, purification and enzymatic assays, which can be used for biochemical characterization of Aro ATs.

Mass spectrometry imaging-based assays for aminotransferase activity reveal a broad substrate spectrum for a previously uncharacterized enzyme.

Aminotransferases (ATs) catalyze pyridoxal 5'-phosphate-dependent transamination reactions between amino donor and keto acceptor substrates and play central roles in nitrogen metabolism of all organisms. ATs are involved in the biosynthesis and degradation of both proteinogenic and nonproteinogenic amino acids and also carry out a wide variety of functions in photorespiration, detoxification, and secondary metabolism. Despite the importance of ATs, their functionality is poorly understood as only a small fraction of putative ATs, predicted from DNA sequences, are associated with experimental data. Even for characterized ATs, the full spectrum of substrate specificity, among many potential substrates, has not been explored in most cases. This is largely due to the lack of suitable high-throughput assays that can screen for AT activity and specificity at scale. Here we present a new high-throughput platform for screening AT activity using bioconjugate chemistry and mass spectrometry imaging-based analysis. Detection of AT reaction products is achieved by forming an oxime linkage between the ketone groups of transaminated amino donors and a probe molecule that facilitates mass spectrometry-based analysis using nanostructure-initiator mass spectrometry or MALDI-mass spectrometry. As a proof-of-principle, we applied the newly established method and found that a previously uncharacterized Arabidopsis thaliana tryptophan AT-related protein 1 is a highly promiscuous enzyme that can utilize 13 amino acid donors and three keto acid acceptors. These results demonstrate that this oxime-mass spectrometry imaging AT assay enables high-throughput discovery and comprehensive characterization of AT enzymes, leading to an accurate understanding of the nitrogen metabolic network.

An automated workflow that generates atom mappings for large-scale metabolic models and its application to Arabidopsis thaliana.

Quantification of reaction fluxes of metabolic networks can help us understand how the integration of different metabolic pathways determines cellular functions. Yet, intracellular fluxes cannot be measured directly but are estimated with metabolic flux analysis (MFA), which relies on the patterns of isotope labeling of metabolites in the network. The application of MFA also requires a stoichiometric model with atom mappings that are currently not available for the majority of large-scale metabolic network models, particularly of plants. While automated approaches such as the Reaction Decoder Toolkit (RDT) can produce atom mappings for individual reactions, tracing the flow of individual atoms of the entire reactions across a metabolic model remains challenging. Here we establish an automated workflow to obtain reliable atom mappings for large-scale metabolic models by refining the outcome of RDT, and apply the workflow to metabolic models of Arabidopsis thaliana. We demonstrate the accuracy of RDT through a comparative analysis with atom mappings from a large database of biochemical reactions, MetaCyc. We further show the utility of our automated workflow by simulating 15 N isotope enrichment and identifying nitrogen (N)-containing metabolites which show enrichment patterns that are informative for flux estimation in future 15 N-MFA studies of A. thaliana. The automated workflow established in this study can be readily expanded to other species for which metabolic models have been established and the resulting atom mappings will facilitate MFA and graph-theoretic structural analyses with large-scale metabolic networks.

Evolutionary origin and functional diversification of aminotransferases.

Aminotransferases (ATs) are pyridoxal 5'-phosphate-dependent enzymes that catalyze the transamination reactions between amino acid donor and keto acid acceptor substrates. Modern AT enzymes constitute ∼2% of all classified enzymatic activities, play central roles in nitrogen metabolism, and generate multitude of primary and secondary metabolites. ATs likely diverged into four distinct AT classes before the appearance of the last universal common ancestor and further expanded to a large and diverse enzyme family. Although the AT family underwent an extensive functional specialization, many AT enzymes retained considerable substrate promiscuity and multifunctionality because of their inherent mechanistic, structural, and functional constraints. This review summarizes the evolutionary history, diverse metabolic roles, reaction mechanisms, and structure-function relationships of the AT family enzymes, with a special emphasis on their substrate promiscuity and multifunctionality. Comprehensive characterization of AT substrate specificity is still needed to reveal their true metabolic functions in interconnecting various branches of the nitrogen metabolic network in different organisms.

Genome-wide association of the metabolic shifts underpinning dark-induced senescence in Arabidopsis.

Dark-induced senescence provokes profound metabolic shifts to recycle nutrients and to guarantee plant survival. To date, research on these processes has largely focused on characterizing mutants deficient in individual pathways. Here, we adopted a time-resolved genome-wide association-based approach to characterize dark-induced senescence by evaluating the photochemical efficiency and content of primary and lipid metabolites at the beginning, or after 3 or 6 days in darkness. We discovered six patterns of metabolic shifts and identified 215 associations with 81 candidate genes being involved in this process. Among these associations, we validated the roles of four genes associated with glycine, galactinol, threonine, and ornithine levels. We also demonstrated the function of threonine and galactinol catabolism during dark-induced senescence. Intriguingly, we determined that the association between tyrosine contents and TYROSINE AMINOTRANSFERASE 1 influences enzyme activity of the encoded protein and transcriptional activity of the gene under normal and dark conditions, respectively. Moreover, the single-nucleotide polymorphisms affecting the expression of THREONINE ALDOLASE 1 and the amino acid transporter gene AVT1B, respectively, only underlie the variation in threonine and glycine levels in the dark. Taken together, these results allow us to present a very detailed model of the metabolic aspects of dark-induced senescence, as well as the process itself.

The Rice Plastidial Phosphorylase Participates Directly in Both Sink and Source Processes.

The plastidial starch phosphorylase (Pho1) functions in starch metabolism. A distinctive structural feature of the higher Pho1 is a 50-82-amino-acid long peptide (L50-L82), which is absent in phosphorylases from non-plant organisms. To study the function of the rice Pho1 L80 peptide, we complemented a pho1- rice mutant (BMF136) with the wild-type Pho1 gene or with a Pho1 gene lacking the L80 region (Pho1ΔL80). While expression of Pho1 in BMF136 restored normal wild-type phenotype, the introduction of Pho1ΔL80 enhanced the growth rate and plant productivity above wild-type levels. Mass spectrometry analysis of proteins captured by anti-Pho1 showed the surprising presence of PsaC, the terminal electron acceptor/donor subunit of photosystem I (PSI). This unexpected interaction was substantiated by reciprocal immobilized protein pull-down assays of seedling extracts and supported by the presence of Pho1 on isolated PSI complexes resolved by blue-native gels. Spectrophotometric studies showed that Pho1ΔL80 plants exhibited modified PSI and enhanced CO2 assimilation properties. Collectively, these findings indicate that the higher plant Pho1 has dual roles as a potential modulator of source and sink processes.

The plastid phosphorylase as a multiple-role player in plant metabolism.

The physiological roles of the plastidial phosphorylase in starch metabolism of higher plants have been debated for decades. While estimated physiological substrate levels favor a degradative role, genetic evidence indicates that the plastidial phosphorylase (Pho1) plays an essential role in starch initiation and maturation of the starch granule in developing rice grains. The plastidial enzyme contains a unique peptide domain, up to 82 residues in length depending on the plant species, not found in its cytosolic counterpart or glycogen phosphorylases. The role of this extra peptide domain is perplexing, as its complete removal does not significantly affect the in vitro catalytic or enzymatic regulatory properties of rice Pho1. This peptide domain may have a regulatory function as it contains potential phosphorylation sites and, in some plant Pho1s, a PEST motif, a substrate for proteasome-mediated degradation. We discuss the potential roles of Pho1 and its L80 domain in starch biosynthesis and photosynthesis.

Re-programming of gene expression in the CS8 rice line over-expressing ADPglucose pyrophosphorylase induces a suppressor of starch biosynthesis.

The CS8 transgenic rice (Oryza sativa L.) lines expressing an up-regulated glgC gene produced higher levels of ADPglucose (ADPglc), the substrate for starch synthases. However, the increase in grain weight was much less than the increase in ADPglc levels suggesting one or more downstream rate-limiting steps. Endosperm starch levels were not further enhanced in double transgenic plants expressing both glgC and the maize brittle-1 gene, the latter responsible for transport of ADPglc into the amyloplast. These studies demonstrate that critical processes within the amyloplast stroma restrict maximum carbon flow into starch. RNA-seq analysis showed extensive re-programming of gene expression in the CS8 with 2073 genes up-regulated and 140 down-regulated. One conspicuous gene, up-regulated ~15-fold, coded for a biochemically uncharacterized starch binding domain-containing protein (SBDCP1) possessing a plastid transit peptide. Confocal microscopy and transmission electron microscopy analysis confirmed that SBDCP1 was located in the amyloplasts. Reciprocal immunoprecipitation and pull-down assays indicated an interaction between SBDCP1 and starch synthase IIIa (SSIIIa), which was down-regulated at the protein level in the CS8 line. Furthermore, binding by SBDCP1 inhibited SSIIIa starch polymerization activity in a non-competitive manner. Surprisingly, artificial microRNA gene suppression of SBDCP1 restored protein expression levels of SSIIIa in the CS8 line resulting in starch with lower amylose content and increased amylopectin chains with a higher degree of polymerization. Collectively, our results support the involvement of additional non-enzymatic factors such as SBDCP in starch biosynthesis.

Glu-370 in the large subunit influences the substrate binding, allosteric, and heat stability properties of potato ADP-glucose pyrophosphorylase.

ADP-glucose pyrophosphorylase (AGPase) is a key allosteric enzyme in plant starch biosynthesis. Plant AGPase is a heterotetrameric enzyme that consists of large (LS) and small subunits (SS), which are encoded by two different genes. In this study, we showed that the conversion of Glu to Gly at position 370 in the LS of AGPase alters the heterotetrameric stability along with the binding properties of substrate and effectors of the enzyme. Kinetic analyses revealed that the affinity of the LSE370GSSWT AGPase for glucose-1-phosphate is 3-fold less than for wild type (WT) AGPase. Additionally, the LSE370GSSWT AGPase requires 3-fold more 3-phosphogyceric acid to be activated. Finally, the LSE370GSSWTAGPase is less heat stable compared with the WT AGPase. Computational analysis of the mutant Gly-370 in the 3D modeled LS AGPase showed that this residue changes charge distribution of the surface and thus affect stability of the LS AGPase and overall heat stability of the heterotetrameric AGPase. In summary, our results show that LSE370 intricately modulate the heat stability and enzymatic activity of potato the AGPase.

Rice Endosperm Starch Phosphorylase (Pho1) Assembles with Disproportionating Enzyme (Dpe1) to Form a Protein Complex That Enhances Synthesis of Malto-oligosaccharides.

Starch synthesis in cereal grain endosperm is dependent on the concerted actions of many enzymes. The starch plastidial phosphorylase (Pho1) plays an important role in the initiation of starch synthesis and in the maturation of starch granule in developing rice seeds. Prior evidence has suggested that the rice enzyme, OsPho1, may have a physical/functional interaction with other starch biosynthetic enzymes. Pulldown experiments showed that OsPho1 as well as OsPho1 devoid of its L80 region, a peptide unique to higher plant phosphorylases, captures disproportionating enzyme (OsDpe1). Interaction of the latter enzyme form with OsDpe1 indicates that the putative regulatory L80 is not responsible for multienzyme assembly. This heterotypic enzyme complex, determined at a molar ratio of 1:1, was validated by reciprocal co-immunoprecipitation studies of native seed proteins and by co-elution chromatographic and co-migration electrophoretic patterns of these enzymes in rice seed extracts. The OsPho1-OsDpe1 complex utilized a broader range of substrates for enhanced synthesis of larger maltooligosaccharides than each individual enzyme and significantly elevated the substrate affinities of OsPho1 at 30 °C. Moreover, the assembly with OsDpe1 enables OsPho1 to utilize products of transglycosylation reactions involving G1 and G3, sugars that it cannot catalyze directly.

Substrate binding properties of potato tuber ADP-glucose pyrophosphorylase as determined by isothermal titration calorimetry.

Substrate binding properties of the large (LS) and small (SS) subunits of potato tuber ADP-glucose pyrophosphorylase were investigated by using isothermal titration calorimetry. Our results clearly show that the wild type heterotetramer (S(WT)L(WT)) possesses two distinct types of ATP binding sites, whereas the homotetrameric LS and SS variant forms only exhibited properties of one of the two binding sites. The wild type enzyme also exhibited significantly increased affinity to this substrate compared to the homotetrameric enzyme forms. No stable binding was evident for the second substrate, glucose-1-phosphate, in the presence or absence of ATPγS suggesting that interaction of glucose-1-phosphate is dependent on hydrolysis of ATP and supports the Theorell-Chance bi bi reaction mechanism.

Enhanced heterotetrameric assembly of potato ADP-glucose pyrophosphorylase using reverse genetics.

ADP-glucose pyrophosphorylase (AGPase) is a key allosteric enzyme in plant starch biosynthesis. Plant AGPase is a heterotetrameric enzyme that consists of large (LS) and small subunits (SS), which are encoded by two different genes. Computational and experimental studies have revealed that the heterotetrameric assembly of AGPase is thermodynamically weak. Modeling studies followed by the mutagenesis of the LS of the potato AGPase identified a heterotetramer-deficient mutant, LS(R88A). To enhance heterotetrameric assembly, LS(R88A) cDNA was subjected to error-prone PCR, and second-site revertants were identified according to their ability to restore glycogen accumulation, as assessed with iodine staining. Selected mutations were introduced into the wild-type (WT) LS and co-expressed with the WT SS in Escherichia coli glgC(-). The biochemical characterization of revertants revealed that LS(I90V)SS(WT), LS(Y378C)SS(WT) and LS(D410G)SS(WT) mutants displayed enhanced heterotetrameric assembly with the WT SS. Among these mutants, LS(Y378C)SS(WT) AGPase displayed increased heat stability compared with the WT enzyme. Kinetic characterization of the mutants indicated that the LS(I90V)SS(WT) and LS(Y378C)SS(WT) AGPases have comparable allosteric and kinetic properties. However, the LS(D410G)SS(WT) mutant exhibited altered allosteric properties of being less responsive and more sensitive to 3-phosphoglyceric acid activation and inorganic phosphate inhibition. This study not only enhances our understanding of the interaction between the SS and the LS of AGPase but also enables protein engineering to obtain enhanced assembled heat-stable variants of AGPase, which can be used for the improvement of plant yields.