13C-Labeled metabolic flux analysis of a fed-batch culture of elutriated Saccharomyces cerevisiae

Comparative study of two Saccharomyces cerevisiae strains with kinetic models at genome-scale

The parameterization of kinetic models requires measurement of fluxes and/or metabolite levels for a base strain and a few genetic perturbations thereof. Unlike stoichiometric models that are mostly invariant to the specific strain, it remains unclear whether kinetic models constructed for different strains of the same species have similar or significantly different kinetic parameters. This important question underpins the applicability range and prediction limits of kinetic reconstructions. To this end, herein we parameterize two separate large-scale kinetic models using K-FIT with genome-wide coverage corresponding to two distinct strains of Saccharomyces cerevisiae: CEN.PK 113-7D strain (model k-sacce306-CENPK), and growth-deficient BY4741 (isogenic to S288c; model k-sacce306-BY4741). The metabolic network for each model contains 306 reactions, 230 metabolites, and 119 substrate-level regulatory interactions. The two models (for CEN.PK and BY4741) recapitulate, within one standard deviation, 77% and 75% of the fitted dataset fluxes, respectively, determined by ¹³C metabolic flux analysis for wild-type and eight single-gene knockout mutants of each strain. Strain-specific kinetic parameterization results indicate that key enzymes in the TCA cycle, glycolysis, and arginine and proline metabolism drive the metabolic differences between these two strains of S. cerevisiae. Our results suggest that although kinetic models cannot be readily used across strains as stoichiometric models, they can capture species-specific information through the kinetic parameterization process.

mfapy: An open-source Python package for 13C-based metabolic flux analysis

