Genealogy profiling through strain improvement by using metabolic network analysis: metabolic flux genealogy of several generations of lysine-producing corynebacteria

A comprehensive approach of metabolite balancing, (13)C tracer studies, gas chromatography-mass spectrometry, matrix-assisted laser desorption ionization-time of flight mass spectrometry, and isotopomer modeling was applied for comparative metabolic network analysis of a genealogy of five successive generations of lysine-producing Corynebacterium glutamicum. The five strains examined (C. glutamicum ATCC 13032, 13287, 21253, 21526, and 21543) were previously obtained by random mutagenesis and selection. Throughout the genealogy, the lysine yield in batch cultures increased markedly from 1.2 to 24.9% relative to the glucose uptake flux. Strain optimization was accompanied by significant changes in intracellular flux distributions. The relative pentose phosphate pathway (PPP) flux successively increased, clearly corresponding to the product yield. Moreover, the anaplerotic net flux increased almost twofold as a consequence of concerted regulation of C(3) carboxylation and C(4) decarboxylation fluxes to cover the increased demand for lysine formation; thus, the overall increase was a consequence of concerted regulation of C(3) carboxylation and C(4) decarboxylation fluxes. The relative flux through isocitrate dehydrogenase dropped from 82.7% in the wild type to 59.9% in the lysine-producing mutants. In contrast to the NADPH demand, which increased from 109 to 172% due to the increasing lysine yield, the overall NADPH supply remained constant between 185 and 196%, resulting in a decrease in the apparent NADPH excess through strain optimization. Extrapolated to industrial lysine producers, the NADPH supply might become a limiting factor. The relative contributions of PPP and the tricarboxylic acid cycle to NADPH generation changed markedly, indicating that C. glutamicum is able to maintain a constant supply of NADPH under completely different flux conditions. Statistical analysis by a Monte Carlo approach revealed high precision for the estimated fluxes, underlining the fact that the observed differences were clearly strain specific.     

The relationships between the metabolic fluxes and 13C-labeled isotopomer distribution for the flux analysis of the main metabolic pathways

It has been known that ¹³C-labeling technique is quite useful in estimating the metabolic fluxes. Although the program-based flux analysis is powerful, it is not easy to be confident with the result obtained without experiences and exhaustive trial and errors based on statistical analysis for the confidence intervals in practice. It is, therefore, quite important to grasp the relationship between the fluxes and the ¹³C-labeled isotopomer distribution to get deeper insight into the metabolic flux analysis. In the present research, it was shown explicitly how the isotopomer distribution changes with respect to the fluxes in relation to the labeling patterns of the substrate, where either labeled glucose, acetate, or pyruvate was used as a carbon source. Some of the analytical expressions were derived based on the matrix representation, and they were utilized for analysis. It was shown that the isotopomer pattern does not necessarily change uniformly with respect to fluxes, but changes in a complicated way in particular for the case of using pyruvate as a carbon source where some isotopomers do not necessarily change monotonically. It was shown to be quite important to grasp how the isotopomer pattern changes with respect to fluxes and the labeling pattern of the substrate for flux determination and the experimental design. It was also shown that the mixture of [1-¹³C] acetate and [2-¹³C] acetate should not be used from the information index point of view. Some of the experimental data were evaluated from the present approach. It was also shown that the isotopomer distribution is less sensitive to the bidirectional fluxes in the reversible pathway.     

Evidence for transketolase-like TKTL1 flux in CHO cells based on parallel labeling experiments and 13C-metabolic flux analysis

