Quantitative analysis of metabolic fluxes in Escherichia coli, using two-dimensional NMR spectroscopy and complete isotopomer models.

Complete isotopomer models that simulate distribution of label in 13C tracer experiments are applied to the quantification of metabolic fluxes in the primary carbon metabolism of E. coli under aerobic and anaerobic conditions. The concept of isotopomer mapping matrices (IMMs) is used to simplify the formulation of isotopomer mass balances by expressing all isotopomer mass balances of a metabolite pool in a single matrix equation. A numerically stable method to calculate the steady-state isotopomer distribution in metabolic networks in introduced. Net values of intracellular fluxes and the degree of reversibility of enzymatic steps are estimated by minimization of the deviations between experimental and simulated measurements. The metabolic model applied includes the Embden-Meyerhof-Parnas and the pentose phosphate pathway, the tricarboxylic acid cycle, anaplerotic reaction sequences and pathways involved in amino acid synthesis. The study clearly demonstrates the value of complete isotopomer models for maximizing the information obtainable from 13C tracer experiments. The approach applied here offers a completely general and comprehensive analysis of carbon tracer experiments where any set of experimental data on the labeling state and extracellular fluxes can be used for the quantification of metabolic fluxes in complex metabolic networks.     

Metabolic isotopomer labeling systems. Part II: structural flux identifiability analysis

Metabolic flux analysis using carbon labeling experiments (CLEs) is an important tool in metabolic engineering where the intracellular fluxes have to be computed from the measured extracellular fluxes and the partially measured distribution of 13C labeling within the intracellular metabolite pools. The relation between unknown fluxes and measurements is described by an isotopomer labeling system (ILS) (see Part I [Math. Biosci. 169 (2001) 173]). Part II deals with the structural flux identifiability of measured ILSs in the steady state. The central question is whether the measured data contains sufficient information to determine the unknown intracellular fluxes. This question has to be decided a priori, i.e. before the CLE is carried out. In structural identifiability analysis the measurements are assumed to be noise-free. A general theory of structural flux identifiability for measured ILSs is presented and several algorithms are developed to solve the identifiability problem. In the particular case of maximal measurement information, a symbolical algorithm is presented that decides the identifiability question by means of linear methods. Several upper bounds of the number of identifiable fluxes are derived, and the influence of the chosen inputs is evaluated. By introducing integer arithmetic this algorithm can even be applied to large networks. For the general case of arbitrary measurement information, identifiability is decided by a local criterion. A new algorithm based on integer arithmetic enables an a priori local identifiability analysis to be performed for networks of arbitrary size. All algorithms have been implemented and flux identifiability is investigated for the network of the central metabolic pathways of a microorganism. Moreover, several small examples are worked out to illustrate the influence of input metabolite labeling and the paradox of information loss due to network simplification.     

Quantification of intracellular metabolic fluxes from fractional enrichment and 13C-13C coupling constraints on the isotopomer distribution in labeled biomass components

A method for the quantification of intracellular metabolic flux distributions from steady-state mass balance constraints and from the constraints posed by the measured 13C labeling state of biomass components is presented. Two-dimensional NMR spectroscopy is used to analyze the labeling state of cell protein hydrolysate and cell wall components. No separation of the biomass hydrolysate is required to measure the degree of 13C-13C coupling and the fractional 13C enrichment in various carbon atom positions. A mixture of [1-13C]glucose and uniformly labeled [13C6]glucose is applied to make fractional 13C enrichment data and measurements of the degree of 13C-13C coupling informative with respect to the intracellular flux distribution. Simulation models that calculate the complete isotopomer distribution in biomass components on the basis of isotopomer mapping matrices are used for the estimation of intracellular fluxes by least-squares minimization. The statistical quality of the estimated intracellular flux distributions is assessed by Monte Carlo methods. Principal component analysis is performed on the outcome of the Monte Carlo procedure to identify groups of fluxes that contribute major parts to the total variance in the multiple flux estimations. The methods described are applied to a steady-state culture of a glucoamylase-producing recombinant Aspergillus niger strain.     

Metabolic flux analysis with a comprehensive isotopomer model in Bacillus subtilis

