On the origin of biochemistry at an alkaline hydrothermal vent

A model for the origin of biochemistry at an alkaline hydrothermal vent has been developed that focuses on the acetyl-CoA (Wood-Ljungdahl) pathway of CO2 fixation and central intermediary metabolism leading to the synthesis of the constituents of purines and pyrimidines. The idea that acetogenesis and methanogenesis were the ancestral forms of energy metabolism among the first free-living eubacteria and archaebacteria, respectively, stands in the foreground. The synthesis of formyl pterins, which are essential intermediates of the Wood-Ljungdahl pathway and purine biosynthesis, is found to confront early metabolic systems with steep bioenergetic demands that would appear to link some, but not all, steps of CO2 reduction to geochemical processes in or on the Earth's crust. Inorganically catalysed prebiotic analogues of the core biochemical reactions involved in pterin-dependent methyl synthesis of the modern acetyl-CoA pathway are considered. The following compounds appear as probable candidates for central involvement in prebiotic chemistry: metal sulphides, formate, carbon monoxide, methyl sulphide, acetate, formyl phosphate, carboxy phosphate, carbamate, carbamoyl phosphate, acetyl thioesters, acetyl phosphate, possibly carbonyl sulphide and eventually pterins. Carbon might have entered early metabolism via reactions hardly different from those in the modern Wood-Ljungdahl pathway, the pyruvate synthase reaction and the incomplete reverse citric acid cycle. The key energy-rich intermediates were perhaps acetyl thioesters, with acetyl phosphate possibly serving as the universal metabolic energy currency prior to the origin of genes. Nitrogen might have entered metabolism as geochemical NH3 via two routes: the synthesis of carbamoyl phosphate and reductive transaminations of alpha-keto acids. Together with intermediates of methyl synthesis, these two routes of nitrogen assimilation would directly supply all intermediates of modern purine and pyrimidine biosynthesis. Thermodynamic considerations related to formyl pterin synthesis suggest that the ability to harness a naturally pre-existing proton gradient at the vent-ocean interface via an ATPase is older than the ability to generate a proton gradient with chemistry that is specified by genes.     

Recreating ancient metabolic pathways before enzymes

The biochemistry of all living organisms uses complex, enzyme-catalyzed metabolic reaction networks. Yet, at life’s origins, enzymes had not yet evolved. Therefore, it has been postulated that non-enzymatic metabolic pathways predated their enzymatic counterparts. In this account article, we describe our recent work to evaluate whether two ancient carbon fixation pathways, the rTCA (reductive tricarboxylic acid) cycle and the reductive AcCoA (Wood-Ljungdahl) pathway, could have operated without enzymes and therefore have originated as prebiotic chemistry. We also describe the discovery of an Fe²⁺-promoted complex reaction network that may represent a prebiotic predecessor to the TCA and glyoxylate cycles. The collective results support the idea that most central metabolic pathways could have roots in prebiotic chemistry.     

Ordering events of biochemical evolution

Metabolic pathways exhibit structures resulting from an evolutionary process. Pathways have been inherited through time with modification, from the earliest periods of life. It is possible to compare the structure of pathways as done in comparative anatomy, i.e. for inferring ancestral pathways or parts of it (ancestral enzymatic functions), using standard phylogenetic reconstruction. Thus a phylogenetic tree of pathways provides a relative ordering of the rise of enzymatic functions. It even becomes possible to order the birth of each complete pathway in time. This particular "DNA-free" conceptual approach to evolutionary biochemistry is reviewed, gathering all the justifications given for it. Then, the method of assigning a given pathway to a time span of biochemical development is revisited. The previous method used an implicit "clock" of metabolic development that is difficult to justify. We develop a new clock-free approach, using functional biochemical arguments. Results of the two methods are not significantly different; our method is just more precise. This suggests that the clock assumed in the first method does not provoke any important artefact in describing the development of biochemical evolution. It is just unnecessary to postulate it. As a result, most of the amino acid metabolic pathways develop forwards, confirming former models of amino acid catabolism evolution, but not those for amino acid anabolism. The order of appearance of sectors of universal cellular metabolism is: (1) amino acid catabolism, (2) amino acid anabolism and closure of the urea cycle, (3) glycolysis and glycogenesis, (4) closure of the pentose-phosphate cycle, (5) closure of the Krebs cycle and fatty acids metabolism, (6) closure of the Calvin cycle.     

