Fluxes through plant metabolic networks: measurements, predictions, insights and challenges

Non-canonical plant metabolism

Metabolism is essential for plant growth and has become a major target for crop improvement by enhancing nutrient use efficiency. Metabolic engineering is also the basis for producing high-value plant products such as pharmaceuticals, biofuels and industrial biochemicals. An inherent problem for such engineering endeavours is the tendency to view metabolism as a series of distinct metabolic pathways-glycolysis, the tricarboxylic acid cycle, the Calvin-Benson cycle and so on. While these canonical pathways may represent a dominant or frequently occurring flux mode, systematic analyses of metabolism via computational modelling have emphasized the inherent flexibility of the metabolic network to carry flux distributions that are distinct from the canonical pathways. Recent experimental estimates of metabolic network fluxes using 13C-labelling approaches have revealed numerous instances in which non-canonical pathways occur under different conditions and in different tissues. In this Review, we bring these non-canonical pathways to the fore, summarizing the evidence for their occurrence and the context in which they operate. We also emphasize the importance of non-canonical pathways for metabolic engineering. We argue that the introduction of a high-flux pathway to a desired metabolic product will, by necessity, require non-canonical supporting fluxes in central metabolism to provide the necessary carbon skeletons, energy and reducing power. We illustrate this using the overproduction of isoprenoids and fatty acids as case studies.

Microalgae growth-promoting bacteria for cultivation strategies: Recent updates and progress

Microalgae growth-promoting bacteria (MGPB), both actinobacteria and non-actinobacteria, have received considerable attention recently because of their potential to develop microalgae-bacteria co-culture strategies for improved efficiency and sustainability of the water-energy-environment nexus. Owing to their diverse metabolic pathways and ability to adapt to diverse conditions, microalgal-MGPB co-cultures could be promising biological systems under uncertain environmental and nutrient conditions. This review proposes the recent updates and progress on MGPB for microalgae cultivation through co-culture strategies. Firstly, potential MGPB strains for microalgae cultivation are introduced. Following, microalgal-MGPB interaction mechanisms and applications of their co-cultures for biomass production and wastewater treatment are reviewed. Moreover, state-of-the-art studies on synthetic biology and metabolic network analysis, along with the challenges and prospects of opting these approaches for microalgal-MGPB co-cultures are presented. It is anticipated that these strategies may significantly improve the sustainability of microalgal-MGPB co-cultures for wastewater treatment, biomass valorization, and bioproducts synthesis in a circular bioeconomy paradigm.

Isotopically nonstationary metabolic flux analysis of plants: recent progress and future opportunities

Metabolic flux analysis (MFA) is a valuable tool for quantifying cellular phenotypes and to guide plant metabolic engineering. By introducing stable isotopic tracers and employing mathematical models, MFA can quantify the rates of metabolic reactions through biochemical pathways. Recent applications of isotopically nonstationary MFA (INST-MFA) to plants have elucidated nonintuitive metabolism in leaves under optimal and stress conditions, described coupled fluxes for fast-growing algae, and produced a synergistic multi-organ flux map that is a first in MFA for any biological system. These insights could not be elucidated through other approaches and show the potential of INST-MFA to correct an oversimplified understanding of plant metabolism.

Engineering Gut Symbionts: A Way to Promote Bee Growth?

Bees play a crucial role as pollinators, contributing significantly to ecosystems. However, the honeybee population faces challenges such as global warming, pesticide use, and pathogenic microorganisms. Promoting bee growth using several approaches is therefore crucial for maintaining their roles. To this end, the bacterial microbiota is well-known for its native role in supporting bee growth in several respects. Maximizing the capabilities of these microorganisms holds the theoretical potential to promote the growth of bees. Recent advancements have made it feasible to achieve this enhancement through the application of genetic engineering. In this review, we present the roles of gut symbionts in promoting bee growth and collectively summarize the engineering approaches that would be needed for future applications. Particularly, as the engineering of bee gut symbionts has not been advanced, the dominant gut symbiotic bacteria Snodgrassella alvi and Gilliamella apicola are the main focus of the paper, along with other dominant species. Moreover, we propose engineering strategies that will allow for the improvement in bee growth with listed gene targets for modification to further encourage the use of engineered gut symbionts to promote bee growth.

Genome-scale metabolic modelling of extremophiles and its applications in astrobiological environments

