High similarity of phylogenetic profiles of rate-limiting enzymes with inhibitory relation in Human, Mouse, Rat, budding Yeast and E. coli
© Zhao and Qu; licensee BioMed Central Ltd. 2011
Published: 30 November 2011
Skip to main content
© Zhao and Qu; licensee BioMed Central Ltd. 2011
Published: 30 November 2011
The phylogenetic profile is widely used to characterize functional linkage and conservation between proteins without amino acid sequence similarity. To survey the conservative regulatory properties of rate-limiting enzymes (RLEs) in metabolic inhibitory network across different species, we define the enzyme inhibiting pair as: where the first enzyme in a pair is the inhibitor provider and the second is the target of the inhibitor. Phylogenetic profiles of enzymes in the inhibiting pairs are further generated to measure the functional linkage of these enzymes during evolutionary history.
We find that the RLEs generate, on average, over half of all in vivo inhibitors in each surveyed model organism. And these inhibitors inhibit on average over 85% targets in metabolic inhibitory network and cover the majority of targets of cross-pathway inhibiting relations. Furthermore, we demonstrate that the phylogenetic profiles of the enzymes in inhibiting pairs in which at least one enzyme is rate-limiting often show higher similarities than those in common inhibiting enzyme pairs. In addition, RLEs, compared to common metabolic enzymes, often tend to produce ADP instead of AMP in conservative inhibitory networks.
Combined with the conservative roles of RLEs in their efficiency in sensing metabolic signals and transmitting regulatory signals to the rest of the metabolic system, the RLEs may be important molecules in balancing energy homeostasis via maintaining the ratio of ATP to ADP in living cells. Furthermore, our results indicate that similarities of phylogenetic profiles of enzymes in the inhibiting enzyme pairs are not only correlated with enzyme topological importance, but also related with roles of the enzymes in metabolic inhibitory network.
The phylogenetic profile of a given protein is a string that encodes the presence or absence of all the homologs of the protein in each organism . The similarity of phylogenetic profiles has been widely used to characterize functional linkage between proteins having no amino acid sequence similarity with each other [2–5]. The phylogenetic profiles of proteins with contextual information such as interactions, pathways, and cellular localizations have been widely studied. It is often supposed that proteins which tend to co-occur across species are more likely to co-evolve. The similarity of phylogenetic profiles is also correlated better with topological properties of metabolic enzymes such as degree and betweenness centrality, but independent of enzyme network importance .
The RLEs are a class of important enzymes playing roles of flux control in cell . According to modern flux control theory, RLEs often catalyze the slowest step in a metabolic pathway with a high flux control coefficient. Thus they are the ideal molecules to link metabolic pathways to regulatory networks. Based on their importance on the flux control and regulation of metabolic pathway, we built the first evidence-based database by systematic collecting RLEs in Human, Mouse, Rat, budding Yeast and E. coli. Our recent study reveals that RLEs in human liver often locate in the branch points and produce nearly half of in vivo metabolic enzyme inhibitors . However, the extent to which RLEs interact with enzyme inhibitors in liver or other tissue is less important in investigating the phylogenetic conservation of the regulatory mechanisms. The crucial question is whether these conservative RLEs involving inhibitory networks also maintain a higher functional linkage as expected.
In this study, enzyme inhibiting pairs and inhibitory network were constructed using enzyme inhibitor data from BRENDA database [10, 11]. We found that the RLEs provided over half of all in vivo inhibitors in five surveyed model organisms. These inhibitors inhibited over 85% targets in vivo in the metabolic inhibitory network. Additionally, we found that the phylogenetic profiles of RLEs in inhibiting pairs often showed higher similarities than those of common enzymes, which implied that the inhibitory relationship between two enzymes in an inhibiting pair was conservative across different species.
In total, 230 RLEs were collected from RLEdb . Based on these 230 RLEs, we extracted species specific RLE datasets using the enzyme organism distribution information from KEGG LIGAND database [12, 13]. In the five most common model organisms including Human, Mouse, Rat, budding Yeast and E. coli, there are 204, 201, 175, 139, 121 RLEs respectively (Additional file 1). According to the six pathway categories from KEGG (i.e., Carbohydrate, Lipid, Nucleotide, Amino acid, Cofactor and vitamin and Others molecular metabolism) as in , we investigated the functional distribution of all RLEs in each organism.
To describe how effectively enzyme inhibitor signals on the RLEs were spread to the metabolic system, all the fractions of RLEs as inhibitor targets to all targets from each pathway category in five model organisms were summarized in Figure 2b. The average fraction in human is over 0.45, which suggests that RLEs tend to be inhibited by enzyme inhibitors. As marked in Figure 2b, the inhibited targets are enriched in RLEs in five sixths pathway categories in four eukaryotic model organisms (hypergeometric tests, p value <0.05).
Besides both important roles as inhibitor providers and targets, further analysis (Figure 2c) shows that most targets of RLEs are inhibited by other RLEs for each pathway category in five model organisms. It implies that the overall targets of RLEs tend to be inhibited by themselves. All fractions of the inhibited RLEs by other RLEs to all inhibited RLEs by all metabolic enzymes in each pathway category in five model organisms were calculated. The average fraction in human is more than 0.86, which means that eighty-six percent of RLEs can be inhibited by other RLEs. Further hypergeometric tests reveal that all the inhibited RLEs are preferably inhibited by RLEs except one blue-marked as shown in Figure 2c (all p value < 0.05).
These global inhibitory relations conserved across several model organisms reflect potential regulatory modules during evolution of these organisms. Based on the number of their occurrence in Human, Mouse, Rat, budding Yeast and E. coli, 1070 conservative inhibiting enzyme pairs across at least three organisms were identified, of which 410 pairs contain RLEs as inhibitor provider. Our further analysis identified 317 and 234 conservative inhibiting pairs initiated by ADP and AMP respectively, among which 171 ADP and 32 AMP were products of RLEs (Additional file 2). Thus, comparing with the conservative inhibiting pairs initiated by ADP and AMP, RLEs are more likely to produce ADP instead of AMP in inhibitory network (hypergeometric tests, p value <0.05).
We conducted a systematic survey of regulatory roles of RLEs in metabolic inhibitory network in Human, Mouse, Rat, budding Yeast and E. coli. Our analysis emphasised their regulatory capacity and their possibility of being regulated in metabolic inhibitory regulation. As over half of in vivo inhibitors are produced by RLEs in five surveyed model organisms, RLEs as a whole can initiate inhibiting regulation and spread flux signals to other enzymes quickly. As inhibitor provider, RLEs cover over 85% in vivo inhibited targets in five model organisms. Viewed from inhibited targets, RLEs also possess their high possibility of being reached by inhibitors. Forty-five percent of inhibited targets are RLEs, which means an enzyme inhibitor has nearly a fifty percent chance to inhibit RLEs. Our recent study also reveals that the RLEs often locate near to branch points . Therefore, the inhibitory signals received by RLEs are easy to be amplified to different branches of the rest metabolic network in rate-limiting mode.
The cellular metabolic network can be thought as a molecular economic system of supply-and-demand relations between all the participating molecules . One of the most important cellular functions is to carry out the overall metabolic activities to adapt dynamically different environmental or in vivo metabolic conditions [17, 18]. As enzyme inhibitor is a short-term regulatory style, the cross pathway inhibiting capacity is useful as a feedback to balance metabolic signals between different pathway categories quickly. Our cross-pathway inhibiting analyses confirm the regulatory capacity of RLEs as precise inhibitor providers to maintain the balance of metabolic flux from different pathways, which may form an ideal mechanism to support a quick self-regulatory system.
To investigate the pathway topological evolution, earlier studies have been conducted on relations between phylogenetic profiles of metabolic enzymes and their topological importance including degree, betweenness centrality and closeness centrality . Our present study on phylogenetic profiles of enzymes in inhibiting pairs highlights the evolution of regulatory mechanisms in pathway. Our results suggest that the inhibiting pairs whose inhibitor providers are RLEs have higher similarity scores of phylogenetic profile than those of all common inhibiting pairs except in Rat. Furthermore the inhibiting pairs whose two members are RLEs have higher similarity scores of phylogenetic profile than all the inhibiting pairs where only inhibitor providers are RLEs except in Mouse. Despite the fact that inhibitor data in some organisms are not yet comprehensive, the higher similarities of phylogenetic profiles of RLEs imply their high functional linkages. Especially for conservative inhibiting pairs constituting two RLEs, their highest functional correlation implies that the RLEs possess modularity in metabolic inhibitory regulation network. In summary, our results indicate that the similarities of phylogenetic profiles of metabolic enzymes are not only correlated with topological importance, but also related with their roles in metabolic inhibitory network.
Our comparison between the conservative inhibiting pairs initiated by ADP and AMP indicates that RLEs are more likely to produce ADP instead of AMP in inhibitory network. Previous studies on regulatory role of small compounds in metabolic network revealed that chemical structures of compounds were related with their suitability for use in regulation . And their results show that ATP, ADP and AMP were the top three regulator compounds. Our results confirmed the important regulatory role of ADP in RLE involved inhibitory network. As molecular unit of currency, ATP is the most important energy source for metabolic reactions. And ADP is the most common and quickest product of hydrolysis of ATP in the cell with release of a large amount of free energy. The ratio of ATP to ADP in living cells is the most important metabolic signal in maintaining energy homeostasis. As the RLEs are often related with conversion between ATP to ADP, they are ideal molecules to sense and maintain energy homeostasis.
In conclusion, our systematic study reveals that over 60% in vivo enzyme inhibitors are products of RLEs in Human, Mouse, Rat and Yeast. And the conservative roles of RLEs in efficiency in sensing metabolic signals and transmitting regulatory signals to the rest of the metabolic network reveal that they are ideal regulatory molecules in balancing energy homeostasis via maintaining the ratio of ATP to ADP in living cells. Furthermore, our results indicate that the similarities of phylogenetic profiles of metabolic enzymes are not only correlated with their topological importance, but also related with their roles in enzyme inhibitory network.
230 potential RLEs from Human, Mouse, Rat, budding Yeast (Saccharomyces cerevisiae) and E. coli were compiled from RLE database (RLEdb) . Although some RLEs have species specific metabolic roles, most of RLEs are overlapped in the three mammalian organisms. To do comparative study for the global RLE roles in different organisms, we assumed that all the RLEs could play similar rate-limiting roles in the same pathway in different organisms. If there was a gene to encode a RLE in an organism, the RLE would be assigned to the organism. KEGG pathway annotations were added to each RLE according to KEGG-LIGAND database as described in .
In vivo enzyme inhibitor data were extracted from the BRENDA database 7.1 [10, 11]. And each inhibitor was assigned a KEGG compound identifier. Then the inhibiting enzyme pairs were constructed as described in . We did not transfer the inhibitor annotation between different organisms as the previous study [12, 13], so there were different amounts of inhibitors between different organisms. For example, the inhibitors in Mouse were far fewer than those in Human and Rat (Figure 1).
To analyze the functional linkage of two enzymes in an inhibiting pair, we assigned a profile to each enzyme involved in inhibiting network. To construct a reliable phylogenetic profile, 23 prokaryotes (Their genomes were completed before the end of 2000) and 59 eukaryotes (Their annotation data were from Refseq database in KEGG database) were selected. Data from these 82 well studied genomes were included to minimize the bias of in-complete enzyme annotation. Then for each given enzyme, a profile including a string with 82 entries were calculated. If there were one or more enzyme genes in a certain organism, the assigned value for the organism of the enzyme would be one, otherwise it would be assigned zero; finally the functional linkages of enzymes in inhibiting pairs were calculated by Jaccard coefficient.
To assess the similarity between two profiles A and B, we used Jaccard coefficient to score the overlapping values between A and B. Each value of A and B could either be 0 or 1.
The Jaccard similarity coefficient J was calculated as:
J = M 11 ÷ (M 01 + M 10 + M 11 )
Let J be Jaccard coefficient. The total number of each combination of values for both A and B were specified as follows:
M11 means the count of values when A and B both were 1.
M01 means the count of values when the value of A was 0 and the value of B was 1.
M10 means the count of values when the value of A was 1 and the value of B was 0.
Throughout the paper, we calculated p values based on hypergeometric test and unequal t test in R package 2.6.2 .
This work was supported by the China National High-tech 863 Programs (2006AA02A312, 2008BAI64B01). And we are deeply grateful to Professor Xiaocheng Gu for her suggestions in grammar, style and content which have greatly improve the manuscript.
This article has been published as part of BMC Genomics Volume 12 Supplement 3, 2011: Tenth International Conference on Bioinformatics – First ISCB Asia Joint Conference 2011 (InCoB/ISCB-Asia 2011): Computational Biology. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2164/12?issue=S3.
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.