Comments on the Role of H2S in the Chemistry of Earth's Early Atmosphere and in Prebiotic Synthesis

Free energy calculations and experimental measurements have been used to show that H₂S/CO₂ mixtures outgassing from a prebiotic Earth's crust would have produced a reducing gas mixture containing CO, H₂, H₂O, and S _x as principal components. Due to rapid recombination of H₂, CO, and S _x to H₂S and CO₂ on cooling from a high temperature to ambient conditions, reducing components would have been retained only if efficient quenching of the reduced gas mixture had been possible. Consequently, subsea vents or vents with efficient infusion of water would have been ideal sites for retention of reduced species and for prebiotic organic synthesis. It is suggested that C/H/O/S ratios are important factors in controlling the degree of prebiotic organic synthesis and, hence, the emergence of life, since if oxygen is abundant, CO₂ and SO₂ would have been dominant species. Received: 5 March 1997 / Accepted: 15 December 1997     

Relics from the RNA World

An RNA world is widely accepted as a probable stage in the early evolution of life. Two implications are that proteins have gradually replaced RNA as the main biological catalysts and that RNA has not taken on any major de novo catalytic function after the evolution of protein synthesis, that is, there is an essentially irreversible series of steps RNA → RNP → protein. This transition, as expected from a consideration of catalytic perfection, is essentially complete for reactions when the substrates are small molecules. Based on these principles we derive criteria for identifying RNAs in modern organisms that are relics from the RNA world and then examine the function and phylogenetic distribution of RNA for such remnants of the RNA world. This allows an estimate of the minimum complexity of the last ribo-organism—the stage just preceding the advent of genetically encoded protein synthesis. Despite the constraints placed on its size by a low fidelity of replication (the Eigen limit), we conclude that the genome of this organism reached a considerable level of complexity that included several RNA-processing steps. It would include a large protoribosome with many smaller RNAs involved in its assembly, pre-tRNAs and tRNA processing, an ability for recombination of RNA, some RNA editing, an ability to copy to the end of each RNA strand, and some transport functions. It is harder to recognize specific metabolic reactions that must have existed but synthetic and bio-energetic functions would be necessary. Overall, this requires that such an organism maintained a multiple copy, double-stranded linear RNA genome capable of recombination and splicing. The genome was most likely fragmented, allowing each ``chromosome' to be replicated with minimum error, that is, within the Eigen limit. The model as developed serves as an outgroup to root the tree of life and is an alternative to using sequence data for inferring properties of the earliest cells. Received: 14 January 1997 / Accepted: 19 May 1997     

Site-Specific Prebiotic Oligomerization Reactions of Glycine on the Surface of Hectorite

Condensation reactions of the amino acid glycine on the surface of Cu(II)-exchanged hectorite are investigated using the technique of scanning force microscopy. Prebiotic conditions are simulated using alternate wetting and heating cycles. Concentration, immobilization, and subsequent polymerization resulting in glycine oligomers are seen to occur primarily at step edges or faults in the topmost layer. Condensation reactions also occur within tiny micropores or defects in the topmost layer. These reactions are facilitated by the availability of intergallery metal cations at the step edges or pores in the surface region. Received: 19 January 1998 / Accepted: 24 April 1998     

Reactive Phosphonic Acids as Prebiotic Carriers of Phosphorus

Phosphonic acids are the only phosphorus-containing organic compounds detected in the Murchison meteorite. We earlier described the synthesis of methyl-, hydroxymethyl-, and 1-hydroxyethyl phosphonic acids using sodium phosphite as a source of phosphite radicals. We now show that ultraviolet irradiation of dilute aqueous solutions of acetylene in the presence of sodium phosphite leads to the synthesis of vinyl phosphonic acid. At neutral to basic pH, vinyl phosphonic acid reacts under photochemical conditions to produce phosphonoacetaldehyde and 2-hydroxyethyl phosphonic acid as the major products, as well as smaller yields of 1-hydroxyethyl phosphonic acid, phosphonoacetic acid, and ethyl phosphonic acid. Of these products, phosphonoacetaldehyde is particularly interesting as a potential precursor of prebiotic carbohydrate derivatives. Received: 22 July 1996 / Accepted: 13 September 1996     

The Path from the RNA World

