Integrative analysis of transgenic alfalfa (Medicago sativa L.) suggests new metabolic control mechanisms for monolignol biosynthesis

The entanglement of lignin polymers with cellulose and hemicellulose in plant cell walls is a major biological barrier to the economically viable production of biofuels from woody biomass. Recent efforts of reducing this recalcitrance with transgenic techniques have been showing promise for ameliorating or even obviating the need for costly pretreatments that are otherwise required to remove lignin from cellulose and hemicelluloses. At the same time, genetic manipulations of lignin biosynthetic enzymes have sometimes yielded unforeseen consequences on lignin composition, thus raising the question of whether the current understanding of the pathway is indeed correct. To address this question systemically, we developed and applied a novel modeling approach that, instead of analyzing the pathway within a single target context, permits a comprehensive, simultaneous investigation of different datasets in wild type and transgenic plants. Specifically, the proposed approach combines static flux-based analysis with a Monte Carlo simulation in which very many randomly chosen sets of parameter values are evaluated against kinetic models of lignin biosynthesis in different stem internodes of wild type and lignin-modified alfalfa plants. In addition to four new postulates that address the reversibility of some key reactions, the modeling effort led to two novel postulates regarding the control of the lignin biosynthetic pathway. The first posits functionally independent pathways toward the synthesis of different lignin monomers, while the second postulate proposes a novel feedforward regulatory mechanism. Subsequent laboratory experiments have identified the signaling molecule salicylic acid as a potential mediator of the postulated control mechanism. Overall, the results demonstrate that mathematical modeling can be a valuable complement to conventional transgenic approaches and that it can provide biological insights that are otherwise difficult to obtain.     

Biosynthesis of monolignols. Genomic and reverse genetic approaches

The biosynthesis of monolignols is one of the most studied pathways of plant natural product biosynthesis. However, the pathway has recently undergone considerable revision, and it would appear that our understanding of the exact routes for synthesis of the building blocks of lignin and lignans is still not fully understood. Early studies of in vitro enzyme specificity failed to appreciate the catalytic promiscuity of some of the enzymes of the monolignol pathway, and the evolving model of a metabolic grid for monolignol biosynthesis may fail to appreciate the possible extent of metabolic channeling within the pathway. New approaches to the study of monolignol biosynthesis include genomics, advanced cellular imaging techniques, and transgenic manipulation. This article summarizes the use of these approaches to gain a better understanding of the operation of a complex metabolic pathway.     

Modeling allosteric regulation of de novo pyrimidine biosynthesis in Escherichia coli

With the emergence of multifaceted bioinformatics-derived data, it is becoming possible to merge biochemical and physiological information to develop a new level of understanding of the metabolic complexity of the cell. The biosynthetic pathway of de novo pyrimidine nucleotide metabolism is an essential capability of all free-living cells, and it occupies a pivotal position relative to metabolic processes that are involved in the macromolecular synthesis of DNA, RNA and proteins, as well as energy production and cell division. This regulatory network in all enteric bacteria involves genetic, allosteric, and physiological control systems that need to be integrated into a coordinated set of metabolic checks and balances. Allosterically regulated pathways constitute an exciting and challenging biosynthetic system to be approached from a mathematical perspective. However, to date, a mathematical model quantifying the contribution of allostery in controlling the dynamics of metabolic pathways has not been proposed. In this study, a direct, rigorous mathematical model of the de novo biosynthesis of pyrimidine nucleotides is presented. We corroborate the simulations with experimental data available in the literature and validate it with derepression experiments done in our laboratory. The model is able to faithfully represent the dynamic changes in the intracellular nucleotide pools that occur during metabolic transitions of the de novo pyrimidine biosynthetic pathway and represents a step forward in understanding the role of allosteric regulation in metabolic control.     

Downregulation of caffeoyl-CoA O-methyltransferase (CCoAOMT) by RNA interference leads to reduced lignin production in maize straw

