Finding sequences for over 270 orphan enzymes

The transferred translocases: An old wine in a new bottle

The role of translocases was underappreciated and was not included as a separate class in the enzyme commission until August 2018. The recent research interests in proteomics of orphan enzymes, ionomics, and metallomics along with high-throughput sequencing technologies generated overwhelming data and revamped this enzyme into a separate class. This offers a great opportunity to understand the role of new or orphan enzymes in general and specifically translocases. The enzymes belonging to translocases regulate/permeate the transfer of ions or molecules across the membranes. These enzyme entries were previously associated with other enzyme classes, which are now transferred to a new enzyme class 7 (EC 7). The entries that are reclassified are important to extend the enzyme list, and it is the need of the hour. Accordingly, there is an upgradation of entries of this class of enzymes in several databases. This review is a concise compilation of translocases with reference to the number of entries currently available in the databases. This review also focuses on function as well as dysfunction of translocases during normal and disordered states, respectively.

A computational workflow for the expansion of heterologous biosynthetic pathways to natural product derivatives

Plant natural products (PNPs) and their derivatives are important but underexplored sources of pharmaceutical molecules. To access this untapped potential, the reconstitution of heterologous PNP biosynthesis pathways in engineered microbes provides a valuable starting point to explore and produce novel PNP derivatives. Here, we introduce a computational workflow to systematically screen the biochemical vicinity of a biosynthetic pathway for pharmaceutical compounds that could be produced by derivatizing pathway intermediates. We apply our workflow to the biosynthetic pathway of noscapine, a benzylisoquinoline alkaloid (BIA) with a long history of medicinal use. Our workflow identifies pathways and enzyme candidates for the production of (S)-tetrahydropalmatine, a known analgesic and anxiolytic, and three additional derivatives. We then construct pathways for these compounds in yeast, resulting in platforms for de novo biosynthesis of BIA derivatives and demonstrating the value of cheminformatic tools to predict reactions, pathways, and enzymes in synthetic biology and metabolic engineering.

Implementation of homology based and non-homology based computational methods for the identification and annotation of orphan enzymes: using Mycobacterium tuberculosis H37Rv as a case study

Background

Homology based methods are one of the most important and widely used approaches for functional annotation of high-throughput microbial genome data. A major limitation of these methods is the absence of well-characterized sequences for certain functions. The non-homology methods based on the context and the interactions of a protein are very useful for identifying missing metabolic activities and functional annotation in the absence of significant sequence similarity. In the current work, we employ both homology and context-based methods, incrementally, to identify local holes and chokepoints, whose presence in the Mycobacterium tuberculosis genome is indicated based on its interaction with known proteins in a metabolic network context, but have not been annotated. We have developed two computational procedures using network theory to identify orphan enzymes (‘Hole finding protocol’) coupled with the identification of candidate proteins for the predicted orphan enzyme (‘Hole filling protocol’). We propose an integrated interaction score based on scores from the STRING database to identify candidate protein sequences for the orphan enzymes from M. tuberculosis, as a case study, which are most likely to perform the missing function.

Results

The application of an automated homology-based enzyme identification protocol, ModEnzA, on M. tuberculosis genome yielded 56 novel enzyme predictions. We further predicted 74 putative local holes, 6 choke points, and 3 high confidence local holes in the genome using ‘Hole finding protocol’. The ‘Hole-filling protocol’ was validated on the E. coli genome using artificial in-silico enzyme knockouts where our method showed 25% increased accuracy, compared to other methods, in assigning the correct sequence for the knocked-out enzyme amongst the top 10 ranks. The method was further validated on 8 additional genomes.

Conclusions

We have developed methods that can be generalized to augment homology-based annotation to identify missing enzyme coding genes and to predict a candidate protein for them. For pathogens such as M. tuberculosis, this work holds significance in terms of increasing the protein repertoire and thereby, the potential for identifying novel drug targets.

Enzyme annotation for orphan and novel reactions using knowledge of substrate reactive sites

