Hybrid optimization for 13C metabolic flux analysis using systems parametrized by compactification

BackgroundThe importance and power of isotope-based metabolic flux analysis and its contribution to understanding the metabolic network is increasingly recognized. Its application is, however, still limited partly due to computational inefficiency. 13C metabolic flux analysis aims to compute in vivo metabolic fluxes in terms of metabolite balancing extended by carbon isotopomer balances and involves a nonlinear least-squares problem. To solve the problem more efficiently, improved numerical optimization techniques are necessary.ResultsFor flux computation, we developed a gradient-based hybrid optimization algorithm. Here, independent flux variables were compactified into [0, 1)-ranged variables using a single transformation rule. The compactified parameters could be discriminated between non-identifiable and identifiable variables after model linearization. The developed hybrid algorithm was applied to the central metabolism of Bacillus subtilis with only succinate and glutamate as carbon sources. This creates difficulties caused by symmetry of succinate leading to limited introduction of 13C labeling information into the system. The algorithm was found to be superior to its parent algorithms and to global optimization methods both in accuracy and speed. The hybrid optimization with tolerance adjustment quickly converged to the minimum with close to zero deviation and exactly re-estimated flux variables. In the metabolic network studied, some fluxes were found to be either non-identifiable or nonlinearly correlated. The non-identifiable fluxes could correctly be predicted a priori using the model identification method applied, whereas the nonlinear flux correlation was revealed only by identification runs using different starting values a posteriori.ConclusionThis fast, robust and accurate optimization method is useful for high-throughput metabolic flux analysis, a posteriori identification of possible parameter correlations, and also for Monte Carlo simulations to obtain statistical qualities for flux estimates. In this way, it contributes to future quantitative studies of central metabolic networks in the framework of systems biology.

[1]  Christoph Wittmann,et al.  Theoretical aspects of 13C metabolic flux analysis with sole quantification of carbon dioxide labeling , 2005, Comput. Biol. Chem..

[2]  T. Coleman,et al.  On the Convergence of Reflective Newton Methods for Large-scale Nonlinear Minimization Subject to Bounds , 1992 .

[3]  Douglas B. Kell,et al.  Non-linear optimization of biochemical pathways: applications to metabolic engineering and parameter estimation , 1998, Bioinform..

[4]  W. Wiechert.  13C metabolic flux analysis. , 2001, Metabolic engineering.

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

[6]  Patricia Brackin,et al.  Using Genetic Algorithms to Set Target Values for Engineering Characteristics in the House of Quality , 2002, J. Comput. Inf. Sci. Eng..

[7]  Arne Elofsson.  Bioinformatics: From nucleic acids and proteins to cell metabolism: Edited by D. Schomburg and U. Lessel, VCH; Weinheim-New York, 1995. viii + 195 pp. DM 148.00 (hb). ISBN 3-527-30072-4 , 1996 .

[8]  K. Shimizu,et al.  Metabolic flux analysis of Escherichia coli K12 grown on 13C-labeled acetate and glucose using GC-MS and powerful flux calculation method. , 2003, Journal of biotechnology.

[9]  G. Stephanopoulos,et al.  Elementary metabolite units (EMU): a novel framework for modeling isotopic distributions. , 2007, Metabolic engineering.

[10]  Paul T. Boggs,et al.  Sequential Quadratic Programming , 1995, Acta Numerica.

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

[12]  Louis Wehenkel,et al.  A hybrid optimization technique coupling an evolutionary and a local search algorithm , 2008 .

[13]  S. D. Jong,et al.  Handbook of Chemometrics and Qualimetrics , 1998 .

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

[15]  Peiliang Xu.  A hybrid global optimization method: the multi-dimensional case , 2003 .

[16]  Christoph Wittmann,et al.  Metabolic network simulation using logical loop algorithm and Jacobian matrix. , 2004, Metabolic engineering.