Metabolic modelling approaches have become the powerful tools in modern biology. These mathematical models are widely used to predict metabolic phenotypes of the organisms or communities of interest, and to identify metabolic targets in metabolic engineering. Apart from a broad range of industrial applications, the possibility of using metabolic modelling in the contexts of astrobiology are poorly explored. In this mini-review, we consolidated the concepts and related applications of applying metabolic modelling in studying organisms in space-related environments, specifically the extremophilic microbes. We recapitulated the current state of the art in metabolic modelling approaches and their advantages in the astrobiological context. Our review encompassed the applications of metabolic modelling in the theoretical investigation of the origin of life within prebiotic environments, as well as the compilation of existing uses of genome-scale metabolic models of extremophiles. Furthermore, we emphasize the current challenges associated with applying this technique in extreme environments, and conclude this review by discussing the potential implementation of metabolic models to explore theoretically optimal metabolic networks under various space conditions. Through this mini-review, our aim is to highlight the potential of metabolic modelling in advancing the study of astrobiology.

Daylength variation affects growth, photosynthesis, leaf metabolism, partitioning, and metabolic fluxes

Daylength, a seasonal and latitudinal variable, exerts a substantial impact on plant growth. However, the relationship between daylength and growth is nonproportional, suggesting the existence of adaptive mechanisms. Thus, our study aimed to comprehensively investigate the adaptive strategies employed by plants in response to daylength variation. We grew false flax (Camelina sativa) plants, a model oilseed crop, under long-day (LD) and short-day (SD) conditions and used growth measurements, gas exchange measurements, and isotopic labeling techniques, including 13C, 14C, and 2H2O, to determine responses to different daylengths. Our findings revealed that daylength influences various growth parameters, photosynthetic physiology, carbon partitioning, metabolic fluxes, and metabolite levels. SD plants employed diverse mechanisms to compensate for reduced CO2 fixation in the shorter photoperiod. These mechanisms included enhanced photosynthetic rates and reduced respiration in the light (RL), leading to increased shoot investment. Additionally, SD plants exhibited reduced rates of the glucose 6-phosphate (G6P) shunt and greater partitioning of sugars into starch, thereby sustaining carbon availability during the longer night. Isotopic labeling results further demonstrated substantial alterations in the partitioning of amino acids and TCA cycle intermediates between rapidly and slowly turning over pools. Overall, the results point to multiple developmental, physiological, and metabolic ways in which plants adapt to different daylengths to maintain growth.

Metabolic network reconstruction of Euglena gracilis: Current state, challenges, and applications

A metabolic model, representing all biochemical reactions in a cell, is a prerequisite for several approaches in systems biology used to explore the metabolic phenotype of an organism. Despite the use of Euglena in diverse industrial applications and as a biological model, there is limited understanding of its metabolic network capacity. The unavailability of the completed genome data and the highly complex evolution of Euglena are significant obstacles to the reconstruction and analysis of its genome-scale metabolic model. In this mini-review, we discuss the current state and challenges of metabolic network reconstruction in Euglena gracilis. We have collated and present the available relevant data for the metabolic network reconstruction of E. gracilis, which could be used to improve the quality of the metabolic model of E. gracilis. Furthermore, we deliver the potential applications of the model in metabolic engineering. Altogether, it is supposed that this mini-review would facilitate the investigation of metabolic networks in Euglena and further lay out a direction for model-assisted metabolic engineering.

Isotopically non-stationary metabolic flux analysis of heterotrophic Arabidopsis thaliana cell cultures

Fluxes are the ultimate phenotype of metabolism and their accurate quantification is fundamental to any understanding of metabolic networks. Steady state metabolic flux analysis has been the method of choice for quantifying fluxes in heterotrophic cells, but it is unable to measure fluxes during short-lived metabolic states, such as a transient oxidative load. Isotopically non-stationary metabolic flux analysis (INST-MFA) can be performed over shorter timescales (minutes - hours) and might overcome this limitation. INST-MFA has recently been applied to photosynthesising leaves, but agriculturally important tissues such as roots and storage organs, or plants during the night are heterotrophic. Here we outline the application of INST-MFA to heterotrophic plant cells. Using INST-MFA we were able to identify changes in the fluxes supported by phosphoenolpyruvate carboxylase and malic enzyme under oxidative load, highlighting the potential of INST-MFA to measure fluxes during short-lived metabolic states. We discuss the challenges in applying INST-MFA, and highlight further development required before it can be routinely used to quantify fluxes in heterotrophic plant cells.

15N-labelling of Leaves Combined with GC-MS Analysis as a Tool for Monitoring the Dynamics of Nitrogen Incorporation into Amino Acids

