A universal framework for 13C metabolic flux analysis

Rethinking 13C-metabolic flux analysis - The Bayesian way of flux inference

Metabolic reaction rates (fluxes) play a crucial role in comprehending cellular phenotypes and are essential in areas such as metabolic engineering, biotechnology, and biomedical research. The state-of-the-art technique for estimating fluxes is metabolic flux analysis using isotopic labelling (13C-MFA), which uses a dataset-model combination to determine the fluxes. Bayesian statistical methods are gaining popularity in the field of life sciences, but the use of 13C-MFA is still dominated by conventional best-fit approaches. The slow take-up of Bayesian approaches is, at least partly, due to the unfamiliarity of Bayesian methods to metabolic engineering researchers. To address this unfamiliarity, we here outline similarities and differences between the two approaches and highlight particular advantages of the Bayesian way of flux analysis. With a real-life example, re-analysing a moderately informative labelling dataset of E. coli, we identify situations in which Bayesian methods are advantageous and more informative, pointing to potential pitfalls of current 13C-MFA evaluation approaches. We propose the use of Bayesian model averaging (BMA) for flux inference as a means of overcoming the problem of model uncertainty through its tendency to assign low probabilities to both, models that are unsupported by data, and models that are overly complex. In this capacity, BMA resembles a tempered Ockham's razor. With the tempered razor as a guide, BMA-based 13C-MFA alleviates the problem of model selection uncertainty and is thereby capable of becoming a game changer for metabolic engineering by uncovering new insights and inspiring novel approaches.

Metabolic regulatory network kinetic modeling with multiple isotopic tracers for iPSCs

The rapidly expanding market for regenerative medicines and cell therapies highlights the need to advance the understanding of cellular metabolisms and improve the prediction of cultivation production process for human induced pluripotent stem cells (iPSCs). In this paper, a metabolic kinetic model was developed to characterize the underlying mechanisms of iPSC culture process, which can predict cell response to environmental perturbation and support process control. This model focuses on the central carbon metabolic network, including glycolysis, pentose phosphate pathway, tricarboxylic acid cycle, and amino acid metabolism, which plays a crucial role to support iPSC proliferation. Heterogeneous measures of extracellular metabolites and multiple isotopic tracers collected under multiple conditions were used to learn metabolic regulatory mechanisms. Systematic cross-validation confirmed the model's performance in terms of providing reliable predictions on cellular metabolism and culture process dynamics under various culture conditions. Thus, the developed mechanistic kinetic model can support process control strategies to strategically select optimal cell culture conditions at different times, ensure cell product functionality, and facilitate large-scale manufacturing of regenerative medicines and cell therapies.

The Carbon Flow Shifts from Primary to Secondary Metabolism during Xylem Vessel Cell Differentiation in Arabidopsis thaliana

Xylem vessel cell differentiation is characterized by the deposition of a secondary cell wall (SCW) containing cellulose, hemicellulose and lignin. VASCULAR-RELATED NAC-DOMAIN7 (VND7), a plant-specific NAC (NAM, ATAF1/2, and CUC2) transcription factor, is a master regulator of xylem vessel cell differentiation in Arabidopsis (Arabidopsis thaliana). Previous metabolome analysis using the VND7-inducible system in tobacco BY-2 cells successfully revealed significant quantitative changes in primary metabolites during xylem vessel cell differentiation. However, the flow of primary metabolites is not yet well understood. Here, we performed a metabolomic analysis of VND7-inducible Arabidopsis T87 suspension cells. Capillary electrophoresis-time-of-flight mass spectrometry quantified 57 metabolites, and subsequent data analysis highlighted active changes in the levels of UDP-glucose and phenylalanine, which are building blocks of cellulose and lignin, respectively. In a metabolic flow analysis using stable carbon 13 (13C) isotope, the 13C-labeling ratio specifically increased in 3-phosphoglycerate after 12 h of VND7 induction, followed by an increase in shikimate after 24 h of induction, while the inflow of 13C into lactate from pyruvate was significantly inhibited, indicating an active shift of carbon flow from glycolysis to the shikimate pathway during xylem vessel cell differentiation. In support of this notion, most glycolytic genes involved in the downstream of glyceraldehyde 3-phosphate were downregulated following the induction of xylem vessel cell differentiation, whereas genes for the shikimate pathway and phenylalanine biosynthesis were upregulated. These findings provide evidence for the active shift of carbon flow from primary metabolic pathways to the SCW polymer biosynthetic pathway at specific points during xylem vessel cell differentiation.

Razor sharp: The role of Occam's razor in science

Occam's razor-the principle of simplicity-has recently been attacked as a cultural bias without rational foundation. Increasingly, belief in pseudoscience and mysticism is growing. I argue that inclusion of Occam's razor is an essential factor that distinguishes science from superstition and pseudoscience. I also describe how the razor is embedded in Bayesian inference and argue that science is primarily the means to discover the simplest descriptions of our world.

Determination of Metabolic Fluxes by Deep Learning of Isotope Labeling Patterns

Fluxomics offers a direct readout of metabolic state but relies on indirect measurement. Stable isotope tracers imprint flux-dependent isotope labeling patterns on metabolites we measure; however, the relationship between labeling patterns and fluxes remains elusive. Here we innovate a two-stage machine learning framework termed ML-Flux that streamlines metabolic flux quantitation from isotope tracing. We train machine learning models by simulating atom transitions across five universal metabolic models starting from 26 13C-glucose, 2H-glucose, and 13C-glutamine tracers within feasible flux space. ML-Flux employs deep-learning-based imputation to take variable measurements of labeling patterns as input and successive neural networks to convert the ensuing comprehensive labeling information into metabolic fluxes. Using ML-Flux with multi-isotope tracing, we obtain fluxes through central carbon metabolism that are comparable to those from a least-squares method but orders-of-magnitude faster. ML-Flux is deployed as a webtool to expand the accessibility of metabolic flux quantitation and afford actionable information on metabolism.

Model-driven insights into the effects of temperature on metabolism

