BrEPS: a flexible and automatic protocol to compute enzyme-specific sequence profiles for functional annotation

Background  

Models for the simulation of metabolic networks require the accurate prediction of enzyme function. Based on a genomic sequence, enzymatic functions of gene products are today mainly predicted by sequence database searching and operon analysis. Other methods can support these techniques: We have developed an automatic method "BrEPS" that creates highly specific sequence patterns for the functional annotation of enzymes.     

Metabolomics- and proteomics-assisted genome annotation and analysis of the draft metabolic network of Chlamydomonas reinhardtii

We present an integrated analysis of the molecular repertoire of Chlamydomonas reinhardtii under reference conditions. Bioinformatics annotation methods combined with GCxGC/MS-based metabolomics and LC/MS-based shotgun proteomics profiling technologies have been applied to characterize abundant proteins and metabolites, resulting in the detection of 1069 proteins and 159 metabolites. Of the measured proteins, 204 currently do not have EST sequence support; thus a significant portion of the proteomics-detected proteins provide evidence for the validity of in silico gene models. Furthermore, the generated peptide data lend support to the validity of a number of proteins currently in the proposed model stage. By integrating genomic annotation information with experimentally identified metabolites and proteins, we constructed a draft metabolic network for Chlamydomonas. Computational metabolic modeling allowed an identification of missing enzymatic links. Some experimentally detected metabolites are not producible by the currently known and annotated enzyme set, thus suggesting entry points for further targeted gene discovery or biochemical pathway research. All data sets are made available as supplementary material as well as web-accessible databases and within the functional context via the Chlamydomonas-adapted MapMan annotation platform. Information of identified peptides is also available directly via the JGI-Chlamydomonas genomic resource database (http://genome.jgi-psf.org/Chlre3/Chlre3.home.html).     

The CanOE strategy: integrating genomic and metabolic contexts across multiple prokaryote genomes to find candidate genes for orphan enzymes

Of all biochemically characterized metabolic reactions formalized by the IUBMB, over one out of four have yet to be associated with a nucleic or protein sequence, i.e. are sequence-orphan enzymatic activities. Few bioinformatics annotation tools are able to propose candidate genes for such activities by exploiting context-dependent rather than sequence-dependent data, and none are readily accessible and propose result integration across multiple genomes. Here, we present CanOE (Candidate genes for Orphan Enzymes), a four-step bioinformatics strategy that proposes ranked candidate genes for sequence-orphan enzymatic activities (or orphan enzymes for short). The first step locates "genomic metabolons", i.e. groups of co-localized genes coding proteins catalyzing reactions linked by shared metabolites, in one genome at a time. These metabolons can be particularly helpful for aiding bioanalysts to visualize relevant metabolic data. In the second step, they are used to generate candidate associations between un-annotated genes and gene-less reactions. The third step integrates these gene-reaction associations over several genomes using gene families, and summarizes the strength of family-reaction associations by several scores. In the final step, these scores are used to rank members of gene families which are proposed for metabolic reactions. These associations are of particular interest when the metabolic reaction is a sequence-orphan enzymatic activity. Our strategy found over 60,000 genomic metabolons in more than 1,000 prokaryote organisms from the MicroScope platform, generating candidate genes for many metabolic reactions, of which more than 70 distinct orphan reactions. A computational validation of the approach is discussed. Finally, we present a case study on the anaerobic allantoin degradation pathway in Escherichia coli K-12.     

Genome-wide metabolic re-annotation of Ashbya gossypii: new insights into its metabolism through a comparative analysis with Saccharomyces cerevisiae and Kluyveromyces lactis

Background

Ashbya gossypii is an industrially relevant microorganism traditionally used for riboflavin production. Despite the high gene homology and gene order conservation comparatively with Saccharomyces cerevisiae, it presents a lower level of genomic complexity. Its type of growth, placing it among filamentous fungi, questions how close it really is from the budding yeast, namely in terms of metabolism, therefore raising the need for an extensive and thorough study of its entire metabolism. This work reports the first manual enzymatic genome-wide re-annotation of A. gossypii as well as the first annotation of membrane transport proteins.

Results

After applying a developed enzymatic re-annotation pipeline, 847 genes were assigned with metabolic functions. Comparatively to KEGG’s annotation, these data corrected the function for 14% of the common genes and increased the information for 52 genes, either completing existing partial EC numbers or adding new ones. Furthermore, 22 unreported enzymatic functions were found, corresponding to a significant increase in the knowledge of the metabolism of this organism. The information retrieved from the metabolic re-annotation and transport annotation was used for a comprehensive analysis of A. gossypii’s metabolism in comparison to the one of S. cerevisiae (post-WGD – whole genome duplication) and Kluyveromyces lactis (pre-WGD), suggesting some relevant differences in several parts of their metabolism, with the majority being found for the metabolism of purines, pyrimidines, nitrogen and lipids. A considerable number of enzymes were found exclusively in A. gossypii comparatively with K. lactis (90) and S. cerevisiae (13). In a similar way, 176 and 123 enzymatic functions were absent on A. gossypii comparatively to K. lactis and S. cerevisiae, respectively, confirming some of the well-known phenotypes of this organism.

Conclusions

This high quality metabolic re-annotation, together with the first membrane transporters annotation and the metabolic comparative analysis, represents a new important tool for the study and better understanding of A. gossypii’s metabolism.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-810) contains supplementary material, which is available to authorized users.     