Lignin is a major cell wall component of vascular plants that provides mechanical strength and hydrophobicity to vascular vessels. However, the presence of lignin limits the effective use of crop straw in many agroindustrial processes. Here, we generated transgenic maize plants in which the expression of a lignin biosynthetic gene encoding CCoAOMT, a key enzyme involved in the lignin biosynthesis pathway was downregulated by RNA interference (RNAi). RNAi of CCoAOMT led to significantly downregulated expression of this gene in transgenic maize compared with WT plants. These transgenic plants exhibited a 22.4% decrease in Klason lignin content and a 23.3% increase in cellulose content compared with WT plants, which may reflect compensatory regulation of lignin and cellulose deposition. We also measured the lignin monomer composition of the RNAi plants by GC-MS and determined that transgenic plants had a 57.08% higher S/G ratio than WT plants. In addition, histological staining of lignin with Wiesner reagent produced slightly more coloration in the xylem and sclerenchyma than WT plants. These results provide a foundation for breeding maize with low-lignin content and reveal novel insights about lignin regulation via genetic manipulation of CCoAOMT expression.     

Modifications in lignin and accumulation of phenolic glucosides in poplar xylem upon down-regulation of caffeoyl-coenzyme A O-methyltransferase, an enzyme involved in lignin biosynthesis

Caffeoyl-coenzyme A O-methyltransferase (CCoAOMT) methylates, in vitro, caffeoyl-CoA and 5-hydroxyferuloyl-CoA, two possible precursors in monolignol biosynthesis in vivo. To clarify the in vivo role of CCoAOMT in lignin biosynthesis, transgenic poplars with 10% residual CCoAOMT protein levels in the stem xylem were generated. Upon analysis of the xylem, the affected transgenic lines had a 12% reduced Klason lignin content, an 11% increased syringyl/guaiacyl ratio in the noncondensed lignin fraction, and an increase in lignin-attached p-hydroxybenzoate but otherwise a lignin composition similar to that of wild type. Stem xylem of the CCoAOMT-down-regulated lines had a pink-red coloration, which coincided with an enhanced fluorescence of mature vessel cell walls. The reduced production of CCoAOMT caused an accumulation of O(3)-beta-d-glucopyranosyl-caffeic acid, O(4)-beta-d-glucopyranosyl-vanillic acid, and O(4)-beta-d-glucopyranosyl-sinapic acid (GSA), as authenticated by (1)H NMR. Feeding experiments showed that O(3)-beta-d-glucopyranosyl-caffeic acid and GSA are storage or detoxification products of caffeic and sinapic acid, respectively. The observation that down-regulation of CCoAOMT decreases lignin amount whereas GSA accumulates to 10% of soluble phenolics indicates that endogenously produced sinapic acid is not a major precursor in syringyl lignin biosynthesis. Our in vivo results support the recently obtained in vitro enzymatic data that suggest that the route from caffeic acid to sinapic acid is not used for lignin biosynthesis.     

Profiling phenolic metabolites in transgenic alfalfa modified in lignin biosynthesis

Soluble phenolics, wall-bound phenolics and soluble and core lignin were analyzed in transgenic alfalfa with genetically down-regulated O-methyltransferase genes involved in lignin biosynthesis. High performance liquid chromatography and principal component analysis were used to distinguish metabolic phenotypes of different transgenic alfalfa genotypes growing under standard greenhouse conditions. Principal component analysis of HPLC chromatograms did not resolve differences in leaf metabolite profiles between wild-type and transgenic plants of the same genetic background, although stem phenolic profiles were clearly different between wild-type and transgenic plants. However, the analytical methods clearly differentiated two non-transgenic alfalfa cultivars based on either leaf or stem profiles. Metabolic profiling provides a useful approach to monitoring the broader biochemical phenotypes of transgenic plants with altered expression of lignin pathway enzymes.     

Dual methylation pathways in lignin biosynthesis

Caffeoyl-coenzyme A (CoA) O-methyltransferase (CCoAOMT) has been proposed to be involved in an alternative methylation pathway of lignin biosynthesis. However, no direct evidence has been available to confirm that CCoAOMT is essential for lignin biosynthesis. To understand further the methylation steps in lignin biosynthesis, we used an antisense approach to alter O-methyltransferase (OMT) gene expression and investigated the consequences of this alteration. We generated transgenic tobacco plants with a substantial reduction in CCoAOMT as well as plants with a simultaneous reduction in both CCoAOMT and caffeic acid O-methyltransferase (CAOMT). Lignin analysis showed that the reduction in CCoAOMT alone resulted in a dramatic decrease in lignin content. The reduction in CCoAOMT also led to a dramatic alteration in lignin composition. Both guaiacyl lignin and syringyl lignin were reduced in the transgenic plants. However, guaiacyl lignin was preferentially reduced, which resulted in an increase in the S/G (syringl/guaiacyl) ratio. We have also analyzed lignin content and composition in transgenic plants having a simultaneous reduction in both CCoAOMT and CAOMT. The reduction in both OMTs resulted in a further decrease in total lignin content. This is in sharp contrast to the effect that resulted from the reduction in CAOMT alone, which only decreased the syringl lignin unit without a reduction in overall lignin content. These results unequivocally demonstrate that methylation reactions in lignin biosynthesis are catalyzed by both CCoAOMT and CAOMT.     