Labeling plant material such as detached leaves with ¹⁵NH₄⁺ is a very instrumental method for the characterization of metabolic pathways of mineral nitrogen assimilation and incorporation into amino acids. A procedure of labeling, followed by amino acid extraction, purification, and derivatization for gas chromatography coupled to mass spectrometry (GC/MS) analysis, is presented. The rationale of heavy isotope abundance calculations and amino acid quantification is detailed. This method is adaptable to various plant species and various kinds of investigations, such as elucidating physiological changes occurring as a result of gene mutations (overexpression or inhibition) in natural variants or genetically modified crops, or characterization of metabolic fluxes in genotypes exhibiting contrasted physiological or developmental adaptive responses to biotic and/or abiotic environmental stresses. Furthermore, the benefit of working on detached organs or pieces of organs is to investigate finely the metabolism of species that are not amenable to laboratory work, such as plants growing in natural environments or under agricultural conditions in the field.

Deuterium in marine organic biomarkers: toward a new tool for quantifying aquatic mixotrophy

The traditional separation between primary producers (autotrophs) and consumers (heterotrophs) at the base of the marine food web is being increasingly replaced by the paradigm that mixoplankton, planktonic protists with the nutritional ability to use both phago(hetero)trophy and photo(auto)trophy to access energy are widespread globally. Thus, many 'phytoplankton' eat, while 50% of 'protozooplankton' also perform photosynthesis. Mixotrophy may enhance primary production, biomass transfer to higher trophic levels and the efficiency of the biological pump to sequester atmospheric CO₂ into the deep ocean. Although this view is gaining traction, science lacks a tool to quantify the relative contributions of autotrophy and heterotrophy in planktonic protists. This hinders our understanding of their impacts on carbon cycling within marine pelagic ecosystems. It has been shown that the hydrogen (H) isotopic signature of lipids is uniquely sensitive to heterotrophy relative to autotrophy in plants and bacteria. Here, we explored whether it is also sensitive to the trophic status in protists. The new understanding of H isotope signature of lipid biomarkers suggests it offers great potential as a novel tool for quantifying the prevalence of mixotrophy in diverse marine microorganisms and thus for investigating the implications of the 'mixoplankton' paradigm.

Integration of relative metabolomics and transcriptomics time-course data in a metabolic model pinpoints effects of ribosome biogenesis defects on Arabidopsis thaliana metabolism

Ribosome biogenesis is tightly associated to plant metabolism due to the usage of ribosomes in the synthesis of proteins necessary to drive metabolic pathways. Given the central role of ribosome biogenesis in cell physiology, it is important to characterize the impact of different components involved in this process on plant metabolism. Double mutants of the Arabidopsis thaliana cytosolic 60S maturation factors REIL1 and REIL2 do not resume growth after shift to moderate 10 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{\circ }\hbox {C}$$\end{document}∘C chilling conditions. To gain mechanistic insights into the metabolic effects of this ribosome biogenesis defect on metabolism, we developed TC-iReMet2, a constraint-based modelling approach that integrates relative metabolomics and transcriptomics time-course data to predict differential fluxes on a genome-scale level. We employed TC-iReMet2 with metabolomics and transcriptomics data from the Arabidopsis Columbia 0 wild type and the reil1-1 reil2-1 double mutant before and after cold shift. We identified reactions and pathways that are highly altered in a mutant relative to the wild type. These pathways include the Calvin–Benson cycle, photorespiration, gluconeogenesis, and glycolysis. Our findings also indicated differential NAD(P)/NAD(P)H ratios after cold shift. TC-iReMet2 allows for mechanistic hypothesis generation and interpretation of system biology experiments related to metabolic fluxes on a genome-scale level.

Modelling Cell Metabolism: A Review on Constraint-Based Steady-State and Kinetic Approaches

