In Equation 1, N is the m × r matrix of the stoichiometric coefficients for the r reactions and m metabolites in the model, and v is the r × 1 vector of the steady-state fluxes through the r reactions in the model. Bounds are placed on the reaction fluxes depending on the reversibility of the reactions:
- (CDW = cell dry weight). When simulating a gene knockout, the bounds on the flux through all reactions associated exclusively with the gene being knocked out (or associated exclusively with a protein complex partially encoded by the gene being knocked out) were reset to zero. When simulating media conditions, only nutrients present in the media were allowed to have a net uptake by the cell. All other transportable nutrients were allowed only to be excreted by the cell. Details on conditions for all FBA simulations performed are provided in Table S8 in Additional data files 1 and 2.