The pentose phosphate pathway (PPP) is a fundamental component of cellular metabolism. It provides precursors for the biosynthesis of nucleotides and contributes to the production of reducing power in the form of NADPH. It has been hypothesized that mammalian cells may contain a hidden reaction in PPP catalyzed by transketolase-like protein 1 (TKTL1) that is closely related to the classical transketolase enzyme; however, until now there has been no direct experimental evidence for this reaction. In this work, we have applied state-of-the-art techniques in ¹³C metabolic flux analysis (¹³C-MFA) based on parallel labeling experiments and integrated flux fitting to estimate the TKTL1 flux in CHO cells. We identified a set of three parallel labeling experiments with [1-¹³C]glucose+[4,5,6-¹³C]glucose, [2-¹³C]glucose+[4,5,6-¹³C]glucose, and [3-¹³C]glucose+[4,5,6-¹³C]glucose and developed a new method to measure ¹³C-labeling of fructose 6-phosphate by GC-MS that allows intuitive interpretation of mass isotopomer distributions to determine key fluxes in the model, including glycolysis, oxidative PPP, non-oxidative PPP, and the TKTL1 flux. Using these tracers we detected a significant TKTL1 flux in CHO cells at the stationary phase. The flux results suggest that the main function of oxidative PPP in CHO cells at the stationary phase is to fuel the TKTL1 reaction. Overall, this study demonstrates for the first time that carbon atoms can be lost in the PPP, by means other than the oxidative PPP, and that this loss of carbon atoms is consistent with the hypothesized TKTL1 reaction in mammalian cells.     

Application of metabolic flux analysis for the identification of metabolic bottlenecks in the biosynthesis of penicillin-G

A detailed stoichiometric model was developed for growth and penicillin-G production in Penicillium chrysogenum. From an a priori metabolic flux analysis using this model it appeared that penicillin production requires significant changes in fluxes through the primary metabolic pathways. This is brought about by the biosynthesis of carbon precursors for the beta-lactan nucleus and an increased demand for NADPH, mainly for sulfate reduction. As a result, significant changes in flux partitioning occur around four principal nodes in primary metabolism. These are located at: (1) glucose-6-phosphate; (2) 3-phosphoglycerate; (3) mitochondrial pyruvate; and (4) mitochondrial isocitrate. These nodes should be regarded as potential bottlenecks for increased productivity. The flexibility of these principal nodes was investigated by experimental manipulation of the fluxes through the central metabolic pathways using a high-producing strain of P. chrysogenum. Metabolic fluxes were manipulated through growth of the cells on different substrates in carbon-limited chemostat culture. Metabolic flux analysis, based on measured input and output fluxes, was used to calculate the fluxes around the principal nodes. It was found that, for growth on glucose, ethanol, and acetate, the flux partitioning around these nodes differed significantly. However, this had hardly any effect on penicillin productivity, showing that primary carbon metabolism is not likely to contain potential bottlenecks. Further experiments were performed to manipulate the total metabolic demand for the cofactor nicotinamide adenine dinucleotide phosphate (NADPH). NADPH demand was increased stepwise by cultivating the cells on glucose or xylose as the carbon source combined with either ammonia or nitrate as the nitrogen source, which resulted in a stepwise decrease of penicillin production. This clearly shows that, in penicillin fermentation, possible limitations in primary metabolism reside in the supply/regeneration of cofactors (NADPH) rather than in the supply of carbon precursors.     

Metabolic flux analysis of CHO cell metabolism in the late non‐growth phase

Chinese hamster ovary (CHO) cell cultures are commonly used for production of recombinant human therapeutic proteins. Often the goal of such a process is to separate the growth phase of the cells, from the non‐growth phase where ideally the cells are diverting resources to produce the protein of interest. Characterizing the way that the cells use nutrients in terms of metabolic fluxes as a function of culture conditions can provide a deeper understanding of the cell biology offering guidance for process improvements. To evaluate the fluxes, metabolic flux analysis of the CHO cell culture in the non‐growth phase was performed by a combination of steady‐state isotopomer balancing and stoichiometric modeling. Analysis of the glycolytic pathway and pentose phosphate pathway (PPP) indicated that almost all of the consumed glucose is diverted towards PPP with a high NADPH production; with even recycle from PPP to G6P in some cases. Almost all of the pyruvate produced from glycolysis entered the TCA cycle with little or no lactate production. Comparison of the non‐growth phase against previously reported fluxes from growth phase cultures indicated marked differences in the fluxes, in terms of the split between glycolysis and PPP, and also around the pyruvate node. Possible reasons for the high NADPH production are also discussed. Evaluation of the fluxes indicated that the medium strength, carbon dioxide level, and temperature with dissolved oxygen have statistically significant impacts on different nodes of the flux network. Biotechnol. Bioeng. 2011; 108:82–92. © 2010 Wiley Periodicals, Inc.     

