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

13CO2-labelling and Sampling in Algae for Flux Analysis of Photosynthetic and Central Carbon Metabolism

The flux in photosynthesis can be studied by performing 13CO2 pulse labelling and analysing the temporal labelling kinetics of metabolic intermediates using gas or liquid chromatography linked to mass spectrometry. Metabolic flux analysis (MFA) is the primary approach for analysing metabolic network function and quantifying intracellular metabolic fluxes. Different MFA approaches differ based on the metabolic state (steady vs. non-steady state) and the use of stable isotope tracers. The main methodology used to investigate metabolic systems is metabolite steady state associated with stable isotope labelling experiments. Specifically, in biological systems like photoautotrophic organisms, isotopic non-stationary 113C metabolic flux analysis at metabolic steady state with transient isotopic labelling (13C-INST-MFA) is required. The common requirement for metabolic steady state, alongside its very short half-timed reactions, complicates robust MFA of photosynthetic metabolism. While custom gas chambers design has addressed these challenges in various model plants, no similar tools were developed for liquid photosynthetic cultures (e.g., algae, cyanobacteria), where diffusion and equilibration of inorganic carbon species in the medium entails a new dimension of complexity. Recently, a novel tailor-made microfluidics labelling system has been introduced, supplying short 13CO2 pulses at steady state, and resolving fluxes across most photosynthetic metabolic pathways in algae. The system involves injecting algal cultures and medium containing pre-equilibrated inorganic 13C into a microfluidic mixer, followed by rapid metabolic quenching, enabling precise seconds-level label pulses. This was complemented by a 13CO2-bubbling-based open labelling system (photobioreactor), allowing long pulses (minutes-hours) required for investigating fluxes into central C metabolism and major products. This combined labelling procedure provides a comprehensive fluxome cover for most algal photosynthetic and central C metabolism pathways, thus allowing comparative flux analyses across algae and plants.

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.

Phosphatidylethanolamine homeostasis under conditions of impaired CDP-ethanolamine pathway or phosphatidylserine decarboxylation

Phosphatidylethanolamine is the major inner-membrane lipid in the plasma and mitochondrial membranes. It is synthesized in the endoplasmic reticulum from ethanolamine and diacylglycerol (DAG) by the CDP-ethanolamine pathway and from phosphatidylserine by decarboxylation in the mitochondria. Recently, multiple genetic disorders that impact these pathways have been identified, including hereditary spastic paraplegia 81 and 82, Liberfarb syndrome, and a new type of childhood-onset neurodegeneration-CONATOC. Individuals with these diseases suffer from multisystem disorders mainly affecting neuronal function. This indicates the importance of maintaining proper phospholipid homeostasis when major biosynthetic pathways are impaired. This study summarizes the current knowledge of phosphatidylethanolamine metabolism in order to identify areas of future research that might lead to the development of treatment options.

Pool size measurements improve precision of flux estimates but increase sensitivity to unmodeled reactions outside the core network in isotopically nonstationary metabolic flux analysis (INST-MFA)

Metabolic flux analysis (MFA) involves model-based estimation of metabolic reaction rates (i.e., fluxes) and, in some cases, metabolite content (i.e., pool sizes) from experimental measurements. Applying MFA to biological data helps determine the fate of substrates and the activity of specific pathways within metabolic networks. However, reliably estimating fluxes by using simplified "core" models to predict the dynamics of larger metabolic networks remains a challenge. One point of uncertainty relates to the advantages and potential pitfalls of including pool size measurements as experimental inputs for isotopically nonstationary MFA (INST-MFA). Here, we directly assessed the role of pool sizes using various core models and simulated datasets. To investigate the effects of pool size measurements on INST-MFA, we assessed the accuracy and precision of flux estimates obtained using different subsets of data (e.g., with or without pool size measurements) and simple network models that either matched or differed from the true network. The inclusion of pool size measurements provided incremental improvements to the precision of the flux estimates. However, adding pool size measurements increased the sensitivity of the flux solution to unmodeled reactions outside the core network. These results were confirmed using a large E. coli model that is representative of realistic metabolic networks examined in MFA studies. Our findings indicate that accurate flux estimates can be obtained in the absence of pool size measurements, even when using core models that lack full network coverage. Addition of pool size measurements to INST-MFA datasets may reveal the activity of non-core reactions that influence the labeling dynamics and therefore necessitate network expansion in order to reconcile all available data to the model. Our findings also emphasize the key role that goodness-of-fit testing plays in assessing the quality of model fits obtained with INST-MFA. This article is protected by copyright. All rights reserved.

Stable Isotopes for Tracing Cardiac Metabolism in Diseases

