Metabolic isotopomer labeling systems. Part I: global dynamic behavior.

In the last few years metabolic flux analysis (MFA) using carbon labeling experiments (CLE) has become a major diagnostic tool in metabolic engineering. The mathematical centerpiece of MFA is the solution of isotopomer labeling systems (ILS). An ILS is a high-dimensional nonlinear differential equation system that describes the distribution of isotopomers over a metabolic network during a carbon labeling experiment. This contribution presents a global analysis of the dynamic behavior of general ILSs. It is proven that an ILS is globally stable under very weak conditions that are always satisfied in practice. In particular it is shown that in some sense ILSs are a nonlinear extension to the classical theory of compartmental systems. The central stability condition for compartmental systems, i.e., the non-existence of traps in linear compartmental networks, is also the major stability condition for ILSs. As an important side result of the proof, it is shown that ILSs can be transformed to a cascade of linear systems with time-dependent inhomogeneous terms. This cascade structure has considerable consequences for the development of efficient numerical algorithms for the solution of ILSs and thus for MFA.

[1]  W. Wiechert,et al.  Bidirectional reaction steps in metabolic networks: I. Modeling and simulation of carbon isotope labeling experiments. , 1997, Biotechnology and bioengineering.

[2]  H Sahm,et al.  Determination of the fluxes in the central metabolism of Corynebacterium glutamicum by nuclear magnetic resonance spectroscopy combined with metabolite balancing , 1996, Biotechnology and bioengineering.

[3]  W. Wiechert,et al.  Bidirectional reaction steps in metabolic networks: III. Explicit solution and analysis of isotopomer labeling systems. , 1999, Biotechnology and bioengineering.

[4]  J. Blum,et al.  Quantitative analysis of intermediary metabolism in hepatocytes incubated in the presence and absence of glucagon with a substrate mixture containing glucose, ribose, fructose, alanine and acetate. , 1985, The Biochemical journal.

[5]  David H. Anderson.  Compartmental Modeling and Tracer Kinetics , 1983 .

[6]  W. Wiechert,et al.  Bidirectional reaction steps in metabolic networks: II. Flux estimation and statistical analysis. , 1997, Biotechnology and bioengineering.

[7]  E M Chance,et al.  Calculation of Absolute Metabolic Flux and the Elucidation of the Pathways of Glutamate Labeling in Perfused Rat Heart by 13C NMR Spectroscopy and Nonlinear Least Squares Analysis (*) , 1995, The Journal of Biological Chemistry.

[8]  G. Stephanopoulos,et al.  Intracellular flux analysis in hybridomas using mass balances and in vitro 13C nmr , 1995, Biotechnology and bioengineering.

[9]  H Sahm,et al.  Response of the central metabolism in Corynebacterium glutamicum to the use of an NADH-dependent glutamate dehydrogenase. , 1999, Metabolic engineering.

[10]  J. Bailey,et al.  Toward a science of metabolic engineering , 1991, Science.

[11]  W Wiechert,et al.  In vivo stationary flux analysis by 13C labeling experiments. , 1996, Advances in biochemical engineering/biotechnology.

[12]  J. Villadsen,et al.  Modeling isotopomer distributions in biochemical networks using isotopomer mapping matrices. , 1997, Biotechnology and bioengineering.

[13]  Stephen Barnett,et al.  Introduction to Mathematical Control Theory , 1975 .

[14]  H. Schwetlick,et al.  Stoer, J. / Bulirsch, R., Einführung in die Numerische Mathematik II, IX, 286 S., 1973. DM 14,80, US $ 5.50. Berlin-Heidelberg-New York. Springer-Verlag , 1978 .

[15]  W Wiechert,et al.  Bidirectional reaction steps in metabolic networks: IV. Optimal design of isotopomer labeling experiments. , 1999, Biotechnology and bioengineering.

[16]  J. Nielsen,et al.  Quantitative analysis of metabolic fluxes in Escherichia coli, using two-dimensional NMR spectroscopy and complete isotopomer models. , 1999, Journal of biotechnology.

[17]  Keith R. Yamamoto,et al.  Biological Regulation and Development , 1982, Springer US.

[18]  A. Sherry,et al.  Evaluation of carbon flux and substrate selection through alternate pathways involving the citric acid cycle of the heart by 13C NMR spectroscopy. , 1988, The Journal of biological chemistry.

[19]  D. Jordan,et al.  Nonlinear Ordinary Differential Equations: An Introduction for Scientists and Engineers , 1979 .

[20]  J. Blum,et al.  On the Analysis of Metabolic Networks , 1982 .

[21]  D. Rittenberg,et al.  DEUTERIUM AS AN INDICATOR IN THE STUDY OF INTERMEDIARY METABOLISM. , 1935, Science.

[22]  D. Fife,et al.  Which linear compartmental systems contain traps , 1972 .

[23]  T Szyperski,et al.  13C-NMR, MS and metabolic flux balancing in biotechnology research , 1998, Quarterly Reviews of Biophysics.

[24]  D Weuster-Botz,et al.  Development and application of a membrane cyclone reactor for in vivo NMR spectroscopy with high microbial cell densities. , 2000, Biotechnology and bioengineering.

[25]  A. D. de Graaf,et al.  Flux partitioning in the split pathway of lysine synthesis in Corynebacterium glutamicum. Quantification by 13C- and 1H-NMR spectroscopy. , 1993, European journal of biochemistry.