The evolution and structural anatomy of the small molecule metabolic pathways in Escherichia coli

The 106 small molecule metabolic (SMM) pathways in Escherichia coli are formed by the protein products of 581 genes. We can define 722 domains, nearly all of which are homologous to proteins of known structure, that form all or part of 510 of these proteins. This information allows us to answer general questions on the structural anatomy of the SMM pathway proteins and to trace family relationships and recruitment events within and across pathways. Half the gene products contain a single domain and half are formed by combinations of between two and six domains. The 722 domains belong to one of 213 families that have between one and 51 members. Family members usually conserve their catalytic or cofactor binding properties; substrate recognition is rarely conserved. Of the 213 families, members of only a quarter occur in isolation, i.e. they form single-domain proteins. Most members of the other families combine with domains from just one or two other families and a few more versatile families can combine with several different partners.Excluding isoenzymes, more than twice as many homologues are distributed across pathways as within pathways. However, serial recruitment, with two consecutive enzymes both being recruited to another pathway, is rare and recruitment of three consecutive enzymes is not observed. Only eight of the 106 pathways have a high number of homologues. Homology between consecutive pairs of enzymes with conservation of the main substrate-binding site but change in catalytic mechanism (which would support a simple model of retrograde pathway evolution) occurs only six times in the whole set of enzymes. Most of the domains that form SMM pathways have homologues in non-SMM pathways. Taken together, these results imply a pervasive "mosaic" model for the formation of protein repertoires and pathways.     

DNA-guided assembly of biosynthetic pathways promotes improved catalytic efficiency

Synthetic scaffolds that permit spatial and temporal organization of enzymes in living cells are a promising post-translational strategy for controlling the flow of information in both metabolic and signaling pathways. Here, we describe the use of plasmid DNA as a stable, robust and configurable scaffold for arranging biosynthetic enzymes in the cytoplasm of Escherichia coli. This involved conversion of individual enzymes into custom DNA-binding proteins by genetic fusion to zinc-finger domains that specifically bind unique DNA sequences. When expressed in cells that carried a rationally designed DNA scaffold comprising corresponding zinc finger binding sites, the titers of diverse metabolic products, including resveratrol, 1,2-propanediol and mevalonate were increased as a function of the scaffold architecture. These results highlight the utility of DNA scaffolds for assembling biosynthetic enzymes into functional metabolic structures. Beyond metabolism, we anticipate that DNA scaffolds may be useful in sequestering different types of enzymes for specifying the output of biological signaling pathways or for coordinating other assembly-line processes such as protein folding, degradation and post-translational modifications.     

Yeast sphingolipid metabolism: clues and connections.

This review of sphingolipid metabolism in the budding yeast Saccharomyces cerevisiae contains information on the enzymes and the genes that encode them, as well as connections to other metabolic pathways. Particular attention is given to yeast homologs, domains, and motifs in the sequence, cellular localization of enzymes, and possible protein-protein interactions. Also included are genetic interactions of special interest that provide clues to the cellular biological roles of particular sphingolipid metabolic pathways and specific sphingolipids.     

Evolution of enzymes in metabolism: a network perspective

Several models have been proposed to explain the origin and evolution of enzymes in metabolic pathways. Initially, the retro-evolution model proposed that, as enzymes at the end of pathways depleted their substrates in the primordial soup, there was a pressure for earlier enzymes in pathways to be created, using the later ones as initial template, in order to replenish the pools of depleted metabolites. Later, the recruitment model proposed that initial templates from other pathways could be used as long as those enzymes were similar in chemistry or substrate specificity. These two models have dominated recent studies of enzyme evolution. These studies are constrained by either the small scale of the study or the artificial restrictions imposed by pathway definitions. Here, a network approach is used to study enzyme evolution in fully sequenced genomes, thus removing both constraints. We find that homologous pairs of enzymes are roughly twice as likely to have evolved from enzymes that are less than three steps away from each other in the reaction network than pairs of non-homologous enzymes. These results, together with the conservation of the type of chemical reaction catalyzed by evolutionarily related enzymes, suggest that functional blocks of similar chemistry have evolved within metabolic networks. One possible explanation for these observations is that this local evolution phenomenon is likely to cause less global physiological disruptions in metabolism than evolution of enzymes from other enzymes that are distant from them in the metabolic network.     

