Yield optimization of regulated metabolic systems using deterministic branch‐and‐reduce methods

This paper addresses product yield optimization in microorganisms grown in continuous culture. Traditional optimization strategies of random mutagenesis and selection will eventually have limited efficacy, thus requiring more focused strategies. The best candidates for such strategies appear to be mathematical models that capture the essence of metabolic systems and permit optimization with computational methods. In the past, models used for this purpose have been stoichiometric, kinetic in the form of S‐systems, or ad hoc. This work presents a deterministic approach based on generalized mass action (GMA) systems. These systems are interesting in that they allow direct merging of stoichiometric and S‐system models. Two illustrations are considered. In the first case, the fermentation pathway in Saccharomyces cerevisiae is optimized for ethanol production under steady‐state conditions. The model of this pathway is relatively small, with five states and eight rate constants. The second example addresses the maximization of citric acid in the mold Aspergillus niger. For the optimization of this larger pathway system with 30 states and 60 reactions, a Mixed Integer Nonlinear Programming (MINLP) is proposed. It is shown that efficient MINLP algorithms, based on convexification, branch‐and‐reduce methods, and binary variable selection, are essential for solving these difficult optimization problems. Biotechnol. Bioeng. 2008;99: 1154–1169. © 2007 Wiley Periodicals, Inc.

[1]  G. Gavalas.  Nonlinear Differential Equations of Chemically Reacting Systems , 1968 .

[2]  P. I. Barton,et al.  Outer approximation algorithms for separable nonconvex mixed-integer nonlinear programs , 2004, Math. Program..

[3]  Eldon Hansen,et al.  Global optimization using interval analysis , 1992, Pure and applied mathematics.

[4]  J. Bailey,et al.  Analysis and design of metabolic reaction networks via mixed‐integer linear optimization , 1996 .

[5]  Eberhard O Voit,et al.  Theoretical Biology and Medical Modelling Identification of Metabolic System Parameters Using Global Optimization Methods , 2022 .

[6]  M. Savageau.  Biochemical Systems Analysis: A Study of Function and Design in Molecular Biology , 1976 .

[7]  J. Bailey,et al.  Fermentation pathway kinetics and metabolic flux control in suspended and immobilized Saccharomyces cerevisiae , 1990 .

[8]  J E Bailey,et al.  MCA has more to say. , 1996, Journal of theoretical biology.

[9]  M. Cánovas,et al.  High-density Escherichia coli cultures for continuous l(−)-carnitine production , 1999, Applied Microbiology and Biotechnology.

[10]  A. M. Geoffrion.  Generalized Benders decomposition , 1972 .

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

[12]  A Sorribas,et al.  Strategies for representing metabolic pathways within biochemical systems theory: reversible pathways. , 1989, Mathematical biosciences.

[13]  E. Voit,et al.  Recasting nonlinear differential equations as S-systems: a canonical nonlinear form , 1987 .

[14]  J. Heijnen.  Approximative kinetic formats used in metabolic network modeling , 2005, Biotechnology and bioengineering.

[15]  M. Savageau.  Biochemical systems analysis. II. The steady-state solutions for an n-pool system using a power-law approximation. , 1969, Journal of theoretical biology.

[16]  J. Bailey,et al.  Optimization of regulatory architectures in metabolic reaction networks , 1996, Biotechnology and bioengineering.

[17]  Eberhard O. Voit,et al.  Biochemical systems theory and metabolic control theory: 1. fundamental similarities and differences , 1987 .

[18]  Manfred Morari,et al.  Propositional logic in control and monitoring problems , 1997, 1997 European Control Conference (ECC).

[19]  R. Horst,et al.  Global Optimization: Deterministic Approaches , 1992 .

[20]  S. Ghosh,et al.  Closing the loop between feasible flux scenario identification for construct evaluation and resolution of realized fluxes via NMR , 2005, Comput. Chem. Eng..

[21]  E. Voit,et al.  Pathway Analysis and Optimization in Metabolic Engineering , 2002 .