Temperature affects cellular processes at different spatiotemporal scales, and identifying the genetic and molecular mechanisms underlying temperature responses paves the way to develop approaches for mitigating the effects of future climate scenarios. A systems view of the effects of temperature on cellular physiology can be obtained by focusing on metabolism since: (i) its functions depend on transcription and translation and (ii) its outcomes support organisms' development, growth, and reproduction. Here we provide a systematic review of modelling efforts directed at investigating temperature effects on properties of single biochemical reactions, system-level traits, metabolic subsystems, and whole-cell metabolism across different prokaryotes and eukaryotes. We compare and contrast computational approaches and theories that facilitate modelling of temperature effects on key properties of enzymes and their consideration in constraint-based as well as kinetic models of metabolism. In addition, we provide a summary of insights from computational approaches, facilitating integration of omics data from temperature-modulated experiments with models of metabolic networks, and review the resulting biotechnological applications. Lastly, we provide a perspective on how different types of metabolic modelling can profit from developments in machine learning and models of different cellular layers to improve model-driven insights into the effects of temperature relevant for biotechnological applications.

Practical sampling of constraint-based models: Optimized thinning boosts CHRR performance

Thinning is a sub-sampling technique to reduce the memory footprint of Markov chain Monte Carlo. Despite being commonly used, thinning is rarely considered efficient. For sampling constraint-based models, a highly relevant use-case in systems biology, we here demonstrate that thinning boosts computational and, thereby, sampling efficiencies of the widely used Coordinate Hit-and-Run with Rounding (CHRR) algorithm. By benchmarking CHRR with thinning with simplices and genome-scale metabolic networks of up to thousands of dimensions, we find a substantial increase in computational efficiency compared to unthinned CHRR, in our examples by orders of magnitude, as measured by the effective sample size per time (ESS/t), with performance gains growing with polytope (effective network) dimension. Using a set of benchmark models we derive a ready-to-apply guideline for tuning thinning to efficient and effective use of compute resources without requiring additional coding effort. Our guideline is validated using three (out-of-sample) large-scale networks and we show that it allows sampling convex polytopes uniformly to convergence in a fraction of time, thereby unlocking the rigorous investigation of hitherto intractable models. The derivation of our guideline is explained in detail, allowing future researchers to update it as needed as new model classes and more training data becomes available. CHRR with deliberate utilization of thinning thereby paves the way to keep pace with progressing model sizes derived with the constraint-based reconstruction and analysis (COBRA) tool set. Sampling and evaluation pipelines are available at https://jugit.fz-juelich.de/IBG-1/ModSim/fluxomics/chrrt.

High throughput 13C-metabolic flux analysis of 3-hydroxypropionic acid producing Pichia pastoris reveals limited availability of acetyl-CoA and ATP due to tight control of the glycolytic flux

BACKGROUND

Production of 3-hydroxypropionic acid (3-HP) through the malonyl-CoA pathway has yielded promising results in Pichia pastoris (Komagataella phaffii), demonstrating the potential of this cell factory to produce this platform chemical and other acetyl-CoA-derived products using glycerol as a carbon source. However, further metabolic engineering of the original P. pastoris 3-HP-producing strains resulted in unexpected outcomes, e.g., significantly lower product yield and/or growth rate. To gain an understanding on the metabolic constraints underlying these observations, the fluxome (metabolic flux phenotype) of ten 3-HP-producing P. pastoris strains has been characterized using a high throughput ¹³C-metabolic flux analysis platform. Such platform enabled the operation of an optimised workflow to obtain comprehensive maps of the carbon flux distribution in the central carbon metabolism in a parallel-automated manner, thereby accelerating the time-consuming strain characterization step in the design-build-test-learn cycle for metabolic engineering of P. pastoris.

RESULTS

We generated detailed maps of the carbon fluxes in the central carbon metabolism of the 3-HP producing strain series, revealing the metabolic consequences of different metabolic engineering strategies aimed at improving NADPH regeneration, enhancing conversion of pyruvate into cytosolic acetyl-CoA, or eliminating by-product (arabitol) formation. Results indicate that the expression of the POS5 NADH kinase leads to a reduction in the fluxes of the pentose phosphate pathway reactions, whereas an increase in the pentose phosphate pathway fluxes was observed when the cytosolic acetyl-CoA synthesis pathway was overexpressed. Results also show that the tight control of the glycolytic flux hampers cell growth due to limited acetyl-CoA biosynthesis. When the cytosolic acetyl-CoA synthesis pathway was overexpressed, the cell growth increased, but the product yield decreased due to higher growth-associated ATP costs. Finally, the six most relevant strains were also cultured at pH 3.5 to assess the effect of a lower pH on their fluxome. Notably, similar metabolic fluxes were observed at pH 3.5 compared to the reference condition at pH 5.

CONCLUSIONS

This study shows that existing fluoxomics workflows for high-throughput analyses of metabolic phenotypes can be adapted to investigate P. pastoris, providing valuable information on the impact of genetic manipulations on the metabolic phenotype of this yeast. Specifically, our results highlight the metabolic robustness of P. pastoris's central carbon metabolism when genetic modifications are made to increase the availability of NADPH and cytosolic acetyl-CoA. Such knowledge can guide further metabolic engineering of these strains. Moreover, insights into the metabolic adaptation of P. pastoris to an acidic pH have also been obtained, showing the capability of the fluoxomics workflow to assess the metabolic impact of environmental changes.

Progress in understanding and improving oil content and quality in seeds

The world's population is projected to increase by two billion by 2050, resulting in food and energy insecurity. Oilseed crops have been identified as key to address these challenges: they produce and store lipids in the seeds as triacylglycerols that can serve as a source of food/feed, renewable fuels, and other industrially-relevant chemicals. Therefore, improving seed oil content and composition has generated immense interest. Research efforts aiming to unravel the regulatory pathways involved in fatty acid synthesis and to identify targets for metabolic engineering have made tremendous progress. This review provides a summary of the current knowledge of oil metabolism and discusses how photochemical activity and unconventional pathways can contribute to high carbon conversion efficiency in seeds. It also highlights the importance of ¹³C-metabolic flux analysis as a tool to gain insights on the pathways that regulate oil biosynthesis in seeds. Finally, a list of key genes and regulators that have been recently targeted to enhance seed oil production are reviewed and additional possible targets in the metabolic pathways are proposed to achieve desirable oil content and quality.

Isolated murine skeletal muscles utilize pyruvate over glucose for oxidation

INTRODUCTION