Studying cell metabolism serves a plethora of objectives such as the enhancement of bioprocess performance, and advancement in the understanding of cell biology, of drug target discovery, and in metabolic therapy. Remarkable successes in these fields emerged from heuristics approaches, for instance, with the introduction of effective strategies for genetic modifications, drug developments and optimization of bioprocess management. However, heuristics approaches have showed significant shortcomings, such as to describe regulation of metabolic pathways and to extrapolate experimental conditions. In the specific case of bioprocess management, such shortcomings limit their capacity to increase product quality, while maintaining desirable productivity and reproducibility levels. For instance, since heuristics approaches are not capable of prediction of the cellular functions under varying experimental conditions, they may lead to sub-optimal processes. Also, such approaches used for bioprocess control often fail in regulating a process under unexpected variations of external conditions. Therefore, methodologies inspired by the systematic mathematical formulation of cell metabolism have been used to address such drawbacks and achieve robust reproducible results. Mathematical modelling approaches are effective for both the characterization of the cell physiology, and the estimation of metabolic pathways utilization, thus allowing to characterize a cell population metabolic behavior. In this article, we present a review on methodology used and promising mathematical modelling approaches, focusing primarily to investigate metabolic events and regulation. Proceeding from a topological representation of the metabolic networks, we first present the metabolic modelling approaches that investigate cell metabolism at steady state, complying to the constraints imposed by mass conservation law and thermodynamics of reactions reversibility. Constraint-based models (CBMs) are reviewed highlighting the set of assumed optimality functions for reaction pathways. We explore models simulating cell growth dynamics, by expanding flux balance models developed at steady state. Then, discussing a change of metabolic modelling paradigm, we describe dynamic kinetic models that are based on the mathematical representation of the mechanistic description of nonlinear enzyme activities. In such approaches metabolic pathway regulations are considered explicitly as a function of the activity of other components of metabolic networks and possibly far from the metabolic steady state. We have also assessed the significance of metabolic model parameterization in kinetic models, summarizing a standard parameter estimation procedure frequently employed in kinetic metabolic modelling literature. Finally, some optimization practices used for the parameter estimation are reviewed.

Multiscale computational models can guide experimentation and targeted measurements for crop improvement

Computational models of plants have identified gaps in our understanding of biological systems, and have revealed ways to optimize cellular processes or organ-level architecture to increase productivity. Thus, computational models are learning tools that help direct experimentation and measurements. Models are simplifications of complex systems, and often simulate specific processes at single scales (e.g. temporal, spatial, organizational, etc.). Consequently, single-scale models are unable to capture the critical cross-scale interactions that result in emergent properties of the system. In this perspective article, we contend that to accurately predict how a plant will respond in an untested environment, it is necessary to integrate mathematical models across biological scales. Computationally mimicking the flow of biological information from the genome to the phenome is an important step in discovering new experimental strategies to improve crops. A key challenge is to connect models across biological, temporal and computational (e.g. CPU versus GPU) scales, and then to visualize and interpret integrated model outputs. We address this challenge by describing the efforts of the international Crops in silico consortium.

Enzymes as Parts in Need of Replacement - and How to Extend Their Working Life

Enzymes catalyze reactions in vivo at different rates and each enzyme molecule has a lifetime limit before it is degraded and replaced to enable catalysis to continue. Considering these rates together as a unitless ratio of catalytic cycles until replacement (CCR) provides a new quantitative tool to assess the replacement schedule of and energy investment into enzymes as they relate to function. Here, we outline the challenges of determining CCRs and new approaches to overcome them and then assess the CCRs of selected enzymes in bacteria and plants to reveal a range of seven orders of magnitude for this ratio. Modifying CCRs in plants holds promise to lower cellular costs, to tailor enzymes for particular environments, and to breed enzyme improvements for crop productivity.

Metabolite database for root, tuber, and banana crops to facilitate modern breeding in understudied crops

Roots, tubers, and bananas (RTB) are vital staples for food security in the world's poorest nations. A major constraint to current RTB breeding programmes is limited knowledge on the available diversity due to lack of efficient germplasm characterization and structure. In recent years large-scale efforts have begun to elucidate the genetic and phenotypic diversity of germplasm collections and populations and, yet, biochemical measurements have often been overlooked despite metabolite composition being directly associated with agronomic and consumer traits. Here we present a compound database and concentration range for metabolites detected in the major RTB crops: banana (Musa spp.), cassava (Manihot esculenta), potato (Solanum tuberosum), sweet potato (Ipomoea batatas), and yam (Dioscorea spp.), following metabolomics-based diversity screening of global collections held within the CGIAR institutes. The dataset including 711 chemical features provides a valuable resource regarding the comparative biochemical composition of each RTB crop and highlights the potential diversity available for incorporation into crop improvement programmes. Particularly, the tropical crops cassava, sweet potato and banana displayed more complex compositional metabolite profiles with representations of up to 22 chemical classes (unknowns excluded) than that of potato, for which only metabolites from 10 chemical classes were detected. Additionally, over 20% of biochemical signatures remained unidentified for every crop analyzed. Integration of metabolomics with the on-going genomic and phenotypic studies will enhance 'omics-wide associations of molecular signatures with agronomic and consumer traits via easily quantifiable biochemical markers to aid gene discovery and functional characterization.

Euglena Central Metabolic Pathways and Their Subcellular Locations

