Optimization of Fed‐Batch Saccharomyces cerevisiae Fermentation Using Dynamic Flux Balance Models

We developed a dynamic flux balance model for fed‐batch Saccharomyces cerevisiae fermentation that couples a detailed steady‐state description of primary carbon metabolism with dynamic mass balances on key extracellular species. Model‐based dynamic optimization is performed to determine fed‐batch operating policies that maximize ethanol productivity and/or ethanol yield on glucose. The initial volume and glucose concentrations, the feed flow rate and dissolved oxygen concentration profiles, and the final batch time are treated as decision variables in the dynamic optimization problem. Optimal solutions are generated to analyze the tradeoff between maximal productivity and yield objectives. We find that for both cases the prediction of a microaerobic region is significant. The optimization results are sensitive to network model parameters for the growth associated maintenance and P/O ratio. The results of our computational study motivate continued development of dynamic flux balance models and further exploration of their application to productivity optimization in biochemical reactors.

[1]  John M. Walker,et al.  Molecular Biology and Biotechnology , 1988 .

[2]  B. Palsson,et al.  Stoichiometric flux balance models quantitatively predict growth and metabolic by-product secretion in wild-type Escherichia coli W3110 , 1994, Applied and environmental microbiology.

[3]  A. Johnson,et al.  The control of fed-batch fermentation processes - A survey , 1987, Autom..

[4]  P A Vanrolleghem,et al.  Validation of a Metabolic Network for Saccharomyces cerevisiae Using Mixed Substrate Studies , 1996, Biotechnology progress.

[5]  A. Burgard,et al.  Optimization-based framework for inferring and testing hypothesized metabolic objective functions. , 2003, Biotechnology and bioengineering.

[6]  L. Biegler,et al.  Advances in simultaneous strategies for dynamic process optimization , 2002 .

[7]  G. Church,et al.  Analysis of optimality in natural and perturbed metabolic networks , 2002 .

[8]  Julio R. Banga,et al.  Stochastic Dynamic Optimization of Batch and Semicontinuous Bioprocesses , 1997 .

[9]  Gerhard K Hoppe,et al.  The effect of micro-aerobic conditions on continuous ethanol production by Saccharomyces cerevisiae , 1984, Biotechnology Letters.

[10]  G. Goma,et al.  Aeration strategy: a need for very high ethanol performance in Saccharomyces cerevisiae fed-batch process , 2004, Applied Microbiology and Biotechnology.

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

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

[13]  Eduardo Agosin,et al.  Modeling of yeast metabolism and process dynamics in batch fermentation , 2003, Biotechnology and bioengineering.

[14]  Jung-Heon Lee Comparison of various optimization approaches for fed-batch ethanol production , 1999, Applied biochemistry and biotechnology.

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

[16]  W. Gropp,et al.  Accepted for publication , 2001 .

[17]  Bernard R. Glick,et al.  Molecular biotechnology : principles and applications ofrecombinant DNA , 2010 .

[18]  Marta Cascante,et al.  Multicriteria optimization of biochemical systems by linear programming: application to production of ethanol by Saccharomyces cerevisiae. , 2003, Biotechnology and bioengineering.

[19]  G. Stephanopoulos CHAPTER 1 – The Essence of Metabolic Engineering , 1998 .

[20]  J. E. Cuthrell,et al.  Simultaneous optimization and solution methods for batch reactor control profiles , 1989 .

[21]  J. Nielsen,et al.  Flux distributions in anaerobic, glucose-limited continuous cultures of Saccharomyces cerevisiae. , 1997, Microbiology.

[22]  R. Mahadevan,et al.  Estimating optimal profiles of genetic alterations using constraint-based models. , 2005, Biotechnology and bioengineering.

[23]  U. Ascher,et al.  A New Basis Implementation for a Mixed Order Boundary Value ODE Solver , 1987 .

[24]  A Lübbert,et al.  Bioreactor performance: a more scientific approach for practice. , 2001, Journal of biotechnology.

[25]  B Sonnleitner,et al.  Growth of Saccharomyces cerevisiae is controlled by its limited respiratory capacity: Formulation and verification of a hypothesis , 1986, Biotechnology and bioengineering.

[26]  Rein Luus,et al.  Application of dynamic programming to differential-algebraic process systems , 1993 .

[27]  C Phalakornkule,et al.  A MILP-based flux alternative generation and NMR experimental design strategy for metabolic engineering. , 2001, Metabolic engineering.

[28]  B. Palsson,et al.  Genome-scale reconstruction of the Saccharomyces cerevisiae metabolic network. , 2003, Genome research.

[29]  Attilio Converti,et al.  Simplified modeling of fed‐batch alcoholic fermentation of sugarcane blackstrap molasses , 2003, Biotechnology and bioengineering.

[30]  A. Burgard,et al.  Exploring the overproduction of amino acids using the bilevel optimization framework OptKnock , 2003, Biotechnology and bioengineering.

[31]  G. Church,et al.  Genome-Scale Metabolic Model of Helicobacter pylori 26695 , 2002, Journal of bacteriology.

[32]  Bethany L. Nicholson,et al.  Mathematical Programs with Equilibrium Constraints , 2021, Pyomo — Optimization Modeling in Python.

[33]  M. Taherzadeh,et al.  Use of dynamic step response for control of fed-batch conversion of lignocellulosic hydrolyzates to ethanol. , 2001, Journal of biotechnology.

[34]  Bernhard Ø Palsson,et al.  Integrated analysis of metabolic phenotypes in Saccharomyces cerevisiae , 2004, BMC Genomics.

[35]  F.-S. Wang,et al.  Optimal feed policy for fed-batch fermentation of ethanol production by Zymomous mobilis , 1997 .

[36]  B. Palsson,et al.  Saccharomyces cerevisiae phenotypes can be predicted by using constraint-based analysis of a genome-scale reconstructed metabolic network , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[37]  Ioannis K Kookos,et al.  Optimization of Batch and Fed‐Batch Bioreactors using Simulated Annealing , 2004, Biotechnology progress.

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

[39]  A. Burgard,et al.  Optknock: A bilevel programming framework for identifying gene knockout strategies for microbial strain optimization , 2003, Biotechnology and bioengineering.

[40]  Ilse Smets,et al.  Optimal adaptive control of (bio)chemical reactors: past, present and future , 2004 .

[41]  G Stephanopoulos,et al.  Physiological, biochemical, and mathematical studies of micro‐aerobic continuous ethanol fermentation by Saccharomyces cerevisiae. I: Hysteresis, oscillations, and maximum specific ethanol productivities in chemostat culture , 1990, Biotechnology and bioengineering.

[42]  Ignacio E. Grossmann,et al.  Retrospective on optimization , 2004, Comput. Chem. Eng..

[43]  Lorenz T. Biegler,et al.  Mathematical programs with equilibrium constraints (MPECs) in process engineering , 2003, Comput. Chem. Eng..

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