Environmental genomics reveals a functional chlorite dismutase in the nitrite-oxidizing bacterium 'Candidatus Nitrospira defluvii'

Nitrite-oxidizing bacteria of the genus Nitrospira are ubiquitous in natural ecosystems and also in wastewater treatment plants. Nitrospira are members of a distinct phylum, not closely related to other nitrifiers, and no genomic sequences from this genus have been available so far. Here we applied an environmental genomics approach to sequence and assemble a 137 kbp-long genome fragment of 'Candidatus Nitrospira defluvii', which had been enriched from activated sludge and belongs to Nitrospira sublineage I without isolated representatives. The annotation of this contig, which carried the 16S rRNA gene of N. defluvii, offered first insight into the genome of Nitrospira. Surprisingly, we found a gene similar to genes encoding chlorite dismutase (CLD), an enzyme degrading chlorite (ClO(2)(-)) to Cl(-) and O(2). To date, CLDs with high catalytic activity have been found only in perchlorate- and chlorate-reducing bacteria but not in nitrifiers. Heterologous expression in E. coli followed by enzymatic tests confirmed that this gene of Nitrospira encodes a highly active CLD, which is also expressed in situ by Nitrospira, indicating that this nitrite oxidizer might be involved in the bioremediation of perchlorate and chlorite. Phylogenetic analyses showed that CLD and related proteins are widely distributed among the Bacteria and Archaea, and indicated that this enzyme family appeared relatively early in evolution, has been subject to functional diversification and might play yet unknown roles in microbial metabolism.     

Metabolomic strategies for the identification of new enzyme functions and metabolic pathways

Recent technological advances in accurate mass spectrometry and data analysis have revolutionized metabolomics experimentation. Activity-based and global metabolomic profiling methods allow simultaneous and rapid screening of hundreds of metabolites from a variety of chemical classes, making them useful tools for the discovery of novel enzymatic activities and metabolic pathways. By using the metabolome of the relevant organism or close species, these methods capitalize on biological relevance, avoiding the assignment of artificial and non-physiological functions. This review discusses state-of-the-art metabolomic approaches and highlights recent examples of their use for enzyme annotation, discovery of new metabolic pathways, and gene assignment of orphan metabolic activities across diverse biological sources.     

Prediction and identification of sequences coding for orphan enzymes using genomic and metagenomic neighbours