Euglenids are a group of algae of great interest for biotechnology, with a large and complex metabolic capability. To study the metabolic network, it is necessary to know where the component enzymes are in the cell, but despite a long history of research into Euglena, the subcellular locations of many major pathways are only poorly defined. Euglena is phylogenetically distant from other commonly studied algae, they have secondary plastids bounded by three membranes, and they can survive after destruction of their plastids. These unusual features make it difficult to assume that the subcellular organization of the metabolic network will be equivalent to that of other photosynthetic organisms. We analysed bioinformatic, biochemical, and proteomic information from a variety of sources to assess the subcellular location of the enzymes of the central metabolic pathways, and we use these assignments to propose a model of the metabolic network of Euglena. Other than photosynthesis, all major pathways present in the chloroplast are also present elsewhere in the cell. Our model demonstrates how Euglena can synthesise all the metabolites required for growth from simple carbon inputs, and can survive in the absence of chloroplasts.

Systems Biology of Plant-Microbiome Interactions

In natural environments, plants are exposed to diverse microbiota that they interact with in complex ways. While plant-pathogen interactions have been intensely studied to understand defense mechanisms in plants, many microbes and microbial communities can have substantial beneficial effects on their plant host. Such beneficial effects include improved acquisition of nutrients, accelerated growth, resilience against pathogens, and improved resistance against abiotic stress conditions such as heat, drought, and salinity. However, the beneficial effects of bacterial strains or consortia on their host are often cultivar and species specific, posing an obstacle to their general application. Remarkably, many of the signals that trigger plant immune responses are molecularly highly similar and often identical in pathogenic and beneficial microbes. Thus, it is unclear what determines the outcome of a particular microbe-host interaction and which factors enable plants to distinguish beneficials from pathogens. To unravel the complex network of genetic, microbial, and metabolic interactions, including the signaling events mediating microbe-host interactions, comprehensive quantitative systems biology approaches will be needed.

Shoot tip culture: a step towards 13C metabolite flux analysis of sink leaf metabolism

BACKGROUND

Better understanding of the physiological and metabolic status of plants can only be obtained when metabolic fluxes are accurately assessed in a growing plant. Steady state ¹³C-MFA has been established as a routine method for analysis of fluxes in plant primary metabolism. However, the experimental system needs to be improved for continuous carbon enrichment from labelled sugars into metabolites for longer periods until complex secondary metabolism reaches steady state.

RESULTS

We developed an in vitro plant culture strategy by using peppermint as a model plant with minimizing unlabelled carbon fixation where growing shoot tip was strongly dependent on labelled glucose for their carbon necessity. We optimized the light condition and detected the satisfactory phenotypical growth under the lower light intensity. Total volatile terpenes were also highest at the same light. Analysis of label incorporation into pulegone monoterpene after continuous U-¹³C₆ glucose feeding revealed nearly 100% ¹³C, even at 15 days after first leaf visibility (DALV). Label enrichment was gradually scrambled with increasing light intensity and leaf age. This study was validated by showing similar levels of label enrichment in proteinogenic amino acids. The efficiency of this method was also verified in oregano.

CONCLUSIONS

Our shoot tip culture depicted a method in achieving long term, stable and a high percentage of label accumulation in secondary metabolites within a fully functional growing plant system. It recommends the potential application for the investigations of various facets of plant metabolism by steady state ¹³C-MFA. The system also provides a greater potential to study sink leaf metabolism. Overall, we propose a system to accurately describe complex metabolic phenotypes in a growing plant.

Advances in metabolic flux analysis toward genome-scale profiling of higher organisms

Methodological and technological advances have recently paved the way for metabolic flux profiling in higher organisms, like plants. However, in comparison with omics technologies, flux profiling has yet to provide comprehensive differential flux maps at a genome-scale and in different cell types, tissues, and organs. Here we highlight the recent advances in technologies to gather metabolic labeling patterns and flux profiling approaches. We provide an opinion of how recent local flux profiling approaches can be used in conjunction with the constraint-based modeling framework to arrive at genome-scale flux maps. In addition, we point at approaches which use metabolomics data without introduction of label to predict either non-steady state fluxes in a time-series experiment or flux changes in different experimental scenarios. The combination of these developments allows an experimentally feasible approach for flux-based large-scale systems biology studies.

Labeling Maize (Zea mays L.) Leaves with 15 NH4+ and Monitoring Nitrogen Incorporation into Amino Acids by GC/MS Analysis

