Dynamic flux balance analysis with nonlinear objective function

Dynamic flux balance analysis (DFBA) extends flux balance analysis and enables the combined simulation of both intracellular and extracellular environments of microbial cultivation processes. A DFBA model contains two coupled parts, a dynamic part at the upper level (extracellular environment) and an optimization part at the lower level (intracellular environment). Both parts are coupled through substrate uptake and product secretion rates. This work proposes a Karush–Kuhn–Tucker condition based solution approach for DFBA models, which have a nonlinear objective function in the lower-level part. To solve this class of DFBA models an extreme-ray-based reformulation is proposed to ensure certain regularity of the lower-level optimization problem. The method is introduced by utilizing two simple example networks and then applied to a realistic model of central carbon metabolism of wild-type Corynebacterium glutamicum.

[1]  J Ricardo Pérez-Correa,et al.  Data reconciliation and parameter estimation in flux‐balance analysis , 2003, Biotechnology and bioengineering.

[2]  P I Barton,et al.  A reliable simulator for dynamic flux balance analysis , 2013, Biotechnology and bioengineering.

[3]  Kenneth J. Kauffman,et al.  Advances in flux balance analysis. , 2003, Current opinion in biotechnology.

[4]  N. Yu. Zolotykh,et al.  New modification of the double description method for constructing the skeleton of a polyhedral cone , 2012 .

[5]  Alexander Schrijver,et al.  Theory of linear and integer programming , 1986, Wiley-Interscience series in discrete mathematics and optimization.

[6]  Chantal Landry,et al.  A Second Order Scheme for Solving Optimization-Constrained Differential Equations with Discontinuities , 2008 .

[7]  W. Wiechert,et al.  Isotopically non-stationary metabolic flux analysis: complex yet highly informative. , 2013, Current opinion in biotechnology.

[8]  B. Palsson Systems Biology: Properties of Reconstructed Networks , 2006 .

[9]  Jerzy Kyparisis,et al.  On uniqueness of Kuhn-Tucker multipliers in nonlinear programming , 1985, Math. Program..

[10]  B O Palsson,et al.  Flux-balance analysis of mitochondrial energy metabolism: consequences of systemic stoichiometric constraints. , 2001, American journal of physiology. Regulatory, integrative and comparative physiology.

[11]  Luca Dieci,et al.  A survey of numerical methods for IVPs of ODEs with discontinuous right-hand side , 2012, J. Comput. Appl. Math..

[12]  D. Sarkar,et al.  Dynamic flux balance analysis of batch fermentation: effect of genetic manipulations on ethanol production , 2014, Bioprocess and Biosystems Engineering.

[13]  B. Palsson,et al.  The Escherichia coli MG1655 in silico metabolic genotype: its definition, characteristics, and capabilities. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[14]  O. Mangasarian,et al.  The Fritz John Necessary Optimality Conditions in the Presence of Equality and Inequality Constraints , 1967 .

[15]  Steffen Klamt,et al.  Two approaches for metabolic pathway analysis? , 2003, Trends in biotechnology.

[16]  D. W. Peterson A REVIEW OF CONSTRAINT QUALIFICATIONS IN FINITE-DIMENSIONAL SPACES* , 1973 .

[17]  G. Stephanopoulos,et al.  Metabolic Engineering: Principles And Methodologies , 1998 .

[18]  Stephen J. Wright,et al.  Numerical Optimization (Springer Series in Operations Research and Financial Engineering) , 2000 .

[19]  Jason A. Papin,et al.  Comparison of network-based pathway analysis methods. , 2004, Trends in biotechnology.

[20]  Paul I. Barton,et al.  DFBAlab: a fast and reliable MATLAB code for dynamic flux balance analysis , 2014, BMC Bioinformatics.

[21]  Chantal Landry,et al.  Optimization Problem Coupled with Differential Equations: A Numerical Algorithm Mixing an Interior-Point Method and Event Detection , 2010, J. Optim. Theory Appl..

[22]  Michael A. Saunders,et al.  SNOPT: An SQP Algorithm for Large-Scale Constrained Optimization , 2005, SIAM Rev..

[23]  F. Doyle,et al.  Dynamic flux balance analysis of diauxic growth in Escherichia coli. , 2002, Biophysical journal.

[24]  Andreas Kremling,et al.  Study of the growth of Escherichia coli on mixed substrates using dynamic flux balance analysis , 2010 .

[25]  Lawrence F. Shampine,et al.  The MATLAB ODE Suite , 1997, SIAM J. Sci. Comput..

[26]  J. Heijnen,et al.  A metabolic network stoichiometry analysis of microbial growth and product formation , 1995, Biotechnology and bioengineering.

[27]  Wolfgang Wiechert,et al.  Visualizing multi-omics data in metabolic networks with the software Omix - A case study , 2011, Biosyst..

[28]  Michael A Henson,et al.  Optimization of Fed‐Batch Saccharomyces cerevisiae Fermentation Using Dynamic Flux Balance Models , 2006, Biotechnology progress.

[29]  Paul I. Barton,et al.  Efficient solution of ordinary differential equations with a parametric lexicographic linear program embedded , 2015, Numerische Mathematik.

[30]  G. T. Tsao,et al.  Cybernetic modeling of microbial growth on multiple substrates , 1984, Biotechnology and bioengineering.

[31]  Ernst Hairer,et al.  Solving Optimization-Constrained Differential Equations with Discontinuity Points, with Application to Atmospheric Chemistry , 2009, SIAM J. Sci. Comput..

[32]  D. Fell,et al.  Reaction routes in biochemical reaction systems: Algebraic properties, validated calculation procedure and example from nucleotide metabolism , 2002, Journal of mathematical biology.

[33]  U. Sauer,et al.  Systematic evaluation of objective functions for predicting intracellular fluxes in Escherichia coli , 2007, Molecular systems biology.