Fuel sources for skeletal muscle tissue include carbohydrates and fatty acids, and utilization depends upon fiber type, workload, and substrate availability. The use of isotopically labeled substrate tracers combined with nuclear magnetic resonance (NMR) enables a deeper examination of not only utilization of substrates by a given tissue, but also their contribution to tricarboxylic acid (TCA) cycle intermediates.

OBJECTIVES

The goal of this study was to determine the differential utilization of substrates in isolated murine skeletal muscle, and to evaluate how isopotomer anlaysis provided insight into skeletal muscle metabolism.

METHODS

Isolated C57BL/6 mouse hind limb muscles were incubated in oxygenated solution containing uniformly labeled 13C6 glucose, 13C3 pyruvate, or 13C2 acetate at room temperature. Isotopomer analysis of 13C labeled glutamate was performed on pooled extracts of isolated soleus and extensor digitorum longus (EDL) muscles.

RESULTS

Pyruvate and acetate were more avidly consumed than glucose with resultant increases in glutamate labeling in both muscle groups. Glucose incubation resulted in glutamate labeling, but with high anaplerotic flux in contrast to the labeling by pyruvate. Muscle fiber type distinctions were evident by differences in lactate enrichment and extent of substrate oxidation.

CONCLUSION

Isotope tracing experiments in isolated muscles reveal that pyruvate and acetate are avidly oxidized by isolated soleus and EDL muscles, whereas glucose labeling of glutamate is accompanied by high anaplerotic flux. We believe our results may set the stage for future examination of metabolic signatures of skeletal muscles from pre-clinical models of aging, type-2 diabetes and neuromuscular disease.

Isotope-Assisted Metabolic Flux Analysis: A Powerful Technique to Gain New Insights into the Human Metabolome in Health and Disease

Cell metabolism represents the coordinated changes in genes, proteins, and metabolites that occur in health and disease. The metabolic fluxome, which includes both intracellular and extracellular metabolic reaction rates (fluxes), therefore provides a powerful, integrated description of cellular phenotype. However, intracellular fluxes cannot be directly measured. Instead, flux quantification requires sophisticated mathematical and computational analysis of data from isotope labeling experiments. In this review, we describe isotope-assisted metabolic flux analysis (iMFA), a rigorous computational approach to fluxome quantification that integrates metabolic network models and experimental data to generate quantitative metabolic flux maps. We highlight practical considerations for implementing iMFA in mammalian models, as well as iMFA applications in in vitro and in vivo studies of physiology and disease. Finally, we identify promising new frontiers in iMFA which may enable us to fully unlock the potential of iMFA in biomedical research.

Metabolic flux analysis: a comprehensive review on sample preparation, analytical techniques, data analysis, computational modelling, and main application areas

Metabolic flux analysis (MFA) quantitatively describes cellular fluxes to understand metabolic phenotypes and functional behaviour after environmental and/or genetic perturbations. In the last decade, the application of stable isotopes became extremely important to determine and integrate in vivo measurements of metabolic reactions in systems biology. 13C-MFA is one of the most informative methods used to study central metabolism of biological systems. This review aims to outline the current experimental procedure adopted in 13C-MFA, starting from the preparation of cell cultures and labelled tracers to the quenching and extraction of metabolites and their subsequent analysis performed with very powerful software. Here, the limitations and advantages of nuclear magnetic resonance spectroscopy and mass spectrometry techniques used in carbon labelled experiments are elucidated by reviewing the most recent published papers. Furthermore, we summarise the most successful approaches used for computational modelling in flux analysis and the main application areas with a particular focus in metabolic engineering.

An optimization method for untargeted MS-based isotopic tracing investigations of metabolism

INTRODUCTION

Stable isotope tracer studies are increasingly applied to explore metabolism from the detailed analysis of tracer incorporation into metabolites. Untargeted LC/MS approaches have recently emerged and provide potent methods for expanding the dimension and complexity of the metabolic networks that can be investigated. A number of software tools have been developed to process the highly complex MS data collected in such studies; however, a method to optimize the extraction of valuable isotopic data is lacking.

OBJECTIVES

To develop and validate a method to optimize automated data processing for untargeted MS-based isotopic tracing investigations of metabolism.

METHODS

The method is based on the application of a suitable reference material to rationally perform parameter optimization throughout the complete data processing workflow. It was applied in the context of ¹³C-labelling experiments and with two different software, namely geoRge and X13CMS. It was illustrated with the study of a E. coli mutant impaired for central metabolism.

RESULTS

The optimization methodology provided significant gain in the number and quality of extracted isotopic data, independently of the software considered. Pascal triangle samples are well suited for such purpose since they allow both the identification of analytical issues and optimization of data processing at the same time.

CONCLUSION

The proposed method maximizes the biological value of untargeted MS-based isotopic tracing investigations by revealing the full metabolic information that is encoded in the labelling patterns of metabolites.

Ex Vivo and In Vitro Analysis Identify a Detrimental Impact of Neutrophil Extracellular Traps on Eye Structures in Equine Recurrent Uveitis

Equine recurrent uveitis (ERU) is a common ocular disease of horses and described as a model for human autoimmune uveitis. This immune-mediated, inflammatory condition progressively destroys the eye, ultimately leading to blindness. Genetic and autoimmune factors, next to infections with Leptospira, are discussed as key factors in the pathogenesis. Furthermore, a release of neutrophil extracellular traps (NETs) by activated neutrophils is involved. NETs are composed of decondensed chromatin and proteins that can immobilize invading pathogens. However, if NETs accumulate, they can contribute to detrimental autoimmune processes. Thus, we aimed to investigate the impact of NETs in ERU patients. Therefore, we quantified several NET-markers (cell-free DNA, nucleosomes, citrullinated histone H3, histone-myeloperoxidase complexes, interleukin-17, equine cathelicidin 1 and DNase I activity) and NET-autoantibodies in sera and vitreous body fluids (VBF) of ERU-diseased horses and correlated the data with the disease status (signalment, ERU scores and Leptospira infection status). NET markers were detected to varying degrees in VBF of diseased horses, and partially correlated to disease severity and the presence of Leptospira spp. Cell-free DNA and nucleosomes as NET markers correlate with ERU severity in total and VBF scores, despite the presence of active DNases. Additionally, a significant correlation between fundus affection in the eye and NET autoantibodies was detectable. Therefore, we further investigated the influence of VBF samples from equine patients and isolated NETs on the blood-retina barrier in a cell culture model. VBF of diseased horses significantly induced cytotoxicity in retinal pigment epithelial cells. Moreover, partially digested NETs also resulted in cytotoxic effects. In the presence of lipopolysaccharide (LPS), the main component of the leptospiral surface, both undigested and completely digested NETs were cytotoxic. Correlations between the ERU-scores and Leptospira were also calculated. Detection of leptospiral DNA, and antibody titers of the serovar Grippotyphosa correlated with disease severity. In addition, a correlation between Leptospira and several NET markers was observed in VBF. Altogether, our findings suggest a positive correlation between NET markers with disease severity and involvement of Leptospira in the VBF of ERU-diseased horses, as well as a cytotoxic effect of NETs in eyes.