Despite the current wealth of sequencing data, one‐third of all biochemically characterized metabolic enzymes lack a corresponding gene or protein sequence, and as such can be considered orphan enzymes. They represent a major gap between our molecular and biochemical knowledge, and consequently are not amenable to modern systemic analyses. As 555 of these orphan enzymes have metabolic pathway neighbours, we developed a global framework that utilizes the pathway and (meta)genomic neighbour information to assign candidate sequences to orphan enzymes. For 131 orphan enzymes (37% of those for which (meta)genomic neighbours are available), we associate sequences to them using scoring parameters with an estimated accuracy of 70%, implying functional annotation of 16 345 gene sequences in numerous (meta)genomes. As a case in point, two of these candidate sequences were experimentally validated to encode the predicted activity. In addition, we augmented the currently available genome‐scale metabolic models with these new sequence–function associations and were able to expand the models by on average 8%, with a considerable change in the flux connectivity patterns and improved essentiality prediction.     

Thio-heterocylic naphthalimides with aminoalkyl side chains: novel alternative tools for photodegradation of genomic DNA without impairment on bioactivities of proteins

Thio-heterocylic naphthalimides (R1-R5) were designed, synthesized and evaluated as nonmetallic and long-wavelength photocleavers. Some of them showed highly efficient abilities in the degradation of plasmid and genomic DNA under the mild conditions without obvious impairment on the proteins' bioactivities, when compared with frequently applied nucleic acids removal reagents or precipitants. Their differences in photodegradation selectivity to DNA rather than proteins were dependent on their photodamage mechanisms and binding modes with bio-macromolecules. When maize genomic DNA was used as substrate, 2.38 x 10(-4) M of R5 exhibited the nuclease activity of 8 Unit DNase I, R5 has some characteristic as a typical catalyst as no consumption after two cycles of the photodegradation for DNA. The experiments of enzymatic activity assay and immunology activity analysis showed that R5 was safe to proteins, suggesting its potential in the removal of transgenic material during the preparation of bioactive proteins or enzyme preparations.     

Surprising arginine biosynthesis: a reappraisal of the enzymology and evolution of the pathway in microorganisms.

Major aspects of the pathway of de novo arginine biosynthesis via acetylated intermediates in microorganisms must be revised in light of recent enzymatic and genomic investigations. The enzyme N-acetylglutamate synthase (NAGS), which used to be considered responsible for the first committed step of the pathway, is present in a limited number of bacterial phyla only and is absent from Archaea. In many Bacteria, shorter proteins related to the Gcn5-related N-acetyltransferase family appear to acetylate l-glutamate; some are clearly similar to the C-terminal, acetyl-coenzyme A (CoA) binding domain of classical NAGS, while others are more distantly related. Short NAGSs can be single gene products, as in Mycobacterium spp. and Thermus spp., or fused to the enzyme catalyzing the last step of the pathway (argininosuccinase), as in members of the Alteromonas-Vibrio group. How these proteins bind glutamate remains to be determined. In some Bacteria, a bifunctional ornithine acetyltransferase (i.e., using both acetylornithine and acetyl-CoA as donors of the acetyl group) accounts for glutamate acetylation. In many Archaea, the enzyme responsible for glutamate acetylation remains elusive, but possible connections with a novel lysine biosynthetic pathway arose recently from genomic investigations. In some Proteobacteria (notably Xanthomonadaceae) and Bacteroidetes, the carbamoylation step of the pathway appears to involve N-acetylornithine or N-succinylornithine rather than ornithine. The product N-acetylcitrulline is deacetylated by an enzyme that is also involved in the provision of ornithine from acetylornithine; this is an important metabolic function, as ornithine itself can become essential as a source of other metabolites. This review insists on the biochemical and evolutionary implications of these findings.     

Transient expression of human neutral alpha-glucosidase AB (glucosidase II) in enzyme-deficient mouse lymphoma cells