Phylogenetic analysis based on genome-scale metabolic pathway reaction content

Phylogenetic classifications based on single genes such as rRNA genes do not provide a complete and accurate picture of evolution because they do not account for evolutionary leaps caused by gene transfer, duplication, deletion and functional replacement. Here, we present a whole-genome-scale phylogeny based on metabolic pathway reaction content. From the genome sequences of 42 microorganisms, we deduced the metabolic pathway reactions and used the relatedness of these contents to construct a phylogenetic tree that represents the similarity of metabolic profiles (relatedness) as well as the extent of metabolic pathway similarity (evolutionary distance). This method accounts for horizontal gene transfer and specific gene loss by comparison of whole metabolic subpathways, and allows evaluation of evolutionary relatedness and changes in metabolic pathways. Thus, a tree based on metabolic pathway content represents both the evolutionary time scale (changes in genetic content) and the evolutionary process (changes in metabolism).     

The evolution of photosynthesis…again?

‘Replaying the tape’ is an intriguing ‘would it happen again?’ exercise. With respect to broad evolutionary innovations, such as photosynthesis, the answers are central to our search for life elsewhere. Photosynthesis permits a large planetary biomass on Earth. Specifically, oxygenic photosynthesis has allowed an oxygenated atmosphere and the evolution of large metabolically demanding creatures, including ourselves. There are at least six prerequisites for the evolution of biological carbon fixation: a carbon-based life form; the presence of inorganic carbon; the availability of reductants; the presence of light; a light-harvesting mechanism to convert the light energy into chemical energy; and carboxylating enzymes. All were present on the early Earth. To provide the evolutionary pressure, organic carbon must be a scarce resource in contrast to inorganic carbon. The probability of evolving a carboxylase is approached by creating an inventory of carbon-fixation enzymes and comparing them, leading to the conclusion that carbon fixation in general is basic to life and has arisen multiple times. Certainly, the evolutionary pressure to evolve new pathways for carbon fixation would have been present early in evolution. From knowledge about planetary systems and extraterrestrial chemistry, if organic carbon-based life occurs elsewhere, photosynthesis—although perhaps not oxygenic photosynthesis—would also have evolved.     

New prospects for deducing the evolutionary history of metabolic pathways in prokaryotes: Aromatic biosynthesis as a case-in-point

Metabolic pathways of prokaryotes are more biochemically diverse than is generally recognized. Distinctive biochemical features are shared by phylogenetic clusters. The hierarchical levels of characterstate clustering depends upon evolutionary events which fortuitously became fixed in the genome of a common ancestor. Prokaryotes can now be ordered on a phylogenetic tree. This allows the evolutionary steps that underlie the construction and regulation of appropriately complex biochemical pathways to be traced in an evolutionary progression of prokaryote types that house these pathways. Essentially the approach is to deduce ancestral character states at ever deeper phylogenetic levels, utilizing logical principles of maximum parsimony. The current perspective on the evolution of the biochemical pathway for biosynthesis of aromatic amino acids is developed as a case-in-point model for analyses that should be feasible with many major metabolic systems. Phenylalanine biosynthesis probably arose prior to the addition of branches leading to tyrosine and tryptophan. An evolutionary scenario is developed that begins with non-enzymatic reactions which may have operated in primitive systems, followed by the evolution of an enzymatic system that pre-dated the divergence of major lineages of modern eubacteria (Gram-positive bacteria, Gram-negative purple bacteria, and cyanobacteria).Florida Agricultural Experiment Station, Journal Series No. 8251.     