Ensemble Modeling of Hepatic Fatty Acid Metabolism with a Synthetic Glyoxylate Shunt

The liver plays a central role in maintaining whole body metabolic and energy homeostasis by consuming and producing glucose and fatty acids. Glucose and fatty acids compete for hepatic substrate oxidation with regulation ensuring glucose is oxidized preferentially. Increasing fatty acid oxidation is expected to decrease lipid storage in the liver and avoid lipid-induced insulin-resistance. To increase hepatic lipid oxidation in the presence of glucose, we previously engineered a synthetic glyoxylate shunt into human hepatocyte cultures and a mouse model and showed that this synthetic pathway increases free fatty acid β-oxidation and confers resistance to diet-induced obesity in the mouse model. Here we used ensemble modeling to decipher the effects of perturbations to the hepatic metabolic network on fatty acid oxidation and glucose uptake. Despite sampling of kinetic parameters using the most fundamental elementary reaction models, the models based on current metabolic regulation did not readily describe the phenotype generated by glyoxylate shunt expression. Although not conclusive, this initial negative result prompted us to probe unknown regulations, and malate was identified as inhibitor of hexokinase 2 expression either through direct or indirect actions. This regulation allows the explanation of observed phenotypes (increased fatty acid degradation and decreased glucose consumption). Moreover, the result is a function of pyruvate-carboxylase, mitochondrial pyruvate transporter, citrate transporter protein, and citrate synthase activities. Some subsets of these flux ratios predict increases in fatty acid and decreases in glucose uptake after glyoxylate expression, whereas others predict no change. Altogether, this work defines the possible biochemical space where the synthetic shunt will produce the desired phenotype and demonstrates the efficacy of ensemble modeling for synthetic pathway design.     

Simultaneous regulation of F5H in COMT‐RNAi transgenic switchgrass alters effects of COMT suppression on syringyl lignin biosynthesis

Ferulate 5‐hydroxylase (F5H) catalyses the hydroxylation of coniferyl alcohol and coniferaldehyde for the biosynthesis of syringyl (S) lignin in angiosperms. However, the coordinated effects of F5H with caffeic acid O‐methyltransferase (COMT) on the metabolic flux towards S units are largely unknown. We concomitantly regulated F5H expression in COMT‐down‐regulated transgenic switchgrass (Panicum virgatum L.) lines and studied the coordination of F5H and COMT in lignin biosynthesis. Down‐regulation of F5H in COMT‐RNAi transgenic switchgrass plants further impeded S lignin biosynthesis and, consequently, increased guaiacyl (G) units and reduced 5‐OH G units. Conversely, overexpression of F5H in COMT‐RNAi transgenic plants reduced G units and increased 5‐OH units, whereas the deficiency of S lignin biosynthesis was partially compensated or fully restored, depending on the extent of COMT down‐regulation in switchgrass. Moreover, simultaneous regulation of F5H and COMT expression had different effects on cell wall digestibility of switchgrass without biomass loss. Our results indicate that up‐regulation and down‐regulation of F5H expression, respectively, have antagonistic and synergistic effects on the reduction in S lignin resulting from COMT suppression. The coordinated effects between lignin genes should be taken into account in future studies aimed at cell wall bioengineering.     

Metabolic ensemble modeling for strain engineers

Previous mathematical modeling efforts have made significant contributions to the development of systems biology for predicting biological behavior quantitatively. However, dynamic metabolic model construction remains challenging due to uncertainties in mechanistic structures and parameters. In addition, parameter estimation and model validation often require designated experiments conducted only for purpose of modeling. Such difficulties have hampered the progress of modeling in biology and biotechnology. To circumvent these problems, ensemble approaches have been used to account for uncertainties in model structure and parameters. Specifically, this review focuses on approaches that utilize readily available fermentation data for parameter screening and model validation. Time course data for metabolite measurements, if available, can further calibrate the model. The basis for this approach is explained in non-mathematical terms accessible to experimentalists. Information gained from such an approach has been shown to be useful in designing Escherichia coli strains for metabolic engineering and synthetic biology.     