Fluxes in central carbon metabolism of a genetically engineered, riboflavin-producing Bacillus subtilis strain were investigated in glucose-limited chemostat cultures at low (0.11 h(-1)) and high (0.44 h(-1)) dilution rates. Using a mixture of 10% [U-(13)C] and 90% glucose labeled at natural abundance, (13)C-labeling experiments were carried out to provide additional information for metabolic flux balancing. The resulting labeling pattern in the proteinogenic amino acids were analyzed by two-dimensional [(13)C, (1)H] nuclear magnetic resonance (NMR) spectroscopy. To account rigorously for all available data from these experiments, we developed a comprehensive isotopomer model of B. subtilis central metabolism. Using this model, intracellular carbon net and exchange fluxes were estimated on the basis of validated physiological data and biomass composition in combination with 2D NMR data from 45 individual carbon atom spectra in the amino acids. Glucose catabolism proceeded primarily via glycolysis but pentose phosphate pathway fluxes increased with increasing growth rate. Moreover, significant back fluxes from the TCA cycle to the lower part of glycolysis via the gluconeogenic PEP carboxykinase were detected. The malic enzyme reaction, in contrast, was found to be inactive. A thorough statistical analysis was performed to prove the reliability of the isotopomer balance model and the obtained results. Specifically, a chi(2) test was applied to validate the model and the chi-square criterion was used to explore the sensitivity of model predictions to the experimental data.     

Using isotopomer path tracing to quantify metabolic fluxes in pathway models containing reversible reactions

As a more complete picture of the genetic and enzymatic composition of cells becomes available, there is a growing need to describe how cellular regulatory elements interact with the cellular environment to affect cell physiology. One means for describing intracellular regulatory mechanisms is concurrent measurement of multiple metabolic pathways and their interactions by metabolic flux analysis. Flux of carbon through a metabolic pathway responds to all cellular regulatory systems, including changes in enzyme and substrate concentrations, enzyme activation or inhibition, and ultimately genetic control. The extent to which metabolic flux analysis can describe cellular physiology depends on the number of pathways in the model and the quality of the data. Intracellular information is obtainable from isotopic tracer experiments, the most extensive being the determination of the isotopomer distribution, or specific labeling pattern, of intracellular metabolites. We present a rapid and novel solution method that determines the flux of carbon through complex pathway models using isotopomer data. This time-consuming problem was solved with the introduction of isotopomer path tracing, which drastically reduces the number of isotopomer variables to the number of isotopomers observed experimentally. We propose a partitioned solution method that takes advantage of the nearly linear relationship between fluxes and isotopomers. Whereas the stoichiometric matrix and the isotopomer matrix are invertible, simulated annealing and the Newton-Raphson method are used for the nonlinear components. Reversible reactions are described by a new parameter, the association factor, which scales hyperbolically with the rate of metabolite exchange. Automating the solution method permits a variety of models to be compared, thus enhancing the accuracy of results. A simplified example that contains all of the complexities of a comprehensive pathway model is presented. Copyright John Wiley & Sons, Inc.     

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.     

Determination of full 13C isotopomer distributions for metabolic flux analysis using heteronuclear spin echo difference NMR spectroscopy

13C-isotopomer labeling experiments play an increasingly important role in the analysis of intracellular metabolic fluxes for genetic engineering purposes. 13C NMR spectroscopy is a key technique in the experimental determination of isotopomer distributions. However, only subsets of isotopomers can be quantitated using this technique due to redundancies in the scalar coupling patterns and due to invisibility of the 12C isotope in NMR. Therefore, we developed and describe in this paper a 1H NMR spectroscopy method that allows to determine the complete isotopomer distribution in metabolites having a backbone consisting of up to at least four carbons. The proposed pulse sequences employ up to three alternately applied frequency-selective inversion pulses in the 13C channel. In a first application study, the complete isotopomer distribution of aspartate isolated from [1-13C]ethanol-grown Ashbya gossypii was determined. A tentative model of the central metabolism of this organism was constructed and used for metabolic flux analysis. The aspartate isotopomer NMR data played a key role in the successful determination of the flux through the peroxisomal glyoxylate pathway. The new NMR method can be highly instrumental in generating the data upon which isotopomer labeling experiments for flux analysis, that are becoming increasingly important, are based.     

Isotopomer analysis using GC-MS