[17]  D K Smith,et al.  Numerical Optimization , 2001, J. Oper. Res. Soc..

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

[19]  Mary C Hill,et al.  Determining Extreme Parameter Correlation in Ground Water Models , 2003, Ground water.

[20]  J J Heijnen,et al.  A priori analysis of metabolic flux identifiability from (13)C-labeling data. , 2001, Biotechnology and bioengineering.

[21]  C. Wittmann,et al.  Application of MALDI-TOF MS to lysine-producing Corynebacterium glutamicum: a novel approach for metabolic flux analysis. , 2001, European journal of biochemistry.

[22]  H. Blanch,et al.  Using isotopomer path tracing to quantify metabolic fluxes in pathway models containing reversible reactions. , 2001, Biotechnology and bioengineering.

[23]  Carmen G. Moles,et al.  Parameter estimation in biochemical pathways: a comparison of global optimization methods. , 2003, Genome research.

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

[25]  K. Schittkowski.  NLPQL: A fortran subroutine solving constrained nonlinear programming problems , 1986 .

[26]  W. Wiechert,et al.  Algebraic Methods for the Analysis of Redundancy and Identifiability in Metabolic 13C‐Labelling Systems , 2007 .

[27]  David G. Luenberger,et al.  Linear and nonlinear programming , 1984 .

[28]  William H. Press,et al.  The Art of Scientific Computing Second Edition , 1998 .

[29]  J L Klepeis,et al.  Hybrid global optimization algorithms for protein structure prediction: alternating hybrids. , 2003, Biophysical journal.

[30]  W Wiechert,et al.  Metabolic isotopomer labeling systems. Part I: global dynamic behavior. , 2001, Mathematical biosciences.

[31]  Steffen Klamt,et al.  Calculability analysis in underdetermined metabolic networks illustrated by a model of the central metabolism in purple nonsulfur bacteria. , 2002, Biotechnology and bioengineering.

[32]  Douglas C. Hittle,et al.  Optimization of autonomous village electrification systems by simulated annealing , 2000 .

[33]  Maria Pia Saccomani,et al.  Some Results on Parameter Identification of Nonlinear Systems , 2004 .

[34]  Z. P. Wu,et al.  A novel hybrid genetic algorithm using local optimizer based on heuristic pattern move , 2001, Appl. Artif. Intell..

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

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

[37]  Christoph Wittmann,et al.  Genealogy Profiling through Strain Improvement by Using Metabolic Network Analysis: Metabolic Flux Genealogy of Several Generations of Lysine-Producing Corynebacteria , 2002, Applied and Environmental Microbiology.

[38]  Panos M. Pardalos,et al.  Recent Advances in Global Optimization , 1991 .

[39]  G. Stephanopoulos,et al.  Metabolic flux analysis in a nonstationary system: fed-batch fermentation of a high yielding strain of E. coli producing 1,3-propanediol. , 2007, Metabolic engineering.

[40]  U. Sauer,et al.  Physiology and metabolic fluxes of wild-type and riboflavin-producing Bacillus subtilis , 1996, Applied and environmental microbiology.

[41]  Edward P. Gatzke,et al.  Nonlinear model predictive control using deterministic global optimization , 2006 .

[42]  G. Mahinthakumar,et al.  Hybrid Genetic Algorithm—Local Search Methods for Solving Groundwater Source Identification Inverse Problems , 2005 .

[43]  F. A. Seiler,et al.  Numerical Recipes in C: The Art of Scientific Computing , 1989 .

[44]  S. Pushpavanam,et al.  Kinetic parameter estimation in hydrocracking using a combination of genetic algorithm and sequential quadratic programming , 2003 .

[45]  Min-Jea Tahk,et al.  A hybrid optimization method of evolutionary and gradient search , 2007 .

[46]  Wolfgang Wiechert,et al.  Metabolic isotopomer labeling systems. Part II: structural flux identifiability analysis. , 2003, Mathematical biosciences.