Production of bulk chemicals via novel metabolic pathways in microorganisms

Metabolic engineering has been playing important roles in developing high performance microorganisms capable of producing various chemicals and materials from renewable biomass in a sustainable manner. Synthetic and systems biology are also contributing significantly to the creation of novel pathways and the whole cell-wide optimization of metabolic performance, respectively. In order to expand the spectrum of chemicals that can be produced biotechnologically, it is necessary to broaden the metabolic capacities of microorganisms. Expanding the metabolic pathways for biosynthesizing the target chemicals requires not only the enumeration of a series of known enzymes, but also the identification of biochemical gaps whose corresponding enzymes might not actually exist in nature; this issue is the focus of this paper. First, pathway prediction tools, effectively combining reactions that lead to the production of a target chemical, are analyzed in terms of logics representing chemical information, and designing and ranking the proposed metabolic pathways. Then, several approaches for potentially filling in the gaps of the novel metabolic pathway are suggested along with relevant examples, including the use of promiscuous enzymes that flexibly utilize different substrates, design of novel enzymes for non-natural reactions, and exploration of hypothetical proteins. Finally, strain optimization by systems metabolic engineering in the context of novel metabolic pathways constructed is briefly described. It is hoped that this review paper will provide logical ways of efficiently utilizing ‘big’ biological data to design and develop novel metabolic pathways for the production of various bulk chemicals that are currently produced from fossil resources.     

Metabolic pathway reconstruction of eugenol to vanillin bioconversion in Aspergillus niger

Identification of missing genes or proteins participating in the metabolic pathways as enzymes are of great interest. One such class of pathway is involved in the eugenol to vanillin bioconversion. Our goal is to develop an integral approach for identifying the topology of a reference or known pathway in other organism. We successfully identify the missing enzymes and then reconstruct the vanillin biosynthetic pathway in Aspergillus niger. The procedure combines enzyme sequence similarity searched through BLAST homology search and orthologs detection through COG & KEGG databases. Conservation of protein domains and motifs was searched through CDD, PFAM & PROSITE databases. Predictions regarding how proteins act in pathway were validated experimentally and also compared with reported data. The bioconversion of vanillin was screened on UV-TLC plates and later confirmed through GC and GC-MS techniques. We applied a procedure for identifying missing enzymes on the basis of conserved functional motifs and later reconstruct the metabolic pathway in target organism. Using the vanillin biosynthetic pathway of Pseudomonas fluorescens as a case study, we indicate how this approach can be used to reconstruct the reference pathway in A. niger and later results were experimentally validated through chromatography and spectroscopy techniques.     

Understanding specificity in metabolic pathways—Structural biology of human nucleotide metabolism

Interactions are the foundation of life at the molecular level. In the plethora of activities in the cell, the evolution of enzyme specificity requires the balancing of appropriate substrate affinity with a negative selection, in order to minimize interactions with other potential substrates in the cell. To understand the structural basis for enzyme specificity, the comparison of structural and biochemical data between enzymes within pathways using similar substrates and effectors is valuable.Nucleotide metabolism is one of the largest metabolic pathways in the human cell and is of outstanding therapeutic importance since it activates and catabolises nucleoside based anti-proliferative drugs and serves as a direct target for anti-proliferative drugs. In recent years the structural coverage of the enzymes involved in human nucleotide metabolism has been dramatically improved and is approaching completion. An important factor has been the contribution from the Structural Genomics Consortium (SGC) at Karolinska Institutet, which recently has solved 33 novel structures of enzymes and enzyme domains in human nucleotide metabolism pathways and homologs thereof. In this review we will discuss some of the principles for substrate specificity of enzymes in human nucleotide metabolism illustrated by a selected set of enzyme families where a detailed understanding of the structural determinants for specificity is now emerging.     