Although metabolic remodeling during cardiovascular diseases has been well-recognized for decades, the recent development of analytical platforms and mathematical tools has driven the emergence of assessing cardiac metabolism using tracers. Metabolism is a critical component of cellular functions and adaptation to stress. The pathogenesis of cardiovascular disease involves metabolic adaptation to maintain cardiac contractile function even in advanced disease stages. Stable-isotope tracer measurements are a powerful tool for measuring flux distributions at the whole organism level and assessing metabolic changes at a systems level in vivo. The goal of this review is to summarize techniques and concepts for in vivo or ex vivo stable isotope labeling in cardiovascular research, to highlight mathematical concepts and their limitations, to describe analytical methods at the tissue and single-cell level, and to discuss opportunities to leverage metabolic models to address important mechanistic questions relevant to all patients with cardiovascular disease.

Integrating molecular markers into metabolic models improves genomic selection for Arabidopsis growth

The current trends of crop yield improvements are not expected to meet the projected rise in demand. Genomic selection uses molecular markers and machine learning to identify superior genotypes with improved traits, such as growth. Plant growth directly depends on rates of metabolic reactions which transform nutrients into the building blocks of biomass. Here, we predict growth of Arabidopsis thaliana accessions by employing genomic prediction of reaction rates estimated from accession-specific metabolic models. We demonstrate that, comparing to classical genomic selection on the available data sets for 67 accessions, our approach improves the prediction accuracy for growth within and across nitrogen environments by 32.6% and 51.4%, respectively, and from optimal nitrogen to low carbon environment by 50.4%. Therefore, integration of molecular markers into metabolic models offers an approach to predict traits directly related to metabolism, and its usefulness in breeding can be examined by gathering matching datasets in crops.

An increase in genomic selection (GS) accuracy can accelerate genetic gain by shortening the breeding cycles. Here, the authors introduce a network-based GS method that uses metabolic models and improves the prediction accuracy of Arabidopsis growth within and across environments.

Sink/Source Balance of Leaves Influences Amino Acid Pools and Their Associated Metabolic Fluxes in Winter Oilseed Rape (Brassica napus L.)

Nitrogen remobilization processes from source to sink tissues in plants are determinant for seed yield and their implementation results in a complete reorganization of the primary metabolism during sink/source transition. Here, we decided to characterize the impact of the sink/source balance on amino acid metabolism in the leaves of winter oilseed rape grown at the vegetative stage. We combined a quantitative metabolomics approach with an instationary ¹⁵N-labeling experiment by using [¹⁵N]_L-glycine as a metabolic probe on leaf ranks with a gradual increase in their source status. We showed that the acquisition of the source status by leaves was specifically accompanied by a decrease in asparagine, glutamine, proline and S-methyl-l-cysteine sulphoxide contents and an increase in valine and threonine contents. Dynamic analysis of ¹⁵N enrichment and concentration of amino acids revealed gradual changes in the dynamics of amino acid metabolism with respect to the sink/source status of leaf ranks. Notably, nitrogen assimilation into valine, threonine and proline were all decreased in source leaves compared to sink leaves. Overall, our results suggested a reduction in de novo amino acid biosynthesis during sink/source transition at the vegetative stage.

Sink/Source Balance of Leaves Influences Amino Acid Pools and Their Associated Metabolic Fluxes in Winter Oilseed Rape (Brassica napus L.)

Nitrogen remobilization processes from source to sink tissues in plants are determinant for seed yield and their implementation results in a complete reorganization of the primary metabolism during sink/source transition. Here, we decided to characterize the impact of the sink/source balance on amino acid metabolism in the leaves of winter oilseed rape grown at the vegetative stage. We combined a quantitative metabolomics approach with an instationary ¹⁵N-labeling experiment by using [¹⁵N]_L-glycine as a metabolic probe on leaf ranks with a gradual increase in their source status. We showed that the acquisition of the source status by leaves was specifically accompanied by a decrease in asparagine, glutamine, proline and S-methyl-l-cysteine sulphoxide contents and an increase in valine and threonine contents. Dynamic analysis of ¹⁵N enrichment and concentration of amino acids revealed gradual changes in the dynamics of amino acid metabolism with respect to the sink/source status of leaf ranks. Notably, nitrogen assimilation into valine, threonine and proline were all decreased in source leaves compared to sink leaves. Overall, our results suggested a reduction in de novo amino acid biosynthesis during sink/source transition at the vegetative stage.

Limitations of Deuterium-Labelled Substrates for Quantifying NADPH Metabolism in Heterotrophic Arabidopsis Cell Cultures

