A dynamic flux balance model based on a genome-scale metabolic network reconstruction is developed for in silico analysis of Saccharomyces cerevisiae metabolism and ethanol production in fed-batch culture. Metabolic engineering strategies previously identified for their enhanced steady-state biomass and/or ethanol yields are evaluated for fed-batch performance in glucose and glucose/xylose media. Dynamic analysis is shown to provide a single quantitative measure of fed-batch ethanol productivity that explicitly handles the possible tradeoff between the biomass and ethanol yields. Productivity optimization conducted to rank achievable fed-batch performance demonstrates that the genetic manipulation strategy and the fed-batch operating policy should be considered simultaneously. A library of candidate gene insertions is assembled and directly screened for their achievable ethanol productivity in fed-batch culture. A number of novel gene insertions with ethanol productivities identical to the best metabolic engineering strategies reported in previous studies are identified, thereby providing additional targets for experimental evaluation. The top performing gene insertions were substrate dependent, with the highest ranked insertions for glucose media yielding suboptimal performance in glucose/xylose media. The analysis results suggest that enhancements in biomass yield are most beneficial for the enhancement of fed-batch ethanol productivity by recombinant xylose utilizing yeast strains. We conclude that steady-state flux balance analysis is not sufficient to predict fed-batch performance and that the media, genetic manipulations, and fed-batch operating policy should be considered simultaneously to achieve optimal metabolite productivity.

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.

Steady-state and dynamic flux balance analysis (DFBA) was used to investigate the effects of metabolic model complexity and parameters on ethanol production predictions for wild-type and engineered Saccharomyces cerevisiae. Three metabolic network models ranging from a single compartment representation of metabolism to a genome-scale reconstruction with seven compartments and detailed charge balancing were studied. Steady-state analysis showed that the models generated similar wild-type predictions for the biomass and ethanol yields, but for ten engineered strains the seven compartment model produced smaller ethanol yield enhancements. Simplification of the seven compartment model to two intracellular compartments produced increased ethanol yields, suggesting that reaction localisation had an impact on mutant phenotype predictions. Further analysis with the seven compartment model demonstrated that steady-state predictions can be sensitive to intracellular model parameters, with the biomass yield exhibiting high sensitivity to ATP utilisation parameters and the biomass composition. The incorporation of gene expression data through the zeroing of metabolic reactions associated with unexpressed genes was shown to produce negligible changes in steady-state predictions when the oxygen uptake rate was suitably constrained. Dynamic extensions of the single and seven compartment models were developed through the addition of glucose and oxygen uptake expressions and transient extracellular balances. While the dynamic models produced similar predictions of the optimal batch ethanol productivity for the wild type, the single compartment model produced significantly different predictions for four implementable gene insertions. A combined deletion/overexpression/insertion mutant with improved ethanol productivity capabilities was computationally identified by dynamically screening multiple combinations of the ten metabolic engineering strategies. The authors concluded that extensive compartmentalisation and detailed charge balancing can be important for reliably screening metabolic engineering strategies that rely on modification of the global redox balance and that DFBA offers the potential to identify novel mutants for enhanced metabolite production in batch and fed-batch cultures.

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.