The human body contains approximately 3.2% nitrogen (N), mainly present as protein and amino acids. Although N exists at a high concentration (78%) in the air, it is not readily available to animals and most plants. Plants are however able to take up both nitrate (NO₃^- ) and ammonium (NH₄⁺ ) ions from the soil and convert them to amino acids and proteins, which are excellent sources for all animals. Most N is available as the stable isotope ¹⁴ N, but a second form, ¹⁵ N, is present in very low concentrations. ¹⁵ N can be detected in extracts of plants by gas chromatography followed by mass spectrometry (GC/MS). In this protocol, the methods are described for tracing the pathway by which plants are able to take up ¹⁵ N-labeled nitrate and ammonium and convert them into amino acids and proteins. A protocol for extracting and quantifying amino acids and ¹⁵ N enrichment in maize (Zea mays L.) leaves labeled with ¹⁵ NH₄⁺ is described. Following amino acid extraction, purification, and separation by GC/MS, a calculation of the ¹⁵ N enrichment of each amino acid is carried out on a relative basis to identify any differences in the dynamics of amino acid accumulation. This will allow a study of the impact of genetic modifications or mutations on key reactions involved in primary nitrogen and carbon metabolism. © 2018 by John Wiley & Sons, Inc.

Cell-Type Specific Metabolic Flux Analysis: A Challenge for Metabolic Phenotyping and a Potential Solution in Plants

Stable isotope labelling experiments are used routinely in metabolic flux analysis (MFA) to determine the metabolic phenotype of cells and tissues. A complication arises in multicellular systems because single cell measurements of transcriptomes, proteomes and metabolomes in multicellular organisms suggest that the metabolic phenotype will differ between cell types. In silico analysis of simulated metabolite isotopomer datasets shows that cellular heterogeneity confounds conventional MFA because labelling data averaged over multiple cell types does not necessarily yield averaged flux values. A potential solution to this problem—the use of cell-type specific reporter proteins as a source of cell-type specific labelling data—is proposed and the practicality of implementing this strategy in the roots of Arabidopsis thaliana seedlings is explored. A protocol for the immunopurification of ectopically expressed green fluorescent protein (GFP) from Arabidopsis thaliana seedlings using a GFP-binding nanobody is developed, and through GC-MS analysis of protein hydrolysates it is established that constitutively expressed GFP reports accurately on the labelling of total protein in root tissues. It is also demonstrated that the constitutive expression of GFP does not perturb metabolism. The principal obstacle to the implementation of the method in tissues with cell-type specific GFP expression is the sensitivity of the GC-MS system.

The best models of metabolism

Biochemical systems are among of the oldest application areas of mathematical modeling. Spanning a time period of over one hundred years, the repertoire of options for structuring a model and for formulating reactions has been constantly growing, and yet, it is still unclear whether or to what degree some models are better than others and how the modeler is to choose among them. In fact, the variety of options has become overwhelming and difficult to maneuver for novices and experts alike. This review outlines the metabolic model design process and discusses the numerous choices for modeling frameworks and mathematical representations. It tries to be inclusive, even though it cannot be complete, and introduces the various modeling options in a manner that is as unbiased as that is feasible. However, the review does end with personal recommendations for the choices of default models. WIREs Syst Biol Med 2017, 9:e1391. doi: 10.1002/wsbm.1391 For further resources related to this article, please visit the WIREs website.

iReMet-flux: constraint-based approach for integrating relative metabolite levels into a stoichiometric metabolic models

MOTIVATION

Understanding the rerouting of metabolic reaction fluxes upon perturbations has the potential to link changes in molecular state of a cellular system to alteration of growth. Yet, differential flux profiling on a genome-scale level remains one of the biggest challenges in systems biology. This is particularly relevant in plants, for which fluxes in autotrophic growth necessitate time-consuming instationary labeling experiments and costly computations, feasible for small-scale networks.

RESULTS

Here we present a computationally and experimentally facile approach, termed iReMet-Flux, which integrates relative metabolomics data in a metabolic model to predict differential fluxes at a genome-scale level. Our approach and its variants complement the flux estimation methods based on radioactive tracer labeling. We employ iReMet-Flux with publically available metabolic profiles to predict reactions and pathways with altered fluxes in photo-autotrophically grown Arabidopsis and four photorespiratory mutants undergoing high-to-low CO2 acclimation. We also provide predictions about reactions and pathways which are most strongly regulated in the investigated experiments. The robustness and variability analyses, tailored to the formulation of iReMet-Flux, demonstrate that the findings provide biologically relevant information that is validated with external measurements of net CO2 exchange and biomass production. Therefore, iReMet-Flux paves the wave for mechanistic dissection of the interplay between pathways of primary and secondary metabolisms at a genome-scale.

AVAILABILITY AND IMPLEMENTATION

The source code is available from the authors upon request.

CONTACT

nikoloski@mpimp-golm.mpg.de

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Fluxomics links cellular functional analyses to whole-plant phenotyping