Knowledge of the complete isotopomer distribution represents the ultimate amount of information on the labeling pattern of a metabolite. One technique for measuring the isotopomer distributions is the analysis of the multiplet intensities arising from the 13C-13C couplings in NMR spectroscopy. While this technique has proven to be very valuable in the elucidation of labeling patterns of C2 and C3 units of various amino acids, fragments larger than C3 are very difficult to measure. Another technique, GC-MS, offers a unique possibility of analyzing fragments larger than C3 and GC-MS is therefore able to give information which is complementary to the information that can be obtained from NMR spectroscopy. In this work we have developed fast, simple, and robust GC-MS methods that can be used to gain information on the labeling patterns of the amino acids in a crude biomass hydrolysate. It is shown that a combination of information obtained from these analyses and information from the NMR spectroscopy is able to yield a much more complete picture of the isotopomer distributions of the amino acids than any of the two techniques alone. The GC-MS method was used for analyzing the labeling patterns of amino acids from a batch cultivation of Penicillium chrysogenum grown on fully labeled glucose. The data from this analysis showed no signs of any significant carbon isotope effects, and the measurements can therefore be used without corrections for metabolic flux analysis.     

Metabolite and isotopomer balancing in the analysis of metabolic cycles: I. Theory

Proper analysis of label distribution in metabolic pathway intermediates is critical for correct interpretation of experimental data and strategic experimental design. While, for example, 13C nuclear magnetic resonance (NMR) spectroscopy is usually limited to the measurement of degrees of 13C enrichment, more information about metabolic fluxes can be extracted from the fine structure of NMR spectra, or molecular weight distributions of isotopomers of metabolic intermediates (measured by gas chromatography-mass spectrometry). For this purpose, rigorous accounting for the contribution of all pathways to label distribution is required, especially contributions resulting from multiple turns of metabolic cycles. In this paper we present a mathematical model developed to analyze isotopomer distributions of tricarboxylic acid cycle (TCA) intermediates following the administration of 13C (or 14C) labeled substrates. The theory presented provides the basis to analyze 13C NMR spectra and molecular weight distributions of metabolites. In a companion paper (Park et al., 1999), the theory is applied to the analysis of several cases of biological significance. Copyright 1999 John Wiley & Sons, 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.     

13C isotopomer analysis of glutamate by J-resolved heteronuclear single quantum coherence spectroscopy

13C NMR isotopomer analysis is a powerful method for measuring metabolic fluxes through pathways intersecting in the tricarboxylic acid cycle. However, the inherent insensitivity of 13C NMR spectroscopy makes application of isotopomer analysis to small tissue samples (mouse tissue, human biopsies, or cells grown in tissue culture) problematic. (1)H NMR is intrinsically more sensitive than 13C NMR and can potentially supply the same information via indirect detection of 13C providing that isotopomer information can be preserved. We report here the use of J-resolved HSQC (J-HSQC) for 13C isotopomer analysis of tissue samples. We show that J-HSQC reports isotopomer multiplet patterns identical to those reported by direct 13C detection but with improved sensitivity.     

Quantification of compartmented metabolic fluxes in developing soybean embryos by employing biosynthetically directed fractional (13)C labeling, two-dimensional [(13)C, (1)H] nuclear magnetic resonance, and comprehensive isotopomer balancing

Metabolic flux quantification in plants is instrumental in the detailed understanding of metabolism but is difficult to perform on a systemic level. Toward this aim, we report the development and application of a computer-aided metabolic flux analysis tool that enables the concurrent evaluation of fluxes in several primary metabolic pathways. Labeling experiments were performed by feeding a mixture of U-(13)C Suc, naturally abundant Suc, and Gln to developing soybean (Glycine max) embryos. Two-dimensional [(13)C, (1)H] NMR spectra of seed storage protein and starch hydrolysates were acquired and yielded a labeling data set consisting of 155 (13)C isotopomer abundances. We developed a computer program to automatically calculate fluxes from this data. This program accepts a user-defined metabolic network model and incorporates recent mathematical advances toward accurate and efficient flux evaluation. Fluxes were calculated and statistical analysis was performed to obtain sds. A high flux was found through the oxidative pentose phosphate pathway (19.99 +/- 4.39 micromol d(-1) cotyledon(-1), or 104.2 carbon mol +/- 23.0 carbon mol per 100 carbon mol of Suc uptake). Separate transketolase and transaldolase fluxes could be distinguished in the plastid and the cytosol, and those in the plastid were found to be at least 6-fold higher. The backflux from triose to hexose phosphate was also found to be substantial in the plastid (21.72 +/- 5.00 micromol d(-1) cotyledon(-1), or 113.2 carbon mol +/-26.0 carbon mol per 100 carbon mol of Suc uptake). Forward and backward directions of anaplerotic fluxes could be distinguished. The glyoxylate shunt flux was found to be negligible. Such a generic flux analysis tool can serve as a quantitative tool for metabolic studies and phenotype comparisons and can be extended to other plant systems.     

