Toward the automated generation of genome-scale metabolic networks in the SEED
© DeJongh et al; licensee BioMed Central Ltd. 2007
Received: 22 November 2006
Accepted: 26 April 2007
Published: 26 April 2007
Current methods for the automated generation of genome-scale metabolic networks focus on genome annotation and preliminary biochemical reaction network assembly, but do not adequately address the process of identifying and filling gaps in the reaction network, and verifying that the network is suitable for systems level analysis. Thus, current methods are only sufficient for generating draft-quality networks, and refinement of the reaction network is still largely a manual, labor-intensive process.
We have developed a method for generating genome-scale metabolic networks that produces substantially complete reaction networks, suitable for systems level analysis. Our method partitions the reaction space of central and intermediary metabolism into discrete, interconnected components that can be assembled and verified in isolation from each other, and then integrated and verified at the level of their interconnectivity. We have developed a database of components that are common across organisms, and have created tools for automatically assembling appropriate components for a particular organism based on the metabolic pathways encoded in the organism's genome. This focuses manual efforts on that portion of an organism's metabolism that is not yet represented in the database. We have demonstrated the efficacy of our method by reverse-engineering and automatically regenerating the reaction network from a published genome-scale metabolic model for Staphylococcus aureus. Additionally, we have verified that our method capitalizes on the database of common reaction network components created for S. aureus, by using these components to generate substantially complete reconstructions of the reaction networks from three other published metabolic models (Escherichia coli, Helicobacter pylori, and Lactococcus lactis). We have implemented our tools and database within the SEED, an open-source software environment for comparative genome annotation and analysis.
Our method sets the stage for the automated generation of substantially complete metabolic networks for over 400 complete genome sequences currently in the SEED. With each genome that is processed using our tools, the database of common components grows to cover more of the diversity of metabolic pathways. This increases the likelihood that components of reaction networks for subsequently processed genomes can be retrieved from the database, rather than assembled and verified manually.
The availability of hundreds of sequenced genomes has ushered in a new era in biology, allowing the study of cellular life at a systems level. One approach to systems level understanding of cellular life is in silico modeling of an organism's metabolic capabilities, as determined by the complement of genes in its genome . For example, flux balance analysis is a widely used modeling technique which enables the prediction of metabolic phenotypes based on the enzymes encoded in an organism's genome [2–5]. The model consists in part of a network of biochemical reactions that represent the activity of these enzymes. The usefulness of the model is subject to the accuracy of the reaction network upon which it is based: the reaction network should be complete, fully covering the metabolic capabilities that are to be modeled, coherent, containing no gaps or dead ends, and correct, faithfully representing the metabolic phenotype of the organism. The accuracy of the reaction network can be tested by comparing the predictions of the model with the known metabolic phenotype of the organism under (pseudo) steady-state conditions.
A major challenge of genome-scale metabolic modeling is to reconstruct an accurate reaction network directly from an annotated genome . There are four steps required to generate an accurate reaction network for an organism:
1. Annotating the genome to identify the enzymes encoded by particular genes.
2. Assembling the network of reactions that correspond to these enzymes.
3. Verifying the completeness and coherence of the assembled reaction network.
4. Testing the correctness of the assembled reaction network using a modeling technique such as flux balance analysis.
Few techniques have been proposed for automating the verification of an assembled reaction network (Fig. 1A). Segrè et al.  describe an algorithm that searches the reaction network to identify essential intermediate metabolites that cannot be synthesized by other reactions in the network. Arakawa et al.  apply a gap-filling heuristic to connect segments of pathways that represent consecutive reactions, and flag the filled gap for subsequent manual curation. Other systems require the user to identify gaps, and assist in this process by providing tools to visualize genome annotations in the context of metabolic pathway maps [13–15]. This information can be used to improve the accuracy of the reaction network, by feeding back into the annotation and assembly steps for iterative refinement.
Only a handful of metabolic reconstructions suitable for modeling have been published to date [20–29]. With the completion of the 1000th microbial genome predicted for late 2007 , there is an urgent need for more substantial automation of the process of generating complete and coherent reaction networks suitable for testing at the systems level.
We have developed an integrated suite of tools that supports iterative annotation, assembly, verification and refinement of genome-scale metabolic reaction networks. Our approach is based on partitioning central and intermediary metabolism into discrete, interconnected components that are shared across organisms (Fig. 1B). Each of these components represents a subnetwork of reactions within an overall metabolic reaction network. In our approach, each reaction subnetwork is assembled and verified independently, and stored in a database so that it can be used for multiple organisms. An organism-specific reaction network is assembled by retrieving appropriate reaction subnetworks from the database and verifying their coherence at the level of their connectivity. Any effort required to assemble and verify reaction subnetworks for components that are not already in the database has a cumulative effect, because those reaction subnetworks subsequently can be stored in the database and used for other organisms.
Software tools developed for genome-scale metabolic reconstruction
Purpose of Tools
Implementation in the SEED
The SEED already provides tools for annotation, based on similarity searching and context-based methods
The SEED already provides a database of high-quality genome annotations organized into subsystems (see )
SUB-ASSEMBLY AND SUB-NETWORK VERIFICATION
Curating associations between functional roles and reactions in a particular metabolic context
Reverse-engineering of published genome-scale metabolic models; analysis of gene-reaction associations in the KEGG database; integrated display of KEGG pathway maps in subsystems, highlighting functional roles and associated reactions
Associations between functional roles and KEGG reactions in subsystems
Assembling and verifying the coherence of reaction subnetworks in subsystems
Petri net representation of KEGG reactions; encoded scenarios in subsystems; finding paths through reaction subnetworks from scenario inputs to scenario outputs
Reuseable coherent reaction subnetworks in subsystems
Assembling and verifying the coherence of connected reaction subnetworks across subsystems
Connections between scenarios in different subsystems; finding paths through connected scenarios, from overall inputs to overall outputs
List of curated subsystems with coherent reaction subnetworks for functional variants that interconnect to cover central and intermediary metabolic pathways
ASSEMBLY AND NETWORK VERIFICATION
Assembling and verifying the coherence and completeness of an organism-specific reaction network
Identifying gaps in the reaction network, by cross-checking inputs and outputs for all paths through implemented scenarios, and checking for paths from minimal substrates to biomass compounds; creating files for FluxAnalyzer 
Organism-specific complete and coherent reaction networks for central and intermediary metabolism
Extending subsystems to represent reaction subnetworks
As mentioned above, subsystems within the SEED already have the capability of associating reaction information with functional roles. The reaction information, when it is present, is specified by links to reaction identifiers in the KEGG database. Prior to our work, only a small percentage of functional roles in the SEED had reaction information associated with them. To address this shortcoming, we have reverse-engineered existing genome-scale metabolic models to determine the correspondence between the reactions in these models and functional roles in subsystems (see Methods).
We have created encodings of segments of metabolic pathways within subsystems, which we call scenarios (Fig. 2B). A scenario represents a set of connected reactions that convert a defined set of substrates (scenario inputs) to a defined set of products (scenario outputs). For example, the Embden-Meyerhof and Gluconeogenesis subsystem contains functional roles and associated reactions for the metabolic pathway of glycolysis. We have encoded a Glycolysis scenario within this subsystem, for which the inputs are alpha- and beta-D-glucose, and the output is pyruvate. The input and output compounds for scenarios in subsystems are encoded using KEGG compound identifiers (e.g., C00267 for alpha-D-glucose). A scenario additionally encodes KEGG pathway map identifiers that provide context for the set of connected reactions that convert the scenario inputs to the scenario outputs. The pathway maps contain information about reaction reversibilities, and identify the "main" compounds in each reaction (e.g., the molecules derived from the breakdown of glucose), as opposed to cofactors. A scenario does not explicitly encode this set of connected reactions – it is assembled dynamically from the reactions associated with functional roles in the subsystem, using a process described below to search for paths through the reactions from the scenario inputs to the scenario outputs. However, a scenario can specify additional reactions that are not associated with functional roles, yet play a part in the scenario (e.g., ionization reactions and spontaneous reactions that are not catalyzed by an enzyme).
Curating associations between functional roles and reactions in subsystems
The first step in ensuring that a subsystem contains a coherent subnetwork of reactions is to associate appropriate reactions with the functional roles in the subsystem. Functional roles that correspond to enzymes usually specify an EC number. As noted above, an EC number often designates a set of related reactions; however, a particular organism's gene product may not implement all of the reactions associated with a given EC number. In the SEED, a functional role can be associated with multiple subsystems, and a different set of reactions can be specified for the functional role in each of these subsystems. Thus an organism's gene product will only be associated with the reactions in a particular subsystem if the gene is linked to the associated functional role by that subsystem's spreadsheet. The metabolic pathways represented by each subsystem provide context for determining which reactions should be associated with a functional role in the subsystem. For example, EC 184.108.40.206 (alcohol dehydrogenase) is associated with 12 reactions in the KEGG database. These 12 reactions have a variety of substrates, including all primary and secondary alcohols. The KEGG entry for E. coli K12 gene b1241  annotates it with EC 220.127.116.11, but does not identify its substrate specificity. The SEED entry for this gene includes annotation with the functional role "Alcohol dehydrogenase (EC 18.104.22.168)," and shows that it is present in several subsystems, including Fermentations: Lactate, for which the associated reaction is R00754 (substrate is ethanol), and Glycerolipid and Glycerophospholipid metabolism, for which the associated reaction is R01036 (substrate is glycerol).
To determine which reactions to associate with a functional role in a particular subsystem, we have integrated the display of KEGG pathway maps into the subsystem environment. KEGG pathway maps display metabolic pathways as graphs, where the nodes represent compounds and the edges represent reactions. Edges that represent enzymatically catalyzed reactions are drawn with a box containing the EC number corresponding to the enzyme. KEGG provides the capability of displaying maps with a particular set of edges highlighted, either by specifying the corresponding EC numbers or the corresponding reactions. We have developed a tool which searches for any KEGG map containing one or more of the EC numbers specified by the functional roles in the subsystem, and displays links to the maps in order of decreasing number of EC matches. For each map, two links are displayed: one that highlights the EC numbers that were matched, and one that highlights the reactions associated with functional roles in the subsystem that were matched. These displays are useful for visually determining the extent to which the functional roles and reactions in the subsystem form a connected and complete subnetwork. We use this tool to identify gaps in the reaction subnetwork, which we fill by adding functional roles to the subsystem corresponding to unmatched EC numbers, and by adding reaction associations to appropriate functional roles for unmatched reactions.
Assembling and verifying reaction subnetworks in subsystems
In order to determine which scenarios, if any, are implemented by a given organism, we filter the set of reactions used by the path-finding tool according to the functional roles annotated to the organism's genome in the subsystem's spreadsheet (Fig. 4B). Gaps in the organism-specific reaction subnetwork indicate that further genome annotation is required to determine whether genes responsible for the missing reactions are present. If these genes are identified and added to the spreadsheet, subsequent reaction subnetwork assembly and verification will result in a coherent organism-specific reaction subnetwork with paths through the scenario. In this case, the existence of paths through the scenario is ensured not only for this organism, but also for any other organism which corresponds to the same functional variant, i.e., whose genome is annotated with the same set of functional roles.
The reaction subnetworks encoded in subsystems form the database of reusable metabolic network components (Fig. 1B). By verifying that the reaction subnetwork in a given subsystem is coherent for all the scenarios implemented by a particular functional variant, we have effectively verified the coherence of the reaction subnetwork for any organism that corresponds to that functional variant. Thus, the task of assembling a coherent reaction subnetwork from the subsystem for a given organism is accomplished by determining the organism's corresponding functional variant in the spreadsheet, and assembling the reactions associated with the corresponding set of functional roles. Then the path-finding tool can be applied to determine which scenarios are implemented by that functional variant.
Assembling and verifying connected reaction subnetworks
As described above, reaction subnetworks across related subsystems are connected by shared scenario inputs and outputs (Fig. 3A). We assemble connected reaction subnetworks by collecting all of the reactions associated with a specified list of subsystems. We use the same path-finding tool to search for paths through connected reaction subnetworks; however, instead of constructing the Petri net using KEGG reactions from the subsystems, we represent each possible path through a given scenario as a distinct reaction. These "higher-order" reactions have the path's input compounds as substrates and the path's output compounds as products, and use the cumulative stoichiometry that was previously determined by the path-finding tool. For example, the path through the Gluconeogenesis scenario shown in Fig. 4A can be represented as a higher-order reaction that converts 2 molecules of pyruvate (plus cofactors) to 1 molecule of alpha-D-Glucose 6-phosphate (plus cofactors). By matching the outputs of one higher-order reaction to the inputs of another, the path-finding tool pushes tokens from a set of overall input compounds to a set of overall output compounds. For example, the PEP (phosphoenolpyruvate) generation scenario in the Embden-Meyerhof and Gluconeogenesis subsystem and the Erythrose4P (erythrose-4-phosphate) generation scenario in the Pentose Phosphate Pathway subsystem generate the scenario inputs for the Chorismate synthesis scenario in the Chorismate Synthesis subsystem. The output of the Chorismate synthesis scenario in turn serves as input to the scenarios for the synthesis of the aromatic amino acids. The Petri net tool can be used to find paths through these scenarios that convert glucose to tryptophan, phenylalanine, and tyrosine on both an organism-independent and organism-specific basis (Fig. 3B). As before, dead end tokens signal gaps in the assembled network, which in this case represent missing scenarios.
Assembling and verifying an organism-specific reaction network
To assemble an organism-specific reaction network, we collect all of the reactions associated with the organism's genome annotation in a specified list of subsystems. To verify that this reaction network is coherent and complete, we use our path-finding tool to determine which of the scenarios encoded in the subsystems are implemented by the organism. We cross-check the input and output compounds for all of these paths to determine which input compounds are not outputs for any of the paths, and which output compounds are not inputs for any of the paths. Ideally, these compounds correspond respectively to potential reaction network inputs (e.g., sugars that are transported into the cell by the organism) and outputs (e.g., amino acids that are synthesized by the organism for biomass growth). If any of these compounds represent intermediate steps in metabolism, the subsystems need to be checked for scenarios that produce or consume them accordingly, and the organism's annotation in these subsystems needs to be checked with respect to the functional roles necessary to implement the missing scenarios.
An organism-specific reaction network that is suitable for flux balance analysis must additionally specify transport and exchange reactions, a biomass equation, and a set of minimal substrates that can be used to generate the biomass compounds for the organism. Currently, these additional components must be supplied by hand for a given organism. We have written a tool that, given these supplied components and subsystem reactions, assembles them into a final reaction network, and outputs the network in a format suitable for loading into the FluxAnalyzer . To verify the completeness of the final reaction network, the tool determines whether there is a path through the network from the set of minimal substrates to each biomass compound. This tool uses a breadth-first search algorithm (similar to the algorithm described in ), starting with the minimal substrates, to determine all the ways that paths through scenarios can be connected to form paths to the biomass compounds. It "borrows" compounds along the way, if necessary, to find complete paths through the scenarios, then reports on the paths it found and the compounds it had to borrow. Any borrowed compounds must be accounted for, either by adding more transport reactions or creating new scenarios to produce them. The final step is to load the network into the FluxAnalyzer, supply minimal substrates, and verify flux through the biomass reaction.
Reverse-engineering and regenerating the S. aureus metabolic reaction network
Regenerated S. aureus network
Reactions Automatically assembled from Subsystems
Reactions Added "As is" from iSB619 Modela
Transport Reactions from iSB619 Model
Exchange Reactions for Transported Compounds
The cumulative effect of the database of coherent reaction subnetworks
Although our initial efforts have focused on generating the metabolic reaction network from the iSB619 S. aureus model, in principle our tools can produce metabolic reaction networks for any organism in the SEED, by reusing components from the database of coherent reaction subnetworks (Fig. 1B). To test this hypothesis, we have quantified the extent to which the metabolic reaction networks of three other published genome-scale metabolic models (the iJR904 E. coli model , the iIT341 H. pylori model , and the Oliveira, et al., L. lactis model ) can be assembled from coherent reaction subnetworks already in the database, based solely on their genome annotations in the 65 subsystems used to assemble the S. aureus metabolic reaction network.
Capitalizing on common aspects of metabolism: reuse of scenarios
S. aureus a
E. coli a
H. pylori a
L. lactis a
Vitamins and Cofactors
The iJR904 E. coli model is the most complete genome-scale metabolic model published to date, and we are currently applying our tools to regenerate the full reaction network in the SEED. To demonstrate further the cumulative effect of our approach to reaction network generation, we have compared the reactions from the three additional models that were not mapped to functional roles in the 65 subsystems, and determined the overlapping set of reactions from the H. pylori and L. lactis models that will be mapped to functional roles in other subsystems in the process of regenerating the E. coli model in the SEED. In the case of the iIT341 H. pylori model, there will only be 20 remaining reactions to map to functional roles in subsystems. In the case of the Oliveira, et al., L. lactis model, there will only be 55 remaining reactions to map to functional roles (Fig. 6B and Additional File 3).
Although our work to date has focused on regenerating reaction networks from published genome-scale metabolic models, our tools can be used to generate a preliminary reaction network for any organism in the SEED based solely on its genome annotation (Fig. 1B). As with other methods for generating genome-scale metabolic reaction networks [12–15, 18, 19], the process of refining the reaction network until it is complete and coherent requires manual effort. The amount of manual effort required depends upon the quality of the preliminary reaction network, which in turn depends upon the quality of the annotation. An underlying goal behind our suite of tools is to narrow the gap between the preliminary reaction network and the complete and coherent reaction network as much as possible, so as to minimize the amount of manual effort required. We anticipate that the database of coherent reaction subnetworks will enable the generation of substantially complete and coherent preliminary reaction networks, focusing manual efforts on resolving gaps that are identified by the path-finding tool (Fig. 4B), and on creating new reaction subnetworks for areas of metabolism not yet represented in the database.
The collection of transport and biomass information from the published physiological literature and addition of corresponding reactions to the reaction network requires a significant manual effort for each organism, and is currently a bottleneck in generating a reaction network suitable for flux balance analysis. To date, this problem has not been adequately addressed in systems designed to facilitate genome-scale reaction network reconstruction [12–15, 18, 19], nor have we addressed this problem with our current set of tools. However, it may be possible to hypothesize computationally which transport reactions an organism uses, and which biomass components an organism synthesizes, based on the existence of paths through scenarios that represent catabolic pathways for common substrates, or anabolic pathways for common biomass components.
Our approach to automated genome-scale metabolic reaction network generation represents two important advances when compared to previously published methods [12–15, 18, 19]. The first advance is a process for verifying the completeness and coherence of an overall reaction network by constructing a database of coherent reaction subnetworks that represent interconnected metabolic components. In contrast, most published methods for network generation do not explicitly provide support for the verification process. IdentiCS  and metaSHARK  provide automated annotation, assembly and visualization of a preliminary reaction network. The GEM System  additionally provides a heuristic for filling gaps in sets of consecutive reactions. However, these three systems do not provide tools for iterative refinement and verification of the preliminary network. The Pathway Tools software  takes an existing annotation and produces a pathway/genome database (PGDB), which includes predicted metabolic pathways and associated reactions for an organism. The software enables the visualization and manual refinement of the preliminary reaction network, but does not include tools for verification of the network. Segrè et al.  describe a process that builds upon Pathway Tools and includes an algorithm for the verification of the completeness and coherence of the overall reaction network, but do not describe a process for resolving gaps that are identified in the network. As discussed above, our approach enables resolving such gaps during the process of creating and assembling coherent reaction subnetworks (Fig. 1B). The AUTOGRAPH-method  is unique among previously published methods in that it incorporates information from published genome-scale reaction networks to produce a preliminary reaction network for a particular organism. However, no process for verifying and refining the preliminary network is described.
The second advance is the tight integration of our approach and our tools with a community-based genome annotation and analysis tool. From its inception, the SEED was designed to serve as a repository and clearinghouse for parallel annotation projects across all sequenced genomes. Because of the tight integration of our approach with the SEED, our tools can be used at all stages of every genome annotation project in the SEED, and reaction subnetworks created for one project are immediately available for all the other projects. This is in contrast with, e.g., the Pathway Tools software , which is downloaded and installed locally by each genome annotation project, and creates a standalone PGDB for each project. Likewise, the published methods discussed above all focus on creating a reaction network for a single organism and, with the exception of the AUTOGRAPH-method, do not describe a process for reusing components of reaction networks previously developed for other organisms.
We have described a method for automating the generation of substantially complete and coherent genome-scale metabolic reaction networks from annotated genomes. Our method builds on the subsystems approach to genome annotation and analysis embodied in the SEED. The SEED already provides well-curated genome annotations for central and intermediary metabolism across many organisms. We have extended the SEED by curating associations between reactions and functional roles in subsystems based on metabolic context. We have created tools for encoding components of reaction networks in subsystems, and verifying their coherence and interconnections. We have created tools for assembling these components into organism-specific complete and coherent reaction networks. We have demonstrated that our process can regenerate the reaction network from a published genome-scale metabolic model, and that it produces a cumulative effect supporting the subsequent generation of other reaction networks from published models. Our future work will focus on applying this process to generating reaction networks for new organisms, eventually extending to all organisms annotated in the SEED, thus producing a repository of organism-specific complete and coherent reaction networks. We envision that this repository will be useful for interpretation of large-scale data sets generated for metabolic genomics .
Reverse-engineering existing genome-scale metabolic models into the SEED
Each of the four genome-scale metabolic models [20, 24, 27, 28] that we reverse-engineered contains a list of compounds, a list of reactions, and a list of gene-reaction associations. The list of reactions can be broken into three categories: (1) reactions that take place within the cell, (2) transport reactions, which move compounds into and out of the cell, and (3) exchange reactions, which simulate compounds entering and leaving the extra-cellular space. Our reverse-engineering process focused on mapping the first category of reactions to KEGG reactions, and associating these with functional roles in subsystems.
Our reverse-engineering process started with mapping the compounds in the existing models to compounds in KEGG. Much of this could be done automatically, since the models for the most part use standard compound names and abbreviations that are also used by KEGG. In some cases, a model compound may map to more than one KEGG compound id, e.g., D-Glucose can be mapped to C00031 (also D-Glucose), as well as C00221 (beta-D-Glucose) and C00267 (alpha-D-Glucose). We retained all possible mappings because different KEGG reactions use different compound ids. The next step was to map model reactions to KEGG reactions. Our criterion was a perfect match on all substrates and products, except for protons, since the models and KEGG sometimes account for them differently (e.g., nitrite in the models is NO2, whereas in KEGG it is HNO2). Again, much of this could be done automatically. For model reactions that did not perfectly match a KEGG reaction, we searched by hand for similar KEGG reactions, looking for cofactor differences, equivalent combinations of reactions, and so on. In these cases we accepted an inexact match, but noted the discrepancy. In a few cases we were unable to map model reactions to any KEGG reaction. For example, the models do not contain reactions representing the individual steps in the synthesis of specific fatty acids, but instead represent the cumulative synthesis of specific fatty acids as one "reaction". In these cases we did not process the reactions any further, but noted that they should be added "as is" to the organism's metabolic reaction network. This required the specification of several more reactions to convert between KEGG compound ids and the model compound ids in these added reactions. Additional Files 4 and 5 document the mapping from compounds and reactions in the iSB619 S. aureus model to KEGG compounds and reactions.
When mapping model reactions from the iSB619 S. aureus model to KEGG reactions, we noted that approximately 40% of the reactions had different reversibilities between these two sources (e.g., a reaction asserted to be one-way in the model is represented as bidirectional in KEGG). Some of these discrepancies are clearly errors in the KEGG database (e.g., the KEGG pathway map agrees with the model whereas the KEGG flat file does not); we have brought these to the attention of the KEGG database curators, and they have been fixed. Efforts are underway to correct and improve upon the compound and reaction information in KEGG (e.g., ). In light of this, we have structured our database so that we have the option of using reaction information from sources other than KEGG.
Once we completed mapping the model reactions to KEGG reactions, we used the gene-reaction associations from each of the models to associate KEGG reactions with functional roles in SEED subsystems. For each gene in a model, we identified all subsystems (if any) where the gene is annotated to a functional role. Most of these functional roles specify an EC number; in some cases, they have also been assigned KEGG reactions by the original subsystem author. In addition, we retrieved the KEGG reaction(s) annotated to the gene from the KEGG GENE database. Thus, for each model, we collected up to three lines of evidence for determining whether to assign a given KEGG reaction to a functional role: the gene-reaction association and corresponding KEGG reaction from a model; KEGG's gene-reaction association for the organism; and the functional role's EC number and KEGG reaction. We recorded this evidence in the SEED database for subsequent manual curation. In addition, if all three lines of evidence matched, we automatically assigned the KEGG reaction to the functional role, and noted in the SEED database that an automatic assignment was made. At the end of this process we had a preliminary set of 805 associations between KEGG reactions and functional roles in 205 subsystems, with links back to the gene-reaction associations from the four models.
Manual curation of associations between functional roles and reactions
After reverse-engineering the four models into the SEED, we manually curated associations between functional roles and reactions by checking every automatic assignment of reactions and all links to reaction evidence. When the three lines of evidence for assigning reactions to functional roles were inconclusive, we examined the evidence in the context of the appropriate KEGG pathway to make a determination. For example, in the Embden-Meyerhof and Gluconeogenesis subsystem, the subsystem author assigned reaction R00299 (ATP + D-Glucose <=> ADP + D-Glucose 6-phosphate) to the Glucokinase (EC 22.214.171.124) functional role. One of the models, L. lactis, matched this reaction, but the other three models and KEGG all matched R01600 (ATP + beta-D-Glucose <=> ADP + beta-D-Glucose 6-phosphate). After consulting the Glycolysis/Gluconeogenesis map in KEGG, which represents the catabolism of both the alpha- and beta- forms of glucose, we decided that it was appropriate to assign both reactions to the functional role, as well as R01786 (ATP + alpha-D-Glucose <=> ADP + alpha-D-Glucose 6-phosphate), so that all three forms of glucose would be represented in the reaction network generated from this subsystem. In some cases, we added reactions to functional roles based on evidence from the functional roles and the KEGG pathways alone, such as R01450 for Lactate racemase (EC 126.96.36.199), which is in the Fermentations: Lactate subsystem and the Pyruvate Metabolism KEGG pathway, but not associated with a gene in any of the four models (it is in the iIT341 model, without a gene-reaction association).
Availability of code and database
We are currently working with the SEED development group to incorporate our extensions to the SEED into the standard open-source SEED distribution, available at . Our associations between reactions and functional roles have already been added to subsystems in this distribution. We anticipate that our tools will be available through the SEED in May of 2007.
The authors would like to thank Ross Overbeek, Rick Stevens, Mike Kubal and Nathan Tintle for helpful discussions and insightful ideas. The authors also thank the Fellowship for Interpretation of Genomes, the SEED annotators, and the SEED development group for providing access to the SEED and support in development within this environment. AB and MD thank the following members of the Hope College Bioinformatics class, Spring 2006, for their contributions: John Gonthier and Tyson VonDerfecht. This research was supported in part by a grant to Hope College from the Howard Hughes Medical Institute through the Undergraduate Science Education Program, and Research Experience for Undergraduates grant number 0353566 from the National Science Foundation to the Hope College Computer Science Department. MD received support from the Argonne National Laboratory Guest Faculty Research Participation Program.
- Di Ventura B, Lemerle C, Michalodimitrakis K, Serrano L: From in vivo to in silico biology and back. Nature 2006, 443(7111):527–533.View ArticlePubMedGoogle Scholar
- Schilling CH, Edwards JS, Palsson BO: Toward metabolic phenomics: analysis of genomic data using flux balances. Biotechnol Prog 1999, 15(3):288–295.View ArticlePubMedGoogle Scholar
- Kauffman KJ, Prakash P, Edwards JS: Advances in flux balance analysis. Curr Opin Biotechnol 2003, 14(5):491–496.View ArticlePubMedGoogle Scholar
- Lee JM, Gianchandani EP, Papin JA: Flux balance analysis in the era of metabolomics. Brief Bioinform 2006, 7(2):140–150.View ArticlePubMedGoogle Scholar
- Price ND, Reed JL, Palsson BO: Genome-scale models of microbial cells: evaluating the consequences of constraints. Nat Rev Microbiol 2004, 2(11):886–897.View ArticlePubMedGoogle Scholar
- Palsson B: Two-dimensional annotation of genomes. Nat Biotechnol 2004, 22(10):1218–1219.View ArticlePubMedGoogle Scholar
- Covert MW, Schilling CH, Famili I, Edwards JS, Goryanin, Selkov E, Palsson BO: Metabolic modeling of microbial strains in silico. Trends Biochem Sci 2001, 26(3):179–186.View ArticlePubMedGoogle Scholar
- Francke C, Siezen RJ, Teusink B: Reconstructing the metabolic network of a bacterium from its genome. Trends Microbiol 2005, 13(11):550–558.View ArticlePubMedGoogle Scholar
- Green ML, Karp PD: A Bayesian method for identifying missing enzymes in predicted metabolic pathway databases. BMC Bioinformatics 2004, 5: 76.PubMed CentralView ArticlePubMedGoogle Scholar
- Osterman A, Overbeek R: Missing genes in metabolic pathways: a comparative genomics approach. Curr Opin Chem Biol 2003, 7(2):238–251.View ArticlePubMedGoogle Scholar
- Overbeek R, Begley T, Butler RM, Choudhuri JV, Chuang HY, Cohoon M, de Crecy-Lagard V, Diaz N, Disz T, Edwards R, Fonstein M, Frank ED, Gerdes S, Glass EM, Goesmann A, Hanson A, Iwata-Reuyl D, Jensen R, Jamshidi N, Krause L, Kubal M, Larsen N, Linke B, McHardy AC, Meyer F, Neuweger H, Olsen G, Olson R, Osterman A, Portnoy V, Pusch GD, Rodionov DA, Ruckert C, Steiner J, Stevens R, Thiele I, Vassieva O, Ye Y, Zagnitko O, Vonstein V: The subsystems approach to genome annotation and its use in the project to annotate 1000 genomes. Nucleic Acids Res 2005, 33(17):5691–5702.PubMed CentralView ArticlePubMedGoogle Scholar
- Arakawa K, Yamada Y, Shinoda K, Nakayama Y, Tomita M: GEM System: automatic prototyping of cell-wide metabolic pathway models from genomes. BMC Bioinformatics 2006, 7: 168.PubMed CentralView ArticlePubMedGoogle Scholar
- Karp PD, Paley S, Romero P: The Pathway Tools software. Bioinformatics 2002, 18 Suppl 1: S225–32.View ArticlePubMedGoogle Scholar
- Pinney JW, Shirley MW, McConkey GA, Westhead DR: metaSHARK: software for automated metabolic network prediction from DNA sequence and its application to the genomes of Plasmodium falciparum and Eimeria tenella . Nucleic Acids Res 2005, 33(4):1399–1409.PubMed CentralView ArticlePubMedGoogle Scholar
- Sun J, Zeng AP: IdentiCS--identification of coding sequence and in silico reconstruction of the metabolic network directly from unannotated low-coverage bacterial genome sequence. BMC Bioinformatics 2004, 5: 112.PubMed CentralView ArticlePubMedGoogle Scholar
- Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M: The KEGG resource for deciphering the genome. Nucleic Acids Res 2004, 32(Database issue):D277–80.PubMed CentralView ArticlePubMedGoogle Scholar
- Krieger CJ, Zhang P, Mueller LA, Wang A, Paley S, Arnaud M, Pick J, Rhee SY, Karp PD: MetaCyc: a multiorganism database of metabolic pathways and enzymes. Nucleic Acids Res 2004, 32(Database issue):D438–42.PubMed CentralView ArticlePubMedGoogle Scholar
- Notebaart RA, van Enckevort FH, Francke C, Siezen RJ, Teusink B: Accelerating the reconstruction of genome-scale metabolic networks. BMC Bioinformatics 2006, 7: 296.PubMed CentralView ArticlePubMedGoogle Scholar
- Segre D, Zucker J, Katz J, Lin X, D'Haeseleer P, Rindone WP, Kharchenko P, Nguyen DH, Wright MA, Church GM: From annotated genomes to metabolic flux models and kinetic parameter fitting. Omics 2003, 7(3):301–316.View ArticlePubMedGoogle Scholar
- Becker SA, Palsson BO: Genome-scale reconstruction of the metabolic network in Staphylococcus aureus N315: an initial draft to the two-dimensional annotation. BMC Microbiol 2005, 5(1):8.PubMed CentralView ArticlePubMedGoogle Scholar
- Duarte NC, Herrgard MJ, Palsson BO: Reconstruction and validation of Saccharomyces cerevisiae iND750, a fully compartmentalized genome-scale metabolic model. Genome Res 2004, 14(7):1298–1309.PubMed CentralView ArticlePubMedGoogle Scholar
- Feist AM, Scholten JC, Palsson BO, Brockman FJ, Ideker T: Modeling methanogenesis with a genome-scale metabolic reconstruction of Methanosarcina barkeri . Mol Syst Biol 2006, 2: 2006.0004.PubMed CentralView ArticlePubMedGoogle Scholar
- Heinemann M, Kummel A, Ruinatscha R, Panke S: In silico genome-scale reconstruction and validation of the Staphylococcus aureus metabolic network. Biotechnol Bioeng 2005, 92(7):850–864.View ArticlePubMedGoogle Scholar
- Reed JL, Vo TD, Schilling CH, Palsson BO: An expanded genome-scale model of Escherichia coli K-12 (iJR904 GSM/GPR). Genome Biol 2003, 4(9):R54.PubMed CentralView ArticlePubMedGoogle Scholar
- Schilling CH, Covert MW, Famili I, Church GM, Edwards JS, Palsson BO: Genome-scale metabolic model of Helicobacter pylori 26695. J Bacteriol 2002, 184(16):4582–4593.PubMed CentralView ArticlePubMedGoogle Scholar
- Teusink B, van Enckevort FH, Francke C, Wiersma A, Wegkamp A, Smid EJ, Siezen RJ: In silico reconstruction of the metabolic pathways of Lactobacillus plantarum : comparing predictions of nutrient requirements with those from growth experiments. Appl Environ Microbiol 2005, 71(11):7253–7262.PubMed CentralView ArticlePubMedGoogle Scholar
- Thiele I, Vo TD, Price ND, Palsson BO: Expanded metabolic reconstruction of Helicobacter pylori (iIT341 GSM/GPR): an in silico genome-scale characterization of single- and double-deletion mutants. J Bacteriol 2005, 187(16):5818–5830.PubMed CentralView ArticlePubMedGoogle Scholar
- Oliveira AP, Nielsen J, Forster J: Modeling Lactococcus lactis using a genome-scale flux model. BMC Microbiol 2005, 5: 39.PubMed CentralView ArticlePubMedGoogle Scholar
- Forster J, Famili I, Fu P, Palsson BO, Nielsen J: Genome-scale reconstruction of the Saccharomyces cerevisiae metabolic network. Genome Res 2003, 13(2):244–253.PubMed CentralView ArticlePubMedGoogle Scholar
- The SEED[http://www.theseed.org]
- Kyoto Encyclopedia of Genes and Genomes[http://www.genome.jp/dbget-bin/www_bget?eco+b1241]
- Hardy S, Robillard PN: Modeling and simulation of molecular biology systems using petri nets: modeling goals of various approaches. J Bioinform Comput Biol 2004, 2(4):595–613.View ArticlePubMedGoogle Scholar
- Koch I, Junker BH, Heiner M: Application of Petri net theory for modelling and validation of the sucrose breakdown pathway in the potato tuber. Bioinformatics 2005, 21(7):1219–1226.View ArticlePubMedGoogle Scholar
- Pinney JW, Westhead DR, McConkey GA: Petri Net representations in systems biology. Biochem Soc Trans 2003, 31(Pt 6):1513–1515.View ArticlePubMedGoogle Scholar
- Reddy VN, Liebman MN, Mavrovouniotis ML: Qualitative analysis of biochemical reaction systems. Comput Biol Med 1996, 26(1):9–24.View ArticlePubMedGoogle Scholar
- Klamt S, Stelling J, Ginkel M, Gilles ED: FluxAnalyzer: exploring structure, pathways, and flux distributions in metabolic networks on interactive flux maps. Bioinformatics 2003, 19(2):261–269.View ArticlePubMedGoogle Scholar
- Systems Biology Research Group[http://gcrg.ucsd.edu/organisms/staph.html]
- Chang DE, Conway T: Metabolic genomics. Adv Microb Physiol 2005, 50: 1–39.View ArticlePubMedGoogle Scholar
- Ott MA, Vriend G: Correcting ligands, metabolites, and pathways. BMC Bioinformatics 2006, 7: 517.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.