Profiling the orphan enzymes
© Sorokina et al.; licensee BioMed Central Ltd. 2014
Received: 27 March 2014
Accepted: 29 May 2014
Published: 6 June 2014
The emergence of Next Generation Sequencing generates an incredible amount of sequence and great potential for new enzyme discovery. Despite this huge amount of data and the profusion of bioinformatic methods for function prediction, a large part of known enzyme activities is still lacking an associated protein sequence. These particular activities are called “orphan enzymes”. The present review proposes an update of previous surveys on orphan enzymes by mining the current content of public databases. While the percentage of orphan enzyme activities has decreased from 38% to 22% in ten years, there are still more than 1,000 orphans among the 5,000 entries of the Enzyme Commission (EC) classification. Taking into account all the reactions present in metabolic databases, this proportion dramatically increases to reach nearly 50% of orphans and many of them are not associated to a known pathway. We extended our survey to “local orphan enzymes” that are activities which have no representative sequence in a given clade, but have at least one in organisms belonging to other clades. We observe an important bias in Archaea and find that in general more than 30% of the EC activities have incomplete sequence information in at least one superkingdom. To estimate if candidate proteins for local orphans could be retrieved by homology search, we applied a simple strategy based on the PRIAM software and noticed that candidates may be proposed for an important fraction of local orphan enzymes. Finally, by studying relation between protein domains and catalyzed activities, it appears that newly discovered enzymes are mostly associated with already known enzyme domains. Thus, the exploration of the promiscuity and the multifunctional aspect of known enzyme families may solve part of the orphan enzyme issue. We conclude this review with a presentation of recent initiatives in finding proteins for orphan enzymes and in extending the enzyme world by the discovery of new activities.
This article was reviewed by Michael Galperin, Daniel Haft and Daniel Kahn.
New progress in sequencing technologies generates thousands of new sequences each day. With the large public sequence databases combined with efficient bioinformatic methods, it is possible to predict the function of some new proteins mainly by comparative genomics approaches. Nevertheless, millions of protein entries are not assigned reliable functions due to the lack of trustworthy annotations and the drawbacks of homology-based predictions . This shortcoming illustrates our limited knowledge of the functional diversity in the protein world and restricts the analyses of an organism starting from its genome. This is particularly the case for enzymatic activities that can be predicted by gene functional assignments and used as a starting point to reconstruct genome-scale metabolic models.
The first enzyme was discovered and isolated in 1833 by Anselme Payen . It was the first time a non-living compound was shown to have properties of an organic catalyst, a discovery which shook the scientific community. This enzyme was named “diastase” (now called α-amylase) and the suffix –‘ase’ will be henceforth used to refer to enzymes. Since then, the number of discovered enzymes has continually increased, thanks to the experimental work of chemists and biologists. In the beginning of enzymology, the naming of enzyme was not systematic. Many different enzymes were given similar names and, on the other hand, the same enzymes had several names. An Enzyme Commission, whose first meeting took place in 1961, was created to give rules and recommendations that could be implemented for the systematic naming of enzymes . Enzyme activities are nowadays classified with EC (Enzyme Commission) numbers, a nomenclature maintained by the IUBMB (International Union of Biochemistry and Molecular Biology) [4–6]. To be integrated into the EC classification, an activity must be observed and biochemically characterized without the necessity to identify the associated protein that catalyzes the reaction.
In 2005, sequence-lacking enzymatic activities were named “orphan enzymes” by Bernard Labedan and Olivier Lespinet in an open letter . They conducted a similar analysis to that of Peter Karp and showed that 42% of the EC numbers were orphan enzymes (1,625 EC numbers among 3,820). One of the main surprises of this study was the fact that 200 organisms had orphan enzymes, despite the availability of their complete genome. They also noticed that, in several cases, the protein catalyzing the enzymatic activity had been identified but not sequenced. The following year they published two exploratory articles on orphan enzymes [10, 11]. The proportion of orphan enzymes was updated, giving a slight decrease of 3% (39% of orphans, 1,525 EC entries among 3,877). They pointed out that a number of pathways (~100) had at least one orphan enzyme. They also made several remarks on the use of EC numbers. Moreover, they created a public database, called ORENZA, listing all orphan enzymes present in the EC nomenclature and allowing users to perform queries by tracking them between organisms and pathways .
In 2007, Lifeng Chen and Dennis Vitkup carried out a very detailed review on the historical accumulation of orphan enzyme activities and a wide range of statistical analyses on their distribution across different classifications . They found 1,360 orphans, representing 34% of the 4,003 valid EC entries. They investigated the number of biochemical characterizations per year of discovery and noticed that it decreased in the 1970s and 1990s. A study of the relation between orphan enzymes and their pathway neighbors was conducted: 39% of network neighbors for orphan activities were orphan themselves, compared with 29% for neighbors of non-orphan activities. They also noticed that a majority of orphan activities were found in the most studied organisms. Finally, they pinpointed a possible bias in the EC classification because many reactions in metabolic databases were not associated with any EC number. Considering this limitation, they estimated that up to 50% of all know biochemical reactions were orphan.
Here, we present an extended review on orphan enzyme activities by updating previously conducted surveys and performing new analyses. We first update the estimation of the number of orphan enzymes and interpret their decrease in the light of past and recent enzyme activity discoveries. As the EC classification does not totally cover all known activities, we briefly introduce two main metabolic databases and analyze their content to estimate orphans at the reaction level. Also, an analysis of their connectivity in metabolic network is made. The concept of orphan enzymes is then extended to local orphans (i.e. activities which have no representative sequence in a given clade, but have one in other organisms) and an analysis is made at the superkingdom level to estimate their number and to evaluate if candidate proteins for local orphans could be retrieved by sequence homology. Finally, we expose the notion of promiscuity and multifunctionality in the enzyme world and explore the relation between protein domains and catalyzed activities. In conclusion, we present some new initiatives and concepts of interest to reduce the number of orphan enzymes but, also, to extend the landscape of enzymes by finding new activities.
An updated view of orphan enzymes
In this study, we estimated the number of orphan enzymes by using EC numbers present in the IntEnz  and UniProt  databases (versions of February 2013). UniProt is a resource of proteins where enzymatic activities are described using the EC classification. Only valid and complete EC entries were considered without taking into account deleted or transferred entries. We also considered as valid entries the nearly 100 provisional EC numbers of IntEnz waiting to be confirmed by the IUBMB. It appears that 22.4% of the enzymatic activities are orphans; among the 5,096 EC numbers, 1,143 entries have no associated protein in UniProt. As noticed previously , the proportion of orphan enzymes is not uniformly distributed across the different classes of the EC nomenclature: in EC class 1 the fraction is 25%, 26% in class 2, 19% in class 3 and 4, 15% in class 5 and 13% in class 6 (Additional file 1: Figure S1.1 and Additional file 2: Table S2.1 for the complete list of orphan EC numbers).
Orphan enzymes in the metabolic world
It is important to distinguish the terms “enzyme” and “enzymatic activity”. The first designates a protein able to catalyze a chemical reaction and the second one the chemical reaction catalyzed by the enzyme. Therefore, an EC number does not represent the enzyme itself, but only the activity. As a consequence, non-homologous isoenzymes (i.e. with different ancestral origin) may share the same EC number as they catalyze the same enzymatic reaction. In the case of substrate promiscuity, different EC numbers may exist to give precision to the nature of transformed compounds. Otherwise, only one EC number may be available and represents a generic transformation that could occur on different substrates (e.g. alcohol dehydrogenase, hexokinase). The promiscuity aspect of enzymes is extensively described below. Besides, a same chemical transformation may be represented by different EC numbers when, for example, different cofactors are used. This multiplicity between related activities and EC numbers may lead to discrepancies in databases and mask some orphan enzymes. Another point, reported by Green et al., is the ambiguity in the use of incomplete EC numbers that could lead to enzyme annotation errors in public databases. This is because incomplete EC numbers don’t distinguish between the lack of knowledge of the exact substrate specificity of an enzyme and the lack of an official EC number to describe the given activity. Consequently, the use of EC numbers may have introduced some biases in our survey. It should be noticed that the UniProt consortium is making improvements in the representation of the enzymatic activities through Rhea  and UniPathway  databases, which are focused on the definition of chemical reactions and metabolic pathways, respectively.
Statistics on orphan reactions in KEGG and MetaCyc metabolic databases
Total number of non-spontaneous reactions
Number of orphan reactions
Number of reactions in a pathway
Number of orphan reactions in a pathway
Number of orphan reactions having a non orphan pathway neighbour
Number of pathways
Average number of reactions per pathway
Number of pathways with only non orphan reactions
Number of pathways with only orphan reactions
Enzymes are classically studied through metabolic pathways, which are groups of activities taking part in a same biological process. In this survey, we studied the orphan enzyme content and their connectivity at the pathway level. As described previously , there are several key differences between the way the databases represent the notion of a pathway: KEGG pathways are a kind of mosaic of similar pathways predicted in different species; in MetaCyc, the overall reactions in a pathway are supposed to occur in a defined group of species. Therefore, there are 12 times more pathways in MetaCyc than in KEGG, as MetaCyc attempts to provide distinct pathway variants for a given metabolic process (Table 1). An important fraction of pathways (87% in KEGG and 36% in MetaCyc) contains at least one orphan activity. There is no pathway in KEGG containing only orphan enzyme activities, whereas it is the case for about a quarter of the MetaCyc pathways. This is explained by the very large number of reactions in KEGG pathways in comparison to MetaCyc (80 on average per pathway versus 4). Considering pathways containing a mix of orphan and non-orphan activities in KEGG and MetaCyc, an average of 26.0% and 39.5% of the reactions per pathway corresponds to orphan enzymes, respectively (Table 1). These statistics show that an important proportion of pathways are still not completely resolved at the gene level, which limits in silico reconstructions of genome-scale metabolic models [28, 29]. To cope with this problem, computational tools were developed to find candidate genes for these missing enzymes by using genome and metabolic context-based methods [30–32]. The concept of these methods and the illustration of integrated tools using genomic and post-genomic data to link gene and function have been reviewed recently . Another illustration is presented through the MicroScope platform as a combination of CanOE and phylogenetic profile methods [32, 34]. Actually, these in silico predictions have not raised a lot of orphan cases despite the sophistication of the methods and their relative independence from classical sequence based methods. As many orphan enzymes (1,223 reactions in KEGG and 915 in MetaCyc) have pathway neighbors that are orphans themselves, one difficulty is the definition of correct genomic contexts including candidate genes and known enzymes. Furthermore, there is some part of the metabolism with a lot of missing knowledge like glycan and lipid pathways. For example, a number of orphan enzymes still exist in ether lipid metabolism, even if some recent progresses were made .
Local orphan enzymes
From a taxonomic point of view, we propose to make the distinction between global and local orphan enzymes. Orphan enzymes were previously defined as activities having no associated gene in any organism, which we called here global orphans. In addition, a local orphan is an experimentally observed activity in at least one organism of a given clade with only associated sequences in organisms from other clades [36, 37]. To illustrate this concept at the superkingdom level, we present here the example of the EC number 18.104.22.168, the aspartate 4-decarboxylase, which catalyzes the transformation of an L-aspartate in an L-alanine by releasing a molecule of CO2. In UniProt, 327 bacterial proteins are annotated with this EC number, including two SwissProt entries, but no eukaryotic or archaeal sequences can be found. Nevertheless, the aspartate 4-decarboxylase activity has been characterized in various mammalians (e.g. rat, pig, chicken) , making the EC number 22.214.171.124 a local orphan activity in eukaryotes. For the Archaea, there is no associated sequence and no literature evidence of its presence in this superkingdom. Thus, the aspartate 4-decarboxylase activity could be considered as absent in the Archaea.
Potential candidates for local orphan enzymes retrieved by PRIAM
local orphan EC
not observed EC
local orphan EC
not observed EC
local orphan EC
not observed EC
Number of predictable
Number of predicted
Percent of predicted
Number of candidate
Enzyme promiscuity and protein families
Multifunctional enzymes are enzymes capable of playing several roles in metabolism by catalyzing different transformations that may occur in different pathways. Several kinds of multifunctionality can be observed. Some enzymes may show broad substrate specificity. This substrate promiscuity is a feature of enzymes able to catalyze the same chemical reaction on a variety of related compounds . Other enzymes may catalyze different chemical transformations. One can observe proteins having two or more functional domains with different active sites . The association of several domains within a protein, which is generally the result of a gene fusion event during evolution, may facilitate substrate conversion and regulation of the metabolic fluxes. Another origin of this catalytic promiscuity is the special case of moonlighting enzymes . These proteins switch between activities under environmental changes according to their cellular localization, expression in a novel cell type, ligand or cofactor concentrations, oligomerization or complex formation with other proteins. A repository of multitasking proteins was recently set up and several examples of moonlighting enzymes may be explored .
The proportion of multifunctional enzymes may be underestimated [47, 48] and only a few enzymes are described as multifunctional in databases: among the ~250,000 enzymes in Swiss-Prot, 5% are associated with two or more EC numbers and 3% with EC numbers having different classification at third-level. This proportion should dramatically increase when we will find a simpler way to detect them. Recently, a bioinformatic method based on reaction molecular signatures was proposed to predict catalytic and substrate promiscuity . Using this method, a complementary study showed that highly promiscuous enzymes are more likely to be widespread in the tree of life . Because multifunctional enzymes are so difficult to discover and annotate, they represent an interesting and relatively unexplored reservoir to find sequences for orphan enzymes. Quite often, biochemists discover a “new” activity performed by an enzyme known to catalyze other type of reactions . The point is that the characterization of a novel protein generally leads to the discovery of only one function, but does not automatically include a search for all possible additional functions. Nevertheless, the characterization of supplementary in vitro activities does not necessarily imply the elucidation of bona fide in vivo functions.
Despite an observed decrease of the number of orphan enzyme activities over the last ten years, the orphan enzyme challenge remains important: more than 30% of the enzymatic activities reported in the EC classification have no or incomplete sequence information. Though NGS, combined with improvements in sequence analysis methods, produces an exponential growth of genomic data, an explosion in the number of newly discovered activities has not occurred unlike the 80’s when the democratization of molecular biology techniques took place. This lack of knowledge is obviously problematic in the overall comprehension of metabolism and in potential biotechnological applications like biocatalysis.
As shown in our survey and as previously reported , a more systematic use of comparative genomics across superkingdoms may help to solve part of the local orphans. For the global ones, a delay of knowledge in databases still exists and could be resolved by intensive bibliographical searches. In this way, the Orphan Enzyme Project initiative  recently conducted a systematic analysis of databases and publications, and found protein sequences for about 270 presumed orphans among an initial list of 1,122 activities established in 2009 . Similarly to what is done for protein structures with the PDB  and nucleic sequences by the INSDC (International Nucleotide Sequence Database Collaboration) , the design of a central and common scientific framework to submit enzymes with their activities is of priority to reduce the loss of knowledge between publications and databases. Indeed, collaborative initiatives were recently established to discover new activities and enzymes: the Enzyme Function Initiative  which addresses the challenge of assigning reliable functions to enzymes discovered in bacterial genome projects, and the COMBREX project , connecting computational and experimental biologists to improve protein annotation and proposing grants to experimentally validate new functions. These kinds of projects combining in silico and wet lab strategies should lead to a breakthrough in the discovery of new enzymes and activities since classical sequence based methods have lost momentum in function prediction. In fact, several recent studies have successfully applied this approach by exploiting mass-spectrometry or high throughput enzymatic assay experiments and computational methods using sequence similarity networks, genomic contexts, structural modeling with metabolite docking and active site comparison [52, 61, 62]. Another field of research concerns non-protein enzymes. The most well-known are ribozymes and all kinds of protein-RNA complexes, like ribosomes, that are a real challenge to study and extremely hard to discover [63, 64]. The existence of active RNA has been known for a long time, but expertize in this area is far from being as exhaustive as in classical biochemistry. More recently, the discovery of a glycolipid playing a “membrane protein integrase” role in Escherichia coli has pushed back the limits of known catalytic activities . After all, not only should we enlarge the limits of potential catalysts, but also enlarge the limits of the known metabolites. Progress in metabolomics will certainly catalyze the discovery of numerous chemical compounds orphan of activities.
We thank the reviewers for their comments. We have revised the manuscript taking into account their remarks.
Reviewer 1 (First Round): Dr. Michael Galperin
The paper by Sorokina et al. addresses an important question and includes some interesting results. However, I think that in order to justify publication in Biology Direct, the paper needs to be much better written. The current version is something intermediate between a review and a regular research paper and does not make for either a good review or a good research paper. As an example, I would suggest moving Figure 1 to Supplementary Materials (it is not a new result) and moving Figure S2 into the main text (it is a new result).
Authors’ response: Our article is not a regular research article but a review paper written in a format similar to previous studies listed in Figure1. It includes updated analyses of existing data from public databanks that substantially enhance our knowledge about orphan enzymes. We thus decided not to move Figure1to Supplementary Materials as it resumes previous studies. Figure S2 (now S1.2) is an estimation of orphan enzymes at the superkingdom level based on non-curated data from FRENDA and AMENDA whereas Figure6was made using manually curated data. Therefore, we prefer not to move Figure S1.2 to the main text.
In addition, I am afraid that the current version of the manuscript does not really benefit the scientific community as it simply enumerates the enzymes in each category without providing the specific lists of these enzymes. I could support publication of this paper only after the authors include (at least as Supplementary Materials) the lists of global and local orphans from Figure S2. Unless this is done, the data in Figures 2, 3 and 4 cannot be independently verified and the entire manuscript cannot be considered acceptable for publication.
Authors’ response: We added the lists of global and local orphans and proteins in Supplementary Materials 2 and 3.
Finally, the entire paper looks like a promotion for the Orphan Enzymes Project [http://www.orphanenzymes.org, ref. 49]. However, according to the Orphan Enzymes web site, this project is also the subject of an upcoming paper “Finding sequences for over 270 orphan enzymes” (currently in press). The reviewers should have been provided the text of that other paper to ensure that there was no significant overlap between the two.
Authors’ response: We have no relation or contact with the Orphan Enzymes Project and had not access to their upcoming paper at the time of writing the present article. This article is now published and sentences were included in the main text to present their work.
The Abstract needs to be revised to clearly explain what are the new results communicated in this work. Right now, the new results seem to start from “Besides, we extended our study”? Please rewrite the first 4 sentences of the Abstract to explain what exactly was the goal of this work and what exactly has been done.
The statement in the Abstract “We developed a simple strategy to rescue these local orphan enzymes” is totally enigmatic and has to be deleted or reformulated.
The last sentence of the Abstract does not seem relevant to the rest of the text. Please either delete or at least reformulate.
Authors’ response: Part of the abstract has been rewritten according to the reviewer suggestions.
The Introduction could (and should) be made more compact and succinct. That said, the last paragraph of the Introduction contains a much better description of the work presented in this paper than the Abstract does.
Authors’ response: We removed the definition of the EC nomenclature but we think that it is important to keep a description of previous analysis reviews on orphan enzymes in the introduction.
Citations of the enzyme and EC number databases in the Introduction and other sections of the paper present are unfortunately biased. The authors should, at the very least acknowledge the official web sites of the EC classification, the IUBMB list (http://www.chem.qmul.ac.uk/iubmb/enzyme/) and/or the ExplorEnz (http://www.enzyme-database.org, PMID: 18776214) as well as the ENZYME database (http://www.expasy.org/enzyme/ PMID: 10592255), That would also make it unnecessary to explain the organization of the EC system in the Introduction section. INSDC should be cited (PMID: 23180798). The section on Enzyme promiscuity should probably mention the availability of the MultiTaskDB (http://wallace.uab.es/multitask/, PMID: 24253302).
Authors’ response: Suggested references have been added.
Reviewer 2 (First Round): Dr. Daniel Haft
The manuscript submission by Sorokina et al., “Profiling the Orphan Enzymes”, functions fairly well as a review article on the chronology of the growth of EC numbers with and without associations with specific sequences. The authors define a problem space - identifying enzymes that have no representative in some superkingdom -. They introduce a strategy for generating lists of candidate sequences to fill the void. The revised form of the manuscript now provides lists of these candidate sequences in supplementary materials, rather than their count only, and it clearly warns that the associations offered by their technique are in no way validated.
The strategy relies on PRIAM, an update from March 2013. But there is no discussion of how PRIAM itself is formed and whether its design could be appropriate to the task. PRIAM was described in 2003, and relies on MKDOM. Therefore, PRIAM requires an unsupervised domain definition algorithm to find signature regions one enzyme has but another enzyme lacks. The domain could be a C-terminal extension with no relevance to enzyme function, and could be eukaryotic only, but PRIAM would make it a signature. Should this method be used to identify probable “local orphan enzymes” in the archaea? Not without validation.
Other homology strategies might do as well PRIAM or better, such as searching for bi-directional best BLAST hit matches that link a known exemplar of enzyme function in one superkingdom to a homolog in another superkingdom. The PRIAM strategy itself could have been benchmarked somewhat be seeing how much its predictions vary from one version to the next. Readers are strongly cautioned that the output from the PRIAM strategy should be viewed only as anecdotal evidence, appropriate to a review article, that simple homology methods could generate lists of sequences that contain candidates to represent the first extension into a new superkingdom of enzymatic activities that have been assigned to sequences in other superkingdoms.
Authors’ response: This strategy is not a methodological development but just a way to estimate if candidate proteins for local orphans could be retrieved by homology search. We agree that PRIAM profiles have limitations but, as far as we know, it is one of the best tools to track potential conserved domains which are enzyme specific and have a wide coverage of Swiss-Prot enzymes. BBH cannot be computed for all the Swiss-Prot enzymes as many of them are not from complete organisms. As mentioned in the manuscript: “these [PRIAM] predictions cannot be transferred directly without supplementary bioinformatics analyses or experiments”.
As a review, the manuscript did not do justice to the methods that might be used to find orphan enzymes in general, or domain orphans. In particular, Yamada et al. (ref 27) struck me as a landmark demonstration of data mining combined with comparative genomics for finding complete sequence orphans. The method would work even better for superkingdom orphans. Because that work followed predictions with validations, it represents a standard that should be discussed in any review article on matching sequences to orphan EC numbers.
Authors’ response: We introduce the main methods of finding candidate genes for global or local orphans and some of their limitations. But, we do not wish to develop more deeply these methods for three reasons: (1) a complete review of these methods would require a dedicated article (2) a methodological review should be done by a third party since authors of the paper are involved in methodological developments on this topic (i.e. the CANOE method was published the same year as Yamada et al. paper) (3) a review has recently been published and presents a practical description of these methods (El Yacoubi et al. 2014, a reference to this paper was added in our article). For information, the two experimentally tested enzymes in Yamada et al. are not supported by enough evidence to validate that they are good candidates for the two orphan activities: (1) the two tested activities are amino acid transaminases, which are known to have in vitro substrate promiscuity (2) the candidate protein (UniProt AC Q8R5Q4) for the histidine transaminase activity has a TIGRFAM result corresponding to HisC protein (TIGR01141), which catalyzes the transamination of imidazole acetol-phosphate in the context of the histidine biosynthesis. Furthermore, the corresponding gene (TTE2137) is in the hisGDCBHAFI operon confirming that this protein should be involved in the histidine biosynthesis and not in the degradation process via the histidine transaminase activity. (3) the candidate protein (UniProt AC Q8DTM1) shares more than 50% of amino acid identity with biochemically characterized aspartate aminotransferases (UniProt ACs P23034, Q59228). This activity is more coherent with the asparaginyl tRNA synthetase genomic context than the asparagine aminotransferase activity proposed by Yamada et al., an activity described only in eukaryotes for asparagine degradation. These two cases are really good examples to illustrate the difficulty in interpreting in vitro activities to elucidate bona fide in vivo functions.
The work introduces a workflow for using PRIAM to find sequences that might resolve numbers of local enzyme orphans. The lack of any testing of the workflow’s results or consideration of whether PRIAM’s design makes it a good choice was a problem. The revision, including author responses to the reviews, helps cement that this work serves as a review article only, and no tested new method is presented. Even in the revised form, the discussion of the PRIAM workflow is a bit troubling. Does the article title, “Profiling the Orphan Enzymes”, refer to PRIAM profiles as used in the untested workflow? If so, a revised title might be more appropriate.
Authors’ response: The title is not related to PRIAM profiles. The aim of our review is to analyze and discuss the orphan enzyme problem in the light of the current knowledge in public databanks.
Reviewer 3 (First Round): Dr.Daniel Kahn
This reviewer provided no comments for publication.
Reviewer 1 (Second Round): Dr. Michael Galperin
Previous authors’ response: We added the lists of global and local orphans and proteins in Supplementary Materials 2 and 3.
Enzymes (EC numbers) that are not associated with any sequences are referred to as “global orphans” even though many (probably most) of these enzymes have been described in a single species, or a group of closely related species, and therefore represent “lineage-specific orphans”, rather than “global orphans”. It would be helpful to explain this in the text to avoid confusion.
Authors’ response: For the definitions of global and local orphans, we use the same as the ones of Orth et al. 2010. These definitions are given in the main text. For global orphans, it is very difficult to estimate if they are mostly associated to specific lineages as experimental data is limited and is far from covering the metabolic diversity of living organisms.
Enzymes (EC numbers) that have not been reported in bacteria are referred to in Table S2.3 as “Missing enzymes in Bacteria”. In all previously published literature, “missing enzymes” referred to the enzymatic activities that are expected - or known - to be present in at least some bacteria but have not yet been assigned to any sequence. Thus, “Missing enzymes in Bacteria” are the ones that have been reported in certain eukaryotes and are not even expected to be encoded in any bacteria. As a result, there are 1521 enzymes “missing in Bacteria” and 3773 enzymes “missing in Archaea”. Again, if the authors choose to keep this - unconventional and counterintuitive - group name, they should explain it in the text to avoid confusion.
Authors’ response: We agree with the reviewer that the term “missing” is confusing. We have replaced “missing” by “not observed” in the additional files and in the main text.
Although the text has been significantly improved, I remain puzzled by the expression “Rescuing the local orphans”. What do the authors mean by “rescuing” here, probably not something that is covered by the existing dictionaries?
Authors’ response: The term “rescuing” has been removed.
Reviewer 2 (Second Round): Dr. Daniel Haft
The revised form of the article makes it clearer that it is a review, not original research, and that a method they introduce produces only a suggestive view, not scientifically validated results. But it is still a little troubling. The title seems to speak of the new method, and there is no peer-reviewed endorsement of that method her.
Authors’ response: These points are discussed in the first round of the review.
We would like to thank Patrick Bowe and Andrew Tolonen for their helpful suggestions on the manuscript, Karine Bastard for her support, presence and constructive comments during all this work and Marcel Salanoubat for reading this manuscript. We thank also François Le Fèvre for helping us with MetaCyc data extraction. This work was not supported by any funding.
- Schnoes AM, Brown SD, Dodevski I, Babbitt PC: Annotation error in public databases: misannotation of molecular function in enzyme superfamilies. PLoS Comput Biol. 2009, 5: e1000605-PubMedPubMed CentralView ArticleGoogle Scholar
- Payen A, Perzoz J: Mémoire sur la diastase, les principaux produits de ses rèactions, et leurs applications aux arts industriels. Annales de la chimie et de la physique. 1833, 53: 73-92.Google Scholar
- Tipton K, Boyce S: History of the enzyme nomenclature system. Bioinformatics. 2000, 16: 34-40.PubMedView ArticleGoogle Scholar
- Enzyme nomenclature. http://www.chem.qmul.ac.uk/iubmb/enzyme/,
- McDonald AG, Boyce S, Tipton KF: ExplorEnz: the primary source of the IUBMB enzyme list. Nucleic Acids Res. 2009, 37 (Database issue): D593-7.PubMedPubMed CentralView ArticleGoogle Scholar
- Bairoch A, Universitaire CM, Servet M: The ENZYME database in 2000. 2000, 28: 304-5.Google Scholar
- Roberts RJ: Identifying protein function–a call for community action. PLoS Biol. 2004, 2: E42-PubMedPubMed CentralView ArticleGoogle Scholar
- Karp PD: Call for an enzyme genomics initiative. Genome Biol. 2004, 5: 401-PubMedPubMed CentralView ArticleGoogle Scholar
- Lespinet O, Labedan B: Orphan enzymes?. Sci. 2005, 307: 42-View ArticleGoogle Scholar
- Lespinet O, Labedan B: ORENZA: a web resource for studying ORphan ENZyme activities. BMC Bioinformatics. 2006, 7: 436-PubMedPubMed CentralView ArticleGoogle Scholar
- Lespinet O, Labedan B: Orphan enzymes could be an unexplored reservoir of new drug targets. Drug Discov Today. 2006, 11: 300-5.PubMedView ArticleGoogle Scholar
- Chen L, Vitkup D: Distribution of orphan metabolic activities. Trends Biotechnol. 2007, 25: 343-8.PubMedView ArticleGoogle Scholar
- Fleischmann A, Darsow M, Degtyarenko K, Fleischmann W, Boyce S, Axelsen KB, Bairoch A, Schomburg D, Tipton KF, Apweiler R: IntEnz, the integrated relational enzyme database. Nucleic Acids Res. 2004, 32 (Database issue): D434-7.PubMedPubMed CentralView ArticleGoogle Scholar
- Apweiler R, Bairoch A, Wu CH, Barker WC, Boeckmann B, Ferro S, Gasteiger E, Huang H, Lopez R, Magrane M, Martin MJ, Natale DA, O’Donovan C, Redaschi N, Yeh L-SL: UniProt: the Universal protein knowledgebase. Nucleic Acids Res. 2004, 32 (Database issue): D115-9.PubMedPubMed CentralView ArticleGoogle Scholar
- Pouliot Y, Karp PD: A survey of orphan enzyme activities. BMC Bioinformatics. 2007, 8: 244-PubMedPubMed CentralView ArticleGoogle Scholar
- Shearer AG, Altman T, Rhee CD: Finding sequences for over 270 orphan enzymes. PLoS One. 2014, 9: e97250-PubMedPubMed CentralView ArticleGoogle Scholar
- Karp PD: What we do not know about sequence analysis and sequence database. Bioinformatics. 1998, 14: 753-4.PubMedView ArticleGoogle Scholar
- Karp PD, Paley S, Zhu J: Database verification studies of SWISS-PROT and GenBank. Bioinformatics. 2001, 17: 526-32.PubMedView ArticleGoogle Scholar
- Lee D, Redfern O, Orengo C: Predicting protein function from sequence and structure. Nat Rev Mol Cell Biol. 2007, 8: 995-1005.PubMedView ArticleGoogle Scholar
- Hanson AD, Pribat A, Waller JC, De Crécy-Lagard V: “Unknown” proteins and “orphan” enzymes: the missing half of the engineering parts list–and how to find it. Biochem J. 2010, 425: 1-11.View ArticleGoogle Scholar
- Green ML, Karp PD: Genome annotation errors in pathway databases due to semantic ambiguity in partial EC numbers. Nucleic Acids Res. 2005, 33: 4035-9.PubMedPubMed CentralView ArticleGoogle Scholar
- Alcántara R, Axelsen KB, Morgat A, Belda E, Coudert E, Bridge A, Cao H, De Matos P, Ennis M, Turner S, Owen G, Bougueleret L, Xenarios I, Steinbeck C: Rhea–a manually curated resource of biochemical reactions. Nucleic Acids Res. 2012, 40 (Database issue): D754-60.PubMedPubMed CentralView ArticleGoogle Scholar
- Morgat A, Coissac E, Coudert E, Axelsen KB, Keller G, Bairoch A, Bridge A, Bougueleret L, Xenarios I, Viari A: UniPathway: a resource for the exploration and annotation of metabolic pathways. Nucleic Acids Res. 2012, 40 (Database issue): D761-9.PubMedPubMed CentralView ArticleGoogle Scholar
- Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28: 27-30.PubMedPubMed CentralView ArticleGoogle Scholar
- Caspi R, Altman T, Dreher K, Fulcher CA, Subhraveti P, Keseler IM, Kothari A, Krummenacker M, Latendresse M, Mueller LA, Ong Q, Paley S, Pujar A, Shearer AG, Travers M, Weerasinghe D, Zhang P, Karp PD: The MetaCyc database of metabolic pathways and enzymes and the BioCyc collection of pathway/genome databases. Nucleic Acids Res. 2012, 40 (Database issue): D742-53.PubMedPubMed CentralView ArticleGoogle Scholar
- Altman T, Travers M, Kothari A, Caspi R, Karp PD: A systematic comparison of the MetaCyc and KEGG pathway databases. BMC Bioinformatics. 2013, 14: 112-PubMedPubMed CentralView ArticleGoogle Scholar
- Lespinet O, Labedan B: Lespinet: Puzzling over orphan enzymes. Cell Mol Life Sci. 2006, 63: 517-23.PubMedView ArticleGoogle Scholar
- Durot M, Bourguignon P-Y, Schachter V: Genome-scale models of bacterial metabolism: reconstruction and applications. FEMS Microbiol Rev. 2009, 33: 164-90.PubMedPubMed CentralView ArticleGoogle Scholar
- Hyduke DR, Lewis NE, Palsson BØ: Analysis of omics data with genome-scale models of metabolism. Mol Biosyst. 2013, 9: 167-74.PubMedPubMed CentralView ArticleGoogle Scholar
- Kharchenko P, Chen L, Freund Y, Vitkup D, Church GM: Identifying metabolic enzymes with multiple types of association evidence. BMC Bioinformatics. 2006, 7: 177-PubMedPubMed CentralView ArticleGoogle Scholar
- Yamada T, Waller AS, Raes J, Zelezniak A, Perchat N, Perret A, Salanoubat M, Patil KR, Weissenbach J, Bork P: Prediction and identification of sequences coding for orphan enzymes using genomic and metagenomic neighbours. Mol Syst Biol. 2012, 8: 581-PubMedPubMed CentralView ArticleGoogle Scholar
- Smith AAT, Belda E, Viari A, Medigue C, Vallenet D: The CanOE strategy: integrating genomic and metabolic contexts across multiple prokaryote genomes to find candidate genes for orphan enzymes. PLoS Comput Biol. 2012, 8: e1002540-PubMedPubMed CentralView ArticleGoogle Scholar
- El Yacoubi B, De Crécy-Lagard V: Integrative data-mining tools to link gene and function. Methods Mol Biol. 2014, 1101: 43-66.PubMedView ArticleGoogle Scholar
- Vallenet D, Belda E, Calteau A, Cruveiller S, Engelen S, Lajus A, Le Fèvre F, Longin C, Mornico D, Roche D, Rouy Z, Salvignol G, Scarpelli C, Thil Smith AA, Weiman M, Médigue C: MicroScope–an integrated microbial resource for the curation and comparative analysis of genomic and metabolic data. Nucleic Acids Res. 2013, 41 (Database issue): D636-47.PubMedPubMed CentralView ArticleGoogle Scholar
- Watschinger K, Werner ER: Orphan enzymes in ether lipid metabolism. Biochimie. 2013, 95: 59-65.PubMedPubMed CentralView ArticleGoogle Scholar
- Orth JD, Palsson BØ: Systematizing the generation of missing metabolic knowledge. Biotechnol Bioeng. 2010, 107: 403-12.PubMedPubMed CentralView ArticleGoogle Scholar
- Chen L, Vitkup D: Predicting genes for orphan metabolic activities using phylogenetic profiles. Genome Biol. 2006, 7: R17-PubMedPubMed CentralView ArticleGoogle Scholar
- Rathod PK, Fellman JH: Identification of mammalian aspartate-4-decarboxylase. Arch Biochem Biophys. 1985, 238: 435-46.PubMedView ArticleGoogle Scholar
- Schomburg I, Chang A, Placzek S, Söhngen C, Rother M, Lang M, Munaretto C, Ulas S, Stelzer M, Grote A, Scheer M, Schomburg D: BRENDA in 2013: integrated reactions, kinetic data, enzyme function data, improved disease classification: new options and contents in BRENDA. Nucleic Acids Res. 2013, 41 (Database issue): D764-72.PubMedPubMed CentralView ArticleGoogle Scholar
- Omelchenko MV, Galperin MY, Wolf YI, Koonin EV: Non-homologous isofunctional enzymes: a systematic analysis of alternative solutions in enzyme evolution. Biol Direct. 2010, 5: 31-PubMedPubMed CentralView ArticleGoogle Scholar
- Claudel-Renard C: Enzyme-specific profiles for genome annotation: PRIAM. Nucleic Acids Res. 2003, 31: 6633-9.PubMedPubMed CentralView ArticleGoogle Scholar
- Plata G, Fuhrer T, Hsiao T-L, Sauer U, Vitkup D: Global probabilistic annotation of metabolic networks enables enzyme discovery. Nat Chem Biol. 2012, 8: 848-54.PubMedPubMed CentralView ArticleGoogle Scholar
- Khersonsky O, Tawfik DS: Enzyme promiscuity: a mechanistic and evolutionary perspective. Annu Rev Biochem. 2010, 79: 471-505.PubMedView ArticleGoogle Scholar
- Hawkins AR, Lamb HK: The molecular biology of multidomain proteins. selected examples. Eur J Biochem. 1995, 232: 7-18.PubMedView ArticleGoogle Scholar
- Jeffery CJ: Moonlighting proteins: old proteins learning new tricks. Trends Genet. 2003, 19: 415-7.PubMedView ArticleGoogle Scholar
- Hernández S, Ferragut G, Amela I, Perez-Pons J, Piñol J, Mozo-Villarias A, Cedano J, Querol E: MultitaskProtDB: a database of multitasking proteins. Nucleic Acids Res. 2014, 42 (Database issue): D517-20.PubMedPubMed CentralView ArticleGoogle Scholar
- Cheng X-Y, Huang W-J, Hu S-C, Zhang H-L, Wang H, Zhang J-X, Lin H-H, Chen Y-Z, Zou Q, Ji Z-L: A global characterization and identification of multifunctional enzymes. PLoS One. 2012, 7: e38979-PubMedPubMed CentralView ArticleGoogle Scholar
- Jia B, Cheong G-W, Zhang S: Multifunctional enzymes in archaea: promiscuity and moonlight. Extremophiles : life under extreme conditions. 2013, 17: 1-11.View ArticleGoogle Scholar
- Carbonell P, Faulon J-L: Molecular signatures-based prediction of enzyme promiscuity. Bioinformatics. 2010, 26: 2012-9.PubMedView ArticleGoogle Scholar
- Carbonell P, Lecointre G, Faulon J-L: Origins of specificity and promiscuity in metabolic networks. J Biol Chem. 2011, 286: 43994-4004.PubMedPubMed CentralView ArticleGoogle Scholar
- Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, Boursnell C, Pang N, Forslund K, Ceric G, Clements J, Heger A, Holm L, Sonnhammer ELL, Eddy SR, Bateman A, Finn RD: The Pfam protein families database. Nucleic Acids Res. 2012, 40 (Database issue): D290-301.PubMedPubMed CentralView ArticleGoogle Scholar
- Bastard K, Smith AAT, Vergne-Vaxelaire C, Perret A, Zaparucha A, De Melo-Minardi R, Mariage A, Boutard M, Debard A, Lechaplais C, Pelle C, Pellouin V, Perchat N, Petit J-L, Kreimeyer A, Medigue C, Weissenbach J, Artiguenave F, De Berardinis V, Vallenet D, Salanoubat M: Revealing the hidden functional diversity of an enzyme family. Nature chemical biology. 2013, 10: 42-49.PubMedView ArticleGoogle Scholar
- Furnham N, Sillitoe I, Holliday GL, Cuff AL, Laskowski RA, Orengo CA, Thornton JM: Exploring the evolution of novel enzyme functions within structurally defined protein superfamilies. PLoS Comput Biol. 2012, 8: e1002403-PubMedPubMed CentralView ArticleGoogle Scholar
- Sillitoe I, Cuff AL, Dessailly BH, Dawson NL, Furnham N, Lee D, Lees JG, Lewis TE, Studer RA, Rentzsch R, Yeats C, Thornton JM, Orengo CA: New functional families (FunFams) in CATH to improve the mapping of conserved functional sites to 3D structures. Nucleic Acids Res. 2013, 41 (Database issue): D490-8.PubMedPubMed CentralView ArticleGoogle Scholar
- Almonacid DE, Yera ER, Mitchell JBO, Babbitt PC: Quantitative comparison of catalytic mechanisms and overall reactions in convergently evolved enzymes: implications for classification of enzyme function. PLoS Comput Biol. 2010, 6: e1000700-PubMedPubMed CentralView ArticleGoogle Scholar
- Orphan enzyme project. http://www.orphanenzymes.org/,
- Berman HM: The protein data bank. Nucleic Acids Res. 2000, 28: 235-42.PubMedPubMed CentralView ArticleGoogle Scholar
- Nakamura Y, Cochrane G, Karsch-Mizrachi I: The International nucleotide sequence database collaboration. Nucleic Acids Res. 2013, 41 (Database issue): D21-4.PubMedPubMed CentralView ArticleGoogle Scholar
- Gerlt JA, Allen KN, Almo SC, Armstrong RN, Babbitt PC, Cronan JE, Dunaway-Mariano D, Imker HJ, Jacobson MP, Minor W, Poulter CD, Raushel FM, Sali A, Shoichet BK, Sweedler JV: The enzyme function initiative. Biochemistry. 2011, 50: 9950-62.PubMedPubMed CentralView ArticleGoogle Scholar
- Anton BP, Chang Y-C, Brown P, Choi H-P, Faller LL, Guleria J, Hu Z, Klitgord N, Levy-Moonshine A, Maksad A, Mazumdar V, McGettrick M, Osmani L, Pokrzywa R, Rachlin J, Swaminathan R, Allen B, Housman G, Monahan C, Rochussen K, Tao K, Bhagwat AS, Brenner SE, Columbus L, De Crécy-Lagard V, Ferguson D, Fomenkov A, Gadda G, Morgan RD, Osterman AL, et al: The COMBREX project: design, methodology, and initial results. PLoS Biol. 2013, 11: e1001638-PubMedPubMed CentralView ArticleGoogle Scholar
- Ramkissoon KR, Miller JK, Ojha S, Watson DS, Bomar MG, Galande AK, Shearer AG: Rapid identification of sequences for orphan enzymes to power accurate protein annotation. PLoS One. 2013, 8: e84508-PubMedPubMed CentralView ArticleGoogle Scholar
- Zhao S, Kumar R, Sakai A, Vetting MW, Wood BM, Brown S, Bonanno JB, Hillerich BS, Seidel RD, Babbitt PC, Almo SC, Sweedler JV, Gerlt JA, Cronan JE, Jacobson MP: Discovery of new enzymes and metabolic pathways by using structure and genome context. Nature. 2013, 502: 698-702.PubMedPubMed CentralView ArticleGoogle Scholar
- Strobel SA, Cochrane JC: RNA catalysis: ribozymes, ribosomes, and riboswitches. Curr Opin Chem Biol. 2007, 11: 636-43.PubMedPubMed CentralView ArticleGoogle Scholar
- Serganov A, Nudler E: A decade of riboswitches. Cell. 2013, 152: 17-24.PubMedPubMed CentralView ArticleGoogle Scholar
- Nishiyama K, Maeda M, Yanagisawa K, Nagase R, Komura H, Iwashita T, Yamagaki T, Kusumoto S, Tokuda H, Shimamoto K: MPIase is a glycolipozyme essential for membrane protein integration. Nat Commun. 2012, 3: 1260-PubMedPubMed CentralView 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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.