We describe a sequential (step by step) Darwinian model for the evolution of life from the late stages of the RNA world through to the emergence of eukaryotes and prokaryotes. The starting point is our model, derived from current RNA activity, of the RNA world just prior to the advent of genetically-encoded protein synthesis. By focusing on the function of the protoribosome we develop a plausible model for the evolution of a protein-synthesizing ribosome from a high-fidelity RNA polymerase that incorporated triplets of oligonucleotides. With the standard assumption that during the evolution of enzymatic activity, catalysis is transferred from RNA → RNP → protein, the first proteins in the ``breakthrough organism' (the first to have encoded protein synthesis) would be nonspecific chaperone-like proteins rather than catalytic. Moreover, because some RNA molecules that pre-date protein synthesis under this model now occur as introns in some of the very earliest proteins, the model predicts these particular introns are older than the exons surrounding them, the ``introns-first' theory. Many features of the model for the genome organization in the final RNA world ribo-organism are more prevalent in the eukaryotic genome and we suggest that the prokaryotic genome organization (a single, circular genome with one center of replication) was derived from a ``eukaryotic-like' genome organization (a fragmented linear genome with multiple centers of replication). The steps from the proposed ribo-organism RNA genome → eukaryotic-like DNA genome → prokaryotic-like DNA genome are all relatively straightforward, whereas the transition prokaryotic-like genome → eukaryotic-like genome appears impossible under a Darwinian mechanism of evolution, given the assumption of the transition RNA → RNP → protein. A likely molecular mechanism, ``plasmid transfer,' is available for the origin of prokaryotic-type genomes from an eukaryotic-like architecture. Under this model prokaryotes are considered specialized and derived with reduced dependence on ssRNA biochemistry. A functional explanation is that prokaryote ancestors underwent selection for thermophily (high temperature) and/or for rapid reproduction (r selection) at least once in their history. Received: 14 January 1997 / Accepted: 19 May 1997     

The Nicotinamide Biosynthetic Pathway Is a By-Product of the RNA World

Many of the biosynthetic pathways, especially those leading to the coenzymes, must have originated very early, perhaps before enzymes were available to catalyze their synthesis. While a number of enzymatic reactions in metabolism are known to proceed nonenzymatically, there are no examples of entire metabolic sequences that can be achieved in this manner. The most primitive pathway for nicotinic acid biosynthesis is the reaction of aspartic acid with dihydroxyacetone phosphate. We report here that nicotinic acid (NAc) and its metabolic precursor, quinolinic acid (QA), are produced in yields as high as 7% in a six-step nonenzymatic sequence from aspartic acid and dihydroxyacetone phosphate (DHAP). The biosynthesis of ribose phosphate could have produced DHAP and other three carbon compounds. Aspartic acid could have been available from prebiotic synthesis or from the ribozyme synthesis of pyrimidines. These results suggest that NAD could have originated in the RNA world and that the nonenzymatic biosynthesis of the cofactor nicotinamide could have been an inevitable consequence of life based on carbohydrates and amino acids. The enzymes of the modern pathway were later added in any order. Received: 22 May 2000 / Accepted: 7 August 2000     

Conserved Gene Clusters in Bacterial Genomes Provide Further Support for the Primacy of RNA

Five complete bacterial genome sequences have been released to the scientific community. These include four (eu)Bacteria, Haemophilus influenzae, Mycoplasma genitalium, M. pneumoniae, and Synechocystis PCC 6803, as well as one Archaeon, Methanococcus jannaschii. Features of organization shared by these genomes are likely to have arisen very early in the history of the bacteria and thus can be expected to provide further insight into the nature of early ancestors. Results of a genome comparison of these five organisms confirm earlier observations that gene order is remarkably unpreserved. There are, nevertheless, at least 16 clusters of two or more genes whose order remains the same among the four (eu)Bacteria and these are presumed to reflect conserved elements of coordinated gene expression that require gene proximity. Eight of these gene orders are essentially conserved in the Archaea as well. Many of these clusters are known to be regulated by RNA-level mechanisms in Escherichia coli, which supports the earlier suggestion that this type of regulation of gene expression may have arisen very early. We conclude that although the last common ancestor may have had a DNA genome, it likely was preceded by progenotes with an RNA genome. Received: 10 March 1996 / Accepted: 20 May 1997     