Recent advances in synthetic biochemistry have resulted in a wealth of novel hypothetical enzymatic reactions that are not matched to protein-encoding genes, deeming them “orphan.” A large number of known metabolic enzymes are also orphan, leaving important gaps in metabolic network maps. Proposing genes for the catalysis of orphan reactions is critical for applications ranging from biotechnology to medicine. In this work, the computational method BridgIT identified potential enzymes of orphan reactions and nearly all theoretically possible biochemical transformations, providing candidate genes to catalyze these reactions to the research community. The BridgIT online tool will allow researchers to fill the knowledge gaps in metabolic networks and will act as a starting point for designing novel enzymes to catalyze nonnatural transformations.

Thousands of biochemical reactions with characterized activities are “orphan,” meaning they cannot be assigned to a specific enzyme, leaving gaps in metabolic pathways. Novel reactions predicted by pathway-generation tools also lack associated sequences, limiting protein engineering applications. Associating orphan and novel reactions with known biochemistry and suggesting enzymes to catalyze them is a daunting problem. We propose the method BridgIT to identify candidate genes and catalyzing proteins for these reactions. This method introduces information about the enzyme binding pocket into reaction-similarity comparisons. BridgIT assesses the similarity of two reactions, one orphan and one well-characterized nonorphan reaction, using their substrate reactive sites, their surrounding structures, and the structures of the generated products to suggest enzymes that catalyze the most-similar nonorphan reactions as candidates for also catalyzing the orphan ones. We performed two large-scale validation studies to test BridgIT predictions against experimental biochemical evidence. For the 234 orphan reactions from the Kyoto Encyclopedia of Genes and Genomes (KEGG) 2011 (a comprehensive enzymatic-reaction database) that became nonorphan in KEGG 2018, BridgIT predicted the exact or a highly related enzyme for 211 of them. Moreover, for 334 of 379 novel reactions in 2014 that were later cataloged in KEGG 2018, BridgIT predicted the exact or highly similar enzymes. BridgIT requires knowledge about only four connecting bonds around the atoms of the reactive sites to correctly annotate proteins for 93% of analyzed enzymatic reactions. Increasing to seven connecting bonds allowed for the accurate identification of a sequence for nearly all known enzymatic reactions.

Functional genetic discovery of enzymes using full-scan mass spectrometry metabolomics

Our understanding of metabolic networks is incomplete, and new enzymatic activities await discovery in well-studied organisms. Mass spectrometric measurement of cellular metabolites reveals compounds inside cells that are unexplained by current maps of metabolic reactions, and existing computational models are unable to account for all activities observed within cells. Additional large-scale genetic and biochemical approaches are required to elucidate metabolic gene function. We have used full-scan mass spectrometry metabolomics of polar small molecules to examine deletion mutants of candidate enzymes in the model yeast Saccharomyces cerevisiae. We report the identification of 25 genes whose deletion results in focal metabolic changes consistent with loss of enzymatic activity and describe the informatic approaches used to enrich for candidate enzymes from uncharacterized open reading frames. Triumphs and pitfalls of metabolic phenotyping screens are discussed, including estimates of the frequency of uncharacterized eukaryotic genes that affect metabolism and key issues to consider when searching for new enzymatic functions in other organisms.

Global untargeted serum metabolomic analyses nominate metabolic pathways responsive to loss of expression of the orphan metallo β-lactamase, MBLAC1

