MANET: tracing evolution of protein architecture in metabolic networks
© Kim et al; licensee BioMed Central Ltd. 2006
Received: 13 April 2006
Accepted: 19 July 2006
Published: 19 July 2006
Cellular metabolism can be characterized by networks of enzymatic reactions and transport processes capable of supporting cellular life. Our aim is to find evolutionary patterns and processes embedded in the architecture and function of modern metabolism, using information derived from structural genomics.
The Molecular Ancestry Network (MANET) project traces evolution of protein architecture in biomolecular networks. We describe metabolic MANET, a database that links information in the Structural Classification of Proteins (SCOP), the Kyoto Encyclopedia of Genes and Genomes (KEGG), and phylogenetic reconstructions depicting the evolution of protein fold architecture. Metabolic MANET literally 'paints' the ancestries of enzymes derived from rooted phylogenomic trees directly onto over one hundred metabolic subnetworks, enabling the study of evolutionary patterns at global and local levels. An initial analysis of painted subnetworks reveals widespread enzymatic recruitment and an early origin of amino acid metabolism.
MANET maps evolutionary relationships directly and globally onto biological networks, and can generate and test hypotheses related to evolution of metabolism. We anticipate its use in the study of other networks, such as signaling and other protein-protein interaction networks.
Cellular metabolism represents a collection of enzymatic reactions and transport processes that convert metabolites into molecules capable of supporting cellular life. It is the best-studied biological network, with highly branched pathways describing the enzymatic processing of metabolites. Though underappreciated, it also represents one of the greatest achievements of science, resulting from almost two centuries of biochemical research.
There is considerable interest in the processes underlying the evolution of cellular metabolism. The existence of a core ensemble of metabolic reactions common to most organisms suggests that the global metabolic structure has been the subject of strong evolutionary constraint. Similarly, network connectivity properties suggest modular components typical of evolved systems [1–3] and emergence of hub metabolites involved in many reactions by enzyme specialization . How metabolic networks function and change as organisms increased in complexity remains an important question, making metabolism an interesting model for the evolution of biomolecular networks.
Metabolism is largely driven by enzymatic specificities. Consequently, the origin and evolution of metabolic networks can be explored advantageously by focusing on protein molecules. However, metabolism is very ancient and parts of the metabolic network probably evolved prior to the origin of cellular life from reactions that could have proceeded without catalysis or with inorganic catalysts . This view is supported to some extent by in vitro experiments that try to simulate pre-biotic chemistry. It is likely that polypeptides became metabolic catalysts through takeover of pre-biotic reactions . The earliest enzymes were probably weakly catalytic and multifunctional with broad specificities. Gradually, more numerous, effective, and specific enzymes evolved from the multifunctional enzymes through gene duplication, mutation and divergence. The only condition necessary for such a scenario appears to be selection for faster growth .
As enzymatic pathways became more complicated, new enzymatic functions and metabolic pathways could have been generated by recruitment of individual enzymes from the same or different pathways, or by enzymatic recruitments en masse from entire pathways. In this regard, several possible scenarios for the evolution of enzymes in metabolic pathways have been proposed . One popular scenario is the "backwards" (or retrograde) evolution hypothesis in which pathways evolve driven by successful production of their end products . Here, biosynthetic pathways undergo retro-evolution, with recruitment of enzymes (from within or outside the pathway) to host sites sequentially more remote from the end product of the pathway. By a symmetrical argument, catabolic pathways could have evolved sequentially from the metabolite being degraded . An alternative scenario is one in which new pathways evolve by "enzyme recruitment" from diverse donor sites throughout metabolism . This hypothesis assumes there is already an active enzymatic core with multifunctional and/or specialized enzymes from which new enzyme recruits are drawn for metabolic innovation. The result is an evolutionary "patchwork" of homologous enzymes that are present in different pathways .
There is considerable evidence supporting the patchwork recruitment scenario . For example, enzymes with α/β barrel fold structure that catalyze similar reactions occur across metabolic pathways . These patterns of structural homology appeared to be pervasive when structural assignments and sequence comparisons were used to analyze the small-molecule metabolism in Escherichia coli [13, 14]. Recruitment occurred with little regularity in these instances. However, proximity of donor and host sites appear to influence the probability of recruitment, with diversification to new host sites occurring mainly from nearby enzymes and varying with metabolite usage and enzyme class . It is noteworthy that sequence comparisons revealed homologous enzyme pairs occurring close to each other in the metabolic network more often than expected by chance . However, these homologous enzyme pairs had similar functions that could be best explained by patchwork recruitment. These observations suggest the retrograde evolution model played a small part in the process of metabolic enzyme evolution. None of these studies however used a phylogenetic approach to establish evolutionary patterns. Because common ancestry is the organizing principle underlying biology, we generate here a database for the evolutionary study of metabolic networks that integrates fragmentary knowledge about molecules and their interaction with phylogenomic information. This database characterizes patterns of evolution in cellular metabolism that will help to resolve the contribution of various plausible evolutionary scenarios.
Proteins consist of domains, compact sections of the protein molecule that have distinct structure, function and evolutionary history  and are used as a basis for structural classification [18, 19]. Protein domains represent a finite number of folding architectures, the so-called protein folds . These folds are highly diverse and are believed to originate from a common ancestor [20–22]. Crystallographic information gathered by structural genomics has enhanced our knowledge of the universe of fold architecture. This effort has been complemented by matching structures defined by a library of folds to genome sequences. For example, "occurrence analysis" methods compare how often a particular fold or fold group occurs in various genomes . This provides insights on the evolution of genomes because proteins with similar sequences have analogous structures, and structures are highly conserved in nature [23, 24]. Using this approach, whole-genome trees were reconstructed based on the occurrence of fold architectures and gene orthologs in genomes. These trees resemble those reconstructed from the sequence of the small subunit of ribosomal RNA (rRNA) and showed the tripartite nature of our organismal world [24–29].
We recently used an approach based on a census of folds to study protein diversification and reconstruct universal phylogenomic trees describing the evolution of protein fold architecture [27, 30, 31]. Our approach is based on two fundamental premises: (1) that protein structure is far more conserved than sequence and consequently carries considerable phylogenetic signal, and (2) that protein folds that are successful and popular in nature are generally more ancestral. Trees reconstructed from global fold-usage statistics showed there were clear evolutionary patterns in the appearance of protein folds. For example, the folds in the α/β protein class appeared at the base of the tree, and were followed by those in the α+β, all-α, all-β, small, and multidomain classes respectively . A similar tendency was recently observed when reconstructing parsimonious scenarios describing occurrence patterns of folds in organismal phylogenies . We also traced the number of enzymatic functions associated with folds in the tree of protein architecture and found that older folds were associated with an increased number of enzymatic functions. These phylogenomic studies suggest enzymatic multifunctionality was replaced by specialized function during evolution. Interestingly, a direct association between protein classes and function was previously revealed in which for example α/β folds were disproportionately associated with enzymes, especially transferases and hydrolases, while all-α and small proteins were associated with non-enzymes . Within α/β folds, five folds were the most functionally versatile (TIM-barrel, Rossmann, ferredoxin, α/β hydrolase, and P-loop NTP hydrolase). These folds were all placed at the base of our phylogenomic tree.
Construction and content
A set of phylogenetic features (characters) was used to describe the frequencies with which individual protein folds occur in an individual genome. This frequency was termed genomic abundance (G). At global levels, it describes how popular is a fold architecture in nature. In order to calculate G, we assigned structural domains to proteins. We did this using the HMM-based SUPERFAMILY database . SUPERFAMILY assigns protein structures to amino acid sequences at the fold superfamily level, i.e. at a hierarchical level capable of pooling proteins for which there is structural and sequence evidence of a common evolutionary ancestor . The HMM searching protocol uses a probability cut-off E of 0.02. Differences in topologies of trees reconstructed with more stringent cut-off values were found negligible , so we did not explore the role of this parameter. We also used the hierarchical scheme of SCOP to assign superfamilies to folds. Version 1.67 classifies 65,122 domains present in 24,037 PDB entries into 1,447 superfamilies and 887 protein folds. We analyzed the genome sequence of 36, 19 and 117 genomes from Eucarya, Archaea, and Bacteria, respectively, and reconstructed phylogenetic trees describing the evolution of 784 protein folds. Phylogenomic characters were coded as previously described . G was normalized using gap-recoding techniques to compensate for differences in genome size and proteome representation, and was then subjected to logarithmic transformation to account for unequal variances. The data was range standardized to a 0–20 scale compatible with most phylogenetic analysis programs, treated as linearly ordered multistate characters using an alphanumeric format with numbers 0–9 and letters A-K, and encoded in the NEXUS format. The ANCSTATES command was used to polarize characters assuming that the number of protein representatives in a genome exhibiting a particular fold increases in the course of evolution. Character argumentation is supported by model and assumptions described and discussed previously [27, 30, 31]. Phylogenies were reconstructed using maximum parsimony (MP) as the optimality criterion in PAUP* , and phylogenetic reliability was evaluated by the bootstrap method . Optimal most-parsimonious trees were obtained from heuristic maximum parsimony searches with tree-bisection-reconnection (TBR) branch swapping and 100 replicates of random addition sequence after exclusion of uninformative phylogenomic characters. To decrease search times during branch swapping of suboptimal trees, not more than one tree was saved in each replicate. The structure of phylogenetic signal in the data was tested by the skewness (g1) of the length distribution of > 104 random trees and permutation tail probability (PTP) tests of cladistic covariation using >103 replicates. Homoplasy and synapomorphy were measured with ensemble consistency (CI) and retention (RI) indices.
Because reconstructed trees were intrinsically rooted, we established the relative age of individual protein folds by measuring a distance in nodes from the hypothetical ancestral fold on a relative 0–1 scale. Concretely, we counted the number of nodes in every lineage from the root to the terminals of the tree and divided this number by the maximum number of nodes in a lineage. This node distance (nd) depicts the number of cladogenic events along a lineage and was used as an indicator of the ancestry of each metabolic enzyme for which a protein structure is known or could be inferred. A Perl script was written to extract ancestry values from phylogenetic trees.
Note that the central assumption that folds are more prevalent and widely shared the more ancient is their origin follows a parsimony rationale supported by patterns of distribution and sharing of protein folds across life, statistical analyses, and phylogenetic considerations . Our model is global and applies to the world of genomes, sampled by phylogenetic characters describing organisms in the three domains of life. Consequently, the model should be relatively insensitive to "genome lineage" specific factors such as selection pressures for genome expansion or reduction, life styles of organisms considered, and horizontal gene transfer events. However, our model is minimalist in that it does not account for differences in evolutionary rates across lineages and changes in the size of the protein world expected to have occurred during evolution. In particular, evolution of individual folds may be influenced by factors affecting rates, induced for example by instrinsic fold properties (e.g. flexibility). Consequently, fold ancestries should be regarded as lower bounds for times of origin.
We retrieved flat files directly from KEGG  and SCOP . The retrieval date for data presented in this paper was December 2004. We used Perl scripts to parse data from the files, store the parsed data into the entities of MANET, and manipulate the collected information in the database system. For example, we parsed data fields such as enzyme entry, pathway, and structural PDB entry from the database file "enzyme" obtained from KEGG, and inserted them into the entity called "MetabolicNetwork-A". Similarly, we parsed PDB entries (structures) and SCOP protein classifiers from database files such as "dir.cla.scop.txt_1.67" obtained from SCOP into "PDBclassification". Ancestry values were recorded into the entity "Ancestry".
Join operations and coloring of enzymes
Join operations necessary to combine information from entities in MANET and coloring of enzymes are described in the supplementary data [see Additional file 1]. The join operation linked 687 enzymes to protein folds. This represents about 35% of total enzymes associated with pathway information in the KEGG database. Coloring depicts graphically the relative age of metabolic enzymes, when these are associated to folds. A 'code generator' was used to output PHP and Visual Basic script files that paint the ancestries of enzymes on metabolic subnetwork diagrams.
Superfamily prediction using HMMs
In order to increase protein fold assignments to enzymes in metabolic pathways, we used a library of HMMs for remote homology detection in SUPERFAMILY. We used the model library, Perl wrapper scripts for sequence alignment, and the Sequence Alignment and Modeling System (SAM), and ran the SUPERFAMILY software package locally in a 15-node dual-processor Xserve cluster using the genes catalog file obtained from KEGG. This file includes amino acid and nucleotide sequence information from complete or partial genome sequences. Details can be found in the supplementary data file.
HMM-based structural prediction increased assignments of protein fold to enzymes in metabolic pathways. However, biases in fold superfamily prediction could affect evolutionary tracings in networks. To test the effect of prediction bias, we selected amino acid sequences associated with enzymes that had structural PDB entries. Statistical tests were performed using the SAS software package . The frequency distributions of ancestries derived directly from structural models (group A) and from HMM prediction (group B) were compared by using the Wilcoxon rank sum test. This test is one of several statistical tests that analyze two groups. It requires equal variances and independence of samples taken at random, but unlike t-based tests, it does not require for groups to have normal distributions . The null hypothesis of no difference between groups was tested using the NPAR1WAY procedure in SAS with grouped data.
We also performed a global statistical analysis of distribution of ancestries in mesonetworks. All average ancestry values of enzymes in the 11 mesonetworks described in KEGG were analyzed by ANOVA and by multiple pair-wise comparisons with the SAS package. There is a high probability of declaring at least one pair of means significantly different when running multiple comparisons unless the per-comparison error rate α is small among sample means . We therefore considered Type I error rates related to multiple comparisons, choosing the Least Squares Means with adjustment for multiple comparisons as the general linear model (GLM) procedure and the Tukey-Kramer method as the post-hoc test. The Tukey-Kramer multiple-comparison procedure controls error rates by testing every pair of means. The mesonetwork distribution data was subjected to logarithmic [log (y+1)] transformation. This reduced variances by about 1% and approximated data to normality.
A substantial body of literature has shown that linking protein structure to proteins in biomolecular networks can be advantageous. Some of these studies involved comprehensive structure-enzyme mapping exercises and explored mechanistic aspects related to enzymatic function. For example, global analysis of small-molecule metabolic pathways in Escherichia coli has shown extended distribution of structural homologues across metabolism [13, 14, 44, 45], sometimes confined to specific subnetworks . This and other evidence suggests the presence of widespread enzymatic recruitment and other evolutionary processes. Linking structure to metabolic function has also shown that in metabolism, catalytic mechanisms and co-factor binding properties are conserved while substrate specificity is variable [14, 40]. It appears it is easier to evolve binding sites than catalytic mechanisms. A recent study also shows for example that only a few fold superfamilies exhibit great substrate diversity, while most do not . Knowledge of structure can therefore help generate hypotheses about possible substrates associated with a protein.
However, proteins with similar structures may or may not be evolutionary related, and other approaches that introduce phylogenetic views are therefore needed. For example, Copley and Bork  used sequence, structure, and function to derive a phylogeny describing the evolution of members of 12 superfamilies with α/β barrel fold structure involved in metabolism. Establishing homologies at these levels provided indications that these fold superfamilies shared a common origin.
Our metabolic MANET database adds an evolutionary component to global assignment of protein structure to enzymes in metabolic subnetworks. Phylogenomic trees that describe the evolutionary relationship between protein fold architectures were used to define ancestries of individual folds, and these were traced onto metabolic subnetwork diagrams representing modern metabolism. Because our phylogenies were reconstructed from fold occurrence in hundreds of organisms with fully sequenced genomes spanning all three domains of life, they represent global phylogenetic statements about the protein world . These statements are therefore appropriately mapped to enzymatic structure in global cross-organismal representations of metabolism, such as those embedded in KEGG.
Tracing evolution of protein architecture in metabolic networks is useful. The exercise can uncover evolutionary patterns of architectural diversification within individual pathways of a subnetwork or between subnetworks and mesonetworks. Because metabolism is highly conserved, with about half of enzymes present in at least one species from the three domains of life , it is therefore possible to identify evolutionary patterns unique to the metabolic core that is universally conserved. MANET can also help identify processes driving the evolution of modern metabolism at local and global levels, including enzyme recruitment. It also allows query of SCOP, PDB, Enzyme Commission (EC) numbers, and subnetwork information that is useful for discovery of links between enzymatic activities and structures.
The web interface of MANET implements server-side scripts with a system of database management that provides visualization, query, and statistical analysis dynamically. All components of MANET are Perl-based and easily updated. Enzymes associated with protein folds are painted directly on subnetwork diagrams. Each enzyme is also hyper-linked to the KEGG database so that the user can retrieve additional information by clicking in the diagram. MANET also provides numerous functionalities, which enable searching for specific protein folds with defined ancestry values, displaying the distribution of enzymes that are painted, and exploring folds in individual subnetworks. Finally, the frequency distribution of ancestry values for each subnetwork can also be visualized.
Results and discussion
Reconstructed trees were based on a genomic census of protein architecture. Consequently, they depend on the accuracy of genomic databases, a balanced genomic sampling of the living world, efficient and accurate assignment of structures to proteins, a structural classification scheme that depicts evolutionary patterns, and methods of phylogenetic tree and character state reconstruction. The influence of these factors has been discussed previously [27, 30, 31]. While there is no possible gold standard that can be used to confirm the validy of phylogenomic statements, the genome census data we use to generate the tree of fold architectures was also used to generate trees of proteomes, and these trees group organisms in the three domains for the most part according to established organismal classification [Wang and Caetano-Anollés, ms. in preparation]. This observation supports the validity of phylogenetic signal embedded in the data.
Our study also rests on the accuracy of SCOP, a robust protein classification scheme [18, 49], and on the monophyletic nature of protein folds and superfamilies. Consequently our inferences should be regarded as rough first approximations. While we do not expect major changes in the operational definition of a protein fold, many folds could be better described by "continuous" rather than "discrete" distributions in structure space . Furthermore, we trust SCOP hierarchies reflect true evolutionary groupings. In SCOP, proteins in families express clear evolutionary relationships. They generally exhibit >30% pairwise residue identities or have functions and/or structures that provide definite evidence of common descent. Similarly, fold superfamilies contain proteins with structural and functional features that are highly suggestive of a common evolutionary origin. However, highly popular folds encompass collections of fold superfamilies that share the same arrangement and topology of secondary structures but may not have a common evolutionary origin. Consequently, the monophyletic nature of protein folds needs to be examined case by case, as has been done for the (βα)8 barrels [12, 51].
We also tested the accuracy of the HMM prediction. We selected PDB entries corresponding to enzymes in KEGG that participated in the join operation and were classified structurally, assigned PDB sequence records downloaded from ASTRAL  to the PDB entries, and analyzed the PDB sequences using the HMM package at E value = 0.02. Superfamily IDs and structural classifications corresponding to the PDB sequences were retrieved. Out of 21,173 PDB entries corresponding to enzymes identified by the join operation, 20,941 PDB entries mapped to ASTRAL. Sequences corresponding to these PDB entries were analyzed further. The HMM-based method rejected 212 sequences, leaving a total of 20,729 PDB entries with an assigned fold structure. Out of these, only 67 PDB entries differed in the expected fold assignment. At the fold superfamily level of classification, 106 PDB entries differed in the expected superfamily assignment. These results indicate that the HMM-based superfamily prediction can be performed at 98% accuracy levels. The details of this analysis can be found in our website.
The assignment of numerical ancestry values to enzymes in cellular metabolism uncovers evolutionary patterns of architectural diversification within the metabolic network. A quick examination of ancestry distributions depicted in each subnetwork and mesonetwork diagram of the MANET database reveals that enzymes of old origin generally coexist with those of recent origin (see example subnetwork; Fig. 4). A more detailed analysis of individual subnetwork paintings reveals the absence of clear patterns in individual pathways. Enzymes of old origin were generally followed haphazardly by enzymes of recent origin, and vice versa, with no apparent pattern along pathways. The patchy appearance of ancestries in subnetworks belonging to all metabolic mesonetworks supports strongly the enzyme recruitment (patchwork) evolutionary scenario as the major evolutionary force responsible of present day metabolism. Metabolic MANET makes visually evident enzymatic recruitment patterns that have been observed previously (e.g. [13, 14]), placing them into a relative evolutionary time frame. This offers the possibility of reconstructing temporal timelines of recruitment episodes in subnetworks and mesonetworks. Other evolutionary alternatives (backward evolution, forward evolution, de novo invention, pathway duplication, etc.) are not readily visible in our evolutionary tracing exercise. A detailed analysis of each subnetwork will be required to reveal the incidence of these possible evolutionary mechanisms. Pathway 'take-over' mechanisms in which new enzymes replace either pre-biotic chemistries or old enzymes, and 'co-option' mechanisms in which old enzymes gain novel functions, are also possible. In this regard, we are currently evaluating possible take-over episodes in metabolic subnetworks that may result from en masse enzymatic recruitment processes occurring in subnetwork pathways. We envision that uncovering take-over patterns in MANET at global levels will require extensive information about possible pre-biotic chemistries and novel phylogenetic approaches.
The evolutionary patterns revealed by MANET have other interesting implications. If we assume that pre-biotic chemistries remain imprinted in modern metabolism as relics of the pre-biotic world, patterns of enzymatic ancestries may reveal fundamental steps in prebiotic evolution. These evolutionary patterns may still manifest in the subnetworks despite obscuring events such as take-overs. Morowitz [53, 54] proposed that metabolism evolved through the sequential addition of shells to an "energy amphiphile" core (shell A), which consisted of the Krebs cycle, glycolysis, and fatty acid biosynthesis. The amination of 2-ketoglutarate was the gateway to shell B, the synthesis of most amino acids. In shell C sulfur was incorporated into cysteine and methionine. The gateways to shell D, ring closure and synthesis of nitrogen and dinitrogen heterocycles, gave access to purines, pyrimidines, and many cofactors, including B12. This scenario suggests that compounds in shell D evolved after enzymes (derived from shell B and C) and were not a part of prebiotic chemistry. The energy amphiphile core is consistent with Wächtershäuser's proposal that life evolve on pyrite (see ). According to this theory of an iron-sulfur world, a reductive citric acid cycle that used thio-organic homologues evolved early and was later coopted for oxidation. The reductive citric acid cycle, an autocatalytic network, expanded by branch reactions into higher homologous cycles. This archaic network included pathways for the synthesis and degradation of phosphorylated sugars, some amino acids (glutamate, aspartate, alanine, lysine), fatty acids and isoprenoids, coenzymes (including tetrapyrroles), and purines.
When ancestry patterns embedded in the subnetworks of MANET were analyzed, sequential evolution of metabolic "shells" was not obvious. However, pervasive enzyme recruitment could have masked the original pre-biotic evolutionary patterns. In fact, we performed a global statistical analysis of the distribution of ancestries of enzymes in metabolism, testing if global evolutionary patterns in metabolism matched possible "shell" scenarios (Fig. 6B). We calculated mean ancestry levels from frequency distribution patterns of ancestry data for mesonetworks, assuming these values were indicative of an average age of the enzymes examined. The statistics of distribution of ancestries in mesonetworks showed that distributions differed significantly in mean ancestry levels (p < 0.0001; ANOVA, F-test). Furthermore, the analysis revealed that amino acid mesonetworks were the oldest and lipid (including steroid) and glycan mesonetworks were relatively recent evolutionary additions (p < 0.05; Tukey-Kramer multiple comparison). The early evolutionary appearance of mesonetworks related to amino acid metabolism suggests that metabolic routes leading to the synthesis of polypeptides (shells B and C of Morowitz) 'internalized' early into the protein-based enzymatic machinery.
While mesonetworks may pool subnetworks of different average ancestry complicating interpretation, our results are nevertheless consistent with the shell hypothesis of Morowitz . In this regard, the early evolution of amino acid metabolic mesonetworks raises an interesting question. Why were the energy amphiphile core pre-biotic functions not the first to be replaced by enzymatic counterparts? These pre-biotic functions were the oldest and probably the most stable. One explanation is that replacement of non-enzymatic amino acid metabolic pathways follows the need to secure amino acid synthesis for protein-based enzymatic activities. It is possible that pre-biotic entities could have competed with each other for environmental resources during this early stage of metabolic evolution. Within this context, the opening of the gateway to amino acid synthesis proposed by Morowitz could have offered the possibility of creating enzymes that would perform pre-biotic functions more effectively.
We constructed a database that links biomolecular networks, protein structure, and phylogenomics. Metabolic MANET traces the evolution of protein structure directly onto metabolic networks defined by KEGG, enabling the study of evolutionary patterns in metabolism at global and local levels. MANET can be a valuable resource and constitutes a discovery tool. Individual pathways within subnetworks and mesonetworks can be examined and evolutionary patterns can be detected by visual inspection or statistical analysis of ancestry distributions. The database has many possible applications. For example, it can be used to search for patterns of fold superfamily sharing between subnetworks, with the aim of displaying the coordinated evolution of the subnetworks. Evolutionary information deposited in metabolic MANET will be enhanced by the exponential increase in the number of genomes that have been sequenced, the number of fold architectures uncovered, and the number of metabolic enzymes with gene assignments. The principles used in the construction of metabolic MANET are general and can be extended to other biomolecular networks. In the near future, we plan to analyze other networks of importance, such as cell signaling and other protein interaction networks.
Availability and requirements
The MANET database can be accessed via the Internet at http://manet.uiuc.edu. Data materials are formatted as EXCEL or column-delimited flat files for parsing with programming tools such as Perl and are openly available at our web site. The use of information obtained from the KEGG and SCOP databases is restricted by licensing conditions specified elsewhere. Contact information: THE MANET PROJECT, Atelier of Plant Bioinformatics, E-mail: email@example.com.
We thank S. Hwang for initial Perl scripts related to phylogenetic tree reconstruction and ancestry value generators. Research was supported in part with funds from UIUC and grants from the Office of Naval Research (TRECC A6538-A76) and NSF (MCB-0343126) to GCA. Any opinions, findings, and conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of the funding agencies.
- Rives AW, Galitski T: Modular organization of cellular networks. Proc Natl Acad Sci USA 2003, 100: 1128–1133. 10.1073/pnas.0237338100PubMed CentralView ArticlePubMedGoogle Scholar
- Ravasz E, Somera AL, Mongru DA, Oltvai ZN, Barabási A-L: Hierarchical organization of modularity in metabolic networks. Science 2002, 297: 1551–1555. 10.1126/science.1073374View ArticlePubMedGoogle Scholar
- Gimera R, Amaral LAN: Functional cartography of complex metabolic networks. Nature 2005, 433: 895–900. 10.1038/nature03288View ArticleGoogle Scholar
- Pfeiffer T, Soyer OS, Bonhoeffer S: The evolution of connectivity in metabolic networks. PLOS Biology 2005, 3: 1269–1275. 10.1371/journal.pbio.0030228View ArticleGoogle Scholar
- Maden BEH: No soup for starters? Autotrophy and the origins of metabolism. Trends Biochem Sci 1995, 20: 337–341. 10.1016/S0968-0004(00)89069-6View ArticlePubMedGoogle Scholar
- Kacser H, Beeby R: On the origin of enzyme species by means of natural selection. J Mol Evol 1984, 20: 38–51. 10.1007/BF02101984View ArticlePubMedGoogle Scholar
- Lazcano A, Miller SL: On the origin of metabolic pathways. J Mol Evol 1999, 49: 424–431. 10.1007/PL00006565View ArticlePubMedGoogle Scholar
- Schmidt S, Sunyaev S, Bork P, Dandekar T: Metabolites: a helping hand for pathway evolution? Trends Biochem Sci 2003, 28: 336–341. 10.1016/S0968-0004(03)00114-2View ArticlePubMedGoogle Scholar
- Horowitz NH: On the evolution of biochemical syntheses. Proc Natl Acad Sci USA 1945, 31: 153–157. 10.1073/pnas.31.6.153PubMed CentralView ArticlePubMedGoogle Scholar
- Cordón F: Tratado Evolucionista de Biologia. Madrid, Spain: Aguilar; 1990.Google Scholar
- Jensen RA: Enzyme recruitment in evolution of new function. Ann Rev Microbiol 1976, 30: 409–425. 10.1146/annurev.mi.30.100176.002205View ArticleGoogle Scholar
- Copley RR, Bork P: Homology among (βα) 8 barrels: implications for the evolution of metabolic pathways. J Mol Biol 2000, 303: 627–640. 10.1006/jmbi.2000.4152View ArticlePubMedGoogle Scholar
- Teichmann SA, Rison SCG, Thornton JM, Riley M, Gough J, Chothia C: Small-molecule metabolism: an enzyme mosaic. Trends Biotech 2001, 19: 482–486. 10.1016/S0167-7799(01)01813-3View ArticleGoogle Scholar
- Teichmann SA, Rison SC, Thornton JM, Riley M, Gough J, Cothia C: The evolution and structural anatomy of the small molecule metabolic pathways in Escherichia coli . J Mol Biol 2001, 311: 693–708. 10.1006/jmbi.2001.4912View ArticlePubMedGoogle Scholar
- Alves R, Chaleil RAG, Sternberg MJE: Evolution of enzymes in metabolism: A network perspective. J Mol Biol 2002, 320: 751–770. 10.1016/S0022-2836(02)00546-6View ArticlePubMedGoogle Scholar
- Light S, Kraulis P: Network analysis of metabolic enzyme evolution in Escherichia coli . BMC Bioinformatics 2004, 5: 15. 10.1186/1471-2105-5-15PubMed CentralView ArticlePubMedGoogle Scholar
- Riley M, Labedan B: Protein evolution viewed through Esherichia coli protein sequences: Introducing the notion of a structural segment of homology, the module. J Mol Biol 1997, 268: 857–868. 10.1006/jmbi.1997.1003View ArticlePubMedGoogle Scholar
- Murzin AG, Brenner SE, Hubbard T, Chothia C: SCOP: A structural classification of proteins database for the investigation of sequences and structures. J Mol Biol 1995, 247: 536–540. 10.1006/jmbi.1995.0159PubMedGoogle Scholar
- Orengo CA, Michie AD, Jones S, Jones DJ, Swindells MB, Thornton JM: CATH: a hierarchic classification of protein domain structures. Structure 1997, 5: 1093–1108. 10.1016/S0969-2126(97)00260-8View ArticlePubMedGoogle Scholar
- Chothia C, Gough J, Vogel C, Teichmann SA: Evolution of the protein repertoire. Science 2003, 300: 1701–1703. 10.1126/science.1085371View ArticlePubMedGoogle Scholar
- Chothia C: Proteins-1000 families for the molecular biologist. Nature 1992, 357: 543–544. 10.1038/357543a0View ArticlePubMedGoogle Scholar
- Orengo CA, Jones DT, Thornton JM: Protein superfamilies and domain superfolds. Nature 1994, 372: 631–634. 10.1038/372631a0View ArticlePubMedGoogle Scholar
- Gerstein M, Hegyi H: Comparing genomes in terms of protein structure: Surveys of a finite parts list. FEMS Microbiol Rev 1998, 22: 277–304. 10.1111/j.1574-6976.1998.tb00371.xView ArticlePubMedGoogle Scholar
- Gerstein M: Patterns of protein-fold usage in eight microbial genomes: A comprehensive structural census. Proteins 1998, 33: 518–534. 10.1002/(SICI)1097-0134(19981201)33:4<518::AID-PROT5>3.0.CO;2-JView ArticlePubMedGoogle Scholar
- Wolf YI, Brenner SE, Bash PA, Koonin EV: Distribution of protein folds in the three superkingdoms of life. Genome Res 1999, 9: 17–26.PubMedGoogle Scholar
- Lin J, Gerstein M: Whole-genome trees based on the folds and orthologs: Implications for comparing genomes on different levels. Genome Res 2000, 10: 808–818. 10.1101/gr.10.6.808PubMed CentralView ArticlePubMedGoogle Scholar
- Caetano-Anollés G, Caetano-Anollés D: An evolutionarily structured universe of protein architecture. Genome Res 2003, 13: 1563–1571. 10.1101/gr.1161903PubMed CentralView ArticlePubMedGoogle Scholar
- Deeds EJ, Hennessey H, Shakhnovich EI: Prokaryotic phylogenies inferred from protein structural domains. Genome Res 2005, 15: 393–402. 10.1101/gr.3033805PubMed CentralView ArticlePubMedGoogle Scholar
- Yang S, Doolittle RF, Bourne PE: Phylogeny determined by protein domain content. Proc Natl Acad Sci USA 2005, 102: 373–378. 10.1073/pnas.0408810102PubMed CentralView ArticlePubMedGoogle Scholar
- Caetano-Anollés G, Caetano-Anollés D: Universal sharing patterns in proteomes and evolution of protein fold architecture and life. J Mol Evol 2005, 60: 484–498. 10.1007/s00239-004-0221-6View ArticlePubMedGoogle Scholar
- Wang M, Boca SM, Kalelkar R, Mittenthal JE, Caetano-Anollés D: A phylogenomic reconstruction of the protein world based on a genomic census of protein fold architecture. Complexity 2006, in press.Google Scholar
- Winstanley HF, Abeln S, Deane CM: How old is your fold? Bioinformatics 2005, 21: i449-i458. 10.1093/bioinformatics/bti1008View ArticlePubMedGoogle Scholar
- Hegyi H, Gerstein M: The relationship between protein structure and function: a comprehensive survey with application to the yeast genome. J Mol Biol 1999, 288: 147–164. 10.1006/jmbi.1999.2661View ArticlePubMedGoogle Scholar
- Kanehisa M: Toward pathway engineering: A new database of genetics and molecular pathways. Institute for Chemical Research, Kyoto University; 1999.Google Scholar
- Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M: The KEGG resources for deciphering the genome. Nucleic Acids Res 2004, 32: D277-D280. 10.1093/nar/gkh063PubMed CentralView ArticlePubMedGoogle Scholar
- Gough J, Karplus K, Hughey R, Chothia C: Assignment of homology to genome sequences using a library of hidden Markov models that represent all proteins of known structure. J Mol Biol 2001, 313: 903–919. 10.1006/jmbi.2001.5080View ArticlePubMedGoogle Scholar
- Swofford DL: Phylogenetic analysis using parsimony and other programs (PAUP*), version 4.0. Sunderland, MA: Sinauer Associates; 1998.Google Scholar
- Felsenstein J: Confidence limits on phylogenies: an approach using the bootstrap. Evolution 1985, 39: 783–791. 10.2307/2408678View ArticleGoogle Scholar
- Kyoto Encyclopedia of Genes and Genomes[ftp://ftp.genome.ad.jp/pub/kegg]
- Structural Classification of Proteins[http://scop.mrc-lmb.cam.ac.uk/scop/parse/index.html]
- SAS: OnlineDoc®, version 8. Cary, NC: SAS Institute Inc; 1999.Google Scholar
- Ott RL, Longnecker M: Statistical methods and data analysis. 5th edition. Pacific Grove: Duxbury Press; 2001.Google Scholar
- Rison SCG, Teichmann SA, Thornton JM: Homology, pathway distance and chromosomal localization of the small molecule metabolism enzymes in Escherichia coli . J Mol Biol 2002, 318: 911–932. 10.1016/S0022-2836(02)00140-7View ArticlePubMedGoogle Scholar
- Saqui MAS, Sternberg JE: A structural census of metabolic networks for E. coli . J Mol Biol 2001, 313: 1195–1206. 10.1006/jmbi.2001.5107View ArticleGoogle Scholar
- Nobeli I, Spriggs RV, George RA, Thornton JML: A ligand-centric analysis of the diversity and evolution of protein-ligand relationships in E. coli . J Mol Biol 2005, 347: 415–436. 10.1016/j.jmb.2005.01.061View ArticlePubMedGoogle Scholar
- Peregrin-Alvarez JM, Tsoka S, Ouzounis CA: The phylogenetic extent of metabolic enzymes and pathways. Genome Res 2003, 13: 422–427. 10.1101/gr.246903PubMed CentralView ArticlePubMedGoogle Scholar
- Bingham J, Sudarsanam S: Visualizing large hierarchical clusters in hyperbolic space. Bioinformatics 2000, 16: 660–661. 10.1093/bioinformatics/16.7.660View ArticlePubMedGoogle Scholar
- Lo Conte L, Brenner SE, Hubbard TJP, Chothia C, Murzin A: SCOP database in 2002: Refinements accommodate structural genomics. Nucleic Acids Res 2002, 30: 264–267. 10.1093/nar/30.1.264PubMed CentralView ArticlePubMedGoogle Scholar
- Harrison A, Pearl F, Mott R, Thornton J, Orengo C: Quantifying the similarities within fold space. J Mol Biol 2002, 323: 909–929. 10.1016/S0022-2836(02)00992-0View ArticlePubMedGoogle Scholar
- Nagano N, Orengo CA, Thornton JM: One fold with many functions: the evolutionary relationships between TIM barrel families based on their sequences, structures and functions. J Mol Biol 2002, 321: 741–765. 10.1016/S0022-2836(02)00649-6View ArticlePubMedGoogle Scholar
- Chandonia JM, Hon G, Walker NS, Lo Conte L, Koehl P, Levitt M, Brenner SE: The ASTRAL compendium in 2004. Nucleic Acids Res 2004, 32: D189-D192. 10.1093/nar/gkh034PubMed CentralView ArticlePubMedGoogle Scholar
- Morowitz H: A theory of biochemical organization, metabolic pathways, and evolution. Complexity 1999, 4: 39–53. Publisher Full Text 10.1002/(SICI)1099-0526(199907/08)4:6<39::AID-CPLX8>3.0.CO;2-2View ArticleGoogle Scholar
- Morowitz HJ, Kostelnik JD, Yang J, Ody GD: The origin of intermediary metabolism. Proc Natl Acad Sci USA 2000, 97: 7704–7708. 10.1073/pnas.110153997PubMed 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.