Fluxomics - New Metabolomics Approaches to Monitor Metabolic Pathways

Fluxomics is an innovative -omics research field that measures the rates of all intracellular fluxes in the central metabolism of biological systems. Fluxomics gathers data from multiple different -omics fields, portraying the whole picture of molecular interactions. Recently, fluxomics has become one of the most relevant approaches to investigate metabolic phenotypes. Metabolic flux using ¹³C-labeled molecules is increasingly used to monitor metabolic pathways, to probe the corresponding gene-RNA and protein-metabolite interaction networks in actual time. Thus, fluxomics reveals the functioning of multi-molecular metabolic pathways and is increasingly applied in biotechnology and pharmacology. Here, we describe the main fluxomics approaches and experimental platforms. Moreover, we summarize recent fluxomic results in different biological systems.

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.

Dissecting Host-Pathogen Interactions in TB Using Systems-Based Omic Approaches

Tuberculosis (TB) is a devastating infectious disease that kills over a million people every year. There is an increasing burden of multi drug resistance (MDR) and extensively drug resistance (XDR) TB. New and improved therapies are urgently needed to overcome the limitations of current treatment. The causative agent, Mycobacterium tuberculosis (Mtb) is one of the most successful pathogens that can manipulate host cell environment for adaptation, evading immune defences, virulence, and pathogenesis of TB infection. Host-pathogen interaction is important to establish infection and it involves a complex set of processes. Metabolic cross talk between the host and pathogen is a facet of TB infection and has been an important topic of research where there is growing interest in developing therapies and drugs that target these interactions and metabolism of the pathogen in the host. Mtb scavenges multiple nutrient sources from the host and has adapted its metabolism to survive in the intracellular niche. Advancements in systems-based omic technologies have been successful to unravel host-pathogen interactions in TB. In this review we discuss the application and usefulness of omics in TB research that provides promising interventions for developing anti-TB therapies.

Recent advances in metabolic engineering-integration of in silico design and experimental analysis of metabolic pathways

Microorganisms are widely used to produce valuable compounds. Because thousands of metabolic reactions occur simultaneously and many metabolic reactions are related to target production and cell growth, the development of a rational design method for metabolic pathway modification to optimize target production is needed. In this paper, recent advances in metabolic engineering are reviewed, specifically considering computational pathway modification design and experimental evaluation of metabolic fluxes by ¹³C-metabolic flux analysis. Computational tools for seeking effective gene deletion targets and recruiting heterologous genes are described in flux balance analysis approaches. A kinetic model and adaptive laboratory evolution were applied to identify and eliminate the rate-limiting step in metabolic pathways. Data science-based approaches for process monitoring and control are described to maximize the performance of engineered cells in bioreactors.

Robustifying Experimental Tracer Design for13C-Metabolic Flux Analysis

¹³C metabolic flux analysis (MFA) has become an indispensable tool to measure metabolic reaction rates (fluxes) in living organisms, having an increasingly diverse range of applications. Here, the choice of the¹³C labeled tracer composition makes the difference between an information-rich experiment and an experiment with only limited insights. To improve the chances for an informative labeling experiment, optimal experimental design approaches have been devised for¹³C-MFA, all relying on some a priori knowledge about the actual fluxes. If such prior knowledge is unavailable, e.g., for research organisms and producer strains, existing methods are left with a chicken-and-egg problem. In this work, we present a general computational method, termed robustified experimental design (R-ED), to guide the decision making about suitable tracer choices when prior knowledge about the fluxes is lacking. Instead of focusing on one mixture, optimal for specific flux values, we pursue a sampling based approach and introduce a new design criterion, which characterizes the extent to which mixtures are informative in view of all possible flux values. The R-ED workflow enables the exploration of suitable tracer mixtures and provides full flexibility to trade off information and cost metrics. The potential of the R-ED workflow is showcased by applying the approach to the industrially relevant antibiotic producer Streptomyces clavuligerus, where we suggest informative, yet economic labeling strategies.

An energetic profile of Corynebacterium glutamicum underpinned by measured biomass yield on ATP

The supply and usage of energetic cofactors in metabolism is a central concern for systems metabolic engineering, particularly in case of energy intensive products. One of the most important parameters for systems wide balancing of energetic cofactors is the ATP requirement for biomass formation Y_ATP/Biomass. Despite its fundamental importance, Y_ATP/Biomass values for non-fermentative organisms are still rough estimates deduced from theoretical considerations. For the first time, we present an approach for the experimental determination of Y_ATP/Biomass using comparative ¹³C metabolic flux analysis (¹³C MFA) of a wild type strain and an ATP synthase knockout mutant. We show that the energetic profile of a cell can then be deduced from a genome wide stoichiometric model and experimental maintenance data. Particularly, the contributions of substrate level phosphorylation (SLP) and electron transport phosphorylation (ETP) to ATP generation become available which enables the overall energetic efficiency of a cell to be characterized. As a model organism, the industrial platform organism Corynebacterium glutamicum is used. C. glutamicum uses a respiratory type of energy metabolism, implying that ATP can be synthesized either by SLP or by ETP with the membrane-bound F₁F_O-ATP synthase using the proton motive force (pmf) as driving force. The presence of two terminal oxidases, which differ in their proton translocation efficiency by a factor of three, further complicates energy balancing for this organism. By integration of experimental data and network models, we show that in the wild type SLP and ETP contribute equally to ATP generation. Thus, the role of ETP in respiring bacteria may have been overrated in the past. Remarkably, in the genome wide setting 65% of the pmf is actually not used for ATP synthesis. However, it turns out that, compared to other organisms C. glutamicum still uses its energy budget rather efficiently.