RNA Folding Argues Against a Hot-Start Origin of Life

Opinion is strongly divided on whether life arose on earth under hot or cold conditions, the hot-start and cold-start scenarios, respectively. The origin of life close to deep thermal vents appears as the majority opinion among biologists, but there is considerable biochemical evidence that high temperatures are incompatible with an RNA world. To be functional, RNA has to fold into a three-dimensional structure. We report both theoretical and experimental results on RNA folding and show that (as expected) hot conditions strongly reduce RNA folding. The theoretical results come from energy-minimization calculations of the average extent of folding of RNA, mainly from 0–90°C, for both random sequences and tRNA sequences. The experimental results are from circular-dichroism measurements of tRNA over a similar range of temperatures. The quantitative agreement between calculations and experiment is remarkable, even to the shape of the curves indicating the cooperative nature of RNA folding and unfolding. These results provide additional evidence for a lower temperature stage being necessary in the origin of life. Received: 1 March 2000 / Accepted: 14 June 2000     

Question 5: On the Chemical Reality of the RNA World

The discovery of catalytic RNA has revolutionised modern molecular biology and bears important implications for the origin of Life research. Catalytic RNA, in particular self-replicating RNA, prompted the hypothesis of an early “RNA world” where RNA molecules played all major roles such information storage and catalysis. The actual role of RNA as primary actor in the origin of life has been under debate for a long time, with a particular emphasis on possible pathways to the prebiotic synthesis of mononucleotides; their polymerization and the possibility of spontaneous emergence of catalytic RNAs synthesised under plausible prebiotic conditions. However, little emphasis has been put on the chemical reality of an RNA world; in particular concerning the chemical constrains that such scenario should have met to be feasible. This paper intends to address those concerns with regard to the achievement of high local RNA molecules concentration and the aetiology of unique sequence under plausible prebiotic conditions. Presented at: International School of Complexity – 4th Course: Basic Questions on the Origins of Life; “Ettore Majorana” Foundation and Centre for Scientific Culture, Erice, Italy, 1–6 October 2006.     

On the Origin of Operons and Their Possible Role in Evolution Toward Thermophily

We present a hypothesis suggesting that close linkage of functionally related anabolic genes and their ultimate integration into operons developed under selective pressure as a molecular strategy which contributed to the viability of ancestral thermophilic cells. Cotranslation of functionally related proteins is viewed as having facilitated the formation of multienzyme complexes channeling thermolabile substrates and the mutual stabilization of inherently thermolabile proteins. In this perspective, the evolutionary scheme considered the most probable is the evolution of both Bacteria and Archaea by thermoreduction (Forterre 1995) from a mesophilic, protoeukaryotic last common ancestor (LCA) endowed with appreciable genetic redundancy.     

Estimating the Contributions of Selection and Self-Organization in RNA Secondary Structure

In addition to characteristic structural properties imposed by evolutionary modification, evolved, single-stranded RNAs also display characteristic structural properties imposed by intrinsic physical constraints on RNA polymer folding. The balance of intrinsic and functionally selected characters in the folded conformation of evolved secondary structures was determined by comparing the predicted secondary structures of evolved and unevolved (random) RNA sequences. Though evolved conformations are significantly more ordered than conformations of random-sequence RNA, this analysis demonstrates that the majority of conformational order within evolved structures results not from evolutionary optimization but from constraints imposed by rules intrinsic to RNA polymer folding. Received: 25 November 1998 / Accepted: 12 February 1999     

Contribution of Mutation and RNA Recombination to the Evolution of a Plant Pathogenic RNA