Fluxes through metabolic pathways reflect the integration of genetic and metabolic regulations. While it is attractive to measure all the mRNAs (transcriptome), all the proteins (proteome), and a large number of the metabolites (metabolome) in a given cellular system, linking and integrating this information remains difficult. Measurement of metabolome-wide fluxes (termed the fluxome) provides an integrated functional output of the cell machinery and a better tool to link functional analyses to plant phenotyping. This review presents and discusses sets of methodologies that have been developed to measure the fluxome. First, the principles of metabolic flux analysis (MFA), its 'short time interval' version Inst-MFA, and of constraints-based methods, such as flux balance analysis and kinetic analysis, are briefly described. The use of these powerful methods for flux characterization at the cellular scale up to the organ (fruits, seeds) and whole-plant level is illustrated. The added value given by fluxomics methods for unravelling how the abiotic environment affects flux, the process, and key metabolic steps are also described. Challenges associated with the development of fluxomics and its integration with 'omics' for thorough plant and organ functional phenotyping are discussed. Taken together, these will ultimately provide crucial clues for identifying appropriate target plant phenotypes for breeding.

Metabolic network modeling with model organisms

Flux balance analysis (FBA) with genome-scale metabolic network models (GSMNM) allows systems level predictions of metabolism in a variety of organisms. Different types of predictions with different accuracy levels can be made depending on the applied experimental constraints ranging from measurement of exchange fluxes to the integration of gene expression data. Metabolic network modeling with model organisms has pioneered method development in this field. In addition, model organism GSMNMs are useful for basic understanding of metabolism, and in the case of animal models, for the study of metabolic human diseases. Here, we discuss GSMNMs of most highly used model organisms with the emphasis on recent reconstructions.

Can stable isotope mass spectrometry replace ‎radiolabelled approaches in metabolic studies?

Metabolic pathways and the key regulatory points thereof can be deduced using isotopically labelled substrates. One prerequisite is the accurate measurement of the labeling pattern of targeted metabolites. The subsequent estimation of metabolic fluxes following incubation in radiolabelled substrates has been extensively used. Radiolabelling is a sensitive approach and allows determination of total label uptake since the total radiolabel content is easy to detect. However, the incubation of cells, tissues or the whole plant in a stable isotope enriched environment and the use of either mass spectrometry or nuclear magnetic resonance techniques to determine label incorporation within specific metabolites offers the possibility to readily obtain metabolic information with higher resolution. It additionally also offers an important complement to other post-genomic strategies such as metabolite profiling providing insights into the regulation of the metabolic network and thus allowing a more thorough description of plant cellular function. Thus, although safety concerns mean that stable isotope feeding is generally preferred, the techniques are in truth highly complementary and application of both approaches in tandem currently probably provides the best route towards a comprehensive understanding of plant cellular metabolism.

Decoding Biosynthetic Pathways in Plants by Pulse-Chase Strategies Using (13)CO₂ as a Universal Tracer †

(13)CO₂ pulse-chase experiments monitored by high-resolution NMR spectroscopy and mass spectrometry can provide (13)C-isotopologue compositions in biosynthetic products. Experiments with a variety of plant species have documented that the isotopologue profiles generated with (13)CO₂ pulse-chase labeling are directly comparable to those that can be generated by the application of [U-(13)C₆]glucose to aseptically growing plants. However, the application of the (13)CO₂ labeling technology is not subject to the experimental limitations that one has to take into account for experiments with [U-(13)C₆]glucose and can be applied to plants growing under physiological conditions, even in the field. In practical terms, the results of biosynthetic studies with (13)CO₂ consist of the detection of pairs, triples and occasionally quadruples of (13)C atoms that have been jointly contributed to the target metabolite, at an abundance that is well above the stochastic occurrence of such multiples. Notably, the connectivities of jointly transferred (13)C multiples can have undergone modification by skeletal rearrangements that can be diagnosed from the isotopologue data. As shown by the examples presented in this review article, the approach turns out to be powerful in decoding the carbon topology of even complex biosynthetic pathways.

Functional centrality as a predictor of shifts in metabolic flux states

Background

The flux phenotype describes the entirety of biochemical conversions in a cell, which renders it a key characteristic of metabolic function. To quantify the functional relevance of individual biochemical reactions, functional centrality has been introduced based on cooperative game theory and structural modeling. It was shown to be capable to determine metabolic control properties utilizing only structural information. Here, we demonstrate the capability of functional centrality to predict changes in the flux phenotype.

Results

We use functional centrality to successfully predict changes of metabolic flux triggered by switches in the environment. The predictions via functional centrality improve upon predictions using control-effective fluxes, another measure aiming at capturing metabolic control using structural information.