[22]  Ignacio E. Grossmann,et al.  An outer-approximation algorithm for a class of mixed-integer nonlinear programs , 1986, Math. Program..

[23]  G. Stephanopoulos,et al.  Strain improvement by metabolic engineering: lysine production as a case study for systems biology. , 2005, Current opinion in biotechnology.

[24]  Nikolaos V. Sahinidis,et al.  Optimization of metabolic pathways under stability considerations , 2005, Comput. Chem. Eng..

[25]  Nikolaos V. Sahinidis,et al.  Global optimization of mixed-integer nonlinear programs: A theoretical and computational study , 2004, Math. Program..

[26]  B. Palsson,et al.  Metabolic Flux Balancing: Basic Concepts, Scientific and Practical Use , 1994, Bio/Technology.

[27]  Eberhard O. Voit,et al.  Computational Analysis of Biochemical Systems: A Practical Guide for Biochemists and Molecular Biologists , 2000 .

[28]  M. Savageau.  Development of fractal kinetic theory for enzyme-catalysed reactions and implications for the design of biochemical pathways. , 1998, Bio Systems.

[29]  Ramon E. Moore.  Methods and applications of interval analysis , 1979, SIAM studies in applied mathematics.

[30]  E O Voit,et al.  Optimization in integrated biochemical systems , 1992, Biotechnology and bioengineering.

[31]  I. D. L. Bogle,et al.  Global optimisation of constrained non-convex programs using reformulation and interval analysis , 1999 .

[32]  Ignacio E. Grossmann,et al.  Global optimization algorithm for heat exchanger networks , 1993 .

[33]  M. Cascante,et al.  Comparative characterization of the fermentation pathway of Saccharomyces cerevisiae using biochemical systems theory and metabolic control analysis: model definition and nomenclature. , 1995, Mathematical biosciences.

[34]  E. Voit,et al.  An indirect optimization method for biochemical systems: description of method and application to the maximization of the rate of ethanol, glycerol, and carbohydrate production in Saccharomyces cerevisiae. , 1997, Biotechnology and bioengineering.

[35]  James N. Currie,et al.  THE CITRIC ACID FERMENTATION OF ASPERGILLUS NIGER , 1917 .

[36]  R. Heinrich,et al.  The Regulation of Cellular Systems , 1996, Springer US.

[37]  Arthur R. Schulz,et al.  Enzyme Kinetics: From Diastase to Multi-enzyme Systems , 1994 .

[38]  J. F. Benders.  Partitioning procedures for solving mixed-variables programming problems , 1962 .

[39]  Garth P. McCormick,et al.  Computability of global solutions to factorable nonconvex programs: Part I — Convex underestimating problems , 1976, Math. Program..

[40]  J. E. Falk,et al.  An Algorithm for Separable Nonconvex Programming Problems , 1969 .

[41]  M A Savageau,et al.  Accuracy of alternative representations for integrated biochemical systems. , 1987, Biochemistry.

[42]  Mahmoud M. El-Halwagi,et al.  Global optimization of nonconvex nonlinear programs via interval analysis , 1994 .

[43]  E O Voit,et al.  Optimization of nonlinear biotechnological processes with linear programming: Application to citric acid production by Aspergillus niger , 2000, Biotechnology and bioengineering.

[44]  Nikolaos V. Sahinidis,et al.  Convexification and Global Optimization in Continuous and Mixed-Integer Nonlinear Programming , 2002 .

[45]  M. Savageau.  Biochemical systems analysis. II. The steady-state solutions for an n-pool system using a power-law approximation. , 1969, Journal of theoretical biology.

[46]  Ignacio E. Grossmann,et al.  Logic-based outer approximation for globally optimal synthesis of process networks , 2005, Comput. Chem. Eng..

[47]  E. Papoutsakis,et al.  Stoichiometric modeling of Clostridium acetobutylicum fermentations with non-linear constraints. , 1999, Journal of biotechnology.