Planning optimal measurements of isotopomer distributions for estimation of metabolic fluxes

MOTIVATION: Flux estimation using isotopomer information of metabolites is currently the most reliable method to obtain quantitative estimates of the activity of metabolic pathways. However, the development of isotopomer measurement techniques for intermediate metabolites is a demanding task. Careful planning of isotopomer measurements is thus needed to maximize the available flux information while minimizing the experimental effort. RESULTS: In this paper we study the question of finding the smallest subset of metabolites to measure that ensure the same level of isotopomer information as the measurement of every metabolite in the metabolic network. We study the computational complexity of this optimization problem in the case of the so-called positional enrichment data, give methods for obtaining exact and fast approximate solutions, and evaluate empirically the efficacy of the proposed methods by analyzing a metabolic network that models the central carbon metabolism of Saccharomyces cerevisiae.     

INCA 2.0: A tool for integrated,dynamic modeling of NMR- and MS-based isotopomer measurements and rigorous metabolic flux analysis

Metabolic flux analysis (MFA) combines experimental measurements and computational modeling to determine biochemical reaction rates in live biological systems. Advancements in analytical instrumentation, such as nuclear magnetic resonance (NMR) spectroscopy and mass spectrometry (MS), have facilitated chemical separation and quantification of isotopically enriched metabolites. However, no software packages have been previously described that can integrate isotopomer measurements from both MS and NMR analytical platforms and have the flexibility to estimate metabolic fluxes from either isotopic steady-state or dynamic labeling experiments. By applying physiologically relevant cardiac and hepatic metabolic models to assess NMR isotopomer measurements, we herein test and validate new modeling capabilities of our enhanced flux analysis software tool, INCA 2.0. We demonstrate that INCA 2.0 can simulate and regress steady-state ¹³C NMR datasets from perfused hearts with an accuracy comparable to other established flux assessment tools. Furthermore, by simulating the infusion of three different ¹³C acetate tracers, we show that MFA based on dynamic ¹³C NMR measurements can more precisely resolve cardiac fluxes compared to isotopically steady-state flux analysis. Finally, we show that estimation of hepatic fluxes using combined ¹³C NMR and MS datasets improves the precision of estimated fluxes by up to 50%. Overall, our results illustrate how the recently added NMR data modeling capabilities of INCA 2.0 can enable entirely new experimental designs that lead to improved flux resolution and can be applied to a wide range of biological systems and measurement time courses.     

Profiling of central metabolism in human cancer cells by two-dimensional NMR, GC-MS analysis, and isotopomer modeling

Tracking metabolic profiles has the potential to reveal crucial enzymatic steps that could be targeted in the drug discovery process. It is of special importance for various types of cancer known to be associated with substantial rewiring of metabolic networks. Here we introduce an integrated approach for the analysis of metabolome that allows us to simultaneously assess pathway activities (fluxes) and concentrations of a large number of the key components involved in central metabolism of human cells. This is accomplished by in vivo labeling with [U-¹³C]glucose followed by two-dimensional nuclear magnetic resonance (NMR) spectroscopy and gas chromatography-mass spectrometry (GC-MS) analysis. A comprehensive isotopomer model was developed, which enabled us to compare fluxes through the key central metabolic pathways including glycolysis, pentose phosphate pathway, tricarboxylic acid cycle, anaplerotic reactions, and biosynthetic pathways of fatty acids and amino acids. The validity and strength of this approach is illustrated by its application to a number of perturbations to breast cancer cells, including exposure to hypoxia, drug treatment, and tumor progression. We observed significant differences in the activities of specific metabolic pathways resulting from these perturbations and providing new mechanistic insights. Based on these findings we conclude that the developed metabolomic approach constitutes a promising analytical tool for revealing specific metabolic phenotypes in a variety of cell types and pathological conditions. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Chen Yang and Adam D. Richardson contributed equally to this work.     

一种中间代谢途径代谢通量的定量分析方法

以¹³C标记的碳源,用二维核磁共振技术(¹H-¹³C,HMQC)测定代谢中产生的氨基酸标记模式,研究对中间代谢途径胞内代谢通量进行定量分析的方法.通过开发软件包,改进同位素分布的数学模型,提出了反应映射矩阵(RMM)等概念.由简化算法,提高程序的执行效率,建立了定量分析胞内代谢通量的平台.代谢模型涉及了糖酵解途径、磷酸戊糖途径、三羧酸循环、几种回补反应、发酵途径和氨基酸合成途径.     