Syringyl lignin is unaltered by severe sinapyl alcohol dehydrogenase suppression in tobacco

The manipulation of lignin could, in principle, facilitate efficient biofuel production from plant biomass. Despite intensive study of the lignin pathway, uncertainty exists about the enzyme catalyzing the last step in syringyl (S) monolignol biosynthesis, the reduction of sinapaldehyde to sinapyl alcohol. Traditional schemes of the pathway suggested that both guaiacyl (G) and S monolignols are produced by a single substrate-versatile enzyme, cinnamyl alcohol dehydrogenase (CAD). This was challenged by the discovery of a novel sinapyl alcohol dehydrogenase (SAD) that preferentially uses sinapaldehyde as a substrate and that was claimed to regulate S lignin biosynthesis in angiosperms. Consequently, most pathway schemes now show SAD (or SAD and CAD) at the sinapaldehyde reduction step, although functional evidence is lacking. We cloned SAD from tobacco (Nicotiana tabacum) and suppressed it in transgenic plants using RNA interference-inducing vectors. Characterization of lignin in the woody stems shows no change to content, composition, or structure, and S lignin is normal. By contrast, plants additionally suppressed in CAD have changes to lignin structure and S:G ratio and have increased sinapaldehyde in lignin, similar to plants suppressed in CAD alone. These data demonstrate that CAD, not SAD, is the enzyme responsible for S lignin biosynthesis in woody angiosperm xylem.     

ZmMYB31 directly represses maize lignin genes and redirects the phenylpropanoid metabolic flux

Few regulators of phenylpropanoids have been identified in monocots having potential as biofuel crops. Here we demonstrate the role of the maize (Zea mays) R2R3-MYB factor ZmMYB31 in the control of the phenylpropanoid pathway. We determined its in vitro consensus DNA-binding sequence as ACC(T)/(A) ACC, and chromatin immunoprecipitation (ChIP) established that it interacts with two lignin gene promoters in vivo. To explore the potential of ZmMYB31 as a regulator of phenylpropanoids in other plants, its role in the regulation of the phenylpropanoid pathway was further investigated in Arabidopsis thaliana. ZmMYB31 downregulates several genes involved in the synthesis of monolignols and transgenic plants are dwarf and show a significantly reduced lignin content with unaltered polymer composition. We demonstrate that these changes increase cell wall degradability of the transgenic plants. In addition, ZmMYB31 represses the synthesis of sinapoylmalate, resulting in plants that are more sensitive to UV irradiation, and induces several stress-related proteins. Our results suggest that, as an indirect effect of repression of lignin biosynthesis, transgenic plants redirect carbon flux towards the biosynthesis of anthocyanins. Thus, ZmMYB31 can be considered a good candidate for the manipulation of lignin biosynthesis in biotechnological applications.     

Modeling and simulation: tools for metabolic engineering.

Mathematical modeling is one of the key methodologies of metabolic engineering. Based on a given metabolic model different computational tools for the simulation, data evaluation, systems analysis, prediction, design and optimization of metabolic systems have been developed. The currently used metabolic modeling approaches can be subdivided into structural models, stoichiometric models, carbon flux models, stationary and nonstationary mechanistic models and models with gene regulation. However, the power of a model strongly depends on its basic modeling assumptions, the simplifications made and the data sources used. Model validation turns out to be particularly difficult for metabolic systems. The different modeling approaches are critically reviewed with respect to their potential and benefits for the metabolic engineering cycle. Several tools that have emerged from the different modeling approaches including structural pathway synthesis, stoichiometric pathway analysis, metabolic flux analysis, metabolic control analysis, optimization of regulatory architectures and the evaluation of rapid sampling experiments are discussed.     

MMT--a pathway modeling tool for data from rapid sampling experiments