Genomic analysis of carbon monoxide utilization and butanol production by Clostridium carboxidivorans strain P7

Increasing demand for the production of renewable fuels has recently generated a particular interest in microbial production of butanol. Anaerobic bacteria, such as Clostridium spp., can naturally convert carbohydrates into a variety of primary products, including alcohols like butanol. The genetics of microorganisms like Clostridium acetobutylicum have been well studied and their solvent-producing metabolic pathways characterized. In contrast, less is known about the genetics of Clostridium spp. capable of converting syngas or its individual components into solvents. In this study, the type of strain of a new solventogenic Clostridium species, C. carboxidivorans, was genetically characterized by genome sequencing. C. carboxidivorans strain P7(T) possessed a complete Wood-Ljungdahl pathway gene cluster, involving CO and CO(2) fixation and conversion to acetyl-CoA. Moreover, with the exception of an acetone production pathway, all the genetic determinants of canonical ABE metabolic pathways for acetate, butyrate, ethanol and butanol production were present in the P7(T) chromosome. The functionality of these pathways was also confirmed by growth of P7(T) on CO and production of CO(2) as well as volatile fatty acids (acetate and butyrate) and solvents (ethanol and butanol). P7(T) was also found to harbour a 19 Kbp plasmid, which did not include essential or butanol production related genes. This study has generated in depth knowledge of the P7(T) genome, which will be helpful in developing metabolic engineering strategies to improve C. carboxidivorans's natural capacity to produce potential biofuels from syngas.     

Probabilistic reconstruction of ancestral protein sequences

Using a maximum-likelihood formalism, we have developed a method with which to reconstruct the sequences of ancestral proteins. Our approach allows the calculation of not only the most probable ancestral sequence but also of the probability of any amino acid at any given node in the evolutionary tree. Because we consider evolution on the amino acid level, we are better able to include effects of evolutionary pressure and take advantage of structural information about the protein through the use of mutation matrices that depend on secondary structure and surface accessibility. The computational complexity of this method scales linearly with the number of homologous proteins used to reconstruct the ancestral sequence.     

Computation of metabolic fluxes and efficiencies for biological carbon dioxide fixation

With rising energy prices and concern over the environmental impact of fossil fuel consumption, the push to develop biomass derived fuels has increased significantly. Although most global carbon fixation occurs via the Calvin Benson Bassham cycle, there are currently five other known pathways for carbon fixation; the goal of this study was to determine the thermodynamic efficiencies of all six carbon fixation pathways for the production of biomass using flux balance analysis. The three chemotrophic pathways, the reductive acetyl-CoA pathway, the 3-hydroxypropionate/4-hydroxybutyrate cycle and the dicarboxylate/4-hydroxybutyrate cycle, were found to be more efficient than photoautotrophic carbon fixation pathways. However, as hydrogen is not freely available, the energetic cost of hydrogen production from sunlight was calculated and included in the overall energy demand, which results in a 5 fold increase in the energy demand of chemoautotrophic carbon fixation. Therefore, when the cost of hydrogen production is included, photoautotrophic pathways are more efficient. However, the energetic cost for the production of 12 metabolic precursors was found to vary widely across the different carbon fixation pathways; therefore, different pathways may be more efficient at producing products from a single precursor than others. The results of this study have significant impact on the selection or design of autotrophic organisms for biofuel or biochemical production. Overall biomass production from solar energy is most efficient in organisms using the reductive TCA cycle, however, products derived from one metabolic precursor may be more efficiently produced using other carbon fixation pathways.     

The enzymatic and metabolic capabilities of early life