The C. elegans gene swip-10 encodes an orphan metallo β-lactamase that genetic studies indicate is vital for limiting neuronal excitability and viability. Sequence analysis indicates that the mammalian gene Mblac1 is the likely ortholog of swip-10, with greatest sequence identity localized to the encoded protein's single metallo β-lactamase domain. The substrate for the SWIP-10 protein remains unknown and to date no functional roles have been ascribed to MBLAC1, though we have shown that the protein binds the neuroprotective β-lactam antibiotic, ceftriaxone. To gain insight into the functional role of MBLAC1 in vivo, we used CRISPR/Cas9 methods to disrupt N-terminal coding sequences of the mouse Mblac1 gene, resulting in a complete loss of protein expression in viable, homozygous knockout (KO) animals. Using serum from both WT and KO mice, we performed global, untargeted metabolomic analyses, resolving small molecules via hydrophilic interaction chromatography (HILIC) based ultra-performance liquid chromatography, coupled to mass spectrometry (UPLC-MS/MS). Unsupervised principal component analysis reliably segregated the metabolomes of MBLAC1 KO and WT mice, with 92 features subsequently nominated as significantly different by ANOVA, and for which we made tentative and putative metabolite assignments. Bioinformatic analyses of these molecules nominate validated pathways subserving bile acid biosynthesis and linoleate metabolism, networks known to be responsive to metabolic and oxidative stress. Our findings lead to hypotheses that can guide future targeted studies seeking to identify the substrate for MBLAC1 and how substrate hydrolysis supports the neuroprotective actions of ceftriaxone.

Exploring biochemical pathways for mono-ethylene glycol (MEG) synthesis from synthesis gas

Mono-ethylene glycol (MEG) is an important petrochemical with widespread use in numerous consumer products. The current industrial MEG-production process relies on non-renewable fossil fuel-based feedstocks, such as petroleum, natural gas, and naphtha; hence, it is useful to explore alternative routes of MEG-synthesis from gases as they might provide a greener and more sustainable alternative to the current production methods. Technologies of synthetic biology and metabolic engineering of microorganisms can be deployed for the expression of new biochemical pathways for MEG-synthesis from gases, provided that such promising alternative routes are first identified. We used the BNICE.ch algorithm to develop novel and previously unknown biological pathways to MEG from synthesis gas by leveraging the Wood-Ljungdahl pathway of carbon fixation of acetogenic bacteria. We developed a set of useful pathway pruning and analysis criteria to systematically assess thousands of pathways generated by BNICE.ch. Published genome-scale models of Moorella thermoacetica and Clostridium ljungdahlii were used to perform the pathway yield calculations and in-depth analyses of seven (7) newly developed biological MEG-producing pathways from gases, including CO₂, CO, and H₂. These analyses helped identify not only better candidate pathways, but also superior chassis organisms that can be used for metabolic engineering of the candidate pathways. The pathway generation, pruning, and detailed analysis procedures described in this study can also be used to develop biochemical pathways for other commodity chemicals from gaseous substrates.

Biochemistry and physiology within the framework of the extended synthesis of evolutionary biology

Functional biologists, like Claude Bernard, ask "How?", meaning that they investigate the mechanisms underlying the emergence of biological functions (proximal causes), while evolutionary biologists, like Charles Darwin, asks "Why?", meaning that they search the causes of adaptation, survival and evolution (remote causes). Are these divergent views on what is life? The epistemological role of functional biology (molecular biology, but also biochemistry, physiology, cell biology and so forth) appears essential, for its capacity to identify several mechanisms of natural selection of new characters, individuals and populations. Nevertheless, several issues remain unsolved, such as orphan metabolic activities, i.e., adaptive functions still missing the identification of the underlying genes and proteins, and orphan genes, i.e., genes that bear no signature of evolutionary history, yet provide an organism with improved adaptation to environmental changes. In the framework of the Extended Synthesis, we suggest that the adaptive roles of any known function/structure are reappraised in terms of their capacity to warrant constancy of the internal environment (homeostasis), a concept that encompasses both proximal and remote causes.

Characterising Complex Enzyme Reaction Data