Metabolic engineering of Zymomonas moblis for ethylene production from straw hydrolysate

Biological ethylene production is a promising sustainable alternative approach for fossil-based ethylene production. The high glucose utilization of Z. mobilis makes it as a promising bioethylene producer. In this study, Zymomonas mobilis has been engineered to produce ethylene through the introduction of the synthetic ethylene-forming enzyme (EFE). We also investigated the effect of systematically knocking out the competitive metabolic pathway of pyruvate in an effort to improve the availability of pyruvate for ethylene production in Z. mobilis expressing EFE. Guided by these results, we tested a number of conjectures that could improve the α-ketoglutarate supply. Optimization of these pathways and different substrate supplies resulted in a greater production of ethylene (from 1.36 to 12.83 nmol/OD₆₀₀/mL), which may guide future engineering work on ethylene production using other organisms. Meanwhile, we achieved an ethylene production of 5.8 nmol/OD₆₀₀/mL in the ZM532-efe strain using enzymatic straw hydrolysate of corn straw as the sole carbon source. As a preferred host in biorefinery technologies using lignocellulosic biomass as feedstock, heterologous expression of EFE in Z. mobilis converts the non-ethylene producing strain into an ethylene-producing one using a metabolic engineering approach, which is of great significance for the utilization of cellulosic biomass in the future. KEY POINTS: • Heterologous expression of EFE in Z. mobilis successfully converted the non-ethylene producing strain into an ethylene producer (1.36 nmol/OD₆₀₀/mL). Targeted modifications of the central carbon metabolism can effectively improve ethylene production (peak production: 8.3 nmol/OD₆₀₀/mL). • The addition of nutrients to the medium can further increase the production of ethylene (peak production: 12.8 nmol/OD₆₀₀/mL). • The ZM532-efe strain achieved an ethylene production of 5.8 nmol/OD₆₀₀/mL when enzymatic hydrolysate of corn straw was used as the sole carbon source.

Exploring the Use of Gas Chromatography Coupled to Chemical Ionization Mass Spectrometry (GC-CI-MS) for Stable Isotope Labeling in Metabolomics

Isotopic-labeling experiments have been valuable to monitor the flux of metabolic reactions in biological systems, which is crucial to understand homeostatic alterations with disease. Experimental determination of metabolic fluxes can be inferred from a characteristic rearrangement of stable isotope tracers (e.g., 13C or 15N) that can be detected by mass spectrometry (MS). Metabolites measured are generally members of well-known metabolic pathways, and most of them can be detected using both gas chromatography (GC)-MS and liquid chromatography (LC)-MS. In here, we show that GC methods coupled to chemical ionization (CI) MS have a clear advantage over alternative methodologies due to GC's superior chromatography separation efficiency and the fact that CI is a soft ionization technique that yields identifiable protonated molecular ion peaks. We tested diverse GC-CI-MS setups, including methane and isobutane reagent gases, triple quadrupole (QqQ) MS in SIM mode, or selected ion clusters using optimized narrow windows (∼10 Da) in scan mode, and standard full scan methods using high resolution GC-(q)TOF and GC-Orbitrap systems. Isobutane as a reagent gas in combination with both low-resolution (LR) and high-resolution (HR) MS showed the best performance, enabling precise detection of isotopologues in most metabolic intermediates of central carbon metabolism. Finally, with the aim of overcoming manual operations, we developed an R-based tool called isoSCAN that automatically quantifies all isotopologues of intermediate metabolites of glycolysis, TCA cycle, amino acids, pentose phosphate pathway, and urea cycle, from LRMS and HRMS data.

Recent Advances in Carbon and Nitrogen Metabolism in C3 Plants

C and N are the most important essential elements constituting organic compounds in plants. The shoots and roots depend on each other by exchanging C and N through the xylem and phloem transport systems. Complex mechanisms regulate C and N metabolism to optimize plant growth, agricultural crop production, and maintenance of the agroecosystem. In this paper, we cover the recent advances in understanding C and N metabolism, regulation, and transport in plants, as well as their underlying molecular mechanisms. Special emphasis is given to the mechanisms of starch metabolism in plastids and the changes in responses to environmental stress that were previously overlooked, since these changes provide an essential store of C that fuels plant metabolism and growth. We present general insights into the system biology approaches that have expanded our understanding of core biological questions related to C and N metabolism. Finally, this review synthesizes recent advances in our understanding of the trade-off concept that links C and N status to the plant's response to microorganisms.

Validity of natural isotope abundance correction for metabolic flux analysis

A pervasive issue in stable isotope tracing and metabolic flux analysis is the presence of naturally occurring isotopes such as ¹³C. For mass isotopomer distributions (MIDs) measured by mass spectrometry, it is common practice to correct for natural occurrence of isotopes within metabolites of interest using a linear transform based on binomial distributions. The resulting corrected MIDs are often used to fit metabolic network models and infer metabolic fluxes, which implicitly assumes that corrected MIDs will yield the same flux solution as the actual observed MIDs. Although this assumption can be empirically verified in special cases by simulation studies, there seems to be no published proof of this important property for the general case. In this paper, we prove that this property holds for the case of noise-free MID data obtained at steady state. On the other hand, for noisy MID data, the flux solution will generally differ between the two representations. These results provide a theoretical foundation for the common practice of MID correction in metabolic flux analysis.

Towards model-driven characterization and manipulation of plant lipid metabolism

Plant lipids have versatile applications and provide essential fatty acids in human diet. Therefore, there has been a growing interest to better characterize the genetic basis, regulatory networks, and metabolic pathways that shape lipid quantity and composition. Addressing these issues is challenging due to context-specificity of lipid metabolism integrating environmental, developmental, and tissue-specific cues. Here we systematically review the known metabolic pathways and regulatory interactions that modulate the levels of storage lipids in oilseeds. We argue that the current understanding of lipid metabolism provides the basis for its study in the context of genome-wide plant metabolic networks with the help of approaches from constraint-based modeling and metabolic flux analysis. The focus is on providing a comprehensive summary of the state-of-the-art of modeling plant lipid metabolic pathways, which we then contrast with the existing modeling efforts in yeast and microalgae. We then point out the gaps in knowledge of lipid metabolism, and enumerate the recent advances of using genome-wide association and quantitative trait loci mapping studies to unravel the genetic regulations of lipid metabolism. Finally, we offer a perspective on how advances in the constraint-based modeling framework can propel further characterization of plant lipid metabolism and its rational manipulation.