An improved method for statistical analysis of metabolic flux analysis using isotopomer mapping matrices with analytical expressions

Analytical expressions were derived for calculating the sensitivities of isotopomer distribution vectors, the weighted output matrix with respect to the fluxes, and the covariance matrix for the metabolic flux analysis based on isotopomers mapping matrices (IMM). These expressions allow us to implement efficient statistical analysis, avoiding the time-consuming Monte Carlo techniques for estimating the confidence interval of the fluxes. The analytical expressions are also useful in implementing a faster design of experiment, which requires repetitive computation of the covariance matrix that is not straightforward to make in practice with the numerical techniques based on the conventional IMM. The proposed method was applied for analyzing the central carbon metabolism of the mixotrophically cultivated Synechocystis sp. PCC6803, and the confidence intervals of all its fluxes were computed based on the isotopomer distribution measured using NMR and GC-MS. It was found that the best feasible mixture for labeling experiment is 70% unlabeled, 10% [U-13C] and 20% [1,2-13C2] labeled glucose to obtain the most reliable metabolic fluxes.     

Simultaneous tracers and a unified model of positional and mass isotopomers for quantification of metabolic flux in liver

Computational models based on the metabolism of stable isotope tracers can yield valuable insight into the metabolic basis of disease. The complexity of these models is limited by the number of tracers and the ability to characterize tracer labeling in downstream metabolites. NMR spectroscopy is ideal for multiple tracer experiments since it precisely detects the position of tracer nuclei in molecules, but it lacks sensitivity for detecting low-concentration metabolites. GC-MS detects stable isotope mass enrichment in low-concentration metabolites, but lacks nuclei and positional specificity. We performed liver perfusions and in vivo infusions of ²H and ¹³C tracers, yielding complex glucose isotopomers that were assigned by NMR and fit to a newly developed metabolic model. Fluxes regressed from ²H and ¹³C NMR positional isotopomer enrichments served to validate GC-MS-based flux estimates obtained from the same experimental samples. NMR-derived fluxes were largely recapitulated by modeling the mass isotopomer distributions of six glucose fragment ions measured by GC-MS. Modest differences related to limited fragmentation coverage of glucose C1–C3 were identified, but fluxes such as gluconeogenesis, glycogenolysis, cataplerosis and TCA cycle flux were tightly correlated between the methods. Most importantly, modeling of GC-MS data could assign fluxes in primary mouse hepatocytes, an experiment that is impractical by ²H or ¹³C NMR.     

Comparative metabolic flux analysis of lysine-producing Corynebacterium glutamicum cultured on glucose or fructose

A comprehensive approach to (13)C tracer studies, labeling measurements by gas chromatography-mass spectrometry, metabolite balancing, and isotopomer modeling, was applied for comparative metabolic network analysis of lysine-producing Corynebacterium glutamicum on glucose or fructose. Significantly reduced yields of lysine and biomass and enhanced formation of dihydroxyacetone, glycerol, and lactate in comparison to those for glucose resulted on fructose. Metabolic flux analysis revealed drastic differences in intracellular flux depending on the carbon source applied. On fructose, flux through the pentose phosphate pathway (PPP) was only 14.4% of the total substrate uptake flux and therefore markedly decreased compared to that for glucose (62.0%). This result is due mainly to (i) the predominance of phosphoenolpyruvate-dependent phosphotransferase systems for fructose uptake (PTS(Fructose)) (92.3%), resulting in a major entry of fructose via fructose 1,6-bisphosphate, and (ii) the inactivity of fructose 1,6-bisphosphatase (0.0%). The uptake of fructose during flux via PTS(Mannose) was only 7.7%. In glucose-grown cells, the flux through pyruvate dehydrogenase (70.9%) was much less than that in fructose-grown cells (95.2%). Accordingly, flux through the tricarboxylic acid cycle was decreased on glucose. Normalized to that for glucose uptake, the supply of NADPH during flux was only 112.4% on fructose compared to 176.9% on glucose, which might explain the substantially lower lysine yield of C. glutamicum on fructose. Balancing NADPH levels even revealed an apparent deficiency of NADPH on fructose, which is probably overcome by in vivo activity of malic enzyme. Based on these results, potential targets could be identified for optimization of lysine production by C. glutamicum on fructose, involving (i) modification of flux through the two PTS for fructose uptake, (ii) amplification of fructose 1,6-bisphosphatase to increase flux through the PPP, and (iii) knockout of a not-yet-annotated gene encoding dihydroxyacetone phosphatase or kinase activity to suppress overflow metabolism. Statistical evaluation revealed high precision of the estimates of flux, so the observed differences for metabolic flux are clearly substrate specific.     

