A dynamic model of lignin biosynthesis in Brachypodium distachyon

Modeling lignin biosynthesis: a pathway to renewable chemicals

Plant biomass contains lignin that can be converted into high-value-added chemicals, fuels, and materials. The precise genetic manipulation of lignin content and composition in plant cells offers substantial environmental and economic benefits. However, the intricate regulatory mechanisms governing lignin formation challenge the development of crops with specific lignin profiles. Mathematical models and computational simulations have recently been employed to gain fundamental understanding of the metabolism of lignin and related phenolic compounds. This review article discusses the strategies used for modeling plant metabolic networks, focusing on the application of mathematical modeling for flux network analysis in monolignol biosynthesis. Furthermore, we highlight how current challenges might be overcome to optimize the use of metabolic modeling approaches for developing lignin-engineered plants.

Hepatic cholesterol biosynthesis and dioxin-induced dysregulation: A multiscale computational approach

Humans are constantly exposed to lipophilic persistent organic pollutants (POPs) that accumulate in fatty foods. Among the numerous POPs, dioxins, in particular 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), can impact several organ systems. While the hazard is clearly recognized, it is still difficult to develop a comprehensive understanding of the overall health impacts of dioxins. As chemical toxicity testing is steadily adopting new approach methodologies (NAMs), it becomes imperative to develop computational models that can bridge the data gaps between in vitro testing and in vivo outcomes. As an effort to address this challenge, we propose a multiscale computational approach using a "template-and-anchor" (T&A) structure. A template is a high-level umbrella model that permits the integration of information from various, detailed anchor models. In the present study, we use this T&A approach to describe the effect of TCDD on cholesterol dynamics. Specifically, we represent hepatic cholesterol biosynthesis as an anchor model that is perturbed by TCDD, leading to steatosis, along with alterations of plasma cholesterol. In the future, incorporating pertinent information from all anchor models into the template model will allow the characterization of the global effects of dioxin, which can subsequently be translated into overall - and ultimately personalized - human health risk assessment.

Developmental changes in lignin composition are driven by both monolignol supply and laccase specificity

The factors controlling lignin composition remain unclear. Catechyl (C)-lignin is a homopolymer of caffeyl alcohol with unique properties as a biomaterial and precursor of industrial chemicals. The lignin synthesized in the seed coat of Cleome hassleriana switches from guaiacyl (G)- to C-lignin at around 12 to 14 days after pollination (DAP), associated with a rerouting of the monolignol pathway. Lack of synthesis of caffeyl alcohol limits C-lignin formation before around 12 DAP, but coniferyl alcohol is still synthesized and highly accumulated after 14 DAP. We propose a model in which, during C-lignin biosynthesis, caffeyl alcohol noncompetitively inhibits oxidation of coniferyl alcohol by cell wall laccases, a process that might limit movement of coniferyl alcohol to the apoplast. Developmental changes in both substrate availability and laccase specificity together account for the metabolic fates of G- and C-monolignols in the Cleome seed coat.

Systematic Analysis and Biochemical Characterization of the Caffeoyl Shikimate Esterase Gene Family in Poplar

Caffeoyl shikimate esterase (CSE) hydrolyzes caffeoyl shikimate into caffeate and shikimate in the phenylpropanoid pathway. In this study, we performed a systematic analysis of the CSE gene family and investigated the possible roles of CSE and CSE-like genes in Populus. We conducted a genome-wide analysis of the CSE gene family, including functional and phylogenetic analyses of CSE and CSE-like genes, using the poplar (Populus trichocarpa) genome. Eighteen CSE and CSE-like genes were identified in the Populus genome, and five phylogenetic groups were identified from phylogenetic analysis. CSEs in Group Ia, which were proposed as bona fide CSEs, have probably been lost in most monocots except Oryza sativa. Primary functional classification showed that PoptrCSE1 and PoptrCSE2 had putative function in lignin biosynthesis. In addition, PoptrCSE2, along with PoptrCSE12, might also respond to stress with a function in cell wall biosynthesis. Enzymatic assay of PoptoCSE1 (Populus tomentosa), -2 and -12 showed that PoptoCSE1 and -2 maintained CSE activity. PoptoCSE1 and 2 had similar biochemical properties, tissue expression patterns and subcellular localization. Most of the PoptrCSE-like genes are homologs of AtMAGL (monoacylglycerol lipase) genes in Arabidopsis and may function as MAG lipase in poplar. Our study provides a systematic understanding of this novel gene family and suggests the function of CSE in monolignol biosynthesis in Populus.

Modeling Plant Metabolism: From Network Reconstruction to Mechanistic Models

Mathematical modeling of plant metabolism enables the plant science community to understand the organization of plant metabolism, obtain quantitative insights into metabolic functions, and derive engineering strategies for manipulation of metabolism. Among the various modeling approaches, metabolic pathway analysis can dissect the basic functional modes of subsections of core metabolism, such as photorespiration, and reveal how classical definitions of metabolic pathways have overlapping functionality. In the many studies using constraint-based modeling in plants, numerous computational tools are currently available to analyze large-scale and genome-scale metabolic networks. For ¹³C-metabolic flux analysis, principles of isotopic steady state have been used to study heterotrophic plant tissues, while nonstationary isotope labeling approaches are amenable to the study of photoautotrophic and secondary metabolism. Enzyme kinetic models explore pathways in mechanistic detail, and we discuss different approaches to determine or estimate kinetic parameters. In this review, we describe recent advances and challenges in modeling plant metabolism.

Modeling cross-regulatory influences on monolignol transcripts and proteins under single and combinatorial gene knockdowns in Populus trichocarpa