Dynamic 13C Flux Analysis Captures the Reorganization of Adipocyte Glucose Metabolism in Response to Insulin

Cellular metabolism is dynamic, but quantifying non-steady metabolic fluxes by stable isotope tracers presents unique computational challenges. Here, we developed an efficient 13C-tracer dynamic metabolic flux analysis (13C-DMFA) framework for modeling central carbon fluxes that vary over time. We used B-splines to generalize the flux parameterization system and to improve the stability of the optimization algorithm. As proof of concept, we investigated how 3T3-L1 cultured adipocytes acutely metabolize glucose in response to insulin. Insulin rapidly stimulates glucose uptake, but intracellular pathways responded with differing speeds and magnitudes. Fluxes in lower glycolysis increased faster than those in upper glycolysis. Glycolysis fluxes rose disproportionally larger and faster than the tricarboxylic acid cycle, with lactate a primary glucose end product. The uncovered array of flux dynamics suggests that glucose catabolism is additionally regulated beyond uptake to help shunt glucose into appropriate pathways. This work demonstrates the value of using dynamic intracellular fluxes to understand metabolic function and pathway regulation.

High-resolution 13C metabolic flux analysis

Precise quantification of metabolic pathway fluxes in biological systems is of major importance in guiding efforts in metabolic engineering, biotechnology, microbiology, human health, and cell culture. ¹³C metabolic flux analysis (¹³C-MFA) is the predominant technique used for determining intracellular fluxes. Here, we present a protocol for ¹³C-MFA that incorporates recent advances in parallel labeling experiments, isotopic labeling measurements, and statistical analysis, as well as best practices developed through decades of experience. Experimental design to ensure that fluxes are estimated with the highest precision is an integral part of the protocol. The protocol is based on growing microbes in two (or more) parallel cultures with ¹³C-labeled glucose tracers, followed by gas chromatography-mass spectrometry (GC-MS) measurements of isotopic labeling of protein-bound amino acids, glycogen-bound glucose, and RNA-bound ribose. Fluxes are then estimated using software for ¹³C-MFA, such as Metran, followed by comprehensive statistical analysis to determine the goodness of fit and calculate confidence intervals of fluxes. The presented protocol can be completed in 4 d and quantifies metabolic fluxes with a standard deviation of ≤2%, a substantial improvement over previous implementations. The presented protocol is exemplified using an Escherichia coli ΔtpiA case study with full supporting data, providing a hands-on opportunity to step through a complex troubleshooting scenario. Although applications to prokaryotic microbial systems are emphasized, this protocol can be easily adjusted for application to eukaryotic organisms.

Position-Specific Metabolic Probing and Metagenomics of Microbial Communities Reveal Conserved Central Carbon Metabolic Network Activities at High Temperatures

Temperature is a primary driver of microbial community composition and taxonomic diversity; however, it is unclear to what extent temperature affects characteristics of central carbon metabolic pathways (CCMPs) at the community level. In this study, 16S rRNA gene amplicon and metagenome sequencing were combined with ¹³C-labeled metabolite probing of the CCMPs to assess community carbon metabolism along a temperature gradient (60–95°C) in Great Boiling Spring, NV. 16S rRNA gene amplicon diversity was inversely proportional to temperature, and Archaea were dominant at higher temperatures. KO richness and diversity were also inversely proportional to temperature, yet CCMP genes were similarly represented across the temperature gradient and many individual metagenome-assembled genomes had complete pathways. In contrast, genes encoding cellulosomes and many genes involved in plant matter degradation and photosynthesis were absent at higher temperatures. In situ¹³C-CO₂ production from labeled isotopomer pairs of glucose, pyruvate, and acetate suggested lower relative oxidative pentose phosphate pathway activity and/or fermentation at 60°C, and a stable or decreased maintenance energy demand at higher temperatures. Catabolism of ¹³C-labeled citrate, succinate, L-alanine, L-serine, and L-cysteine was observed at 85°C, demonstrating broad heterotrophic activity and confirming functioning of the TCA cycle. Together, these results suggest that temperature-driven losses in biodiversity and gene content in geothermal systems may not alter CCMP function or maintenance energy demands at a community level.

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.

Ex vivo and in vivo stable isotope labelling of central carbon metabolism and related pathways with analysis by LC-MS/MS

Targeted tandem mass spectrometry (LC-MS/MS) has been extremely useful for profiling small molecules extracted from biological sources, such as cells, bodily fluids and tissues. Here, we present a protocol for analysing incorporation of the non-radioactive stable isotopes carbon-13 (¹³C) and nitrogen-15 (¹⁵N) into polar metabolites in central carbon metabolism and related pathways. Our platform utilizes selected reaction monitoring (SRM) with polarity switching and amide hydrophilic interaction liquid chromatography (HILIC) to capture transitions for carbon and nitrogen incorporation into selected metabolites using a hybrid triple quadrupole (QQQ) mass spectrometer. This protocol represents an extension of a previously published protocol for targeted metabolomics of unlabeled species and has been used extensively in tracing the metabolism of nutrients such as ¹³C-labeled glucose, ¹³C-glutamine and ¹⁵N-glutamine in a variety of biological settings (e.g., cell culture experiments and in vivo mouse labelling via i.p. injection). SRM signals are integrated to produce an array of peak areas for each labelling form that serve as the output for further analysis. The processed data are then used to obtain the degree and distribution of labelling of the targeted molecules (termed fluxomics). Each method can be customized on the basis of known unlabeled Q1/Q3 SRM transitions and adjusted to account for the corresponding ¹³C or ¹⁵N incorporation. The entire procedure takes ~6-7 h for a single sample from experimental labelling and metabolite extraction to peak integration.

Therapeutic Implication of SOCS1 Modulation in the Treatment of Autoimmunity and Cancer