To define new methods for gene isolation exploiting mutant mammalian cells we transformed a mutant mouse cell line deficient in glucosidase II with total human genomic DNA and detected transient expression of the human glucosidase II gene. Maximum gene expression was detected 48 h after addition of DNA as a 2.5-fold increase in neutral alpha-glucosidase activity (2.47 +/- 0.15, n = 4). When mutant mouse DNA was used for transformation, no increase in enzyme activity was seen. The increased enzyme activity was due to expression of the human gene product. Thus, by rocket immunoelectrophoresis, cells transformed with human DNA yielded a "rocket" which reacted with antibody to human but not to mouse glucosidase II and which hydrolyzed substrate in situ. Specific DNA sequences were required for expression of the enzyme activity, since digestion of DNA with EcoRI and SstI rendered the DNA ineffective for eliciting expression of the enzyme, while digestion of DNA with BamHI and XhoI did not affect the increase. Transfection with intact phage from a human genomic DNA library also resulted in transient expression of the human gene. These results demonstrate the feasibility of detecting, by enzymatic assay, transient expression of a human gene for an intracellular enzyme following DNA-mediated transformation both with total human DNA and with intact phage from a human recombinant library. This system could be used as an assay for isolation of a gene from a genomic library by sibling selection.     

KEGG: kyoto encyclopedia of genes and genomes

