Metabolic trade-offs between biomass synthesis and photosynthate export at different light intensities in a genome-scale metabolic model of rice

Multiscale plant modeling: from genome to phenome and beyond

Plants are complex organisms that adapt to changes in their environment using an array of regulatory mechanisms that span across multiple levels of biological organization. Due to this complexity, it is difficult to predict emergent properties using conventional approaches that focus on single levels of biology such as the genome, transcriptome, or metabolome. Mathematical models of biological systems have emerged as useful tools for exploring pathways and identifying gaps in our current knowledge of biological processes. Identification of emergent properties, however, requires their vertical integration across biological scales through multiscale modeling. Multiscale models that capture and predict these emergent properties will allow us to predict how plants will respond to a changing climate and explore strategies for plant engineering. In this review, we (1) summarize the recent developments in plant multiscale modeling; (2) examine multiscale models of microbial systems that offer insight to potential future directions for the modeling of plant systems; (3) discuss computational tools and resources for developing multiscale models; and (4) examine future directions of the field.

Environment-coupled models of leaf metabolism

The plant leaf is the main site of photosynthesis. This process converts light energy and inorganic nutrients into chemical energy and organic building blocks for the biosynthesis and maintenance of cellular components and to support the growth of the rest of the plant. The leaf is also the site of gas–water exchange and due to its large surface, it is particularly vulnerable to pathogen attacks. Therefore, the leaf's performance and metabolic modes are inherently determined by its interaction with the environment. Mathematical models of plant metabolism have been successfully applied to study various aspects of photosynthesis, carbon and nitrogen assimilation and metabolism, aided suggesting metabolic intervention strategies for optimized leaf performance, and gave us insights into evolutionary drivers of plant metabolism in various environments. With the increasing pressure to improve agricultural performance in current and future climates, these models have become important tools to improve our understanding of plant–environment interactions and to propel plant breeders efforts. This overview article reviews applications of large-scale metabolic models of leaf metabolism to study plant–environment interactions by means of flux-balance analysis. The presented studies are organized in two ways — by the way the environment interactions are modelled — via external constraints or data-integration and by the studied environmental interactions — abiotic or biotic.

Biochemical Characteristics and a Genome-Scale Metabolic Model of an Indian Euryhaline Cyanobacterium with High Polyglucan Content

Marine cyanobacteria are promising microbes to capture and convert atmospheric CO₂ and light into biomass and valuable industrial bio-products. Yet, reports on metabolic characteristics of non-model cyanobacteria are scarce. In this report, we show that an Indian euryhaline Synechococcus sp. BDU 130192 has biomass accumulation comparable to a model marine cyanobacterium and contains approximately double the amount of total carbohydrates, but significantly lower protein levels compared to Synechococcus sp. PCC 7002 cells. Based on its annotated chromosomal genome sequence, we present a genome scale metabolic model (GSMM) of this cyanobacterium, which we have named as iSyn706. The model includes 706 genes, 908 reactions, and 900 metabolites. The difference in the flux balance analysis (FBA) predicted flux distributions between Synechococcus sp. PCC 7002 and Synechococcus sp. BDU130192 strains mimicked the differences in their biomass compositions. Model-predicted oxygen evolution rate for Synechococcus sp. BDU130192 was found to be close to the experimentally-measured value. The model was analyzed to determine the potential of the strain for the production of various industrially-useful products without affecting growth significantly. This model will be helpful to researchers interested in understanding the metabolism as well as to design metabolic engineering strategies for the production of industrially-relevant compounds.

Leaf Energy Balance Requires Mitochondrial Respiration and Export of Chloroplast NADPH in the Light

Key aspects of leaf mitochondrial metabolism in the light remain unresolved. For example, there is debate about the relative importance of exporting reducing equivalents from mitochondria for the peroxisomal steps of photorespiration versus oxidation of NADH to generate ATP by oxidative phosphorylation. Here, we address this and explore energetic coupling between organelles in the light using a diel flux balance analysis model. The model included more than 600 reactions of central metabolism with full stoichiometric accounting of energy production and consumption. Different scenarios of energy availability (light intensity) and demand (source leaf versus a growing leaf) were considered, and the model was constrained by the nonlinear relationship between light and CO2 assimilation rate. The analysis demonstrated that the chloroplast can theoretically generate sufficient ATP to satisfy the energy requirements of the rest of the cell in addition to its own. However, this requires unrealistic high light use efficiency and, in practice, the availability of chloroplast-derived ATP is limited by chloroplast energy dissipation systems, such as nonphotochemical quenching, and the capacity of the chloroplast ATP export shuttles. Given these limitations, substantial mitochondrial ATP synthesis is required to fulfill cytosolic ATP requirements, with only minimal, or zero, export of mitochondrial reducing equivalents. The analysis also revealed the importance of exporting reducing equivalents from chloroplasts to sustain photorespiration. Hence, the chloroplast malate valve and triose phosphate-3-phosphoglycerate shuttle are predicted to have important metabolic roles, in addition to their more commonly discussed contribution to the avoidance of photooxidative stress.