The suppressor of cytokine signaling (SOCS) family of intracellular proteins has a vital role in the regulation of the immune system and resolution of inflammatory cascades. SOCS1, also called STAT-induced STAT inhibitor (SSI) or JAK-binding protein (JAB), is a member of the SOCS family with actions ranging from immune modulation to cell cycle regulation. Knockout of SOCS1 leads to perinatal lethality in mice and increased vulnerability to cancer, while several SNPs associated with the SOCS1 gene have been implicated in human inflammation-mediated diseases. In this review, we describe the mechanism of action of SOCS1 and its potential therapeutic role in the prevention and treatment of autoimmunity and cancer. We also provide a brief outline of the other JAK inhibitors, both FDA-approved and under investigation.

13C-Fingerprinting and Metabolic Flux Analysis of Bacterial Metabolisms

¹³C-assisted metabolism analysis provides rigorous calculations of the intracellular reaction rates (i.e., fluxes) within the central metabolism of microbial hosts. This mapping of the intracellular network within microbes has proven to be essential for understanding the cell physiology. The approach is also a key to identifying central metabolic nodes, probing the rigidity of a metabolic network, revealing cofactor balances, and delineating hidden pathways. Here we present the methodology of using stable isotopic carbon substrates for both qualitative (¹³C-fingerprinting of functional pathways) and quantitative (Metabolic Flux Analysis) metabolism studies on bacterial species. In this methodology, we include step-by-step instructions to use the open source WUflux software for the steady-state flux calculations based on labeling information of amino acids or free metabolites.

Analytical Considerations of Stable Isotope Labelling in Lipidomics

Over the last two decades, lipids have come to be understood as far more than merely components of cellular membranes and forms of energy storage, and are now also being implicated to play important roles in a variety of diseases, with lipid biomarker research one of the most widespread applications of lipidomic techniques both in research and in clinical settings. Stable isotope labelling has become a staple technique in the analysis of small molecule metabolism and dynamics, as it is the only experimental setup by which biosynthesis, remodelling and degradation of biomolecules can be directly measured. Using state-of-the-art analytical technologies such as chromatography-coupled high resolution tandem mass spectrometry, the stable isotope label can be precisely localized and quantified within the biomolecules. The application of stable isotope labelling to lipidomics is however complicated by the diversity of lipids and the complexity of the necessary data analysis. This article discusses key experimental aspects of stable isotope labelling in the field of mass spectrometry-based lipidomics, summarizes current applications and provides an outlook on future developments and potential.

FluxPyt: a Python-based free and open-source software for 13C-metabolic flux analyses

¹³C-Metabolic flux analysis (MFA) is a powerful approach to estimate intracellular reaction rates which could be used in strain analysis and design. Processing and analysis of labeling data for calculation of fluxes and associated statistics is an essential part of MFA. However, various software currently available for data analysis employ proprietary platforms and thus limit accessibility. We developed FluxPyt, a Python-based truly open-source software package for conducting stationary ¹³C-MFA data analysis. The software is based on the efficient elementary metabolite unit framework. The standard deviations in the calculated fluxes are estimated using the Monte-Carlo analysis. FluxPyt also automatically creates flux maps based on a template for visualization of the MFA results. The flux distributions calculated by FluxPyt for two separate models: a small tricarboxylic acid cycle model and a larger Corynebacterium glutamicum model, were found to be in good agreement with those calculated by a previously published software. FluxPyt was tested in Microsoft™ Windows 7 and 10, as well as in Linux Mint 18.2. The availability of a free and open ¹³C-MFA software that works in various operating systems will enable more researchers to perform ¹³C-MFA and to further modify and develop the package.

Comprehensive assessment of measurement uncertainty in 13C-based metabolic flux experiments

In the field of metabolic engineering ¹³C-based metabolic flux analysis experiments have proven successful in indicating points of action. As every step of this approach is affected by an inherent error, the aim of the present work is the comprehensive evaluation of factors contributing to the uncertainty of nonnaturally distributed C-isotopologue abundances as well as to the absolute flux value calculation. For this purpose, a previously published data set, analyzed in the course of a ¹³C labeling experiment studying glycolysis and the pentose phosphate pathway in a yeast cell factory, was used. Here, for isotopologue pattern analysis of these highly polar metabolites that occur in multiple isomeric forms, a gas chromatographic separation approach with preceding derivatization was used. This rendered a natural isotope interference correction step essential. Uncertainty estimation of the resulting C-isotopologue distribution was performed according to the EURACHEM guidelines with Monte Carlo simulation. It revealed a significant increase for low-abundance isotopologue fractions after application of the necessary correction step. For absolute flux value estimation, isotopologue fractions of various sugar phosphates, together with the assessed uncertainties, were used in a metabolic model describing the upper part of the central carbon metabolism. The findings pinpointed the influence of small isotopologue fractions as sources of error and highlight the need for improved model curation.

Graphical abstract

Prediction of metabolic fluxes from gene expression data with Huber penalty convex optimization function

As one of the critical parameters of a metabolic pathway, the metabolic flux in a metabolic network serves as an essential role in physiology and pathology. Constraint-based metabolic models are the widely used frameworks for predicting metabolic fluxes in genome-scale metabolic networks. Integrating the transcriptomic data into the constraint-based metabolic models can effectively predict context-specific fluxes across different conditions. However, these methods always need user-defined thresholds to identify the expression levels of metabolic genes or restrain the rate of biomass production, and the predictive results are sensitive to the thresholds. In this work, we present the Huber penalty convex optimization function (HPCOF) combined with the flux minimization principle to predict metabolic fluxes. Our HPCOF method integrates gene expression profiles into the genome-scale metabolic models (GEMs) to reduce the sensitivity to outliers, and uses continuous expression data to avoid selection of arbitrary threshold parameters. In the case studies of Saccharomyces cerevisiae (S. cerevisiae) and Escherichia coli (E. coli) strains under different conditions, the results show that our HPCOF method has a better fit to the experimentally measured values, and has a higher Pearson correlation coefficient, a smaller P-value and a lower sum of squared error than other methods. The HPCOF code can be freely downloaded from for academic users.

Computational Approaches on Stoichiometric and Kinetic Modeling for Efficient Strain Design

Engineering biological systems that are capable of overproducing products of interest is the ultimate goal of any biotechnology application. To this end, stoichiometric (or steady state) and kinetic models are increasingly becoming available for a variety of organisms including prokaryotes, eukaryotes, and microbial communities. This ever-accelerating pace of such model reconstructions has also spurred the development of optimization-based strain design techniques. This chapter highlights a number of such frameworks developed in recent years in order to generate testable hypotheses (in terms of genetic interventions), thus addressing the challenges in metabolic engineering. In particular, three major methods are covered in detail including two methods for designing strains (i.e., one stoichiometric model-based and the other by integrating kinetic information into a stoichiometric model) and one method for analyzing microbial communities.