MRAD: Metabolic reaction analysis database--an entity-relationship approach

The Metabolic Reaction Analysis Database (MRAD) is a relational database based on the Entity-Relationship (ER) model which combines information about organisms, biochemical pathways, reactions, enzymes, substrates, products and genes. It describes 244,596 genes in 79 organisms, 6,552 enzymes, and 3,552 reactions, 3,100 substrates, 2,866 products and 118 metabolic pathways. The MRAD graphical user interface allows for the identification of metabolic reactions which are similar and dissimilar in multiple organisms, reactions in a pathway which are missing in an organism and using any combination between one to six of the biological entities of organisms, genes, pathways, enzymes, substrates and products to determine metabolic reactions. MRAD provides a powerful and efficient tool for the construction of flux balance models for metabolic engineering applications.     

Adaptive evolution of metabolic pathways in Drosophila

The adaptive significance of enzyme variation has been of central interest in population genetics. Yet, how natural selection operates on enzymes in the larger context of biochemical pathways has not been broadly explored. A basic expectation is that natural selection on metabolic phenotypes will target enzymes that control metabolic flux, but how adaptive variation is distributed among enzymes in metabolic networks is poorly understood. Here, we use population genetic methods to identify enzymes responding to adaptive selection in the pathways of central metabolism in Drosophila melanogaster and Drosophila simulans. We report polymorphism and divergence data for 17 genes that encode enzymes of 5 metabolic pathways that converge at glucose-6-phosphate (G6P). Deviations from neutral expectations were observed at five loci. Of the 10 genes that encode the enzymes of glycolysis, only aldolase (Ald) deviated from neutrality. The other 4 genes that were inconsistent with neutral evolution (glucose-6-phosphate dehydrogenase [G6pd]), phosphoglucomutase [Pgm], trehalose-6-phosphate synthetase [Tps1], and glucose-6phosphatase [G6pase] encode G6P branch point enzymes that catalyze reactions at the entry point to the pentose-phosphate, glycogenic, trehalose synthesis, and gluconeogenic pathways. We reconcile these results with population genetics theory and existing arguments on metabolic regulation and propose that the incidence of adaptive selection in this system is related to the distribution of flux control. The data suggest that adaptive evolution of G6P branch point enzymes may have special significance in metabolic adaptation.     

The University of Minnesota Biocatalysis/Biodegradation Database: emphasizing enzymes