Accurate manipulation of metabolites in monolignol biosynthesis is a key step for controlling lignin content, structure, and other wood properties important to the bioenergy and biomaterial industries. A crucial component of this strategy is predicting how single and combinatorial knockdowns of monolignol specific gene transcripts influence the abundance of monolignol proteins, which are the driving mechanisms of monolignol biosynthesis. Computational models have been developed to estimate protein abundances from transcript perturbations of monolignol specific genes. The accuracy of these models, however, is hindered by their inability to capture indirect regulatory influences on other pathway genes. Here, we examine the manifestation of these indirect influences on transgenic transcript and protein abundances, identifying putative indirect regulatory influences that occur when one or more specific monolignol pathway genes are perturbed. We created a computational model using sparse maximum likelihood to estimate the resulting monolignol transcript and protein abundances in transgenic Populus trichocarpa based on targeted knockdowns of specific monolignol genes. Using in-silico simulations of this model and root mean square error, we showed that our model more accurately estimated transcript and protein abundances, in comparison to previous models, when individual and families of monolignol genes were perturbed. We leveraged insight from the inferred network structure obtained from our model to identify potential genes, including PtrHCT, PtrCAD, and Ptr4CL, involved in post-transcriptional and/or post-translational regulation. Our model provides a useful computational tool for exploring the cascaded impact of single and combinatorial modifications of monolignol specific genes on lignin and other wood properties.

Lignin biosynthesis: old roads revisited and new roads explored

Lignin is a major component of secondarily thickened plant cell walls and is considered to be the second most abundant biopolymer on the planet. At one point believed to be the product of a highly controlled polymerization procedure involving just three potential monomeric components (monolignols), it is becoming increasingly clear that the composition of lignin is quite flexible. Furthermore, the biosynthetic pathways to the major monolignols also appear to exhibit flexibility, particularly as regards the early reactions leading to the formation of caffeic acid from coumaric acid. The operation of parallel pathways to caffeic acid occurring at the level of shikimate esters or free acids may help provide robustness to the pathway under different physiological conditions. Several features of the pathway also appear to link monolignol biosynthesis to both generation and detoxification of hydrogen peroxide, one of the oxidants responsible for creating monolignol radicals for polymerization in the apoplast. Monolignol transport to the apoplast is not well understood. It may involve passive diffusion, although this may be targeted to sites of lignin initiation/polymerization by ordered complexes of both biosynthetic enzymes on the cytosolic side of the plasma membrane and structural anchoring of proteins for monolignol oxidation and polymerization on the apoplastic side. We present several hypothetical models to illustrate these ideas and stimulate further research. These are based primarily on studies in model systems, which may or may not reflect the major lignification process in forest trees.

Low Lignin Mutants and Reduction of Lignin Content in Grasses for Increased Utilisation of Lignocellulose

Biomass rich in lignocellulose from grasses is a major source for biofuel production and animal feed. However, the presence of lignin in cell walls limits its efficient utilisation such as in its bioconversion to biofuel. Reduction of the lignin content or alteration of its structure in crop plants have been pursued, either by regulating genes encoding enzymes in the lignin biosynthetic pathway using biotechnological techniques or by breeding naturally-occurring low lignin mutant lines. The aim of this review is to provide a summary of these studies, focusing on lignin (monolignol) biosynthesis and composition in grasses and, where possible, the impact on recalcitrance to bioconversion. An overview of transgenic crops of the grass family with regulated gene expression in lignin biosynthesis is presented, including the effect on lignin content and changes in the ratio of p-hydroxyphenyl (H), guaiacyl (G) and syringyl (S) units. Furthermore, a survey is provided of low-lignin mutants in grasses, including cereals in particular, summarising their origin and phenotypic traits together with genetics and the molecular function of the various genes identified.

The regulatory role of shikimate in plant phenylalanine metabolism

In higher plants, the amino acid phenylalanine is a substrate of both primary and secondary metabolic pathways. The primary pathway that consumes phenylalanine, protein biosynthesis, is essential for the viability of all cells. Meanwhile, the secondary pathways are not necessary for the survival of individual cells, but benefit of the plant as a whole. Here we focus on the monolignol pathway, a secondary metabolic pathway in the cytosol that rapidly consumes phenylalanine to produce the precursors of lignin during wood formation. In planta monolignol biosynthesis involves a series of seemingly redundant steps wherein shikimate, a precursor of phenylalanine synthesized in the plastid, is transiently ligated to the main substrate of the pathway. However, shikimate is not catalytically involved in the reactions of the monolignol pathway, and is only needed for pathway enzymes to recognize their main substrates. After some steps the shikimate moiety is removed unaltered, and the main substrate continues along the pathway. It has been suggested that this portion of the monolignol pathway fulfills a regulatory role in the following way. Low phenylalanine concentrations (viz. availability) correlate with low shikimate concentrations. When shikimate concentratios are low, flux into the monolignol pathway will be limited by means of the steps requiring shikimate. Thus, when the concentration of phenylalanine is low it will be reserved for protein biosynthesis. Here we employ a theoretical approach to test this hypothesis. Simplified versions of plant phenylalanine metabolism are modelled as systems of ordinary differential equations. Our analysis shows that the seemingly redundant steps can be sufficient for the prioritization of protein biosynthesis over the monolignol pathway when the availability of phenylalanine is low, depending on system parameters. Thus, the phenylalanine precursor shikimate may signal low phenylalanine availability to secondary pathways. Because our models have been abstracted from plant phenylalanine metabolism, this mechanism of metabolic signalling, which we call the Precursor Shutoff Valve (PSV), may also be present in other biochemical networks comprised of two pathways that share a common substrate.