The nucleotide sequence of 17 variants of the satellite RNA of cucumber mosaic virus (CMV-satRNA) isolated from field-infected tomato plants in the springs of 1989, 1990, and 1991 was determined. The sequence of each of the 17 satRNAs was unique and was between 334 and 340 nucleotides in length; 57 positions were polymorphic. There was much genetic divergence, ranging from 0.006 to 0.141 nucleotide substitutions per site for pairwise comparisons, and averaging 0.074 for any pair. When the polymorphic positions were analyzed relative to a secondary structure model proposed for CMV-satRNAs, it was found that there were significantly different numbers of changes in base-paired and non–base-paired positions, and that mutations that did not disrupt base pairing were preferred at the putatively paired sites. This supports the concept that the need to maintain a functional structure may limit genetic divergence of CMV-satRNA. Phylogenetic analyses showed that the 17 CMV-satRNA variants clustered into two subgroups, I and II, and evolutionary lines proceeding by the sequential accumulation of mutations were apparent. Three satRNA variants were outliers for these two phylogenetic groups. They were shown to be recombinants of subgroup I and II satRNAs by calculating phylogenies for different molecular regions and by using Sawyer's test for gene conversion. At least two recombination events were required to produce these three recombinant satRNAs. Thus, recombinants were found to be frequent (∼17%) in natural populations of CMV-satRNA, and recombination may make an important contribution to the generation of new variants. To our knowledge this is the first report of data allowing the frequency of recombinant isolates in natural populations of an RNA replicon to be estimated. Received: 14 May 1996 / Accepted: 17 July 1996     

The Vertebrate 7S K RNA Separates Hagfish (Myxine glutinosa) and Lamprey (Lampetra fluviatilis)

7S RNA sequences from the hagfish (Myxiniformes) and lamprey (Petromyzontiformes) were cloned and analyzed. In both species, 7S L RNA (also designated SRP RNA, since it represents the RNA constituent of the signal recognition particle) was clearly detectable. The sequence similarity between the two species was 86%, compared with about 75% similarity between either of these species and mammals. 7S K RNA was also cloned from the lamprey. The similarity between the 7S K RNA of the lamprey and that of mammals was 68%. Interestingly, several interspersed elements were found with nearly 100% similarity compared with mammals. In contrast to the lamprey, no 7S K RNA-related sequences were detectable among hagfish RNA, neither in northern blots nor with the PCR assay. In view of the significant conservation between the 7S K RNA of lamprey and that of mammals (human), this unexpected result clearly separates lamprey and hagfish. In addition, the lack of detectable 7S K RNA sequences in an outgroup, such as amphioxus, indicates that these results do not reflect an autapomorphy of hagfish. Therefore, our data provide additional support to the notion of a sister group relationship between Petromyzontiformes and gnathostomous vertebrates to the exclusion of Myxiniformes. Received: 24 September 1999 / Accepted: 9 February 2000     

Use of RNA Secondary Structure for Studying the Evolution of RNase P and RNase MRP

Secondary structure is evaluated for determining evolutionary relationships between catalytic RNA molecules that are so distantly related they are scarcely alignable. The ribonucleoproteins RNase P (P) and RNase MRP (MRP) have been suggested to be evolutionarily related because of similarities in both function and secondary structure. However, their RNA sequences cannot be aligned with any confidence, and this leads to uncertainty in any trees inferred from sequences. We report several approaches to using secondary structures for inferring evolutionary trees and emphasize quantitative tests to demonstrate that evolutionary information can be recovered. For P and MRP, three hypotheses for the relatedness are considered. The first is that MRP is derived from P in early eukaryotes. The next is that MRP is derived from P from an early endosymbiont. The third is that both P and MRP evolved in the RNA-world (and the need for MRP has since been lost in prokaryotes). Quantitative comparisons of the pRNA and mrpRNA secondary structures have found that the possibility of an organellar origin of MRP is unlikely. In addition, comparison of secondary structures support the identity of an RNase P–like sequence in the maize chloroplast genome. Overall, it is concluded that RNA secondary structure is useful for evaluating evolutionary relatedness, even with sequences that cannot be aligned with confidence. Received: 19 July 1999 / Accepted: 3 May 2000     

Mineral Catalysis of a Potentially Prebiotic Aldol Condensation

Minerals may have played a significant role in chemical evolution. In the course of investigating the chemistry of phosphonoacetaldehyde (PAL), an analogue of glycolaldehyde phosphate, we have observed a striking case of catalysis by the layered hydroxide mineral hydrotalcite ([Mg₂Al(OH)₆][Cl.nH₂O]). In neutral or moderately basic aqueous solutions, PAL is unreactive even at a concentration of 0.1 M. In the presence of a large excess of NaOH (2 M), the compound undergoes aldol condensation to produce a dimer containing a C3–C4 double-bond. In dilute neutral solutions and in the presence of the mineral, however, condensation takes place rapidly, to produce a dimer which is almost exclusively the C2–C3 unsaturated product. Received: 11 February 1998 / Accepted: 12 May 1998     