Application of MALDI-TOF MS to lysine-producing Corynebacterium glutamicum: a novel approach for metabolic flux analysis.

In the present work, a novel comprehensive approach of (13)C-tracer studies with labeling measurements by MALDI-TOF MS, and metabolite balancing was developed to elucidate key fluxes in the central metabolism of lysine producing Corynebacterium glutamicum during batch culture. MALDI-TOF MS methods established allow the direct quantification of labeling patterns of low molecular mass Corynebacterium products from 1 microL of diluted culture supernatant. A mathematical model of the central Corynebacterium metabolism was developed, that describes the carbon transfer through the network via matrix calculations in a generally applicable way and calculates steady state mass isotopomer distributions of the involved metabolites. The model was applied for both experimental planning of tracer experiments and parameter estimation. Metabolic fluxes were calculated from stoichiometric data and from selected mass intensity ratios of lysine, alanine, and trehalose measured by MALDI-TOF MS in tracer experiments either with 1-(13)C glucose or with mixtures of (13)C6/(12)C6 glucose. During the phase of maximum lysine production C. glutamicum ATCC 21253 exhibited high relative fluxes into the pentose phosphate pathway of 71%, a highly reversible glucose-6-phosphate isomerase, significant backfluxes from the tricarboxylic acid cycle to the pyruvate node consuming the lysine precursor oxaloacetate, 36% net flux of anaplerotic carboxylation and 63% contribution of the dehydrogenase branch in the lysine biosynthetic pathway. Due to the straightforward and simple measurements of selected labeling patterns by MALDI-TOF MS sensitively reflecting the flux parameters of interest, the presented approach has an excellent potential to extend metabolic flux analysis from single experiments with enormous experimental effort to a broadly applied technique.     

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.     

Determination of metabolic fluxes in a non-steady-state system

Estimation of fluxes through metabolic networks from redistribution patterns of (13)C has become a well developed technique in recent years. However, the approach is currently limited to systems at metabolic steady-state; dynamic changes in metabolic fluxes cannot be assessed. This is a major impediment to understanding the behaviour of metabolic networks, because steady-state is not always experimentally achievable and a great deal of information about the control hierarchy of the network can be derived from the analysis of flux dynamics. To address this issue, we have developed a method for estimating non-steady-state fluxes based on the mass-balance of mass isotopomers. This approach allows multiple mass-balance equations to be written for the change in labelling of a given metabolite pool and thereby permits over-determination of fluxes. We demonstrate how linear regression methods can be used to estimate non-steady-state fluxes from these mass balance equations. The approach can be used to calculate fluxes from both mass isotopomer and positional isotopomer labelling information and thus has general applicability to data generated from common spectrometry- or NMR-based analytical platforms. The approach is applied to a GC-MS time-series dataset of (13)C-labelling of metabolites in a heterotrophic Arabidopsis cell suspension culture. Threonine biosynthesis is used to demonstrate that non-steady-state fluxes can be successfully estimated from such data while organic acid metabolism is used to highlight some common issues that can complicate flux estimation. These include multiple pools of the same metabolite that label at different rates and carbon skeleton rearrangements.     

Bidirectional reaction steps in metabolic networks: III. Explicit solution and analysis of isotopomer labeling systems

The last few years have brought tremendous progress in experimental methods for metabolic flux determination by carbon-labeling experiments. A significant enlargement of the available measurement data set has been achieved, especially when isotopomer fractions within intracellular metabolite pools are quantitated. This information can be used to improve the statistical quality of flux estimates. Furthermore, several assumptions on bidirectional intracellular reaction steps that were hitherto indispensable may now become obsolete. To make full use of the complete measurement information a general mathematical model for isotopomer systems is established in this contribution. Then, by introducing the important new concept of cumomers and cumomer fractions, it is shown that the arising nonlinear isotopomer balance equations can be solved analytically in all cases. In particular, the solution of the metabolite flux balances and the positional carbon-labeling balances presented in part I of this series turn out to be just the first two steps of the general solution procedure for isotopomer balances. A detailed analysis of the isotopomer network structure then opens up new insights into the intrinsic structure of isotopomer systems. In particular, it turns out that isotopomer systems are not as complex as they appear at first glance. This enables some far-reaching conclusions to be drawn on the information potential of isotopomer experiments with respect to flux identification. Finally, some illustrative examples are examined to show that an information increase is not guaranteed when isotopomer measurements are used in addition to positional enrichment data.