NADPH is the primary source of cellular reductant for biosynthesis, and strategies for increasing productivity via metabolic engineering need to take account of the requirement for reducing power. In plants, while the oxidative pentose phosphate pathway is the most direct route for NADPH production in heterotrophic tissues, there is increasing evidence that other pathways make significant contributions to redox balance. Deuterium-based isotopic labelling strategies have recently been developed to quantify the relative production of NADPH from different pathways in mammalian cells, but the application of these methods to plants has not been critically evaluated. In this study, LC-MS was used to measure deuterium incorporation into metabolites extracted from heterotrophic Arabidopsis cell cultures grown on [1-²H]glucose or D₂O. The results show that a high rate of flavin-enzyme-catalysed water exchange obscures labelling of NADPH from deuterated substrates and that this exchange cannot be accurately accounted for due to exchange between triose- and hexose-phosphates. In addition, the duplication of NADPH generating reactions between subcellular compartments can confound analysis based on whole cell extracts. Understanding how the structure of the metabolic network affects the applicability of deuterium labelling methods is a prerequisite for development of more effective flux determination strategies, ensuring data are both quantitative and representative of endogenous biological processes.

The Design of FluxML: A Universal Modeling Language for 13C Metabolic Flux Analysis

¹³C metabolic flux analysis (MFA) is the method of choice when a detailed inference of intracellular metabolic fluxes in living organisms under metabolic quasi-steady state conditions is desired. Being continuously developed since two decades, the technology made major contributions to the quantitative characterization of organisms in all fields of biotechnology and health-related research. ¹³C MFA, however, stands out from other "-omics sciences," in that it requires not only experimental-analytical data, but also mathematical models and a computational toolset to infer the quantities of interest, i.e., the metabolic fluxes. At present, these models cannot be conveniently exchanged between different labs. Here, we present the implementation-independent model description language FluxML for specifying ¹³C MFA models. The core of FluxML captures the metabolic reaction network together with atom mappings, constraints on the model parameters, and the wealth of data configurations. In particular, we describe the governing design processes that shaped the FluxML language. We demonstrate the utility of FluxML to represent many contemporary experimental-analytical requirements in the field of ¹³C MFA. The major aim of FluxML is to offer a sound, open, and future-proof language to unambiguously express and conserve all the necessary information for model re-use, exchange, and comparison. Along with FluxML, several powerful computational tools are supplied for easy handling, but also to maintain a maximum of flexibility. Altogether, the FluxML collection is an "all-around carefree package" for ¹³C MFA modelers. We believe that FluxML improves scientific productivity as well as transparency and therewith contributes to the efficiency and reproducibility of computational modeling efforts in the field of ¹³C MFA.

Cyanobacterial carboxysome mutant analysis reveals the influence of enzyme compartmentalization on cellular metabolism and metabolic network rigidity

Cyanobacterial carboxysomes encapsulate carbonic anhydrase and ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO). Genetic deletion of the major structural proteins encoded within the ccm operon in Synechococcus sp. PCC 7002 (ΔccmKLMN) disrupts carboxysome formation and significantly affects cellular physiology. Here we employed both metabolite pool size analysis and isotopically nonstationary metabolic flux analysis (INST-MFA) to examine metabolic regulation in cells lacking carboxysomes. Under high CO₂ environments (1%), the ΔccmKLMN mutant could recover growth and had a similar central flux distribution as the control strain, with the exceptions of moderately decreased photosynthesis and elevated biomass protein content and photorespiration activity. Metabolite analyses indicated that the ΔccmKLMN strain had significantly larger pool sizes of pyruvate (>18 folds), UDPG (uridine diphosphate glucose), and aspartate as well as higher levels of secreted organic acids (e.g., malate and succinate). These results suggest that the ΔccmKLMN mutant is able to largely maintain a fluxome similar to the control strain by changing in intracellular metabolite concentrations and metabolite overflows under optimal growth conditions. When CO₂ was insufficient (0.2%), provision of acetate moderately promoted mutant growth. Interestingly, the removal of microcompartments may loosen the flux network and promote RuBisCO side-reactions, facilitating redirection of central metabolites to competing pathways (i.e., pyruvate to heterologous lactate production). This study provides important insights into metabolic regulation via enzyme compartmentation and cyanobacterial compensatory responses.

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.

Isotopically nonstationary metabolic flux analysis (INST-MFA): putting theory into practice