Metabolic flux analysis in Escherichia coli by integrating isotopic dynamic and isotopic stationary 13C labeling data

The novel concept of isotopic dynamic 13C metabolic flux analysis (ID-13C MFA) enables integrated analysis of isotopomer data from isotopic transient and/or isotopic stationary phase of a 13C labeling experiment, short-time experiments, and an extended range of applications of 13C MFA. In the presented work, an experimental and computational framework consisting of short-time 13C labeling, an integrated rapid sampling procedure, a LC-MS analytical method, numerical integration of the system of isotopomer differential equations, and estimation of metabolic fluxes was developed and applied to determine intracellular fluxes in glycolysis, pentose phosphate pathway (PPP), and citric acid cycle (TCA) in Escherichia coli grown in aerobic, glucose-limited chemostat culture at a dilution rate of D = 0.10 h(-1). Intracellular steady state concentrations were quantified for 12 metabolic intermediates. A total of 90 LC-MS mass isotopomers were quantified at sampling times t = 0, 91, 226, 346, 589 s and at isotopic stationary conditions. Isotopic stationarity was reached within 10 min in glycolytic and PPP metabolites. Consistent flux solutions were obtained by ID-13C MFA using isotopic dynamic and isotopic stationary 13C labeling data and by isotopic stationary 13C MFA (IS-13C MFA) using solely isotopic stationary data. It is demonstrated that integration of dynamic 13C labeling data increases the sensitivity of flux estimation, particularly at the glucose-6-phosphate branch point. The identified split ratio between glycolysis and PPP was 55%:44%. These results were confirmed by IS-13C MFA additionally using labeling data in proteinogenic amino acids (GC-MS) obtained after 5 h from sampled biomass.     

An optimized algorithm for flux estimation from isotopomer distribution in glucose metabolites

MOTIVATION: Analysis of the conversion of (13)C glucose within the metabolic network allows the evaluation of the biochemical fluxes in interconnecting metabolic pathways. Such analyses require solving hundreds of equations with respect to individual isotopomer concentrations, and this assumes applying special software even for constructing the equations. The algorithm, proposed by others could be improved. METHOD: A C-code linked to the program written in Mathematica (Wolfram Research Inc.), constructs and solves differential equations for all isotopomer concentrations, using the general enzyme characteristics (K(m), equilibrium constant, etc.). This code uses innovative algorithm of determination for the isotopomers-products, thus essentially decreasing the computation time. Feasible metabolic fluxes are provided by the parameters of enzyme kinetics found from the data fitting. RESULTS: The software effectively evaluates metabolic fluxes based on the measured isotopomer distribution, as was illustrated by the analysis of glycolysis and pentose phosphate cycle. The mechanism of transketolase and transaldolase catalysis was shown to induce a specific kind of isotopomer re-distribution, which, despite the significance of its effect, usually is not taken into account. AVAILABILITY: The software could be freely downloaded from the site: http://bq.ub.es/bioqint/label_distribution/.     

Metabolic flux analysis of postburn hepatic hypermetabolism

The hepatic response to severe injury is characterized by a marked upregulation of glucose, fatty acid, and amino acid turnover, which, if persistent, predisposes the patient to progressive organ dysfunction. To study the effect of injury on liver intermediary metabolism, metabolic flux analysis was applied to isolated perfused livers of burned and sham-burned rats. Intracellular fluxes were calculated using metabolite measurements and a stoichiometric balance model. Significant flux increases were found for multiple pathways, including mitochondrial electron transport, the TCA and urea cycles, gluconeogenesis, and pentose phosphate pathway (PPP). The burn-induced increase in gluconeogenesis did not significantly increase glucose output. Instead, glucose-6-phosphate was diverted into the PPP. These changes were paralleled by increases in glucose-6-phosphate dehydrogenase (G6PDH) and glutathione reductase (GR) activities. Given that G6PDH and GR are the most significant NADPH producers and consumers in the liver, respectively, and that GR is responsible for recycling the free radical scavenger glutathione, these data are consistent with the notion that hepatic metabolic changes are in part due to the induction of liver antioxidant defenses.     

