Improved phenotypic predictions of metabolic models by integrating regulatory constraints

Genome-scale metabolic networks (GSMs) are mathematic representation of a set of stoichiometrically balanced reactions. However, such static GSMs do not reflect or incorporate functional organization of genes and their dynamic regulation (e.g., operons and regulons). Specifically, there are numerous topologically coupled local reactions through which fluxes are coordinated; and downstream metabolites often dynamically regulate the gene expression of their reactions via feedback. Here, we present a method which reconstructs GSMs with locally coupled reactions and transcriptional regulation of metabolism by key metabolites. The proposed method has outstanding performance in phenotype prediction of wild-type and mutants in Escherichia coli (E. coli), Saccharomyces cerevisiae (S. cerevisiae) and Bacillus subtilis (B. subtilis) growing in various conditions, outperforming existing methods. The predicted growth rate and metabolic fluxes are highly correlated with those experimentally measured. More importantly, our method can also explain the observed growth rates by capturing the real (experimentally measured) changes in flux between the wild-types and mutants. Overall, by identifying and incorporating locally organized and regulated functional modules into GSMs, Decrem achieves accurate predictions of phenotypes and has broad applications in bioengineering, synthetic biology and microbial pathology.