Determination of confidence intervals of metabolic fluxes estimated from stable isotope measurements.

Metabolic fluxes, estimated from stable isotope studies, provide a key to quantifying physiology in fields ranging from metabolic engineering to the analysis of human metabolic diseases. A serious drawback of the flux estimation method in current use is that it does not produce confidence limits for the estimated fluxes. Without this information it is difficult to interpret flux results and expand the physiological significance of flux studies. To address this shortcoming we derived analytical expressions of flux sensitivities with respect to isotope measurements and measurement errors. These tools allow the determination of local statistical properties of fluxes and relative importance of measurements. Furthermore, we developed an efficient algorithm to determine accurate flux confidence intervals and demonstrated that confidence intervals obtained with this method closely approximate true flux uncertainty. In contrast, confidence intervals approximated from local estimates of standard deviations are inappropriate due to inherent system nonlinearities. We applied these methods to analyze the statistical significance and confidence of estimated gluconeogenesis fluxes from human studies with [U-13C]glucose as tracer and found true limits for flux estimation in specific human isotopic protocols.

[1]  Marc K Hellerstein,et al.  In vivo measurement of fluxes through metabolic pathways: the missing link in functional genomics and pharmaceutical research. , 2003, Annual review of nutrition.

[2]  J. Katz,et al.  Glucose production, recycling, and gluconeogenesis in normals and diabetics: a mass isotopomer [U-13C]glucose study. , 1996, The American journal of physiology.

[3]  J E Bailey,et al.  Metabolic flux analysis with a comprehensive isotopomer model in Bacillus subtilis. , 2001, Biotechnology and bioengineering.

[4]  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.

[5]  Philip E. Gill,et al.  Numerical Linear Algebra and Optimization , 1991 .

[6]  Kazuyuki Shimizu,et al.  An improved method for statistical analysis of metabolic flux analysis using isotopomer mapping matrices with analytical expressions. , 2003, Journal of biotechnology.

[7]  G. Stephanopoulos,et al.  Systematic quantification of complex metabolic flux networks using stable isotopes and mass spectrometry. , 2003, European journal of biochemistry.

[8]  Wolfgang M. Hartmann,et al.  Computing the Moore-Penrose Inverse for the Covariance Matrix in Constrained Nonlinear Estimation , 1996, SIAM J. Optim..

[9]  J. Wahren,et al.  Limitations in estimating gluconeogenesis and Cori cycling from mass isotopomer distributions using [U-13C6]glucose. , 1998, American journal of physiology. Endocrinology and metabolism.

[10]  A. Sunehag The Role of Parenteral Lipids in Supporting Gluconeogenesis in Very Premature Infants , 2003, Pediatric Research.

[11]  G. Stephanopoulos Metabolic fluxes and metabolic engineering. , 1999, Metabolic engineering.

[12]  W Wiechert,et al.  A universal framework for 13C metabolic flux analysis. , 2001, Metabolic engineering.

[13]  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.

[14]  J. Kelleher Estimating gluconeogenesis with [U-13C]glucose: molecular condensation requires a molecular approach. , 1999, American journal of physiology. Endocrinology and metabolism.

[15]  J. Nielsen,et al.  Metabolic network analysis of Penicillium chrysogenum using (13)C-labeled glucose. , 2000, Biotechnology and bioengineering.

[16]  T. Szyperski Biosynthetically directed fractional 13C-labeling of proteinogenic amino acids. An efficient analytical tool to investigate intermediary metabolism. , 1995, European journal of biochemistry.

[17]  J. Radziuk,et al.  Measurement of gluconeogenesis and mass isotopomer analysis based on [U-13C]glucose. , 1999, American journal of physiology. Endocrinology and metabolism.

[18]  Boris Kalitventzeff,et al.  Sensitivity calculations and variance analysis in plant measurement reconciliation , 1996 .

[19]  U. Sauer,et al.  Metabolic fluxes in riboflavin-producing Bacillus subtilis , 1997, Nature Biotechnology.

[20]  L. Foster Rank and null space calculations using matrix decomposition without column interchanges , 1986 .

[21]  J. Wahren,et al.  Contributions of gluconeogenesis to glucose production in the fasted state. , 1996, The Journal of clinical investigation.

[22]  J K Kelleher,et al.  Applications of mass isotopomer analysis to nutrition research. , 1997, Annual review of nutrition.

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

[24]  J. Heijnen,et al.  Metabolic-flux analysis of Saccharomyces cerevisiae CEN.PK113-7D based on mass isotopomer measurements of (13)C-labeled primary metabolites. , 2005, FEMS yeast research.

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

[26]  Jorge Nocedal,et al.  An Interior Point Algorithm for Large-Scale Nonlinear Programming , 1999, SIAM J. Optim..

[27]  T. Szyperski Biosynthetically Directed Fractional 13C‐labeling of Proteinogenic Amino Acids , 1995 .

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

[29]  Sahm,et al.  13C tracer experiments and metabolite balancing for metabolic flux analysis: comparing two approaches , 1998, Biotechnology and bioengineering.

[30]  W. J. Conover,et al.  Practical Nonparametric Statistics , 1972 .

[31]  P. Reeds,et al.  Gluconeogenesis measured with [U-13C]glucose and mass isotopomer analysis of apoB-100 amino acids in pigs. , 1998, American journal of physiology. Endocrinology and metabolism.

[32]  J. Kelleher,et al.  Probing metabolic pathways with isotopic tracers: insights from mammalian metabolic physiology. , 2004, Metabolic engineering.

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

[34]  A J Sinskey,et al.  Metabolite and isotopomer balancing in the analysis of metabolic cycles: II. Applications. , 1999, Biotechnology and bioengineering.

[35]  P. Verheijen,et al.  Cumulative bondomers: a new concept in flux analysis from 2D [13C,1H] COSY NMR data. , 2002, Biotechnology and bioengineering.

[36]  L. Groop,et al.  Glucose and free fatty acid metabolism in non-insulin-dependent diabetes mellitus. Evidence for multiple sites of insulin resistance. , 1989, The Journal of clinical investigation.