Evolved RNA Secondary Structure and the Rooting of the Universal Tree of Life

The origin and diversification of RNA secondary structure were traced using cladistic methods. Structural components were coded as polarized and ordered multi-state characters, following a model of character state transformation outlined by considerations in statistical mechanics. Several classes of functional RNA were analyzed, including ribosomal RNA (rRNA). Considerable phylogenetic signal was present in their secondary structure. The intrinsically rooted phylogenies reconstructed from evolved RNA structure depicted those derived from nucleic acid sequence at all taxonomical levels, and grouped organisms in concordance with traditional classification, especially in the archaeal and eukaryal domains. Natural selection appears therefore to operate early in the information flow that originates in sequence and ends in an adapted phenotype. When examining the hierarchical classification of the living world, phylogenetic analysis of secondary structure of the small and large rRNA subunits reconstructed a universal tree of life that branched in three monophyletic groups corresponding to Eucarya, Archaea, and Bacteria, and was rooted in the eukaryotic branch. Ribosomal characters involved in the translational cycle could be easily traced and showed that transfer RNA (tRNA) binding domains in the large rRNA subunit evolved concurrently with the rest of the rRNA molecule. Results suggest it is equally parsimonious to consider that ancestral unicellular eukaryotes or prokaryotes gave rise to all extant life forms and provide a rare insight into the early evolution of nucleic acid and protein biosynthesis. Received: 13 September 2000 / Accepted: 27 August 2001     

Energy from Redox Disproportionation of Sugar Carbon Drives Biotic and Abiotic Synthesis

To identify the energy source that drives the biosynthesis of amino acids, lipids, and nucleotides from glucose, we calculated the free energy change due to redox disproportionation of the substrate carbon of (1) 26-carbon fermentation reactions and (2) the biosynthesis of amino acids and lipids of E. coli from glucose. The free energy (cal/mmol of carbon) of these reactions was plotted as a function of the degree of redox disproportionation of carbon (disproportionative electron transfers (mmol)/mmol of carbon). The zero intercept and proportionality between energy yield and degree of redox disporportionation exhibited by this plot demonstrate that redox disproportionation is the principal energy source of these redox reactions (slope of linear fit =−10.4 cal/mmol of disproportionative electron transfers). The energy and disproportionation values of E. coli amino acid and lipid biosynthesis from glucose lie near this linear curve fit with redox disproportionation accounting for 84% and 96% (and ATP only 6% and 1%) of the total energy of amino acid and lipid biosynthesis, respectively. These observations establish that redox disproportionation of carbon, and not ATP, is the primary energy source driving amino acid and lipid biosynthesis from glucose. In contrast, we found that nucleotide biosynthesis involves very little redox disproportionation, and consequently depends almost entirely on ATP for energy. The function of sugar redox disproportionation as the major source of free energy for the biosynthesis of amino acids and lipids suggests that sugar disproportionation played a central role in the origin of metabolism, and probably the origin of life. Received: 18 April 1996 / Accepted: 31 October 1996     

Selective Pressures on RNA Hairpins In Vivo and In Vitro

Comparison of the most stable potential hairpins in the sequences of natural ribozymes with those in the randomized sequences has revealed that the hairpin loop energies are lower than expected by chance. Although these hairpins are not necessarily parts of functional structures, there is a selective pressure to diminish the destabilizing free energies of the hairpin loops. In contrast, no significant bias is observed in the stacking values of the most stable stems. In the ribozymes isolated in vitro the loops of potential hairpins are closer to random values, which can result in less efficient folding rates. Furthermore, the effects of kinetic traps seem to be more significant in the folding pathways of the in vitro isolates due to a potential to form stable stacks incompatible with the functional folds. Similarly to natural ribozyme sequences, the untranslated regions of viral RNAs also form hairpins with relatively low loop free energies. These evolutionary trends suggest ways for efficient engineering of improved RNA constructs on the basis of analysis of in vitro isolates and approaches for the search of regions coding for functional RNA structures in large genome sequences. Received: 12 January 2001 / Accepted: 21 May 2001