We introduce the concept of metaconsensus and employ it to make high confidence predictions of early enzyme functions and the metabolic properties that they may have produced. Several independent studies have used comparative bioinformatics methods to identify taxonomically broad features of genomic sequence data, protein structure data, and metabolic pathway data in order to predict physiological features that were present in early, ancestral life forms. But all such methods carry with them some level of technical bias. Here, we cross-reference the results of these previous studies to determine enzyme functions predicted to be ancient by multiple methods. We survey modern metabolic pathways to identify those that maintain the highest frequency of metaconsensus enzymes. Using the full set of modern reactions catalyzed by these metaconsensus enzyme functions, we reconstruct a representative metabolic network that may reflect the core metabolism of early life forms. Our results show that ten enzyme functions, four hydrolases, three transferases, one oxidoreductase, one lyase, and one ligase, are determined by metaconsensus to be present at least as late as the last universal common ancestor. Subnetworks within central metabolic processes related to sugar and starch metabolism, amino acid biosynthesis, phospholipid metabolism, and CoA biosynthesis, have high frequencies of these enzyme functions. We demonstrate that a large metabolic network can be generated from this small number of enzyme functions.     

Rapid Pathway Evolution Facilitated by Horizontal Gene Transfers across Prokaryotic Lineages

The evolutionary history of biological pathways is of general interest, especially in this post-genomic era, because it may provide clues for understanding how complex systems encoded on genomes have been organized. To explain how pathways can evolve de novo, some noteworthy models have been proposed. However, direct reconstruction of pathway evolutionary history both on a genomic scale and at the depth of the tree of life has suffered from artificial effects in estimating the gene content of ancestral species. Recently, we developed an algorithm that effectively reconstructs gene-content evolution without these artificial effects, and we applied it to this problem. The carefully reconstructed history, which was based on the metabolic pathways of 160 prokaryotic species, confirmed that pathways have grown beyond the random acquisition of individual genes. Pathway acquisition took place quickly, probably eliminating the difficulty in holding genes during the course of the pathway evolution. This rapid evolution was due to massive horizontal gene transfers as gene groups, some of which were possibly operon transfers, which would convey existing pathways but not be able to generate novel pathways. To this end, we analyzed how these pathways originally appeared and found that the original acquisition of pathways occurred more contemporaneously than expected across different phylogenetic clades. As a possible model to explain this observation, we propose that novel pathway evolution may be facilitated by bidirectional horizontal gene transfers in prokaryotic communities. Such a model would complement existing pathway evolution models.     

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.     

The complete genome sequence of Moorella thermoacetica (f. Clostridium thermoaceticum)

This paper describes the genome sequence of Moorella thermoacetica (f. Clostridium thermoaceticum), which is the model acetogenic bacterium that has been widely used for elucidating the Wood-Ljungdahl pathway of CO and CO(2) fixation. This pathway, which is also known as the reductive acetyl-CoA pathway, allows acetogenic (often called homoacetogenic) bacteria to convert glucose stoichiometrically into 3 mol of acetate and to grow autotrophically using H(2) and CO as electron donors and CO(2) as an electron acceptor. Methanogenic archaea use this pathway in reverse to grow by converting acetate into methane and CO(2). Acetogenic bacteria also couple the Wood-Ljungdahl pathway to a variety of other pathways to allow the metabolism of a wide variety of carbon sources and electron donors (sugars, carboxylic acids, alcohols and aromatic compounds) and electron acceptors (CO(2), nitrate, nitrite, thiosulfate, dimethylsulfoxide and aromatic carboxyl groups). The genome consists of a single circular 2 628 784 bp chromosome encoding 2615 open reading frames (ORFs), which includes 2523 predicted protein-encoding genes. Of these, 1834 genes (70.13%) have been assigned tentative functions, 665 (25.43%) matched genes of unknown function, and the remaining 24 (0.92%) had no database match. A total of 2384 (91.17%) of the ORFs in the M. thermoacetica genome can be grouped in orthologue clusters. This first genome sequence of an acetogenic bacterium provides important information related to how acetogens engage their extreme metabolic diversity by switching among different carbon substrates and electron donors/acceptors and how they conserve energy by anaerobic respiration. Our genome analysis indicates that the key genetic trait for homoacetogenesis is the core acs gene cluster of the Wood-Ljungdahl pathway.     

