AbstractConstraint-based metabolic flux analysis of knockout strategies is an efficient method to simulate the production of useful metabolites in microbes. Owing to the recent development of technologies for artificial DNA synthesis, it may become important in the near future to mathematically design minimum metabolic networks to simulate metabolite production. Accordingly, we have developed a computational method where parsimonious metabolic flux distribution is computed for designated constraints on growth and production rates which are represented by grids. When the growth rate of this obtained parsimonious metabolic network is maximized, higher production rates compared to those noted using existing methods are observed for many target metabolites. The set of reactions used in this parsimonious flux distribution consists of reactions included in the original genome scale model iAF1260. The computational experiments show that the grid size affects the obtained production rates. Under the conditions that the growth rate is maximized and the minimum cases of flux variability analysis are considered, the developed method produced more than 90% of metabolites, while the existing methods produced less than 50%. Mathematical explanations using examples are provided to demonstrate potential reasons for the ability of the proposed algorithm to identify design strategies that the existing methods could not identify. The source code is freely available, and is implemented in MATLAB and COBRA toolbox.Author summaryMetabolic networks represent the relationships between biochemical reactions and compounds in living cells. By computationally modifying a given metabolic network of microbes, we can simulate the effect of knockouts and estimate the production of valuable metabolites. A common mathematical model of metabolic networks is the constraint-based flux model. In constraint-based flux balance analysis, a pseudo-steady state is assumed to predict the metabolic profile where the sum of all incoming fluxes is equal to the sum of all outgoing fluxes for each internal metabolite. Based on these constraints, the biomass objective function, written as a linear combination of fluxes, is maximized. In this study, we developed an efficient method for computing the design of minimum metabolic networks by using constraint-based flux balance analysis to simulate the production of useful metabolites.
Grid-based computational methods for the design of constraint-based parsimonious chemical reaction networks to simulate metabolite production: GridProd