The identification of metabolic regulation is a major concern in metabolic engineering. Metabolic regulation phenomena depend on intracellular compounds such as enzymes, metabolites and cofactors. A complete understanding of metabolic regulation requires quantitative information about these compounds under in vivo conditions. This quantitative knowledge in combination with the known network of metabolic pathways allows the construction of mathematical models that describe the dynamic changes in metabolite concentrations over time. Rapid sampling combined with pulse experiments is a useful tool for the identification of metabolic regulation owing to the transient data they provide. Enzymatic tests in combination with ESI-LC-MS (Electrospray Ionization Liquid Chromatographic Tandem Mass Spectrometry) and HPLC measurements have been used to identify up to 30 metabolites and nucleotides from rapid sampling experiments. A metabolic modeling tool (MMT) that is built on a relational database was developed specifically for analysis of rapid sampling experiments. The tool allows to construct complex pathway models with information stored in the relational database. Parameter fitting and simulation algorithms for the resulting system of Ordinary Differential Equations (ODEs) are part of MMT. Additionally explicit sensitivity functions are calculated. The integration of all necessary algorithms in one tool allows fast model analysis and comparison. Complex models have been developed to describe the central metabolic pathways of Escherichia coli during a glucose pulse experiment.     

Proteomic and metabolic disturbances in lignin-modified Brachypodium distachyon

Lignin biosynthesis begins with the deamination of phenylalanine and tyrosine (Tyr) as a key branch point between primary and secondary metabolism in land plants. Here, we used a systems biology approach to investigate the global metabolic responses to lignin pathway perturbations in the model grass Brachypodium distachyon. We identified the lignin biosynthetic protein families and found that ammonia-lyases (ALs) are among the most abundant proteins in lignifying tissues in grasses. Integrated metabolomic and proteomic data support a link between lignin biosynthesis and primary metabolism mediated by the ammonia released from ALs that is recycled for the synthesis of amino acids via glutamine. RNA interference knockdown of lignin genes confirmed that the route of the canonical pathway using shikimate ester intermediates is not essential for lignin formation in Brachypodium, and there is an alternative pathway from Tyr via sinapic acid for the synthesis of syringyl lignin involving yet uncharacterized enzymatic steps. Our findings support a model in which plant ALs play a central role in coordinating the allocation of carbon for lignin synthesis and the nitrogen available for plant growth. Collectively, these data also emphasize the value of integrative multiomic analyses to advance our understanding of plant metabolism.

Ammonia-lyases play a key role in coordinating the allocation of carbon for lignin synthesis and nitrogen availability for plant growth.     

Application of a generalized MWC model for the mathematical simulation of metabolic pathways regulated by allosteric enzymes

In our effort to elucidate the systems biology of the model organism, Escherichia coli, we have developed a mathematical model that simulates the allosteric regulation for threonine biosynthesis pathway starting from aspartate. To achieve this goal, we used kMech, a Cellerator language extension that describes enzyme mechanisms for the mathematical modeling of metabolic pathways. These mechanisms are converted by Cellerator into ordinary differential equations (ODEs) solvable by Mathematica. In this paper, we describe a more flexible model in Cellerator, which generalizes the Monod, Wyman, Changeux (MWC) model for enzyme allosteric regulation to allow for multiple substrate, activator and inhibitor binding sites. Furthermore, we have developed a model that describes the behavior of the bifunctional allosteric enzyme aspartate kinase I-homoserine dehydrogenase I (AKI-HDHI). This model predicts the partition of enzyme activities in the steady state which paves the way for a more generalized prediction of the behavior of bifunctional enzymes.     

O-Methylation of benzaldehyde derivatives by "lignin specific" caffeic acid 3-O-methyltransferase

Although S-adenosyl-l-methionine (SAM) dependent caffeic acid/5-hydroxyferulic acid 3/5-O-methyltransferase (COMT) is one of the key enzymes in lignin biosynthesis, the present work demonstrates that alfalfa COMT methylates benzaldehyde derivatives more efficiently than lignin pathway intermediates. 3,4-Dihydroxy, 5-methoxybenzaldehyde and protocatechuic aldehyde were the best in vitro substrates for OMT activity in extracts from developing alfalfa stems, and these compounds were preferred over lignin pathway intermediates for 3-O-methylation by recombinant alfalfa COMT expressed in Escherichia coli. OMT activity with benzaldehydes was strongly reduced in extracts from stems of transgenic alfalfa down-regulated in COMT. However, although COMT down-regulation drastically affects lignin composition, it does not appear to significantly impact metabolism of benzaldehyde derivatives in alfalfa. Structurally designed site-directed mutants of COMT showed altered relative substrate preferences for lignin precursors and benzaldehyde derivatives. Taken together, these results indicate that COMT may have more than one role in phenylpropanoid metabolism (but probably not in alfalfa), and that engineered COMT enzymes could be useful for metabolic engineering of both lignin and benzaldehyde-derived flavors and fragrances.