Techniques for determining protein sequence evolution applied to primates

Each amino acid in a protein is considered to be an individual, mutable characteristic of the species from which the protein is extracted. For a branching tree representing the evolutionary history of the known sequences in different species, our computer programs use majority logic and parsimony of mutations to determine the most likely ancestral amino acid for each position of the protein at each node of the tree. The number of mutations necessary between the ancestral and present species is summed for each branch and the entire tree. The programs then move branches to make many different configurations, from which we select the one with the minimum number of mutations as the most likely evolutionary history. We used this method to elucidate primate phylogeny from sequences of fibrinopeptides, carbonic anhydrase, and the hemoglobin beta, delta and alpha chains. All available sequences indicate that the early Pongidae had diverged into two lines before the divergence of an ancestor for the human line alone. We have constructed some probable ancestral sequences at major points during primate evolution and have developed tentative trees showing the order of divergences and evolutionary distances among primate groups. Further questions on primate evolution could be answered in the future by the detemination of the appropriate sequences.     

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.     

The method of parsimony: an experimental test and theoretical analysis of the adequacy of molecular restoration studies

The ability of the principle of parsimony to accurately reconstruct molecular evolutionary pathways from an analysis of amino acid or nucleic acid sequences from extant organisms is tested by direct comparison with a known pathway. Topological errors occur under specified conditions. Importantly, given no errors in the topology, and error-free experimental sequences, the ancestral sequences inferred by the parsimony principle err significantly, the magnitude of the error increasing with the distance of the nodal sequence from the present. These errors are irreducible as an inherent consequence of any evolutionary process in which chance processes operate within the constraints imposed by Darwinian selection. Formulae are derived which predict the errors in the ancestral sequences from a knowledge of only the internodal distances. The parsimony solution is not a reliably good solution. It is necessary to develop a detailed understanding of the interaction between chance processes and natural selection to further advance our understanding of molecular change in proteins and nucleic acids.     

The primordial metabolism: an ancestral interconnection between leucine,arginine, and lysine biosynthesis

Background

It is generally assumed that primordial cells had small genomes with simple genes coding for enzymes able to react with a wide range of chemically related substrates, interconnecting different metabolic routes. New genes coding for enzymes with a narrowed substrate specificity arose by paralogous duplication(s) of ancestral ones and evolutionary divergence. In this way new metabolic pathways were built up by primordial cells. Useful hints to disclose the origin and evolution of ancestral metabolic routes and their interconnections can be obtained by comparing sequences of enzymes involved in the same or different metabolic routes. From this viewpoint, the lysine, arginine, and leucine biosynthetic routes represent very interesting study-models. Some of the lys, arg and leu genes are paralogs; this led to the suggestion that their ancestor genes might interconnect the three pathways. The aim of this work was to trace the evolutionary pathway leading to the appearance of the extant biosynthetic routes and to try to disclose the interrelationships existing between them and other pathways in the early stages of cellular evolution.

Results

The comparative analysis of the genes involved in the biosynthesis of lysine, leucine, and arginine, their phylogenetic distribution and analysis revealed that the extant metabolic "grids" and their interrelationships might be the outcome of a cascade of duplication of ancestral genes that, according to the patchwork hypothesis, coded for unspecific enzymes able to react with a wide range of substrates. These genes belonged to a single common pathway in which the three biosynthetic routes were highly interconnected between them and also to methionine, threonine, and cell wall biosynthesis. A possible evolutionary model leading to the extant metabolic scenarios was also depicted.

Conclusion

The whole body of data obtained in this work suggests that primordial cells synthesized leucine, lysine, and arginine through a single common metabolic pathway, whose genes underwent a set of duplication events, most of which can have predated the appearance of the last common universal ancestor of the three cell domains (Archaea, Bacteria, and Eucaryotes). The model proposes a relative timing for the appearance of the three routes and also suggests a possible evolutionary pathway for the assembly of bacterial cell-wall.