The University of Minnesota Biocatalysis/Biodegradation Database (UM-BBD, http://umbbd.ahc.umn.edu/) provides curated information on microbial catabolic enzymes and their organization into metabolic pathways. Currently, it contains information on over 400 enzymes. In the last year the enzyme page was enhanced to contain more internal and external links; it also displays the different metabolic pathways in which each enzyme participates. In collaboration with the Nomenclature Commission of the International Union of Biochemistry and Molecular Biology, 35 UM-BBD enzymes were assigned complete EC codes during 2000. Bacterial oxygenases are heavily represented in the UM-BBD; they are known to have broad substrate specificity. A compilation of known reactions of naphthalene and toluene dioxygenases were recently added to the UM-BBD; 73 and 108 were listed respectively. In 2000 the UM-BBD is mirrored by two prestigious groups: the European Bioinformatics Institute and KEGG (the Kyoto Encyclopedia of Genes and Genomes). Collaborations with other groups are being developed. The increased emphasis on UM-BBD enzymes is important for predicting novel metabolic pathways that might exist in nature or could be engineered. It also is important for current efforts in microbial genome annotation.     

Metabolons involving plant cytochrome P450s

Arranging biological processes into “compartments” is a key feature of all eukaryotic cells. Through this mechanism, cells can drastically increase metabolic efficiency and manage complex cellular processes more efficiently, saving space and energy. Compartmentation at the molecular level is mediated by metabolons. A metabolon is an ordered protein complex of sequential metabolic enzymes and associated cellular structural elements. The sub-cellular organization of enzymes involved in the synthesis and storage of plant natural products appears to involve the anchoring of metabolons by cytochrome P450 monooxygenases (P450s) to specific domains of the endoplasmic reticulum (ER) membrane. This review focuses on the current evidence supporting the organization of metabolons around P450s on the surface of the ER. We␣outline direct and indirect experimental data that describes P450 enzymes in the phenylpropanoid, flavonoid, cyanogenic glucoside, and other biosynthetic pathways. We also discuss the limitations and future directions of metabolon research and the potential for application to metabolic engineering endeavors.     

The ATP-grasp enzymes

The ATP-grasp enzymes consist of a superfamily of 21 proteins that contain an atypical ATP-binding site, called the ATP-grasp fold. The ATP-grasp fold is comprised of two α+β domains that "grasp" a molecule of ATP between them and members of the family typically have an overall structural design containing three common conserved focal domains. The founding members of the family consist of biotin carboxylase, d-ala-d-ala ligase and glutathione synthetase, all of which catalyze the ATP-assisted reaction of a carboxylic acid with a nucleophile via the formation of an acylphosphate intermediate. While most members of the superfamily follow this mechanistic pathway, studies have demonstrated that two enzymes catalyze only the phosphoryl transfer step and thus are kinases instead of ligases. Members of the ATP-grasp superfamily are found in several metabolic pathways including de novo purine biosynthesis, gluconeogenesis, and fatty acid synthesis. Given the critical nature of these enzymes, researchers have actively sought the development of potent inhibitors of several members of the superfamily as antibacterial and anti-obseity agents. In this review, we will discuss the structure, function, mechanism, and inhibition of the ATP-grasp enzymes.     

Interaction interfaces of protein domains are not topologically equivalent across families within superfamilies: Implications for metabolic and signaling pathways

Using a data set of aligned protein domain superfamilies of known three-dimensional structure, we compared the location of interdomain interfaces on the tertiary folds between members of distantly related protein domain superfamilies. The data set analyzed is comprised of interdomain interfaces, with domains occurring within a polypeptide chain and those between two polypeptide chains. We observe that, in general, the interfaces between protein domains are formed entirely in different locations on the tertiary folds in such pairs. This variation in the location of interface happens in protein domains involved in a wide range of functions, such as enzymes, adapters, and domains that bind protein ligands, or cofactors. While basic biochemical functionality is preserved at the domain superfamily level, the effect of biochemical function on protein assemblies is different in these protein domains related by superfamily. The divergence between proteins, in most cases, is coupled with domain recruitment, with different modes of interaction with the recruited domain. This is in complete contrast to the observation that in closely related homologous protein domains, almost always the interaction interfaces are topologically equivalent. In a small subset of interacting domains within proteins related by remote homology, we observe that the relative positioning of domains with respect to one another is preserved. Based on the analysis of multidomain proteins of known or unknown structure, we suggest that variation in protein-protein interactions in members within a superfamily could serve as diverging points in otherwise parallel metabolic or signaling pathways. We discuss a few representative cases of diverging pathways involving domains in a superfamily.     

Insight into trichomonas vaginalis genome evolution through metabolic pathways comparison

Trichomonas vaginalis causes the trichomoniasis, in women and urethritis and prostate cancer in men. Its genome draft published by TIGR in 2007 presents many unusual genomic and biochemical features like, exceptionally large genome size, the presence of hydrogenosome, gene duplication, lateral gene transfer mechanism and the presence of miRNA. To understand some of genomic features we have performed a comparative analysis of metabolic pathways of the T. vaginalis with other 22 significant common organisms. Enzymes from the biochemical pathways of T. vaginalis and other selected organisms were retrieved from the KEGG metabolic pathway database. The metabolic pathways of T. vaginalis common in other selected organisms were identified. Total 101 enzymes present in different metabolic pathways of T. vaginalis were found to be orthologous by using BLASTP program against the selected organisms. Except two enzymes all identified orthologous enzymes were also identified as paralogous enzymes. Seventy-five of identified enzymes were also identified as essential for the survival of T. vaginalis, while 26 as non-essential. The identified essential enzymes also represent as good candidate for novel drug targets. Interestingly, some of the identified orthologous and paralogous enzymes were found playing significant role in the key metabolic activities while others were found playing active role in the process of pathogenesis. The N-acetylneuraminate lyase was analyzed as the candidate of lateral genes transfer. These findings clearly suggest the active participation of lateral gene transfer and gene duplication during evolution of T. vaginalis from the enteric to the pathogenic urogenital environment.     

The protein acetylome and the regulation of metabolism

Acetyl-coenzyme A (CoA) is a central metabolite involved in numerous anabolic and catabolic pathways, as well as in protein acetylation. Beyond histones, a large number of metabolic enzymes are acetylated in both animal and bacteria, and the protein acetylome is now emerging in plants. Protein acetylation is influenced by the cellular level of both acetyl-CoA and NAD(+), and regulates the activity of several enzymes. Acetyl-CoA is thus ideally placed to act as a key molecule linking the energy balance of the cell to the regulation of gene expression and metabolic pathways via the control of protein acetylation. Better knowledge over how to influence acetyl-CoA levels and the acetylation process promises to be an invaluable tool to control metabolic pathways.     

Genome-Wide Analysis of Biotin Biosynthesis in Eukaryotic Photosynthetic Algae

Biotin is a cofactor responsible for carbon dioxide transfer in several carboxylase enzymes, which play a significant role in various metabolic reactions such as fatty acid synthesis, branched chain amino acid catabolism, and gluconeogenesis. Biotin is also involved in citric acid cycle, which is the process of biochemical energy generation during aerobic respiration. Though the function of biotin in the growth of algae has been extensively investigated, little is known about the biosynthetic routes of biotin in the algal kingdom. In the present study, 44 biotin biosynthesis-related genes were identified from 14 eukaryotic photosynthetic algal genomes by BLASTP and TBLASN programs. A comprehensive analysis was performed to characterize distribution, phylogeny, structure domains, and coevolution patterns of those genes. Forty-four biotin biosynthesis-related enzymes (BBREs) were found to be distributed in three groups: 7-keto-8-aminopelargonic acid synthase, diaminopelargonic acid synthase/dethiobiotin synthetase, and biotin synthase. Structure domains were considerably conserved among the subfamilies of BBREs. The intramolecular coevolutionary sites are widely distributed in biotin synthase. The present study provides new insights into the origin and evolution of biotin biosynthetic pathways in eukaryotic photosynthetic algae. Furthermore, the characterization of biotin biosynthesis-related genes from algae will promote the identification and functional studies of BBREs.     

Discovery of novel pathways for carbohydrate metabolism

Closing the gap between the increasing availability of complete genome sequences and the discovery of novel enzymes in novel metabolic pathways is a significant challenge. Here, we review recent examples of assignment of in vitro enzymatic activities and in vivo metabolic functions to uncharacterized proteins, with a focus on enzymes and metabolic pathways involved in the catabolism and biosynthesis of monosaccharides and polysaccharides. The most effective approaches are based on analyses of sequence-function space in protein families that provide clues for the predictions of the functions of the uncharacterized enzymes. As summarized in this Opinion, this approach allows the discovery of the catabolism of new molecules, new pathways for common molecules, and new enzymatic chemistries.