KEGG (Kyoto Encyclopedia of Genes and Genomes) is a knowledge base for systematic analysis of gene functions, linking genomic information with higher order functional information. The genomic information is stored in the GENES database, which is a collection of gene catalogs for all the completely sequenced genomes and some partial genomes with up-to-date annotation of gene functions. The higher order functional information is stored in the PATHWAY database, which contains graphical representations of cellular processes, such as metabolism, membrane transport, signal transduction and cell cycle. The PATHWAY database is supplemented by a set of ortholog group tables for the information about conserved subpathways (pathway motifs), which are often encoded by positionally coupled genes on the chromosome and which are especially useful in predicting gene functions. A third database in KEGG is LIGAND for the information about chemical compounds, enzyme molecules and enzymatic reactions. KEGG provides Java graphics tools for browsing genome maps, comparing two genome maps and manipulating expression maps, as well as computational tools for sequence comparison, graph comparison and path computation. The KEGG databases are daily updated and made freely available (http://www. genome.ad.jp/kegg/).     

The enzymatic nature of an anonymous protein sequence cannot reliably be inferred from superfamily level structural information alone

As the largest fraction of any proteome does not carry out enzymatic functions, and in order to leverage 3D structural data for the annotation of increasingly higher volumes of sequence data, we wanted to assess the strength of the link between coarse grained structural data (i.e., homologous superfamily level) and the enzymatic versus non-enzymatic nature of protein sequences. To probe this relationship, we took advantage of 41 phylogenetically diverse (encompassing 11 distinct phyla) genomes recently sequenced within the GEBA initiative, for which we integrated structural information, as defined by CATH, with enzyme level information, as defined by Enzyme Commission (EC) numbers. This analysis revealed that only a very small fraction (about 1%) of domain sequences occurring in the analyzed genomes was found to be associated with homologous superfamilies strongly indicative of enzymatic function. Resorting to less stringent criteria to define enzyme versus non-enzyme biased structural classes or excluding highly prevalent folds from the analysis had only modest effect on this proportion. Thus, the low genomic coverage by structurally anchored protein domains strongly associated to catalytic activities indicates that, on its own, the power of coarse grained structural information to infer the general property of being an enzyme is rather limited.     

Computational identification of significantly regulated metabolic reactions by integration of data on enzyme activity and gene expression

The concentrations and catalytic activities of enzymes control metabolic rates. Previous studies have focused on enzyme concentrations because there are no genome-wide techniques used for the measurement of enzyme activity. We propose a method for evaluating the significance of enzyme activity by integrating metabolic network topologies and genome-wide microarray gene expression profiles. We quantified the enzymatic activity of reactions and report the 388 significant reactions in five perturbation datasets. For the 388 enzymatic reactions, we identified 70 that were significantly regulated (P-value < 0.001). Thirty-one of these reactions were part of anaerobic metabolism, 23 were part of low-pH aerobic metabolism, 8 were part of high-pH anaerobic metabolism, 3 were part of low-pH aerobic reactions, and 5 were part of high-pH anaerobic metabolism.     

油料作物EST资源的生物信息学分析

利用生物信息学方法,收集整理GenBank数据库中截至2008年5月收录的油料作物油菜、花生、芝麻、大豆、向日葵、蓖麻、亚麻、棕榈等八种油料作物的表达序列标签(EST)序列信息,共获得1,185,911条EST序列,使用Crosmatch、RepeatMask-er、Phrap、CAP3、EMBOSS、Blast、EST-pipeline、ORF finder、Interproscan、blast2go、IdentiCS等软件,基于Linux操作系统,进行了综合及分类分析。共获得289,892条UniEST序列,通过以上对EST序列信息的基因注释信息,筛选出与油脂代谢相关的基因信息,并以此为基础构建了油料作物油脂代谢途径比较结构图。本研究为油料作物油脂代谢相关基因数据库的构建和不同油料作物油脂代谢异同的比较打下基础。     

Cloning, nucleotide sequence, and expression of the Escherichia coli gene encoding carnitine dehydratase.

Carnitine dehydratase from Escherichia coli O44 K74 is an inducible enzyme detectable in cells grown anaerobically in the presence of L-(-)-carnitine or crotonobetaine. The purified enzyme catalyzes the dehydration of L-(-)-carnitine to crotonobetaine (H. Jung, K. Jung, and H.-P. Kleber, Biochim. Biophys. Acta 1003:270-276, 1989). The caiB gene, encoding carnitine dehydratase, was isolated by oligonucleotide screening from a genomic library of E. coli O44 K74. The caiB gene is 1,215 bp long, and it encodes a protein of 405 amino acids with a predicted M(r) of 45,074. The identity of the gene product was first assessed by its comigration in sodium dodecyl sulfate-polyacrylamide gels with the purified enzyme after overexpression in the pT7 system and by its enzymatic activity. Moreover, the N-terminal amino acid sequence of the purified protein was found to be identical to that predicted from the gene sequence. Northern (RNA) analysis showed that caiB is likely to be cotranscribed with at least one other gene. This other gene could be the gene encoding a 47-kDa protein, which was overexpressed upstream of caiB.     

Automated metabolic reconstruction for Methanococcus jannaschii

We present the computational prediction and synthesis of the metabolic pathways in Methanococcus jannaschii from its genomic sequence using the PathoLogic software. Metabolic reconstruction is based on a reference knowledge base of metabolic pathways and is performed with minimal manual intervention. We predict the existence of 609 metabolic reactions that are assembled in 113 metabolic pathways and an additional 17 super-pathways consisting of one or more component pathways. These assignments represent significantly improved enzyme and pathway predictions compared with previous metabolic reconstructions, and some key metabolic reactions, previously missing, have been identified. Our results, in the form of enzymatic assignments and metabolic pathway predictions, form a database (MJCyc) that is accessible over the World Wide Web for further dissemination among members of the scientific community.     

Caveat emptor: limitations of the automated reconstruction of metabolic pathways in Plasmodium

The functional reconstruction of metabolic pathways from an annotated genome is a tedious and demanding enterprise. Automation of this endeavor using bioinformatics algorithms could cope with the ever-increasing number of sequenced genomes and accelerate the process. Here, the manual reconstruction of metabolic pathways in the functional genomic database of Plasmodium falciparum--Malaria Parasite Metabolic Pathways--is described and compared with pathways generated automatically as they appear in PlasmoCyc, metaSHARK and the Kyoto Encyclopedia for Genes and Genomes. A critical evaluation of this comparison discloses that the automatic reconstruction of pathways generates manifold paths that need an expert manual verification to accept some and reject most others based on manually curated gene annotation.     

Covalent enzyme-substrate intermediates in transferase reactions

Enzymatic transfer reactions are often depicted as occurring on the surface of the enzyme by direct transfer of a chemical group from a donor substrate molecule to an acceptor, in a single-displacement reaction without covalent participation of the enzyme. This picture of enzyme action is purely speculative, because no positive evidence in support of it exists. The kinetic argument, which formerly afforded the sole support for the single-displacement mechanism, is now known to be in error. Contrasting sharply with this dearth of proof is the substantial and growing body of evidence which upholds the double-displacement mechanism of enzymatic transfer. In this mechanism, a catalytic group (or atom) within the active site of the enzyme links covalently with one of the substrates, or some fragment of it, at some stage of the reaction. Through use of the covalent enzyme-substrate intermediate it may be that the enzyme is enabled to overcome certain entropic difficulties, thus accounting in part for the well-known speed of enzymatic reactions. Covalent participation by the enzyme also brings enzymatic catalysis into closer accord with homogenous and heterogeneous catalysis, in which the key reaction is the transient formation of a covalent bond between substrate and catalyst.     

Structure-function relationship of arginyl-tRNA synthetase from Escherichia coli: isolation and characterization of the argS mutation MA5002.

The Escherichia coli K12 argS MA5002 mutant appears to have a functionally altered arginyl-tRNA synthetase (ArgRS). The gene coding for this enzyme was isolated from E. coli genomic DNA using the PCR procedure and inserted into a pUC18 multicopy vector. Sequencing revealed that it differs from the wildtype ArgRS structural gene only by one mutation: a replacement of a C by an A residue which results in substitution of an arginine by a serine at position 134, located two residues downstream from the HVGH consensus sequence. As compared to the genomic enzyme level, this recombinant vector, containing the mutated gene, produces in E. coli JM103, about 100 times as much modified ArgRS. This enzyme was obtained nearly pure after only two chromatographic steps; it exhibits a 4-6 times as low activity and a 5 times as high Km value for ATP as the wildtype enzyme in the aminoacylation and ATP-PPi reactions; Km values for arginine and tRNAArg remained unaltered. The position of this mutation and its effect on enzymatic properties suggest the implication of arginine 134 in ATP binding as well as in the activation catalytic process.     

metaSHARK: software for automated metabolic network prediction from DNA sequence and its application to the genomes of Plasmodium falciparum and Eimeria tenella

The metabolic SearcH And Reconstruction Kit (metaSHARK) is a new fully automated software package for the detection of enzyme-encoding genes within unannotated genome data and their visualization in the context of the surrounding metabolic network. The gene detection package (SHARKhunt) runs on a Linux system and requires only a set of raw DNA sequences (genomic, expressed sequence tag and/or genome survey sequence) as input. Its output may be uploaded to our web-based visualization tool (SHARKview) for exploring and comparing data from different organisms. We first demonstrate the utility of the software by comparing its results for the raw Plasmodium falciparum genome with the manual annotations available at the PlasmoDB and PlasmoCyc websites. We then apply SHARKhunt to the unannotated genome sequences of the coccidian parasite Eimeria tenella and observe that, at an E-value cut-off of 10⁻²⁰, our software makes 142 additional assertions of enzymatic function compared with a recent annotation package working with translated open reading frame sequences. The ability of the software to cope with low levels of sequence coverage is investigated by analyzing assemblies of the E.tenella genome at estimated coverages from 0.5× to 7.5×. Lastly, as an example of how metaSHARK can be used to evaluate the genomic evidence for specific metabolic pathways, we present a study of coenzyme A biosynthesis in P.falciparum and E.tenella.