Identification of metabolic pathways using pathfinding approaches: a systematic review

Metabolic pathways have become increasingly available for various microorganisms. Such pathways have spurred the development of a wide array of computational tools, in particular, mathematical pathfinding approaches. This article can facilitate the understanding of computational analysis of metabolic pathways in genomics. Moreover, stoichiometric and pathfinding approaches in metabolic pathway analysis are discussed. Three major types of studies are elaborated: stoichiometric identification models, pathway-based graph analysis and pathfinding approaches in cellular metabolism. Furthermore, evaluation of the outcomes of the pathways with mathematical benchmarking metrics is provided. This review would lead to better comprehension of metabolism behaviors in living cells, in terms of computed pathfinding approaches.

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.

Simultaneous tracing of carbon and nitrogen isotopes in human cells

Stable isotope tracing is a powerful method for interrogating metabolic enzyme activities across the metabolic network of living cells. However, most studies of mammalian cells have used (13)C-labeled tracers only and focused on reactions in central carbon metabolism. Cellular metabolism, however, involves other biologically important elements, including nitrogen, hydrogen, oxygen, phosphate and sulfur. Tracing stable isotopes of such elements may help shed light on poorly understood metabolic pathways. Here, we demonstrate the use of high-resolution mass spectrometry to simultaneously trace carbon and nitrogen metabolism in human cells cultured with (13)C- and (15)N-labeled glucose and glutamine. To facilitate interpretation of the complex isotopomer data generated, we extend current methods for metabolic flux analysis to handle multivariate mass isotopomer distributions (MMIDs). We find that observed MMIDs are broadly consistent with known biochemical pathways. Whereas measured (13)C MIDs were informative for central carbon metabolism, (15)N isotopes provided evidence for nitrogen-carrying reactions in amino acid and nucleotide metabolism. This computational and experimental methodology expands the scope of metabolic flux analysis beyond carbon metabolism, and may prove important to understanding metabolic phenotypes in health and disease.

A comprehensive evaluation of constraining amino acid biosynthesis in compartmented models for metabolic flux analysis

Recent advances in the availability and applicability of genetic tools for non-conventional yeasts have raised high hopes regarding the industrial applications of such yeasts; however, quantitative physiological data on these yeasts, including intracellular flux distributions, are scarce and have rarely aided in the development of novel yeast applications. The compartmentation of eukaryotic cells adds to model complexity. Model constraints are ideally based on biochemical evidence, which is rarely available for non-conventional yeast and eukaryotic cells. A small-scale model for ¹³C-based metabolic flux analysis of central yeast carbon metabolism was developed that is universally valid and does not depend on localization information regarding amino acid anabolism. The variable compartmental origin of traced metabolites is a feature that allows application of the model to yeasts with uncertain genomic and transcriptional backgrounds. The presented test case includes the baker's yeast Saccharomyces cerevisiae and the methylotrophic yeast Hansenula polymorpha. Highly similar flux solutions were computed using either a model with undefined pathway localization or a model with constraints based on curated (S. cerevisiae) or computationally predicted (H. polymorpha) localization information, while false solutions were found with incorrect localization constraints. These results indicate a potentially adverse effect of universally assuming Saccharomyces-like constraints on amino acid biosynthesis for non-conventional yeasts and verify the validity of neglecting compartmentation constraints using a small-scale metabolic model. The model was specifically designed to investigate the intracellular metabolism of wild-type yeasts under various growth conditions but is also expected to be useful for computing fluxes of other eukaryotic cells.

•

Compartmentation influences computed intracellular fluxes.

•

Improper localization constraints potentially produce false flux solutions.

•

Minimal compartmentation constraints result in high-quality flux computations.

Metabolic flux analysis of heterotrophic growth in Chlamydomonas reinhardtii

Despite the wealth of knowledge available for C. reinhardtii, the central metabolic fluxes of growth on acetate have not yet been determined. In this study, ¹³C-metabolic flux analysis (¹³C-MFA) was used to determine and quantify the metabolic pathways of primary metabolism in C. reinhardtii cells grown under heterotrophic conditions with acetate as the sole carbon source. Isotopic labeling patterns of compartment specific biomass derived metabolites were used to calculate the fluxes. It was found that acetate is ligated with coenzyme A in the three subcellular compartments (cytosol, mitochondria and plastid) included in the model. Two citrate synthases were found to potentially be involved in acetyl-coA metabolism; one localized in the mitochondria and the other acting outside the mitochondria. Labeling patterns demonstrate that Acetyl-coA synthesized in the plastid is directly incorporated in synthesis of fatty acids. Despite having a complete TCA cycle in the mitochondria, it was also found that a majority of the malate flux is shuttled to the cytosol and plastid where it is converted to oxaloacetate providing reducing equivalents to these compartments. When compared to predictions by flux balance analysis, fluxes measured with ¹³C-MFA were found to be suboptimal with respect to biomass yield; C. reinhardtii sacrifices biomass yield to produce ATP and reducing equivalents.

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.

Synergizing 13C Metabolic Flux Analysis and Metabolic Engineering for Biochemical Production

Metabolic engineering of industrial microorganisms to produce chemicals, fuels, and drugs has attracted increasing interest as it provides an environment-friendly and renewable route that does not depend on depleting petroleum sources. However, the microbial metabolism is so complex that metabolic engineering efforts often have difficulty in achieving a satisfactory yield, titer, or productivity of the target chemical. To overcome this challenge, ¹³C Metabolic Flux Analysis (¹³C-MFA) has been developed to investigate rigorously the cell metabolism and quantify the carbon flux distribution in central metabolic pathways. In the past decade, ¹³C-MFA has been widely used in academic labs and the biotechnology industry to pinpoint the key issues related to microbial-based chemical production and to guide the development of the appropriate metabolic engineering strategies for improving the biochemical production. In this chapter we introduce the basics of ¹³C-MFA and illustrate how ¹³C-MFA has been applied to synergize with metabolic engineering to identify and tackle the rate-limiting steps in biochemical production.