Cytosolic NADPH metabolism in penicillin-G producing and non-producing chemostat cultures of Penicillium chrysogenum

This study addresses the relation between NADPH supply and penicillin synthesis, by comparing the flux through the oxidative branch of the pentose phosphate pathway (PPP; the main source of cytosolic NADPH) in penicillin-G producing and non-producing chemostat cultures of Penicillium chrysogenum. The fluxes through the oxidative part of the PPP were determined using the recently introduced gluconate-tracer method. Significantly higher oxidative PPP fluxes were observed in penicillin-G producing chemostat cultures, indicating that penicillin production puts a major burden on the supply of cytosolic NADPH. To our knowledge this is the first time direct experimental proof is presented for the causal relationship between penicillin production and NADPH supply. Additional insight in the metabolism of P. chrysogenum was obtained by comparing the PPP fluxes from the gluconate-tracer experiment to oxidative PPP fluxes derived via metabolic flux analysis, using different assumptions for the stoichiometry of NADPH consumption and production.     

Gas chromatography-mass spectrometry analysis of 13C labeling in sugars for metabolic flux analysis

Metabolic flux analysis, using 13C labeled substrates, has become a powerful methodology for quantifying intracellular fluxes. Most often, analysis is restricted to nuclear magnetic resonance or mass spectrometry measurement of 13C label incorporation into protein amino acids. However, amino acid isotopomer distribution insufficiently covers the entire network of central metabolism, especially in plant cells with highly compartmented metabolism, and analysis of other metabolites is required. Analysis of label in saccharides provides complementary data to better define fluxes around hexose, pentose, and triose phosphate pools. Here, we propose a gas chromatography-mass spectrometry (GC-MS) method to analyze 13C labeling in glucose and fructose moieties of sucrose, free glucose, fructose, maltose, inositol, and starch. Our results show that saccharide labeling for isotopomer quantification is better analyzed by chemical ionization than by electron ionization. The structure of the generated fragments was simulated and validated using labeled standards. The method is illustrated by analysis of saccharides extracted from developing rapeseed (Brassica napus L.) embryos. It is shown that glucose 6-phosphate isomerase and plastidial glucose 6-phosphate transport reactions are not at equilibrium, and light is shed on the pathways leading to fructose, maltose, and inositol synthesis.     

GC-MS analysis of amino acids rapidly provides rich information for isotopomer balancing

Gas chromatography-mass spectrometry (GC-MS) is a rapid method that provides rich information on isotopomer distributions for metabolic flux analysis. First, we established a fast and reliable experimental protocol for GC-MS analysis of amino acids from total biomass hydrolyzates, and common experimental pitfalls are discussed. Second, a suitable interface for the use of mass isotopomer data is presented. For this purpose, a general, matrix-based correction procedure that accounts for naturally occurring isotopes is introduced. Simulated and experimentally determined mass distributions of unlabeled amino acids showed a deviation of less than 3% for 89% of the analyzed fragments. Third, to investigate general properties of GC-MS-based isotopomer balancing, altered flux ratios through glycolysis and pentose phosphate pathway and/or exchange fluxes were simulated. Different fluxes were found to exert specific and significant influence on the mass isotopomer distributions, thus indicating that GC-MS data contain valuable information for metabolic flux analysis. Fourth, comparison of different methods revealed that GC-MS analysis provides the largest number of independent constraints on amino acid isotopomer abundance, followed by correlation spectroscopy and fractional enrichment analysis by nuclear magnetic resonance.     

High-throughput metabolic flux analysis based on gas chromatography-mass spectrometry derived 13C constraints

