Compartmentalization of the Edinburgh Human Metabolic Network
© Hao et al; licensee BioMed Central Ltd. 2010
Received: 30 October 2009
Accepted: 22 July 2010
Published: 22 July 2010
Direct in vivo investigation of human metabolism is complicated by the distinct metabolic functions of various sub-cellular organelles. Diverse micro-environments in different organelles may lead to distinct functions of the same protein and the use of different enzymes for the same metabolic reaction. To better understand the complexity in the human metabolism, a compartmentalized human metabolic network with integrated sub-cellular location information is required.
We extended the previously reconstructed Edinburgh Human Metabolic Network (EHMN) [Ma, et al. Molecular Systems Biology, 3:135, 2007] by integrating the sub-cellular location information for the reactions, adding transport reactions and refining the protein-reaction relationships based on the location information. Firstly, protein location information was obtained from Gene Ontology and complemented by a Swiss-Prot location keywords search. Then all the reactions in EHMN were assigned to a location based on the protein-reaction relationships to get a preliminary compartmentalized network. We investigated the localized sub-networks in each pathway to identify gaps and isolated reactions by connectivity analysis and refined the location information based on information from literature. As a result, location information for hundreds of reactions was revised and hundreds of incorrect protein-reaction relationships were corrected. Over 1400 transport reactions were added to link the location specific metabolic network. To validate the network, we have done pathway analysis to examine the capability of the network to synthesize or degrade certain key metabolites. Compared with a previously published human metabolic network (Human Recon 1), our network contains over 1000 more reactions assigned to clear cellular compartments.
By combining protein location information, network connectivity analysis and manual literature search, we have reconstructed a more complete compartmentalized human metabolic network. The whole network is available at http://www.ehmn.bioinformatics.ed.ac.uk and free for academic use.
Direct in vivo investigation of human metabolism is complicated by the distinct metabolic functions of different sub-cellular locations. For example, lysosomes are organelles containing digestive enzymes that break down polymeric macromolecules into their smaller building blocks. The lysosome membrane enables an acidic internal environment (pH less than 5 rather than around 7 in the cytosol) to maximize the enzyme activities . Due to localization of metabolic enzymes, many metabolic processes involve coordinated interactions between different organelles, and one metabolic step may be dependent upon the successful completion of the previous step. For example, the decomposition of very long chain fatty acids (VLCFAs) is a process shared by peroxisomes and mitochondria. Likewise, the final steps in the synthesis of plasmologens occur in the endoplasmic reticulum, but the process depends on precursors which are synthesized in peroxisomes . At the regulation level, the efficacy of many cellular processes is dependent on proper regulation of proteins trafficking to and from their site(s) of action. The endoplasmic reticulum (ER) and Golgi apparatus (GA) are known as the main organelles for protein targeting or protein sorting which transport proteins to the appropriate locations inside a cell or outside of it . Therefore it is quite usual that an enzyme synthesized in ER may be active only in another sub-cellular location. In addition, diverse micro-environments in different organelles may lead to distinct functions of the same enzyme. For example, acid ceramidase (EC 22.214.171.124; AC) is the lipid hydrolase responsible for the degradation of ceramide into sphingosine and free fatty acids within lysosomes. However, at higher pH in the cytosol AC can also synthesize ceramide from sphingosine and free fatty acids . Reduced lysosomal AC activity causes Farber disease, which is a member of a group of diseases called lysosomal storage diseases (LSDs) which results from defects in lysosomal enzyme function . Therefore, determining the location of enzymes and reactions is important for the investigation of the mechanism of a metabolic process and its related diseases. Currently there are two high quality literature based human metabolic networks available, the Edinburgh Human Metabolic Network (EHMN) reconstructed by our group  and the Human Recon 1 reconstructed by Palsson's group . Human Recon 1 contains eight sub-cellular locations while EHMN did not include location information . Localization of reactions in Human Recon 1 was determined from "protein localization data, sequence targeting signals, and indirect physiological evidence". If these data were unavailable, reactions were assigned to cytoplasm (cytosol in the data downloaded from BiGG) . As mentioned in our previous paper, EHMN is a more complete network with 1028 more reactions and 1202 more metabolites . Therefore, it is valuable to EHMN to include the information on sub-cellular location distribution of proteins and reactions. In this work, we chose 8 locations for enzyme proteins mainly based on information from GO (Gene Ontology) . The localization of reactions was initially determined by protein locations and then the gaps and isolated reactions in the preliminary compartmentalized network were examined and revised based on information from literature. By integrating the protein location information and the reaction (metabolite) location information, we also corrected hundreds of wrong protein-reaction relationships in EHMN. Furthermore, over 1400 transport reactions were added in order to obtain a connected compartmentalized network. The metabolic pathways for the synthesis and degradation of dozens of metabolites were investigated to validate the reconstructed network.
Results and discussion
1. Localization based on protein location
1.1 Protein location information
Protein distribution in different locations
Number of proteins
Proteins in EHMN
1.2 Sub-cellular localization of reactions
The distribution of reactions in different locations
Number of reactions
after revision of type c reaction
after gap filling
after IRLR revision
after literature based revision
We paid special attention to type c in the four types. Proteins catalyzing the same reaction in different locations are called "complementary proteins" in the study. We identified 16 such complementary protein groups (See additional file 1: Complementary groups and deleted protein-reaction relationships in revision of type c reaction). The two most common cases are: (a) the proteins are in the cytosol and mitochondria respectively (as for R00342 mentioned above); (b) the proteins are in peroxisomes and mitochondria respectively. Most of the protein groups in case b are involved in lipid metabolism and often catalyze a large class of reactions. As reported in literature, peroxisomes are the main places for the oxidation of VLCFAs . Therefore the protein location information may be used to identify the correct reactions catalyzed and thus improve the annotation of the protein-reaction relationships in EHMN. For example, RE3139 (CoA + 3-oxo-all-cis-6,9,12,15,18,21-tetracosahexaenoyl-CoA = acetyl-CoA + docosa-4,7,10,13,16,19-all-cis-hexaenoyl-CoA) is a step in beta-oxidation of Omega-3 fatty acid, which is a family of long chain and very long chain unsaturated fatty acids and their oxidation was known to take place in peroxisomes . RE3139 is associated with four proteins (ACAA1, ACAA2, HADHB and HADHA) in EHMN. However, only ACAA1 is in peroxisomes and the other three are mitochondria proteins. Therefore, the relationships between RE3139 and the three mitochondrial proteins in EHMN are wrong and should be corrected. Using this method, we have removed 43 incorrect protein-reaction relationships (See additional file 1: Complementary groups and deleted protein-reaction relationships in revision of type c reaction).
2. Location revision based on pathway connectivity analysis
According to the location of the reactions, we obtained the reaction location distribution for all the pathways in EHMN (See additional file 2: The reaction location distribution in the pathways in EHMN). We found that there are many poorly-studied pathways in which the number of reactions in any one specific location (excluding "uncertain") is less than one quarter of the total number of reactions in the pathway. This rather patchy feature of the localized network indicates that the location information from GO and SwissProt is still quite limited and more information is needed to improve the network for functional analysis. As the first step, we tried to improve the network connectivity by identifying the gaps or isolated reactions in the localized network by network structure analysis.
2.1 Gap filling
As depicted in Figure 3, a gap is a reaction which is in a different location with its two neighbour reactions in a pathway. Details about the methods for gap identification can be found in the Methods section. Altogether we found 510 gaps from 360 reactions in the 69 pathways ("isolated" pathway was excluded). Interestingly, a major part of the gap reactions (316 of 360) are in "uncertain". We directly assigned these reactions to the corresponding location to fill the gaps. This gap filling method is similar to previous approaches to identify the gaps in a reconstructed metabolic network for searching the missing genes/enzymes in pathways [13–15]. Now we applied this connectivity-based gap filling method to the localized network, thus complementing the patchy location data obtained solely from protein location information.
For the 49 gaps whose gap reactions are not in "uncertain", we used the information from the literature to decide whether we should assign the gap reaction to the corresponding location or not. Only those with positive evidences in the literature to support the existence of a gap reaction in the missing location were filled. Due to the limited research on reaction location, we only found 8 gaps with positive evidence in the literature. For example, R03617 (Galactosylceramide + H2O = beta-D-Galactose + N-Acylsphingosine) was identified as a gap in extracellular whereas Dolcetta et al  reported that the corresponding protein GALC actually catalyzed the reaction in the extracellular space. For the other 7 reactions, we found some new proteins that can catalyze the reaction in the missing locations. For example, one of the gaps in the di-unsaturated fatty acid beta-oxidation pathway in peroxisomes can be filled by adding reaction RE1516 (FAD + linoleoyl-CoA = FADH2 + trans,cis,cis-2,9,12-octadecatrienoyl-CoA). In EHMN, this reaction is catalyzed only by protein ACADL which is active in mitochondria. However, Bode and Couee  reported that two peroxisomal proteins: ACOX1 and ACOX3 have the same catalytic activity as ACADL. Therefore, RE1516 is catalyzed by ACOX1 and ACOX3 in peroxisomes. In this case, we updated EHMN by adding the protein-reaction relationships to the database. Altogether we added 14 protein-reaction relationships from literature and 469 gaps were filled, including 461 gaps from reactions with "uncertain" location and 8 gaps from reactions with other clear locations. The distribution of the reactions in each location after gap filling is shown in Table 2.
2.2 Isolated reactions
For some isolated reactions, we found positive evidence in the literature. For example, the L-Pipecolic acid oxidation reaction R02204 (L-Pipecolate + Oxygen = 2,3,4,5-Tetrahydropyridine-2-carboxylate + H2O2) should occur in human peroxisomes. Its deficiency in this location will cause a disease belonging to a group of disorders called peroxisome biogenesis disorders (PBDs) . Even though it is the only reaction in Lysine metabolism which occurs in peroxisomes, its location should not be changed, especially considering that oxygen is a reactant in this reaction. There are 19 such cases in our network. It is very likely that these reactions may be linked with other parts of the network through transport reactions.
We did not find any related literature for more than 80% of the IRLRs. In this case, we decided to keep an IRLR if the corresponding protein exists only in that location and it does not catalyze other reactions. As discussed above, if the protein exists in several locations, it may be inactive in the isolated location. Alternatively, if a protein catalyzes several reactions, it may imply that the protein may prefer to catalyze the other reaction(s) rather than the isolated one. For example, reaction R05989 (2 UDP-D-galactose + (GlcNAc)4 (LFuc)1 (Man)3 (Asn)1 = 2 UDP + (Gal)2 (GlcNAc)4 (LFuc)1 (Man)3 (Asn)1) in the N-Glycan biosynthesis pathway is an isolated reaction in the cytosol. Protein B4GALT1 is the protein related to this reaction in the cytosol but this protein also functions in extracellular and GA where it is connected with other reactions. In addition, B4GALT1 is responsible for the synthesis of complex-type N-linked oligosaccharides in many glycoproteins. Therefore, even if B4GALT1 has activity in the cytosol, its function may not be for reaction R05989 but for some other reactions. Accordingly, the location annotation of R05989 in the cytosol was deleted. However, further experimental studies would be needed to validate such predictions. In all the 177 isolated reaction-location relationships, 153 of them were revised, including 9 revised because of incorrect protein-reaction relationships and 144 without any related literature citations. The distribution of reactions in each location is shown in Table 2.
3. Literature-based revision
The objective of gap filling and IRLR revision is to make the compartmentalized network more amenable for further network analysis and functional analysis. However, even after such structure-based revision, there are still many pathways in which the number of reactions in any one specific location (excluding "uncertain") is less than one third of the total number of reactions in the pathway. To further improve the usability of the network, we manually examined the location distribution of all the pathways in EHMN, with a focus on pathway level, to check the literature for location information about a biological process including a series of reactions.
Although in total such information is still limited, we did find evidences in textbooks or in the literature for a number of pathways. In addition to the well known mitochondrial occurrence of the TCA cycle pathway, many studies are focused on the location of fatty acid beta-oxidation pathways. Three main criteria were used to decide the location of the beta-oxidation reaction. Firstly, it was reported that very long chain (greater than C-22) fatty acids underwent initial oxidation in peroxisomes which ceases at octanyl CoA and was then followed by mitochondrial oxidation . Thus, the length of fatty acid chain can be considered as a rule to distinguish between beta-oxidation in mitochondria and in peroxisomes. Secondly, oxidation reactions coupled with ATP should occur in mitochondria while those coupled with H2O2 should be in peroxisomes . The third criterion is the knowledge on the location of some specific metabolites. It has been reported that stearic acid and palmitic acid are oxidized firstly in peroxisomes to octanyl CoA and then oxidized in mitochondria. The unsaturated fatty acids (except very long chain acids) can be oxidized in both peroxisomes and mitochondria . According to the above criteria, the locations of 103 reactions in the beta-oxidation pathways were revised. Due to the revised reaction location, the incorrect protein-reaction relationships were also corrected. For example, the location of R03856 (Lauroyl-CoA + NADP+ = 2-trans-Dodecenoyl-CoA + NADPH + H+) was revised from both peroxisomes and mitochondria to peroxisomes only because it is a step in the degradation of palmitic acid in peroxisomes before ceasing at octanyl CoA. In EHMN, R03856 is linked with two proteins: MECR in mitochondria and PECR in peroxisomes. The MECR-R03856 relationship can therefore be deleted because the reaction only occurs in peroxisomes. In this way, a total of 129 protein-reaction relationships were deleted in the beta-oxidation pathways. We also revised the location information for 337 reactions in other pathways (e.g. Pentose phosphate pathway, Valine, leucine and isoleucine degradation pathway) based on evidences from textbooks or the literature and more than 120 incorrect protein-reaction relationships were deleted (See additional file 3: Result of literature-based revision). The reaction distribution in each location is shown in Table 2.
4. Transport reactions
It is important to include transport reactions between different compartments in a localized network in order to make the network connected. It is unreasonable to add a transport reaction involving the "uncertain" location. Therefore we assigned all the reactions with "uncertain" location to cytosol as that used in Recon 1 . As Recon 1 already included transport reactions, we first took it as a reference to add the transport reactions to our network. There are 1190 transport reactions in Recon 1 and 857 of them are for metabolites that are also in EHMN. Among the 857 reactions, only 378 reactions that are based on biochemical or genetic experiments (with confidence score 3 in Recon 1) were directly added into EHMN. For the 310 reactions based on physiological data or biochemical/genetic evidence from a nonhuman mammalian cell (with confidence score 2 in Recon 1), a transport reaction was added only if the corresponding metabolite exists in both compartments involved in the transport reaction. The reactions with confidence score less than 2 in Recon 1 were ignored because those reactions were added merely based on in sillico simulation. Finally 546 transport reactions were added from Recon 1 and 429 of them are related with 212 transport proteins.
At the second step, we searched human transport proteins in TransportDB  and all the proteins with gene ontology annotation GO: 0005215 (transporter activity) in Uniprot. There are over 1000 transport proteins obtained from the two resources. However, a large number of the proteins are ion channel related proteins or proteins annotated based on sequence similarity but without clear substrate information. Excluding the 212 transport proteins already obtained from Recon 1, only a small proportion of the transport proteins were investigated to add corresponding transport reactions based on information from the databases and literature. Altogether 109 new transport reactions from 83 transport proteins were added to EHMN. Some of them are actually also in Recon 1 but with a low confidence score and without assigned transport protein.
After adding 655 reactions based on Recon 1 and transporter information from databases and literature, there are still many gaps between different compartments in the network mainly because of our limited knowledge on the function of transport proteins and the existence of non-protein required transport reactions. To make the network connected and thus amenable for pathway analysis, we calculated the "dead ends" in each location for each pathway (excluding the isolated pathway) and tried to link the dead ends in different locations by adding new transport reactions. A dead end was defined as a metabolite that is only produced or consumed in a location. If a dead end metabolite is also appeared in other locations then transport reactions are added for it. By default, cytosol was used as a bridge to connect other locations and transport reaction to/from cytosol were added for the dead end metabolites. Altogether 961 transport reactions were identified for the dead end metabolites. Among them, 193 reactions were already added in the first two steps (from Recon 1 and transport proteins). Therefore 768 new transport reactions were needed to make the network connected. It should be noted that the reactions added solely based on dead end analysis are of low quality and the existence are subjected to further experimental investigation. Among these reactions, 73 are between cytosol and nucleus which actually can be carried out by the nuclear pore complexes formed by the 49 nucleoporin proteins in human .
As a summary 1423 transport reactions were added to the network and 611 of them are associated with transport proteins. The whole network containing 4804 reactions is available at http://www.ehmn.bioinformatics.ed.ac.uk and free for academic use.
5. Comparison of EHMN with other human metabolic networks
Sub-cellular location comparison of Human Recon 1, EHMN, and Reactome database
Human Recon 1a
Non-repeat metabolic reactionse
6. Function analysis of the compartmentalized EHMN
A compartmentalized human metabolic network has been reconstructed by adding the sub-cellular location information to enzymes and reactions in EHMN. Network structure analysis and manual curation of literature were used to improve the quality of the localized network which was mainly based on protein location information. Transport reactions were added to obtain a connected sub-cellular location specific metabolic network. The refined network was validated by in silico investigation of the metabolic functions for the production/degradation of certain important metabolites.
Localization of human proteins
The protein location information in this work is extracted mainly from Gene Ontology (GO). Files containing gene-GO association in human and the hierarchically organized GO terms ("OBO" file) can be easily downloaded from the GO website. This file is maintained by the GOA group at EBI which aims to provide high-quality GO annotations to proteins in the UniProtKB. GO database provides a controlled vocabulary of hierarchically organized terms to describe the location in the "cellular component" category. There are over one thousand GO terms in the "cellular component" category. We chose eight top locations to compartmentalize the EHMN based on the cellular structure. The selected locations and their hierarchical structure are shown in Figure 1. The organelles considered as individual compartments are nucleus, Endoplasmic Reticulum (ER), Golgi apparatus (GA), peroxisomes, lysosomes and mitochondria. Other organelles such as microsome, cytoskeleton and small vesicles were not considered because their functions are often not related with cellular metabolism. The term "plasma membrane" in GO was merged with "cytosol" in this work, so all the proteins which are annotated in "plasma membrane" in GO were classified to "cytosol" because they usually catalyze reactions occurring in the cytosol. From the gene-GO association file, many proteins are actually associated with a subordinate term such as mitochondrial membrane (GO: 0031966). In these cases, we traced the GO tree back to a higher-level term within the eight chosen locations (mitochondria in this case). The proteins with uncertain locations or other locations not included in the eight specific locations were classified to an "uncertain" location. In this way, all the human proteins coming from GO were classified into the chosen locations.
To obtain more protein location information and cross-validate information from different sources, we also used the Swiss-Prot keywords to extract all the human proteins related with the keywords. There are more than eighty locations listed as keywords in the "cellular component" category in Swiss-Prot keywords, in which six locations, including nucleus, ER, GA, peroxisomes, lysosomes and mitochondria, are same with the GO terms used in our study. We extracted all the human proteins related with these location keywords to get the protein-location relationships.
We developed a method to automatically identify gaps in the network. For each pathway (excluding isolated pathway), we extracted a subset of reactions which are in a specific location and then converted it into a graph. The group of reactions in a location which can be connected were classified into one component, so the reactions in a location may be divided into several separated components. It should be noted that only the main compounds in a reaction were used in linking reactions and the currency metabolites were ignored . Then we calculated the connected components in the reaction graph  and obtained a list of main metabolites in each component. In the next step, we checked every reaction which is not in this location to see if it contains main metabolites in two different components. If so, this reaction will be identified as a gap for that location. This process was repeated for each location (excluding "uncertain") and one reaction may be identified as gaps in different locations. A gap is a missing reaction-location relationship rather than a reaction.
isolated reaction-location relationship
Very long chain fatty acid.
This work was supported by the National Natural Science Foundation of China (NSFC-20806055, NSFC-20875068), the National Project of Key Fundamental Research (2007CB707802), the Development Project of Science and Technology of Tianjin (05YFGZGX04500), Programme of Introducing Talents of Discipline to Universities (B06006) and State Scholarship Fund of China Scholarship Council (2008625087).
- De Duve C, Wattiaux R: Functions of lysosomes. Annu Rev Physiol 1966, 28: 435–492. 10.1146/annurev.ph.28.030166.002251View ArticlePubMedGoogle Scholar
- Colin J, Masters DC: The peroxisome:a vital organelle. Cambridge University Press; 2005:101–105.Google Scholar
- Presley JF, Cole NB, Schroer TA, Hirschberg K, Zaal KJ, Lippincott-Schwartz J: ER-to-Golgi transport visualized in living cells. Nature 1997, 389(6646):81–85. 10.1038/38001View ArticlePubMedGoogle Scholar
- Park JH, Schuchman EH: Acid ceramidase and human disease. Biochim Biophys Acta 2006, 1758(12):2133–2138. 10.1016/j.bbamem.2006.08.019View ArticlePubMedGoogle Scholar
- Winchester B, Vellodi A, Young E: The molecular basis of lysosomal storage diseases and their treatment. Biochem Soc Trans 2000, 28(2):150–154.View ArticlePubMedGoogle Scholar
- Ma H, Sorokin A, Mazein A, Selkov A, Selkov E, Demin O, Goryanin I: The Edinburgh human metabolic network reconstruction and its functional analysis. Mol Syst Biol 2007, 3: 135. 10.1038/msb4100177View ArticlePubMedPubMed CentralGoogle Scholar
- Duarte NC, Becker SA, Jamshidi N, Thiele I, Mo ML, Vo TD, Srivas R, Palsson BO: Global reconstruction of the human metabolic network based on genomic and bibliomic data. Proceedings of the National Academy of Sciences of the United States of America 2007, 104(6):1777–1782. 10.1073/pnas.0610772104View ArticlePubMedPubMed CentralGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al.: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet 2000, 25(1):25–29. 10.1038/75556View ArticlePubMedPubMed CentralGoogle Scholar
- Gomez-Fabre PM, Aledo JC, Del Castillo-Olivares A, Alonso FJ, Nunez De Castro I, Campos JA, Marquez J: Molecular cloning, sequencing and expression studies of the human breast cancer cell glutaminase. Biochemical Journal 2000, 345(2):365–375. 10.1042/0264-6021:3450365View ArticlePubMedPubMed CentralGoogle Scholar
- Modi WS, Pollock DD, Mock BA, Banner C, Renauld JC, Van Snick J: Regional localization of the human glutaminase (GLS) and interleukin-9 (IL9) genes by in situ hybridization. Cytogenetics and cell genetics 1991, 57:(2–3):114–116.View ArticleGoogle Scholar
- Ide Y, Miyoshi E, Nakagawa T, Gu J, Tanemura M, Nishida T, Ito T, Yamamoto H, Kozutsumi Y, Taniguchi N: Aberrant expression of N-acetylglucosaminyltransferase-IVa and IVb (GnT-IVa and b) in pancreatic cancer. Biochemical and Biophysical Research Communications 2006, 341(2):478–482. 10.1016/j.bbrc.2005.12.208View ArticlePubMedGoogle Scholar
- Ferdinandusse S, Denis S, Dacremont G, Wanders RJ: Studies on the metabolic fate of n-3 polyunsaturated fatty acids. Journal of lipid research 2003, 44(10):1992–1997. 10.1194/jlr.M300223-JLR200View ArticlePubMedGoogle Scholar
- Claudel-Renard C, Chevalet C, Faraut T, Kahn D: Enzyme-specific profiles for genome annotation: PRIAM. Nucleic Acids Res 2003, 31(22):6633–6639. 10.1093/nar/gkg847View ArticlePubMedPubMed CentralGoogle Scholar
- Osterman A, Overbeek R: Missing genes in metabolic pathways: a comparative genomics approach. Curr Opin Chem Biol 2003, 7(2):238–251. 10.1016/S1367-5931(03)00027-9View ArticlePubMedGoogle Scholar
- Cordwell SJ: Microbial genomes and "missing" enzymes: redefining biochemical pathways. Arch Microbiol 1999, 172(5):269–279. 10.1007/s002030050780View ArticlePubMedGoogle Scholar
- Dolcetta D, Perani L, Givogri MI, Galbiati F, Amadio S, Del Carro U, Finocchiaro G, Fanzani A, Marchesini S, Naldini L, et al.: Design and optimization of lentiviral vectors for transfer of GALC expression in Twitcher brain. The journal of gene medicine 2006, 8(8):962–971. 10.1002/jgm.924View ArticlePubMedGoogle Scholar
- Bode K, Hooks MA, Couee II: Identification, separation, and characterization of acyl-coenzyme A dehydrogenases involved in mitochondrial beta-oxidation in higher plants. Plant physiology 1999, 119(4):1305–1314. 10.1104/pp.119.4.1305View ArticlePubMedPubMed CentralGoogle Scholar
- Steinberg SJ, Wang SJ, Kim DG, Mihalik SJ, Watkins PA: Human very-long-chain acyl-CoA synthetase: cloning, topography, and relevance to branched-chain fatty acid metabolism. Biochemical and biophysical research communications 1999, 257(2):615–621. 10.1006/bbrc.1999.0510View ArticlePubMedGoogle Scholar
- Wanders RJ, Heymans HS, Schutgens RB, Barth PG, van den Bosch H, Tager JM: Peroxisomal disorders in neurology. Journal of the neurological sciences 1988, 88(1–3):1–39. 10.1016/0022-510X(88)90203-1View ArticlePubMedGoogle Scholar
- Singh I, Pahan K, Dhaunsi GS, Lazo O, Ozand P: Phytanic acid alpha-oxidation. Differential subcellular localization in rat and human tissues and its inhibition by nycodenz. J Biol Chem 1993, 268(14):9972–9979.PubMedGoogle Scholar
- Verhoeven NM, Roe DS, Kok RM, Wanders RJ, Jakobs C, Roe CR: Phytanic acid and pristanic acid are oxidized by sequential peroxisomal and mitochondrial reactions in cultured fibroblasts. J Lipid Res 1998, 39(1):66–74.PubMedGoogle Scholar
- Mihalik SJ, Steinberg SJ, Pei Z, Park J, Kim DG, Heinzer AK, Dacremont G, Wanders RJ, Cuebas DA, Smith KD, et al.: Participation of two members of the very long-chain acyl-CoA synthetase family in bile acid synthesis and recycling. The Journal of biological chemistry 2002, 277(27):24771–24779. 10.1074/jbc.M203295200View ArticlePubMedGoogle Scholar
- Dodt G, Kim DG, Reimann SA, Reuber BE, McCabe K, Gould SJ, Mihalik SJ: L-Pipecolic acid oxidase, a human enzyme essential for the degradation of L-pipecolic acid, is most similar to the monomeric sarcosine oxidases. Biochemical Journal 2000, 345(3):487–494. 10.1042/0264-6021:3450487View ArticlePubMedPubMed CentralGoogle Scholar
- Wilcke M, Alexson SE: Characterization of acyl-CoA thioesterase activity in isolated rat liver peroxisomes. Partial purification and characterization of a long-chain acyl-CoA thioesterase. European journal of biochemistry/FEBS 1994, 222(3):803–811. 10.1111/j.1432-1033.1994.tb18927.xView ArticlePubMedGoogle Scholar
- Hutson SM, Fenstermacher D, Mahar C: Role of mitochondrial transamination in branched chain amino acid metabolism. J Biol Chem 1988, 263(8):3618–3625.PubMedGoogle Scholar
- Suryawan A, Hawes JW, Harris RA, Shimomura Y, Jenkins AE, Hutson SM: A molecular model of human branched-chain amino acid metabolism. Am J Clin Nutr 1998, 68(1):72–81.PubMedGoogle Scholar
- Hutson SM, Sweatt AJ, Lanoue KF: Branched-chain [corrected] amino acid metabolism: implications for establishing safe intakes. J Nutr 2005, 135(6 Suppl):1557S-1564S.PubMedGoogle 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. 10.1101/gr.2250904View ArticlePubMedPubMed CentralGoogle Scholar
- Ren Q, Chen K, Paulsen IT: TransportDB: a comprehensive database resource for cytoplasmic membrane transport systems and outer membrane channels. Nucleic Acids Res 2007, (35 Database):D274–279. 10.1093/nar/gkl925
- McMorrow I, Bastos R, Horton H, Burke B: Sequence analysis of a cDNA encoding a human nuclear pore complex protein, hnup153. Biochim Biophys Acta 1994, 1217(2):219–223.View ArticlePubMedGoogle Scholar
- Matthews L, Gopinath G, Gillespie M, Caudy M, Croft D, de Bono B, Garapati P, Hemish J, Hermjakob H, Jassal B, et al.: Reactome knowledgebase of human biological pathways and processes. Nucleic Acids Res 2009, (37 Database):D619–622. 10.1093/nar/gkn863
- Romero P, Wagg J, Green ML, Kaiser D, Krummenacker M, Karp PD: Computational prediction of human metabolic pathways from the complete human genome. Genome Biol 2005, 6(1):R2. 10.1186/gb-2004-6-1-r2View ArticlePubMedPubMed CentralGoogle Scholar
- Patton SM, Pinero DJ, Surguladze N, Beard J, Connor JR: Subcellular localization of iron regulatory proteins to Golgi and ER membranes. J Cell Sci 2005, 118(Pt 19):4365–4373. 10.1242/jcs.02570View ArticlePubMedGoogle Scholar
- Ponka P: Tissue-specific regulation of iron metabolism and heme synthesis: distinct control mechanisms in erythroid cells. Blood 1997, 89(1):1–25.PubMedGoogle Scholar
- Straka JG, Rank JM, Bloomer JR: Porphyria and porphyrin metabolism. Annu Rev Med 1990, 41: 457–469. 10.1146/annurev.me.41.020190.002325View ArticlePubMedGoogle Scholar
- Ma H, Zeng AP: Reconstruction of metabolic networks from genome data and analysis of their global structure for various organisms. Bioinformatics (Oxford, England) 2003, 19(2):270–277. 10.1093/bioinformatics/19.2.270View ArticleGoogle Scholar
- Ma HW, Zeng AP: The connectivity structure, giant strong component and centrality of metabolic networks. Bioinformatics (Oxford, England) 2003, 19(11):1423–1430. 10.1093/bioinformatics/btg177View ArticleGoogle 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.