The relationship between enzyme-catalysed reactions and the Enzyme Commission (EC) number, the widely accepted classification scheme used to characterise enzyme activity, is complex and with the rapid increase in our knowledge of the reactions catalysed by enzymes needs revisiting. We present a manual and computational analysis to investigate this complexity and found that almost one-third of all known EC numbers are linked to more than one reaction in the secondary reaction databases (e.g., KEGG). Although this complexity is often resolved by defining generic, alternative and partial reactions, we have also found individual EC numbers with more than one reaction catalysing different types of bond changes. This analysis adds a new dimension to our understanding of enzyme function and might be useful for the accurate annotation of the function of enzymes and to study the changes in enzyme function during evolution.

Mesaconase Activity of Class I Fumarase Contributes to Mesaconate Utilization by Burkholderia xenovorans

Pseudomonas aeruginosa, Yersinia pestis, and many other bacteria are able to utilize the C5-dicarboxylic acid itaconate (methylenesuccinate). Itaconate degradation starts with its activation to itaconyl coenzyme A (itaconyl-CoA), which is further hydrated to (S)-citramalyl-CoA, and citramalyl-CoA is finally cleaved into acetyl-CoA and pyruvate. The xenobiotic-degrading betaproteobacterium Burkholderia xenovorans possesses a P. aeruginosa-like itaconate degradation gene cluster and is able to grow on itaconate and its isomer mesaconate (methylfumarate). Although itaconate degradation proceeds in B. xenovorans in the same way as in P. aeruginosa, the pathway of mesaconate utilization is not known. Here, we show that mesaconate is metabolized through its hydration to (S)-citramalate. The latter compound is then metabolized to acetyl-CoA and pyruvate with the participation of two enzymes of the itaconate degradation pathway, a promiscuous itaconate-CoA transferase able to activate (S)-citramalate in addition to itaconate and (S)-citramalyl-CoA lyase. The first reaction of the pathway, the mesaconate hydratase (mesaconase) reaction, is catalyzed by a class I fumarase. As this enzyme (Bxe_A3136) has similar efficiencies (kcat/Km) for both fumarate and mesaconate hydration, we conclude that B. xenovorans class I fumarase is in fact a promiscuous fumarase/mesaconase. This promiscuity is physiologically relevant, as it allows the growth of this bacterium on mesaconate as a sole carbon and energy source.

UniProt: a hub for protein information

UniProt is an important collection of protein sequences and their annotations, which has doubled in size to 80 million sequences during the past year. This growth in sequences has prompted an extension of UniProt accession number space from 6 to 10 characters. An increasing fraction of new sequences are identical to a sequence that already exists in the database with the majority of sequences coming from genome sequencing projects. We have created a new proteome identifier that uniquely identifies a particular assembly of a species and strain or subspecies to help users track the provenance of sequences. We present a new website that has been designed using a user-experience design process. We have introduced an annotation score for all entries in UniProt to represent the relative amount of knowledge known about each protein. These scores will be helpful in identifying which proteins are the best characterized and most informative for comparative analysis. All UniProt data is provided freely and is available on the web at http://www.uniprot.org/.

Updates in Rhea--a manually curated resource of biochemical reactions

Rhea (http://www.ebi.ac.uk/rhea) is a comprehensive and non-redundant resource of expert-curated biochemical reactions described using species from the ChEBI (Chemical Entities of Biological Interest) ontology of small molecules. Rhea has been designed for the functional annotation of enzymes and the description of genome-scale metabolic networks, providing stoichiometrically balanced enzyme-catalyzed reactions (covering the IUBMB Enzyme Nomenclature list and additional reactions), transport reactions and spontaneously occurring reactions. Rhea reactions are extensively curated with links to source literature and are mapped to other publicly available enzyme and pathway databases such as Reactome, BioCyc, KEGG and UniPathway, through manual curation and computational methods. Here we describe developments in Rhea since our last report in the 2012 database issue of Nucleic Acids Research. These include significant growth in the number of Rhea reactions and the inclusion of reactions involving complex macromolecules such as proteins, nucleic acids and other polymers that lie outside the scope of ChEBI. Together these developments will significantly increase the utility of Rhea as a tool for the description, analysis and reconciliation of genome-scale metabolic models.

Profiling the orphan enzymes

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.