Stoichiometric analysis of the energetics and metabolic impact of photorespiration in C3 plants

Analysis of the impact of photorespiration on plant metabolism is usually based on manual inspection of small network diagrams. Here we create a structural metabolic model that contains the reactions that participate in photorespiration in the plastid, peroxisome, mitochondrion and cytosol, and the metabolite exchanges between them. This model was subjected to elementary flux modes analysis, a technique that enumerates all the component, minimal pathways of a network. Any feasible photorespiratory metabolism in the plant will be some combination of the elementary flux modes (EFMs) that contain the Rubisco oxygenase reaction. Amongst the EFMs we obtained was the classic photorespiratory cycle, but there were also modes that involve photorespiration coupled with mitochondrial metabolism and ATP production, the glutathione-ascorbate cycle and nitrate reduction to ammonia. The modes analysis demonstrated the underlying basis of the metabolic linkages with photorespiration that have been inferred experimentally. The set of reactions common to all the elementary modes showed good agreement with the gene products of mutants that have been reported to have a defective phenotype in photorespiratory conditions. Finally, the set of modes provided a formal demonstration that photorespiration itself does not impact on the CO₂ :O₂ ratio (assimilation quotient), except in those modes associated with concomitant nitrate reduction.

Using metabolome data for mathematical modeling of plant metabolic systems

Plant metabolism is characterized by a wide diversity of metabolites, with systems far more complicated than those of microorganisms. Mathematical modeling is useful for understanding dynamic behaviors of plant metabolic systems for metabolic engineering. Time-series metabolome data has great potential for estimating kinetic model parameters to construct a genome-wide metabolic network model. However, data obtained by current metabolomics techniques does not meet the requirement for constructing accurate models. In this article, we highlight novel strategies and algorithms to handle the underlying difficulties and construct dynamic in vivo models for large-scale plant metabolic systems. The coarse but efficient modeling enables the prediction of unknown mechanisms regulating plant metabolism.

Modeling Rice Metabolism: From Elucidating Environmental Effects on Cellular Phenotype to Guiding Crop Improvement

Crop productivity is severely limited by various biotic and abiotic stresses. Thus, it is highly needed to understand the underlying mechanisms of environmental stress response and tolerance in plants, which could be addressed by systems biology approach. To this end, high-throughput omics profiling and in silico modeling can be considered to explore the environmental effects on phenotypic states and metabolic behaviors of rice crops at the systems level. Especially, the advent of constraint-based metabolic reconstruction and analysis paves a way to characterize the plant cellular physiology under various stresses by combining the mathematical network models with multi-omics data. Rice metabolic networks have been reconstructed since 2013 and currently six such networks are available, where five are at genome-scale. Since their publication, these models have been utilized to systematically elucidate the rice abiotic stress responses and identify agronomic traits for crop improvement. In this review, we summarize the current status of the existing rice metabolic networks and models with their applications. Furthermore, we also highlight future directions of rice modeling studies, particularly stressing how these models can be used to contextualize the affluent multi-omics data that are readily available in the public domain. Overall, we envisage a number of studies in the future, exploiting the available metabolic models to enhance the yield and quality of rice and other food crops.

A Method of Accounting for Enzyme Costs in Flux Balance Analysis Reveals Alternative Pathways and Metabolite Stores in an Illuminated Arabidopsis Leaf

Flux balance analysis of plant metabolism is an established method for predicting metabolic flux phenotypes and for exploring the way in which the plant metabolic network delivers specific outcomes in different cell types, tissues, and temporal phases. A recurring theme is the need to explore the flexibility of the network in meeting its objectives and, in particular, to establish the extent to which alternative pathways can contribute to achieving specific outcomes. Unfortunately, predictions from conventional flux balance analysis minimize the simultaneous operation of alternative pathways, but by introducing flux-weighting factors to allow for the variable intrinsic cost of supporting each flux, it is possible to activate different pathways in individual simulations and, thus, to explore alternative pathways by averaging thousands of simulations. This new method has been applied to a diel genome-scale model of Arabidopsis (Arabidopsis thaliana) leaf metabolism to explore the flexibility of the network in meeting the metabolic requirements of the leaf in the light. This identified alternative flux modes in the Calvin-Benson cycle revealed the potential for alternative transitory carbon stores in leaves and led to predictions about the light-dependent contribution of alternative electron flow pathways and futile cycles in energy rebalancing. Notable features of the analysis include the light-dependent tradeoff between the use of carbohydrates and four-carbon organic acids as transitory storage forms and the way in which multiple pathways for the consumption of ATP and NADPH can contribute to the balancing of the requirements of photosynthetic metabolism with the energy available from photon capture.