Construction and completion of flux balance models from pathway databases

Abstract

Motivation: Flux balance analysis (FBA) is a well-known technique for genome-scale modeling of metabolic flux. Typically, an FBA formulation requires the accurate specification of four sets: biochemical reactions, biomass metabolites, nutrients and secreted metabolites. The development of FBA models can be time consuming and tedious because of the difficulty in assembling completely accurate descriptions of these sets, and in identifying errors in the composition of these sets. For example, the presence of a single non-producible metabolite in the biomass will make the entire model infeasible. Other difficulties in FBA modeling are that model distributions, and predicted fluxes, can be cryptic and difficult to understand.Results: We present a multiple gap-filling method to accelerate the development of FBA models using a new tool, called MetaFlux, based on mixed integer linear programming (MILP). The method suggests corrections to the sets of reactions, biomass metabolites, nutrients and secretions. The method generates FBA models directly from Pathway/Genome Databases. Thus, FBA models developed in this framework are easily queried and visualized using the Pathway Tools software. Predicted fluxes are more easily comprehended by visualizing them on diagrams of individual metabolic pathways or of metabolic maps. MetaFlux can also remove redundant high-flux loops, solve FBA models once they are generated and model the effects of gene knockouts. MetaFlux has been validated through construction of FBA models for Escherichia coli and Homo sapiens.Availability: Pathway Tools with MetaFlux is freely available to academic users, and for a fee to commercial users. Download from: biocyc.org/download.shtml.Contact: mario.latendresse@sri.comSupplementary information: Supplementary data are available at Bioinformatics online.

Bioenergetics of Monoterpenoid Essential Oil Biosynthesis in Nonphotosynthetic Glandular Trichomes

“…Additional metabolic reactions were extracted from the AraCyc and MetaCyc databases (Mueller et al, 2003;Caspi et al, 2016). A process termed instantiation (Latendresse et al, 2012) was employed, using an implementation of the relevant algorithm in YASMEnv, to ensure that only mass-balanced reactions were included in the model. Each reaction was associated with a subcellular compartment based on the records in the Arabidopsis Core Model, AraCyc, and MetaCyc, and these associations were confirmed by searches against the TargetP and SUBAcon databases (Emanuelsson et al, 2000;Hooper et al, 2014).…”

High-throughput comparison, functional annotation, and metabolic modeling of plant genomes using the PlantSEED resource

“…Through a process called "gapfilling," models rapidly identify when these criteria are not satisfied and propose reactions, annotations, and genes that may be associated with missing pathway steps. This type of gapfilling is offered in other annotation and modeling environments, including the Constraints-Based Reconstruction and Analysis (COBRA) toolbox (45) and Pathway Tools (46), although PlantSEED is distinct in offering model-driven gapfilling of plant genomes based on a refined and curated database of plant biochemistry.…”

Bioinformatics computation of metabolic models from sequenced genomes

“…An important component of the model generation step is gap filling (Latendresse, et al, 2012;Satish Kumar, et al, 2007). Genome-based metabolic network reconstructions are usually missing reactions because of incompleteness and inaccuracy in all of the preceding steps; a single missing reaction can prevent the network from producing one or more components of the organism's biomass.…”