13C-constrained flux balancing analysis based on gas chromatography-mass spectrometry data is presented here as a simple and robust method for the estimation of intracellular carbon fluxes. In this approach, the underdetermined system of metabolite balances deduced from stoichiometric relations and measured extracellular rates is complemented with 13C constraints from metabolic flux ratio analysis. Fluxes in central carbon metabolism of exponentially growing Escherichia coli were estimated by 13C-constrained flux balancing from three different 13C-labeled glucose experiments. The best resolution of the network was achieved using 13C constraints derived from [U-13C]glucose and [1-13C]glucose experiments. The corresponding flux estimate was in excellent agreement with a solution that was independently obtained with a comprehensive isotopomer model. This new methodology was also demonstrated to faithfully capture the intracellular flux distribution in E. coli shake flasks and 1-ml deep-well microtiter plates. Due to its simplicity, speed, and robustness, 13C-constrained metabolic flux balancing is promising for routine and high-throughput analysis on a miniaturized scale.     

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 ¹³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. © 2010 American Institute of Chemical Engineers Biotechnol. Prog., 2010     

Ion-trap mass spectrometry used in combination with gas chromatography for high-resolution metabolic flux determination

Metabolic fluxes provide a detailed metric of the cellular metabolic state. Fluxes are estimated indirectly from available measurements, and various methods have been developed for this purpose. Of particular interest are methods that make use of stable isotopic tracers because they enable flux estimation at a fine resolution. In this report, we present a protocol for the use of ion-trap mass spectrometry (MS) in combination with gas chromatography to measure the mass isotopomer distribution of biomass hydrolysates. At physiological steady-state, these measurements directly reflect the isotopic tracer distribution in the amino acid central carbon metabolism precursors. Because the extent to which a metabolic flux network can be accurately resolved strongly depends on the reliability and precision of the MS measurements and, in light of the current need for quantitative high-throughput biological analysis at the microscale, we discuss every step of the measurement process, indicate possible sources of error, and suggest solutions to avoid them. Potential advantages to using ion-trap versus quadrupole MS are also addressed. The final protocol requires 0.5 mg of dry biomass to detect the mass isotopomer distribution of 2-4 fragments of 13 amino acids, with a relative variance less than 1% for the most abundant peaks.     

Isotopically nonstationary 13C metabolic flux analysis in resting and activated human platelets

Platelet metabolism is linked to platelet hyper- and hypoactivity in numerous human diseases. Developing a detailed understanding of the link between metabolic shifts and platelet activation state is integral to improving human health. Here, we show the first application of isotopically nonstationary ¹³C metabolic flux analysis to quantitatively measure carbon fluxes in both resting and thrombin activated platelets. Metabolic flux analysis results show that resting platelets primarily metabolize glucose to lactate via glycolysis, while acetate is oxidized to fuel the tricarboxylic acid cycle. Upon activation with thrombin, a potent platelet agonist, platelets increase their uptake of glucose 3-fold. This results in an absolute increase in flux throughout central metabolism, but when compared to resting platelets they redistribute carbon dramatically. Activated platelets decrease relative flux to the oxidative pentose phosphate pathway and TCA cycle from glucose and increase relative flux to lactate. These results provide the first report of reaction-level carbon fluxes in platelets and allow us to distinguish metabolic fluxes with much higher resolution than previous studies.     

Current status of 13C-metabolic flux analysis and future perspectives

Current ¹³C-metabolic flux analysis methods were reviewed as well as the weakness of the conventional metabolic flux analysis without ¹³C-labeled experiments. Although it has been recognized that ¹³C-labeling technique is powerful in estimating the metabolic fluxes, and the program-based flux analysis is necessary, one may not be confident with the result obtained without experiences and exhaustive trial and errors in practice due to its black box nature. In the present article, we call attention to the importance of investigating the relationships between fluxes and isotopomer or mass isotopomer distributions to understand the mechanism of generating specific isotopomers. Then, the experimental design for the preferred mixture of the specific ¹³C-labeled substrate was discussed. The effect of the reversibility in the bidirectional flux on the isotopomer distribution was also mentioned, and it was shown why the reliability of the bidirectional fluxes becomes lower. Moreover, by noting that recent development of measurement techniques enables us to measure the isotopomer patterns of intracellular metabolites instead of proteinogenic amino acids, it is mentioned that this enables us to estimate the flux changes during time-variant batch culture. Some future perspectives are discussed in relation to the integration of different levels of information in the cell.