Typically, ¹³C flux analysis relies on assumptions of both metabolic and isotopic steady state. If metabolism is steady but isotope labeling is not allowed to fully equilibrate, isotopically nonstationary metabolic flux analysis (INST-MFA) can be used to estimate fluxes. This requires solution of differential equations that describe the time-dependent labeling of network metabolites, while iteratively adjusting the flux and pool size parameters to match the transient labeling measurements. INST-MFA holds a number of unique advantages over approaches that rely solely upon steady-state isotope enrichments. First, INST-MFA can be applied to estimate fluxes in autotrophic systems, which consume only single-carbon substrates. Second, INST-MFA is ideally suited to systems that label slowly due to the presence of large intermediate pools or pathway bottlenecks. Finally, INST-MFA provides increased measurement sensitivity to estimate reversible exchange fluxes and metabolite pool sizes, which represents a potential framework for integrating metabolomic analysis with ¹³C flux analysis. This review highlights the unique capabilities of INST-MFA, describes newly available software tools that automate INST-MFA calculations, presents several practical examples of recent INST-MFA applications, and discusses the technical challenges that lie ahead.

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.

Elucidating dynamic metabolic physiology through network integration of quantitative time-course metabolomics

The increasing availability of metabolomics data necessitates novel methods for deeper data analysis and interpretation. We present a flux balance analysis method that allows for the computation of dynamic intracellular metabolic changes at the cellular scale through integration of time-course absolute quantitative metabolomics. This approach, termed "unsteady-state flux balance analysis" (uFBA), is applied to four cellular systems: three dynamic and one steady-state as a negative control. uFBA and FBA predictions are contrasted, and uFBA is found to be more accurate in predicting dynamic metabolic flux states for red blood cells, platelets, and Saccharomyces cerevisiae. Notably, only uFBA predicts that stored red blood cells metabolize TCA intermediates to regenerate important cofactors, such as ATP, NADH, and NADPH. These pathway usage predictions were subsequently validated through ¹³C isotopic labeling and metabolic flux analysis in stored red blood cells. Utilizing time-course metabolomics data, uFBA provides an accurate method to predict metabolic physiology at the cellular scale for dynamic systems.

Photorespiratory Bypasses Lead to Increased Growth in Arabidopsis thaliana: Are Predictions Consistent with Experimental Evidence?

Arguably, the biggest challenge of modern plant systems biology lies in predicting the performance of plant species, and crops in particular, upon different intracellular and external perturbations. Recently, an increased growth of Arabidopsis thaliana plants was achieved by introducing two different photorespiratory bypasses via metabolic engineering. Here, we investigate the extent to which these findings match the predictions from constraint-based modeling. To determine the effect of the employed metabolic network model on the predictions, we perform a comparative analysis involving three state-of-the-art metabolic reconstructions of A. thaliana. In addition, we investigate three scenarios with respect to experimental findings on the ratios of the carboxylation and oxygenation reactions of Ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO). We demonstrate that the condition-dependent growth phenotypes of one of the engineered bypasses can be qualitatively reproduced by each reconstruction, particularly upon considering the additional constraints with respect to the ratio of fluxes for the RuBisCO reactions. Moreover, our results lend support for the hypothesis of a reduced photorespiration in the engineered plants, and indicate that specific changes in CO₂ exchange as well as in the proxies for co-factor turnover are associated with the predicted growth increase in the engineered plants. We discuss our findings with respect to the structure of the used models, the modeling approaches taken, and the available experimental evidence. Our study sets the ground for investigating other strategies for increase of plant biomass by insertion of synthetic reactions.

Omics data input for metabolic modeling

Recent advancements in high-throughput large-scale analytical methods to sequence genomes of organisms, and to quantify gene expression, proteins, lipids and metabolites have changed the paradigm of metabolic modeling. The cost of data generation and analysis has decreased significantly, which has allowed exponential increase in the amount of omics data being generated for an organism in a very short time. Compared to progress made in microbial metabolic modeling, plant metabolic modeling still remains limited due to its complex genomes and compartmentalization of metabolic reactions. Herein, we review and discuss different omics-datasets with potential application in the functional genomics. In particular, this review focuses on the application of omics-datasets towards construction and reconstruction of plant metabolic models.

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.

Inference and Prediction of Metabolic Network Fluxes

In this Update, we cover the basic principles of the estimation and prediction of the rates of the many interconnected biochemical reactions that constitute plant metabolic networks. This includes metabolic flux analysis approaches that utilize the rates or patterns of redistribution of stable isotopes of carbon and other atoms to estimate fluxes, as well as constraints-based optimization approaches such as flux balance analysis. Some of the major insights that have been gained from analysis of fluxes in plants are discussed, including the functioning of metabolic pathways in a network context, the robustness of the metabolic phenotype, the importance of cell maintenance costs, and the mechanisms that enable energy and redox balancing at steady state. We also discuss methodologies to exploit 'omic data sets for the construction of tissue-specific metabolic network models and to constrain the range of permissible fluxes in such models. Finally, we consider the future directions and challenges faced by the field of metabolic network flux phenotyping.