Conclusions

The predictions of flux changes via functional centrality corroborate the capability of the measure to gain a mechanistic understanding of metabolic control from the structure of metabolic networks.

Quantifying plant phenotypes with isotopic labeling & metabolic flux analysis

Analyses of metabolic flux using stable isotopes in plants have traditionally been restricted to tissues with presumed homogeneous cell populations and long metabolic steady states such as developing seeds, cell suspensions, or cultured roots and root tips. It is now possible to describe these and other metabolically more dynamic tissues such as leaves in greater detail using novel methods in mass spectrometry, isotope labeling strategies, and transient labeling-based flux analyses. Such studies are necessary for a systems level description of plant function that more closely represents biological reality, and provides insights into the genes that will need to be modified as natural resources become ever more limited and environments change.

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.

Plant phenomics and the need for physiological phenotyping across scales to narrow the genotype-to-phenotype knowledge gap

Plants are affected by complex genome×environment×management interactions which determine phenotypic plasticity as a result of the variability of genetic components. Whereas great advances have been made in the cost-efficient and high-throughput analyses of genetic information and non-invasive phenotyping, the large-scale analyses of the underlying physiological mechanisms lag behind. The external phenotype is determined by the sum of the complex interactions of metabolic pathways and intracellular regulatory networks that is reflected in an internal, physiological, and biochemical phenotype. These various scales of dynamic physiological responses need to be considered, and genotyping and external phenotyping should be linked to the physiology at the cellular and tissue level. A high-dimensional physiological phenotyping across scales is needed that integrates the precise characterization of the internal phenotype into high-throughput phenotyping of whole plants and canopies. By this means, complex traits can be broken down into individual components of physiological traits. Since the higher resolution of physiological phenotyping by 'wet chemistry' is inherently limited in throughput, high-throughput non-invasive phenotyping needs to be validated and verified across scales to be used as proxy for the underlying processes. Armed with this interdisciplinary and multidimensional phenomics approach, plant physiology, non-invasive phenotyping, and functional genomics will complement each other, ultimately enabling the in silico assessment of responses under defined environments with advanced crop models. This will allow generation of robust physiological predictors also for complex traits to bridge the knowledge gap between genotype and phenotype for applications in breeding, precision farming, and basic research.

Parallel labeling experiments for pathway elucidation and (13)C metabolic flux analysis

Metabolic pathway models provide the foundation for quantitative studies of cellular physiology through the measurement of intracellular metabolic fluxes. For model organisms metabolic models are well established, with many manually curated genome-scale model reconstructions, gene knockout studies and stable-isotope tracing studies. However, for non-model organisms a similar level of knowledge is often lacking. Compartmentation of cellular metabolism in eukaryotic systems also presents significant challenges for quantitative (13)C-metabolic flux analysis ((13)C-MFA). Recently, innovative (13)C-MFA approaches have been developed based on parallel labeling experiments, the use of multiple isotopic tracers and integrated data analysis, that allow more rigorous validation of pathway models and improved quantification of metabolic fluxes. Applications of these approaches open new research directions in metabolic engineering, biotechnology and medicine.

Pool size measurements facilitate the determination of fluxes at branching points in non-stationary metabolic flux analysis: the case of Arabidopsis thaliana

Pool size measurements are important for the estimation of absolute intracellular fluxes in particular scenarios based on data from heavy carbon isotope experiments. Recently, steady-state fluxes estimates were obtained for central carbon metabolism in an intact illuminated rosette of Arabidopsis thaliana grown photoautotrophically (Szecowka et al., 2013; Heise et al., 2014). Fluxes were estimated therein by integrating mass-spectrometric data of the dynamics of the unlabeled metabolic fraction, data on metabolic pool sizes, partitioning of metabolic pools between cellular compartments and estimates of photosynthetically inactive pools, with a simplified model of plant central carbon metabolism. However, the fluxes were determined by treating the pool sizes as fixed parameters. Here we investigated whether and, if so, to what extent the treatment of pool sizes as parameters to be optimized in three scenarios may affect the flux estimates. The results are discussed in terms of benchmark values for canonical pathways and reactions, including starch and sucrose synthesis as well as the ribulose-1,5-bisphosphate carboxylation and oxygenation reactions. In addition, we discuss pathways emerging from a divergent branch point for which pool sizes are required for flux estimation, irrespective of the computational approach used for the simulation of the observable labeling pattern. Therefore, our findings indicate the necessity for development of techniques for accurate pool size measurements to improve the quality of flux estimates from non-stationary flux estimates in intact plant cells in the absence of alternative flux measurements.