¹³C-based metabolic flux analysis (¹³C-MFA) is an essential tool for estimating intracellular metabolic flux levels in metabolic engineering and biology. In ¹³C-MFA, a metabolic flux distribution that explains the observed isotope labeling data was computationally estimated using a non-linear optimization method. Herein, we report the development of mfapy, an open-source Python package developed for more flexibility and extensibility for ¹³C-MFA. mfapy compels users to write a customized Python code by describing each step in the data analysis procedures of the isotope labeling experiments. The flexibility and extensibility provided by mfapy can support trial-and-error performance in the routine estimation of metabolic flux distributions, experimental design by computer simulations of ¹³C-MFA experiments, and development of new data analysis techniques for stable isotope labeling experiments. mfapy is available to the public from the Github repository (https://github.com/fumiomatsuda/mfapy).

•

An open-source Python package, mfapy, is developed for ¹³C-MFA.

•

mfapy enables users to write Python codes for data analysis procedures of ¹³C-MFA.

•

mfapy has a flexibility and extensibility to support various data analysis procedures.

•

Computer simulations of ¹³C-MFA experiments is supported for experimental design.

Physicochemical and metabolic constraints for thermodynamics-based stoichiometric modelling under mesophilic growth conditions

Metabolic engineering in the post-genomic era is characterised by the development of new methods for metabolomics and fluxomics, supported by the integration of genetic engineering tools and mathematical modelling. Particularly, constraint-based stoichiometric models have been widely studied: (i) flux balance analysis (FBA) (in silico), and (ii) metabolic flux analysis (MFA) (in vivo). Recent studies have enabled the incorporation of thermodynamics and metabolomics data to improve the predictive capabilities of these approaches. However, an in-depth comparison and evaluation of these methods is lacking. This study presents a thorough analysis of two different in silico methods tested against experimental data (metabolomics and ¹³C-MFA) for the mesophile Escherichia coli. In particular, a modified version of the recently published matTFA toolbox was created, providing a broader range of physicochemical parameters. Validating against experimental data allowed the determination of the best physicochemical parameters to perform the TFA (Thermodynamics-based Flux Analysis). An analysis of flux pattern changes in the central carbon metabolism between ¹³C-MFA and TFA highlighted the limited capabilities of both approaches for elucidating the anaplerotic fluxes. In addition, a method based on centrality measures was suggested to identify important metabolites that (if quantified) would allow to further constrain the TFA. Finally, this study emphasised the need for standardisation in the fluxomics community: novel approaches are frequently released but a thorough comparison with currently accepted methods is not always performed.

Biotechnology has benefitted from the development of high throughput methods characterising living systems at different levels (e.g. concerning genes or proteins), allowing the industrial production of chemical commodities. Recently, focus has been placed on determining reaction rates (or metabolic fluxes) in the metabolic network of certain microorganisms, in order to identify bottlenecks hindering their exploitation. Two main approaches are commonly used, termed metabolic flux analysis (MFA) and flux balance analysis (FBA), based on measuring and estimating fluxes, respectively. While the influence of thermodynamics in living systems was accepted several decades ago, its application to study biochemical networks has only recently been enabled. In this sense, a multitude of different approaches constraining well-established modelling methods with thermodynamics has been suggested. However, physicochemical parameters are generally not properly adjusted to the experimental conditions, which might affect their predictive capabilities. In this study, we have explored the reliability of currently available tools by investigating the impact of varying said parameters in the simulation of metabolic fluxes and metabolite concentration values. Additionally, our in-depth analysis allowed us to highlight limitations and potential solutions that should be considered in future studies.

Metabolic switches from quiescence to growth in synchronized Saccharomyces cerevisiae

INTRODUCTION

The switch from quiescence (G0) into G1 and cell cycle progression critically depends on specific nutrients and metabolic capabilities. Conversely, metabolic networks are regulated by enzyme-metabolite interaction and transcriptional regulation that lead to flux modifications to support cell growth. How cells process and integrate environmental information into coordinated responses is challenging to analyse and not yet described quantitatively.

OBJECTIVES

To quantitatively monitor the central carbon metabolism during G0 exit and the first 2 h after reentering the cell cycle from synchronized Saccharomyces cerevisiae.

METHODS

Dynamic tailored ¹³C metabolic flux analysis was used to observe the intracellular metabolite flux changes, and the metabolome and proteome were observed to identify regulatory mechanisms.

RESULTS

G0 cells responded immediately to an extracellular increase of glucose. The intracellular metabolic flux changed in time and specific events were observed. High fluxes into trehalose and glycogen synthesis were observed during the G0 exit. Both fluxes then decreased, reaching a minimum at t = 65 min. Here, storage degradation contributed significantly (i.e. 21%) to the glycolytic flux. In contrast to these changes, the glucose uptake rate remained constant after the G0 exit. The flux into the oxidative pentose phosphate pathway was highest (29-fold increase, 36.4% of the glucose uptake) at t = 65 min, while it was very low at other time points. The maximum flux seems to correlate with a late G1 state preparing for the S phase transition. In the G1/S phase (t = 87 min), anaplerotic reactions such as glyoxylate shunt increased. Protein results show that during this transition, proteins belonging to clusters related with ribosome biogenesis and assembly, and initiation transcription factors clusters were continuously synthetised.

CONCLUSION

The intracellular flux distribution changes dynamically and these major rearrangements highlight the coordinate reorganization of metabolic flux to meet requirements for growth during different cell state.

Temporal fluxomics reveals oscillations in TCA cycle flux throughout the mammalian cell cycle

Cellular metabolic demands change throughout the cell cycle. Nevertheless, a characterization of how metabolic fluxes adapt to the changing demands throughout the cell cycle is lacking. Here, we developed a temporal‐fluxomics approach to derive a comprehensive and quantitative view of alterations in metabolic fluxes throughout the mammalian cell cycle. This is achieved by combining pulse‐chase LC‐MS‐based isotope tracing in synchronized cell populations with computational deconvolution and metabolic flux modeling. We find that TCA cycle fluxes are rewired as cells progress through the cell cycle with complementary oscillations of glucose versus glutamine‐derived fluxes: Oxidation of glucose‐derived flux peaks in late G1 phase, while oxidative and reductive glutamine metabolism dominates S phase. These complementary flux oscillations maintain a constant production rate of reducing equivalents and oxidative phosphorylation flux throughout the cell cycle. The shift from glucose to glutamine oxidation in S phase plays an important role in cell cycle progression and cell proliferation.

Investigation of useful carbon tracers for 13C-metabolic flux analysis of Escherichia coli by considering five experimentally determined flux distributions

The ¹³C-MFA experiments require an optimal design since the precision or confidence intervals of the estimated flux levels depends on factors such as the composition of ¹³C-labeled carbon sources, as well as the metabolic flux distribution of interest. In this study, useful compositions of ¹³C-labeled glucose for ¹³C-metabolic flux analysis (¹³C-MFA) of Escherichia coli are investigated using a computer simulation of the stable isotope labeling experiment. Following the generation of artificial mass spectra datasets of amino acid fragments using five literature-reported flux distributions of E. coli, the best fitted flux distribution and the 95% confidence interval were estimated by the ¹³C-MFA procedure. A comparison of the precision scores showed that [1, 2-¹³C]glucose and a mixture of [1-¹³C] and [U-¹³C]glucose at 8:2 are one of the best carbon sources for a precise estimation of flux levels of the pentose phosphate pathway, glycolysis and the TCA cycle. Although the precision scores of the anaplerotic and glyoxylate pathway reactions were affected by both the carbon source and flux distribution, it was also shown that the mixture of non-labeled, [1-¹³C], and [U-¹³C]glucose at 4:1:5 was specifically effective for the flux estimation of the glyoxylate pathway reaction. These findings were confirmed by wet ¹³C-MFA experiments.

•

Useful compositions of ¹³C-labeled glucose are investigated for ¹³C-MFA of E. coli.

•

Computer simulations revealed that [1,2-¹³C] was one of the best first choices.

•

Mixture of non-labeled, [1-¹³C] and [U-¹³C] at 0:8:2 was also suitable for ¹³C-MFA.

•

Mixture at 4:1:5 was specifically effective for estimation of glyoxylate pathway.

•

The wet ¹³C-MFA experiments of E. coli confirmed the findings.

Achieving Metabolic Flux Analysis for S. cerevisiae at a Genome-Scale: Challenges, Requirements, and Considerations

Recent advances in 13C-Metabolic flux analysis (13C-MFA) have increased its capability to accurately resolve fluxes using a genome-scale model with narrow confidence intervals without pre-judging the activity or inactivity of alternate metabolic pathways. However, the necessary precautions, computational challenges, and minimum data requirements for successful analysis remain poorly established. This review aims to establish the necessary guidelines for performing 13C-MFA at the genome-scale for a compartmentalized eukaryotic system such as yeast in terms of model and data requirements, while addressing key issues such as statistical analysis and network complexity. We describe the various approaches used to simplify the genome-scale model in the absence of sufficient experimental flux measurements, the availability and generation of reaction atom mapping information, and the experimental flux and metabolite labeling distribution measurements to ensure statistical validity of the obtained flux distribution. Organism-specific challenges such as the impact of compartmentalization of metabolism, variability of biomass composition, and the cell-cycle dependence of metabolism are discussed. Identification of errors arising from incorrect gene annotation and suggested alternate routes using MFA are also highlighted.

¹³C-based metabolic flux analysis of Saccharomyces cerevisiae with a reduced Crabtree effect

Saccharomyces cerevisiae shows a Crabtree effect that produces ethanol in a high glucose concentration even under fully aerobic condition. For efficient production of cake yeast or compressed yeast for baking, ethanol by-production is not desired since glucose limited chemostat or fed-batch cultivations are performed to suppress the Crabtree effect. In this study, the (13)C-based metabolic flux analysis ((13)C-MFA) was performed for the S288C derived S. cerevisiae strain to characterize a metabolic state under the reduced Crabtree effect. S. cerevisiae cells were cultured at a low dilution rate (0.1 h(-1)) under the glucose-limited chemostat condition. The estimated metabolic flux distribution showed that the acetyl-CoA in mitochondria was mainly produced from pyruvate by pyruvate dehydrogenase (PDH) reaction and that the level of the metabolic flux through the pentose phosphate pathway was much higher than that of the Embden-Meyerhof-Parnas pathway, which contributes to high biomass yield at low dilution rate by supplying NADPH required for cell growth.

(13)C metabolic flux analysis of recombinant expression hosts

Identifying host cell metabolic phenotypes that promote high recombinant protein titer is a major goal of the biotech industry. (13)C metabolic flux analysis (MFA) provides a rigorous approach to quantify these metabolic phenotypes by applying isotope tracers to map the flow of carbon through intracellular metabolic pathways. Recent advances in tracer theory and measurements are enabling more information to be extracted from (13)C labeling experiments. Sustained development of publicly available software tools and standardization of experimental workflows is simultaneously encouraging increased adoption of (13)C MFA within the biotech research community. A number of recent (13)C MFA studies have identified increased citric acid cycle and pentose phosphate pathway fluxes as consistent markers of high recombinant protein expression, both in mammalian and microbial hosts. Further work is needed to determine whether redirecting flux into these pathways can effectively enhance protein titers while maintaining acceptable glycan profiles.

OpenMebius: an open source software for isotopically nonstationary 13C-based metabolic flux analysis

The in vivo measurement of metabolic flux by ¹³C-based metabolic flux analysis (¹³C-MFA) provides valuable information regarding cell physiology. Bioinformatics tools have been developed to estimate metabolic flux distributions from the results of tracer isotopic labeling experiments using a ¹³C-labeled carbon source. Metabolic flux is determined by nonlinear fitting of a metabolic model to the isotopic labeling enrichment of intracellular metabolites measured by mass spectrometry. Whereas ¹³C-MFA is conventionally performed under isotopically constant conditions, isotopically nonstationary ¹³C metabolic flux analysis (INST-¹³C-MFA) has recently been developed for flux analysis of cells with photosynthetic activity and cells at a quasi-steady metabolic state (e.g., primary cells or microorganisms under stationary phase). Here, the development of a novel open source software for INST-¹³C-MFA on the Windows platform is reported. OpenMebius (Open source software for Metabolic flux analysis) provides the function of autogenerating metabolic models for simulating isotopic labeling enrichment from a user-defined configuration worksheet. Analysis using simulated data demonstrated the applicability of OpenMebius for INST-¹³C-MFA. Confidence intervals determined by INST-¹³C-MFA were less than those determined by conventional methods, indicating the potential of INST-¹³C-MFA for precise metabolic flux analysis. OpenMebius is the open source software for the general application of INST-¹³C-MFA.

Reliable Metabolic Flux Estimation in Escherichia coli Central Carbon Metabolism Using Intracellular Free Amino Acids

¹³C metabolic flux analysis (MFA) is a tool of metabolic engineering for investigation of in vivo flux distribution. A direct ¹³C enrichment analysis of intracellular free amino acids (FAAs) is expected to reduce time for labeling experiments of the MFA. Measurable FAAs should, however, vary among the MFA experiments since the pool sizes of intracellular free metabolites depend on cellular metabolic conditions. In this study, minimal ¹³C enrichment data of FAAs was investigated to perform the FAAs-based MFA. An examination of a continuous culture of Escherichia coli using ¹³C-labeled glucose showed that the time required to reach an isotopically steady state for FAAs is rather faster than that for conventional method using proteinogenic amino acids (PAAs). Considering 95% confidence intervals, it was found that the metabolic flux distribution estimated using FAAs has a similar reliability to that of the PAAs-based method. The comparative analysis identified glutamate, aspartate, alanine and phenylalanine as the common amino acids observed in E. coli under different culture conditions. The results of MFA also demonstrated that the ¹³C enrichment data of the four amino acids is required for a reliable analysis of the flux distribution.

Metabolic flux analysis using ¹³C peptide label measurements

¹³C metabolic flux analysis (MFA) has become the experimental method of choice to investigate the cellular metabolism of microbes, cell cultures and plant seeds. Conventional steady-state MFA utilizes isotopic labeling measurements of amino acids obtained from protein hydrolysates. To retain spatial information in conventional steady-state MFA, tissues or subcellular fractions must be dissected or biochemically purified. In contrast, peptides retain their identity in complex protein extracts, and may therefore be associated with a specific time of expression, tissue type and subcellular compartment. To enable 'single-sample' spatially and temporally resolved steady-state flux analysis, we investigated the suitability of peptide mass distributions (PMDs) as an alternative to amino acid label measurements. PMDs are the discrete convolution of the mass distributions of the constituent amino acids of a peptide. We investigated the requirements for the unique deconvolution of PMDs into amino acid mass distributions (AAMDs), the influence of peptide sequence length on parameter sensitivity, and how AAMD and flux estimates that are determined through deconvolution compare to estimates from a conventional GC-MS measurement-based approach. Deconvolution of PMDs of the storage protein β-conglycinin of soybean (Glycine max) resulted in good AAMD and flux estimates if fluxes were directly fitted to PMDs. Unconstrained deconvolution resulted in inferior AAMD and flux estimates. PMD measurements do not include amino acid backbone fragments, which increase the information content in GC-MS-derived analyses. Nonetheless, the resulting flux maps were of comparable quality due to the precision of Orbitrap quantification and the larger number of peptide measurements.

Labelling analysis for ¹³C MFA using NMR spectroscopy

NMR spectroscopy is an efficient method for analyzing (13)C labelling of cellular metabolites. The strength of it is especially the ability to provide direct quantitative positional information on the (13)C labelling status of carbon atoms in metabolites. NMR spectroscopic methods allow also for detection of contiguously (13)C-labelled fragments in the carbon backbones of the metabolites. Furthermore, the recent developments of NMR spectroscopy hardware have substantially improved the sensitivity of the methods. In this chapter we describe a method for analyzing the (13)C labelling of the biomass amino acids for metabolic flux analysis, sample preparation for NMR spectroscopy, acquiring and processing the NMR spectra, and extracting the (13)C labelling information from the NMR data. Different NMR methods are applied depending on the (13)C labelling strategy chosen. These strategies include uniform (13)C labelling, positional (13)C labelling, or a combination of both. Not only the preparation of sample for analysis of (13)C labelling in proteinogenic amino acids in biomass is described, but also the necessary modifications to the method when analysis of (13)C labelling in free metabolic intermediates is of interest. Finally the strategies for using the different NMR-detected (13)C labelling data in (13)C MFA are discussed.

Using multiple tracers for 13C metabolic flux analysis

(13)C-Metabolic flux analysis ((13)C-MFA) is a powerful technique for quantifying intracellular metabolic fluxes in living cells. These in vivo fluxes provide important information on the physiology of cells in culture, which can be used for metabolic engineering purposes and serve as inputs for systems biology modeling. The (13)C-MFA technique consists of several steps: (1) selecting appropriate tracers for a given system of interest, (2) performing isotopic labeling experiments, (3) measuring isotopic labeling distributions in metabolic products, (4) estimating metabolic fluxes using least-squares regression, and (5) evaluating the goodness of fit and computing confidence intervals for estimated fluxes. In this chapter, we provide guidelines for performing (13)C-MFA studies using multiple isotopic tracers, a technique that is especially useful for elucidating fluxes in complex biological systems where multiple carbon sources are present. Here, as an example, we describe key steps and decision points for designing (13)C-MFA studies for microbes grown on mixtures of glucose and xylose. The general concepts described in this chapter are applicable to many other biological systems. For example, the same procedures can be applied to design (13)C-MFA studies in mammalian cells, which are generally grown in complex media containing multiple substrates such as glucose and amino acids.

Systematic applications of metabolomics in metabolic engineering

The goals of metabolic engineering are well-served by the biological information provided by metabolomics: information on how the cell is currently using its biochemical resources is perhaps one of the best ways to inform strategies to engineer a cell to produce a target compound. Using the analysis of extracellular or intracellular levels of the target compound (or a few closely related molecules) to drive metabolic engineering is quite common. However, there is surprisingly little systematic use of metabolomics datasets, which simultaneously measure hundreds of metabolites rather than just a few, for that same purpose. Here, we review the most common systematic approaches to integrating metabolite data with metabolic engineering, with emphasis on existing efforts to use whole-metabolome datasets. We then review some of the most common approaches for computational modeling of cell-wide metabolism, including constraint-based models, and discuss current computational approaches that explicitly use metabolomics data. We conclude with discussion of the broader potential of computational approaches that systematically use metabolomics data to drive metabolic engineering.

Metabolic modelling in the development of cell factories by synthetic biology

Cell factories are commonly microbial organisms utilized for bioconversion of renewable resources to bulk or high value chemicals. Introduction of novel production pathways in chassis strains is the core of the development of cell factories by synthetic biology. Synthetic biology aims to create novel biological functions and systems not found in nature by combining biology with engineering. The workflow of the development of novel cell factories with synthetic biology is ideally linear which will be attainable with the quantitative engineering approach, high-quality predictive models, and libraries of well-characterized parts. Different types of metabolic models, mathematical representations of metabolism and its components, enzymes and metabolites, are useful in particular phases of the synthetic biology workflow. In this minireview, the role of metabolic modelling in synthetic biology will be discussed with a review of current status of compatible methods and models for the in silico design and quantitative evaluation of a cell factory.

Metabolic flux analysis and visualization

One of the ultimate goals of systems biology research is to obtain a comprehensive understanding of the control mechanisms of complex cellular metabolisms. Metabolic Flux Analysis (MFA) is a important method for the quantitative estimation of intracellular metabolic flows through metabolic pathways and the elucidation of cellular physiology. The primary challenge in the use of MFA is that many biological networks are underdetermined systems; it is therefore difficult to narrow down the solution space from the stoichiometric constraints alone. In this tutorial, we present an overview of Flux Balance Analysis (FBA) and (13)C-Metabolic Flux Analysis ((13)C-MFA), both of which are frequently used to solve such underdetermined systems, and we demonstrate FBA and (13)C-MFA using the genome-scale model and the central carbon metabolism model, respectively. Furthermore, because such comprehensive study of intracellular fluxes is inherently complex, we subsequently introduce various pathway mapping and visualization tools to facilitate understanding of these data in the context of the pathways. Specific visualization of MFA results using the BioCyc Omics Viewer and Pathway Projector are shown as illustrative examples.

Integration of metabolic reactions and gene regulation

Metabolic reactions and gene regulation are two primary processes of cells. In response to environmental changes cells often adjust the regulatory programs and shift the metabolic states. An integrative investigation and modeling of these two processes would improve our understanding about the cellular systems and may generate substantial impacts in medicine, agriculture, environmental protection, and energy production. We review the studies of the various aspects of the crosstalk between metabolic reactions and gene regulation, including models, empirical evidence, and available databases.

13C-metabolic flux analysis for batch culture of Escherichia coli and its Pyk and Pgi gene knockout mutants based on mass isotopomer distribution of intracellular metabolites

Since most bio-production processes are conducted in a batch or fed-batch manner, the evaluation of metabolism with respect to time is highly desirable. Toward this aim, we applied (13)C-metabolic flux analysis to nonstationary conditions by measuring the mass isotopomer distribution of intracellular metabolites. We performed our analysis on batch cultures of wild-type Escherichia coli, as well as on Pyk and Pgi mutants, obtained the fluxes and metabolite concentrations as a function of time. Our results for the wild-type indicated that the TCA cycle flux tended to increase during growth on glucose. Following glucose exhaustion, cells controlled the branch ratio between the glyoxylate pathway and the TCA cycle, depending on the availability of acetate. In the Pyk mutant, the concentrations of glycolytic intermediates changed drastically over time due to the dumping and feedback inhibition caused by PEP accumulation. Nevertheless, the flux distribution and free amino acid concentrations changed little. The growth rate and the fluxes remained constant in the Pgi mutant and the glucose-6-phosphate dehydrogenase reaction was the rate-limiting step. The measured fluxes were compared with those predicted by flux balance analysis using maximization of biomass yield or ATP production. Our findings indicate that the objective function of biosynthesis became less important as time proceeds on glucose in the wild-type, while it remained highly important in the Pyk mutant. Furthermore, ATP production was the primary objective function in the Pgi mutant. This study demonstrates how cells adjust their metabolism in response to environmental changes and/or genetic perturbations in the batch cultivation.

Yeast dynamic metabolic flux measurement in nutrient-rich media by HPLC and accelerator mass spectrometry

Metabolic flux, the flow of metabolites through networks of enzymes, represents the dynamic productive output of cells. Improved understanding of intracellular metabolic fluxes will enable targeted manipulation of metabolic pathways of medical and industrial importance to a greater degree than is currently possible. Flux balance analysis (FBA) is a constraint-based approach to modeling metabolic fluxes, but its utility is limited by a lack of experimental measurements. Incorporation of experimentally measured fluxes as system constraints will significantly improve the overall accuracy of FBA. We applied a novel, two-tiered approach in the yeast Saccharomyces cerevisiae to measure nutrient consumption rates (extracellular fluxes) and a targeted intracellular flux using a (14)C-labeled precursor with HPLC separation and flux quantitation by accelerator mass spectrometry (AMS). The use of AMS to trace the intracellular fate of (14)C-glutamine allowed the calculation of intracellular metabolic flux through this pathway, with glutathione as the metabolic end point. Measured flux values provided global constraints for the yeast FBA model which reduced model uncertainty by more than 20%, proving the importance of additional constraints in improving the accuracy of model predictions and demonstrating the use of AMS to measure intracellular metabolic fluxes. Our results highlight the need to use intracellular fluxes to constrain the models. We show that inclusion of just one such measurement alone can reduce the average variability of model predicted fluxes by 10%.

Dynamic flux responses in riboflavin overproducing Bacillus subtilis to increasing glucose limitation in fed-batch culture

How do intracellular fluxes respond to dynamically increasing glucose limitation when the physiology changes from strong overflow metabolism near to exclusively maintenance metabolism? Here we investigate this question in a typical industrial, glucose-limited fed-batch cultivation with a riboflavin overproducing Bacillus subtilis strain. To resolve dynamic flux changes, a novel approach to (13)C flux analysis was developed that is based on recording (13)C labeling patterns in free intracellular amino acids. Fluxes are then estimated with stationary flux ratio and iterative isotopomer balancing methods, for which a decomposition of the process into quasi-steady states and estimation of isotopic steady state (13)C labeling patterns was necessary. By this approach, we achieve a temporal resolution of 30-60 min that allows us to resolve the slow metabolic transients that typically occur in such cultivations. In the late process phase we found, most prominently, almost exclusive respiratory metabolism, significantly increased pentose phosphate pathway contribution and a strongly decreased futile cycle through the PEP carboxykinase. As a consequence, higher catabolic NADPH formation occurred than was necessary to satisfy the anabolic demands, suggesting a transhydrogenase-like mechanism to close the balance of reducing equivalents.

Metabolic flux analysis in eukaryotes

Metabolic flux analysis (MFA) represents a powerful tool for systems biology research on eukaryotic cells. This review describes recent advances, the challenges as well as applications of metabolic flux analysis comprising fungi, mammalian cells and plants. While MFA is widely established and applied in microorganisms, it remains still a challenge to adapt these methods to eukaryotic cell systems having a higher complexity particularly concerning compartmentation or media composition. In fungi MFA was used in the past few years to analyze a variety of conditions and factors and their effects on cellular metabolism. In mammalian cells MFA was applied mainly in cell culture technology and in medical and toxicological research. (13)C metabolic studies on native whole plants are additionally challenging by the fact that CO(2) is usually the only carbon source.