The virtual circular genome model for primordial RNA replication

Nonenzymatic RNA copying with a potentially primordial genetic alphabet

Nonenzymatic RNA copying is thought to have been responsible for the replication of genetic information during the origin of life. However, chemical copying with the canonical nucleotides (A, U, G, and C) strongly favors the incorporation of G and C and disfavors the incorporation of A and especially U because of the stronger G:C vs. A:U base pair and the weaker stacking interactions of U. Recent advances in prebiotic chemistry suggest that the 2-thiopyrimidines were precursors to the canonical pyrimidines, raising the possibility that they may have played an important early role in RNA copying chemistry. Furthermore, 2-thiouridine (s2U) and inosine (I) form by deamination of 2-thiocytidine (s2C) and A, respectively. We used thermodynamic and crystallographic analyses to compare the I:s2C and A:s2U base pairs. We find that the I:s2C base pair is isomorphic and isoenergetic with the A:s2U base pair. The I:s2C base pair is weaker than a canonical G:C base pair, while the A:s2U base pair is stronger than the canonical A:U base pair, so that a genetic alphabet consisting of s2U, s2C, I, and A generates RNA duplexes with uniform base pairing energies. Consistent with these results, kinetic analysis of nonenzymatic template-directed primer extension reactions reveals that s2C and s2U substrates bind similarly to I and A in the template, and vice versa. Our work supports the plausibility of a potentially primordial genetic alphabet consisting of s2U, s2C, I, and A and offers a potential solution to the long-standing problem of biased nucleotide incorporation during nonenzymatic template copying.

Autocatalytic assembly of a chimeric aminoacyl-RNA synthetase ribozyme

Autocatalytic reactions driving the self-assembly of biological polymers are important for the origin of life, yet few experimental examples of such reactions exist. Here we report an autocatalytic assembly pathway that generates a chimeric, amino acid-bridged aminoacyl-RNA synthetase ribozyme. The noncovalent complex of ribozyme fragments initiates low-level aminoacylation of one of the fragments, which, after loop-closing ligation, generates a highly active covalently linked chimeric ribozyme. The generation of this ribozyme is increasingly efficient over time due to the autocatalytic assembly cycle that sustains the ribozyme over indefinite cycles of serial dilution. Because of its trans activity, this ribozyme also assembles ribozymes distinct from itself, such as the hammerhead, suggesting that RNA aminoacylation, coupled with nonenzymatic ligation, could have facilitated the emergence and propagation of ribozymes.

Prebiotic RNA self-assembling and the origin of life: Mechanistic and molecular modeling rationale for explaining the prebiotic origin and replication of RNA

In recent years, a great interest has been focused on the prebiotic origin of nucleic acids and life on Earth. An attractive idea is that life was initially based on an autocatalytic and autoreplicative RNA (the RNA-world). RNA duplexes are right-handed helical chains with antiparallel orientation, but the rationale for these features is not yet known. An antiparallel (inverted) stacking of purine nucleosides was reported in crystallographic studies. Molecular modeling also supports the inverted orientation of nucleosides. This preferential stacking can also appear when nucleosides are included in a montmorillonite clay matrix. Free-energy values and geometrical parameters show that D-ribose chirality is preferred for the formation of right-handed RNA molecules. Thus, a "zipper" model with antiparallel and auto-intercalated nucleosides linked by phosphate groups can be proposed to form single RNA chains. Unstacking with strand separation and base pairing by H-bonding, results in shortening and inclination of ribose-phosphate chains, leading to right-handed helicity and antiparallel duplexes. Incorporation of complementary precursors on the major groove template by a self-assembly mechanism provides a prebiotic (non-enzymatic) "tetris" replication model by formation of a transient RNA tetrad and tetraplex. Original hairpin motifs appear as simple building units that form typical RNA structures such as hammerheads, cloverleaves and dumbbells. They occur today in the circular viroids and virusoids, as well as in highly branched and complex rRNA molecules.

Maintaining an autocatalytic reaction system in a protocell: Nonenzymatic RNA templating and the link between replication and metabolism

The first protocells must have been driven by a reaction system in which autocatalysis is maintained inside the cell while food molecules outside the cell are unreactive. We show that if there is a second-order autocatalytic mechanism, then an active state can be stable inside the cell with a high catalyst concentration while the environment remains stable in an unreactive state with no catalyst. Addition of a small amount of catalyst to the environment does not cause the initiation of the autocatalytic cycle outside the cell. In contrast, for a first-order mechanism, addition of a small amount of catalyst initiates the reaction outside the cell unless there is continual removal of the catalyst from the environment. Hence, a second-order reaction mechanism maintains the difference between the inside and outside of a protocell much better than a first-order mechanism. The formose reaction, although a prebiotically plausible autocatalytic system, is first order and therefore is unlikely as a means of support for the first protocells. We give other theoretical examples of first- and second-order reaction networks but note there are few known real-world chemical systems that fit these schemes. However, we show that nonenzymatic RNA templating constitutes a second-order autocatalytic system with the necessary properties to support a protocell. Templating is maintained inside the cell but is not initiated outside the cell. If the reaction is driven by an external supply of activated nucleotides, then templating is itself a metabolic cycle. It is not necessary to have an additional separate metabolic cycle before templating reactions can occur. In this view, templating reactions, which are usually considered as the origin of replication and heredity, are also the origin of metabolism.

Making, Using, and Understanding Molecular Systems: The 57th Bürgenstock Conference

The 57th SCS Conference on Stereochemistry, better known as the Bürgenstock Conference, brought together a diverse range of chemistry expertise in Brunnen, Switzerland.

Overcoming nucleotide bias in the nonenzymatic copying of RNA templates

The RNA World hypothesis posits that RNA was the molecule of both heredity and function during the emergence of life. This hypothesis implies that RNA templates can be copied, and ultimately replicated, without the catalytic aid of evolved enzymes. A major problem with nonenzymatic template-directed polymerization has been the very poor copying of sequences containing rA and rU. Here, we overcome that problem by using a prebiotically plausible mixture of RNA mononucleotides and random-sequence oligonucleotides, all activated by methyl isocyanide chemistry, that direct the uniform copying of arbitrary-sequence templates, including those harboring rA and rU. We further show that the use of this mixture in copying reactions suppresses copying errors while also generating a more uniform distribution of mismatches than observed for simpler systems. We find that oligonucleotide competition for template binding sites, oligonucleotide ligation and the template binding properties of reactant intermediates work together to reduce product sequence bias and errors. Finally, we show that iterative cycling of templated polymerization and activation chemistry improves the yields of random-sequence products. These results for random-sequence template copying are a significant advance in the pursuit of nonenzymatic RNA replication.

A potential role for RNA aminoacylation prior to its role in peptide synthesis

Coded ribosomal peptide synthesis could not have evolved unless its sequence and amino acid-specific aminoacylated tRNA substrates already existed. We therefore wondered whether aminoacylated RNAs might have served some primordial function prior to their role in protein synthesis. Here, we show that specific RNA sequences can be nonenzymatically aminoacylated and ligated to produce amino acid-bridged stem-loop RNAs. We used deep sequencing to identify RNAs that undergo highly efficient glycine aminoacylation followed by loop-closing ligation. The crystal structure of one such glycine-bridged RNA hairpin reveals a compact internally stabilized structure with the same eponymous T-loop architecture that is found in many noncoding RNAs, including the modern tRNA. We demonstrate that the T-loop-assisted amino acid bridging of RNA oligonucleotides enables the rapid template-free assembly of a chimeric version of an aminoacyl-RNA synthetase ribozyme. We suggest that the primordial assembly of amino acid-bridged chimeric ribozymes provides a direct and facile route for the covalent incorporation of amino acids into RNA. A greater functionality of covalently incorporated amino acids could contribute to enhanced ribozyme catalysis, providing a driving force for the evolution of sequence and amino acid-specific aminoacyl-RNA synthetase ribozymes in the RNA World. The synthesis of specifically aminoacylated RNAs, an unlikely prospect for nonenzymatic reactions but a likely one for ribozymes, could have set the stage for the subsequent evolution of coded protein synthesis.

Diaminopurine in Nonenzymatic RNA Template Copying

In the RNA World before the emergence of an RNA polymerase, nonenzymatic template copying would have been essential for the transmission of genetic information. However, the products of chemical copying with the canonical nucleotides (A, U, C, and G) are heavily biased toward the incorporation of G and C, which form a more stable base pair than A and U. We therefore asked whether replacing adenine (A) with diaminopurine (D) might lead to more efficient and less biased nonenzymatic template copying by making a stronger version of the A:U pair. As expected, primer extension substrates containing D bound to U in the template more tightly than substrates containing A. However, primer extension with D exhibited elevated reaction rates on a C template, leading to concerns about fidelity. Our crystallographic studies revealed the nature of the D:C mismatch by showing that D can form a wobble-type base pair with C. We then asked whether competition with G would decrease the mismatched primer extension. We performed nonenzymatic primer extension with all four activated nucleotides on randomized RNA templates containing all four letters and used deep sequencing to analyze the products. We found that the DUCG genetic system exhibited a more even product distribution and a lower mismatch frequency than the canonical AUCG system. Furthermore, primer extension is greatly reduced following all mismatches, including the D:C mismatch. Our study suggests that D deserves further attention for its possible role in the RNA World and as a potentially useful component of artificial nonenzymatic RNA replication systems.

Constraints on the emergence of RNA through non-templated primer extension with mixtures of potentially prebiotic nucleotides

The emergence of RNA on the early Earth is likely to have been influenced by chemical and physical processes that acted to filter out various alternative nucleic acids. For example, UV photostability is thought to have favored the survival of the canonical nucleotides. In a recent proposal for the prebiotic synthesis of the building blocks of RNA, ribonucleotides share a common pathway with arabino- and threo-nucleotides. We have therefore investigated non-templated primer extension with 2-aminoimidazole-activated forms of these alternative nucleotides to see if the synthesis of the first oligonucleotides might have been biased in favor of RNA. We show that non-templated primer extension occurs predominantly through 5'-5' imidazolium-bridged dinucleotides, echoing the mechanism of template-directed primer extension. Ribo- and arabino-nucleotides exhibited comparable rates and yields of non-templated primer extension, whereas threo-nucleotides showed lower reactivity. Competition experiments confirmed the bias against the incorporation of threo-nucleotides. The incorporation of an arabino-nucleotide at the end of the primer acts as a chain terminator and blocks subsequent extension. These biases, coupled with potentially selective prebiotic synthesis, and the templated copying that is known to favour the incorporation of ribonucleotides, provide a plausible model for the effective exclusion of arabino- and threo-nucleotides from primordial oligonucleotides.

Emergence of catalytic function in prebiotic information-coding polymers

Life as we know it relies on the interplay between catalytic activity and information processing carried out by biological polymers. Here we present a plausible pathway by which a pool of prebiotic information-coding oligomers could acquire an early catalytic function, namely sequence-specific cleavage activity. Starting with a system capable of non-enzymatic templated replication, we demonstrate that even non-catalyzed spontaneous cleavage would promote proliferation by generating short fragments that act as primers. Furthermore, we show that catalytic cleavage function can naturally emerge and proliferate in this system. Specifically, a cooperative catalytic network with four subpopulations of oligomers is selected by the evolution in competition with chains lacking catalytic activity. The cooperative system emerges through the functional differentiation of oligomers into catalysts and their substrates. The model is inspired by the structure of the hammerhead RNA enzyme as well as other DNA- and RNA-based enzymes with cleavage activity that readily emerge through natural or artificial selection. We identify the conditions necessary for the emergence of the cooperative catalytic network. In particular, we show that it requires the catalytic rate enhancement over the spontaneous cleavage rate to be at least 102-103, a factor consistent with the existing experiments. The evolutionary pressure leads to a further increase in catalytic efficiency. The presented mechanism provides an escape route from a relatively simple pairwise replication of oligomers toward a more complex behavior involving catalytic function. This provides a bridge between the information-first origin of life scenarios and the paradigm of autocatalytic sets and hypercycles, albeit based on cleavage rather than synthesis of reactants.

Transient states during the annealing of mismatched and bulged oligonucleotides

Oligonucleotide hybridization is crucial in various biological, prebiotic and nanotechnological processes, including gene regulation, non-enzymatic primer extension and DNA nanodevice assembly. Although extensive research has focused on the thermodynamics and kinetics of nucleic acid hybridization, the behavior of complex mixtures and the outcome of competition for target binding remain less well understood. In this study, we investigate the impact of mismatches and bulges in a 12 bp DNA or RNA duplex on its association (kon) and dissociation (koff) kinetics. We find that such defects have relatively small effects on the association kinetics, while the dissociation kinetics vary in a position-dependent manner by up to 6 orders of magnitude. Building upon this observation, we explored a competition scenario involving multiple oligonucleotides, and observed a transient low specificity of probe hybridization to fully versus partially complementary targets in solution. We characterize these long-lived metastable states and their evolution toward equilibrium, and show that sufficiently long-lived mis-paired duplexes can serve as substrates for prebiotically relevant chemical copying reactions. Our results suggest that transient low accuracy states may spontaneously emerge within all complex nucleic acid systems comprising a large enough number of competing strands, with potential repercussions for gene regulation in the realm of modern biology and the prebiotic preservation of genetic information.

Structure-guided aminoacylation and assembly of chimeric RNAs

Coded ribosomal peptide synthesis could not have evolved unless its sequence and amino acid specific aminoacylated tRNA substrates already existed. We therefore wondered whether aminoacylated RNAs might have served some primordial function prior to their role in protein synthesis. Here we show that specific RNA sequences can be nonenzymatically aminoacylated and ligated to produce amino acid-bridged stem-loop RNAs. We used deep sequencing to identify RNAs that undergo highly efficient glycine aminoacylation followed by loop-closing ligation. The crystal structure of one such glycine-bridged RNA hairpin reveals a compact internally stabilized structure with the same eponymous T-loop architecture found in modern tRNA. We demonstrate that the T-loop assisted amino acid bridging of RNA oligonucleotides enables the rapid template-free assembly of a chimeric version of an aminoacyl-RNA synthetase ribozyme. We suggest that the primordial assembly of such chimeric ribozymes would have allowed the greater functionality of amino acids to contribute to enhanced ribozyme catalysis, providing a driving force for the evolution of sequence and amino acid specific aminoacyl-RNA synthetase enzymes prior to their role in protein synthesis.

Natural soda lakes provide compatible conditions for RNA and membrane function that could have enabled the origin of life

The origin of life likely occurred within environments that concentrated cellular precursors and enabled their co-assembly into cells. Soda lakes (those dominated by Na+ ions and carbonate species) can concentrate precursors of RNA and membranes, such as phosphate, cyanide, and fatty acids. Subsequent assembly of RNA and membranes into cells is a long-standing problem because RNA function requires divalent cations, e.g. Mg2+, but Mg2+ disrupts fatty acid membranes. The low solubility of Mg-containing carbonates limits soda lakes to moderate Mg2+ concentrations (∼1 mM), so we investigated whether both RNAs and membranes function within these lakes. We collected water from Last Chance Lake and Goodenough Lake in Canada. Because we sampled after seasonal evaporation, the lake water contained ∼1 M Na+ and ∼1 mM Mg2+ near pH 10. In the laboratory, nonenzymatic, RNA-templated polymerization of 2-aminoimidazole-activated ribonucleotides occurred at comparable rates in lake water and standard laboratory conditions (50 mM MgCl2, pH 8). Additionally, we found that a ligase ribozyme that uses oligonucleotide substrates activated with 2-aminoimidazole was active in lake water after adjusting pH from ∼10 to 9. We also observed that decanoic acid and decanol assembled into vesicles in a dilute solution that resembled lake water after seasonal rains, and that those vesicles retained encapsulated solutes despite salt-induced flocculation when the external solution was replaced with dry-season lake water. By identifying compatible conditions for nonenzymatic and ribozyme-catalyzed RNA assembly, and for encapsulation by membranes, our results suggest that soda lakes could have enabled cellular life to emerge on Earth, and perhaps elsewhere.

Combinatorial and frequency properties of the ribosome ancestors

BACKGROUND

The current ribosome has evolved from the primitive stages of life on Earth. Its function is to build proteins and on the basis of this role, we are looking for a universal common ancestor to the ribosome which could: i) present optimal combinatorial properties, and ii) have left vestiges in the current molecules composing the ribosome (rRNA or r-proteins) or helping in its construction and functioning.

METHODS

Genomic public databases are used for finding the nucleotide sequences of rRNAs and mRNA of r-proteins and statistical calculations are performed on the occurrence in these genes of some pentamers belonging to the RNA proposed as optimal ribosome ancestor.

RESULTS

After having exhibited a possible solution to the problem of an RNA capable of catalyzing peptide genesis, traces of this RNA are found in many rRNAs and mRNA of r-proteins, as well as in factors contributing to the construction of the current ribosome.

CONCLUSIONS

The existence of an optimal primordial RNA whose function is to facilitate the creation of peptide bonds between amino acids may have contributed to accelerate the emergence of the first vital processes. Its traces should be found in many living species inside structures structurally and functionally close to the ribosome, which is already the case in the species studied in this article.

Cyclophospholipids Enable a Protocellular Life Cycle

There is currently no plausible path for the emergence of a self-replicating protocell, because prevalent formulations of model protocells are built with fatty acid vesicles that cannot withstand the concentrations of Mg2+ needed for the function and replication of nucleic acids. Although prebiotic chelates increase the survivability of fatty acid vesicles, the resulting model protocells are incapable of growth and division. Here, we show that protocells made of mixtures of cyclophospholipids and fatty acids can grow and divide in the presence of Mg2+-citrate. Importantly, these protocells retain encapsulated nucleic acids during growth and division, can acquire nucleotides from their surroundings, and are compatible with the nonenzymatic extension of an RNA oligonucleotide, chemistry needed for the replication of a primitive genome. Our work shows that prebiotically plausible mixtures of lipids form protocells that are active under the conditions necessary for the emergence of Darwinian evolution.

Thermodynamics and kinetics of DNA and RNA dinucleotide hybridization to gaps and overhangs

Hybridization of short nucleic acid segments (<4 nt) to single-strand templates occurs as a critical intermediate in processes such as nonenzymatic nucleic acid replication and toehold-mediated strand displacement. These templates often contain adjacent duplex segments that stabilize base pairing with single-strand gaps or overhangs, but the thermodynamics and kinetics of hybridization in such contexts are poorly understood because of the experimental challenges of probing weak binding and rapid structural dynamics. Here we develop an approach to directly measure the thermodynamics and kinetics of DNA and RNA dinucleotide dehybridization using steady-state and temperature-jump infrared spectroscopy. Our results suggest that dinucleotide binding is stabilized through coaxial stacking interactions with the adjacent duplex segments as well as from potential noncanonical base-pairing configurations and structural dynamics of gap and overhang templates revealed using molecular dynamics simulations. We measure timescales for dissociation ranging from 0.2-40 μs depending on the template and temperature. Dinucleotide hybridization and dehybridization involve a significant free energy barrier with characteristics resembling that of canonical oligonucleotides. Together, our work provides an initial step for predicting the stability and kinetics of hybridization between short nucleic acid segments and various templates.

Enhanced nonenzymatic RNA copying with in-situ activation of short oligonucleotides

The nonenzymatic copying of RNA is thought to have been necessary for the transition between prebiotic chemistry and ribozyme-catalyzed RNA replication in the RNA World. We have previously shown that a potentially prebiotic nucleotide activation pathway based on phospho-Passerini chemistry can lead to the efficient synthesis of 2-aminoimidazole activated mononucleotides when carried out under freeze-thaw cycling conditions. Such activated nucleotides react with each other to form 5'-5' 2-aminoimidazolium bridged dinucleotides, enabling template-directed primer extension to occur within the same reaction mixture. However, mononucleotides linked to oligonucleotides by a 5'-5' 2-aminoimidazolium bridge are superior substrates for nonenzymatic primer extension; their higher intrinsic reactivity and their higher template affinity enable faster template copying at lower substrate concentrations. Here we show that eutectic phase phospho-Passerini chemistry efficiently activates short oligonucleotides and promotes the formation of monomer-bridged-oligonucleotide species during freeze-thaw cycles. We then demonstrate that in-situ generated monomer-bridged-oligonucleotides lead to efficient nonenzymatic template copying in the same reaction mixture. Our demonstration that multiple steps in the pathway from activation chemistry to RNA copying can occur together in a single complex environment simplifies this aspect of the origin of life.

Direct monitoring of the thermodynamics and kinetics of DNA and RNA dinucleotide dehybridization from gaps and overhangs

Hybridization of short nucleic acid segments (<4 nucleotides) to single-strand templates occurs as a critical intermediate in processes such as non-enzymatic nucleic acid replication and toehold-mediated strand displacement. These templates often contain adjacent duplex segments that stabilize base pairing with single-strand gaps or overhangs, but the thermodynamics and kinetics of hybridization in such contexts are poorly understood due to experimental challenges of probing weak binding and rapid structural dynamics. Here we develop an approach to directly measure the thermodynamics and kinetics of DNA and RNA dinucleotide dehybridization using steady-state and temperature-jump infrared spectroscopy. Our results suggest that dinucleotide binding is stabilized through coaxial stacking interactions with the adjacent duplex segments as well as from potential non-canonical base pairing configurations and structural dynamics of gap and overhang templates revealed using molecular dynamics simulations. We measure timescales for dissociation ranging from 0.2 to 40 µs depending on the template and temperature. Dinucleotide hybridization and dehybridization involves a significant free energy barrier with characteristics resembling that of canonical oligonucleotides. Together, our work provides an initial step for predicting the stability and kinetics of hybridization between short nucleic acid segments and various templates.

Experimental Tests of the Virtual Circular Genome Model for Nonenzymatic RNA Replication

The virtual circular genome (VCG) model was proposed as a means of going beyond template copying to indefinite cycles of nonenzymatic RNA replication during the origin of life. In the VCG model, the protocellular genome is a collection of short oligonucleotides that map to both strands of a virtual circular sequence. Replication is driven by templated nonenzymatic primer extensions on a subset of kinetically trapped partially base-paired configurations, followed by the shuffling of these configurations to enable continued oligonucleotide elongation. Here, we describe initial experimental studies of the feasibility of the VCG model for replication. We designed a small 12-nucleotide model VCG and synthesized all 247 oligonucleotides of lengths 2 to 12 corresponding to this genome. We experimentally monitored the fate of individual labeled primers in the pool of VCG oligonucleotides following the addition of activated nucleotides and investigated the effect of factors such as oligonucleotide length, concentration, composition, and temperature on the extent of primer extension. We observe a surprisingly prolonged equilibration process in the VCG system that enables a considerable extent of reaction. We find that environmental fluctuations would be essential for continuous templated extension of the entire VCG system since the shortest oligonucleotides can only bind to templates at low temperatures, while the longest oligonucleotides require high-temperature spikes to escape from inactive configurations. Finally, we demonstrate that primer extension is significantly enhanced when the mix of VCG oligonucleotides is preactivated. We discuss the necessity of ongoing in situ activation chemistry for continuous and accurate VCG replication.

Hybridization kinetics of out-of-equilibrium mixtures of short RNA oligonucleotides

Hybridization and strand displacement kinetics determine the evolution of the base paired configurations of mixtures of oligonucleotides over time. Although much attention has been focused on the thermodynamics of DNA and RNA base pairing in the scientific literature, much less work has been done on the time dependence of interactions involving multiple strands, especially in RNA. Here we provide a study of oligoribonucleotide interaction kinetics and show that it is possible to calculate the association, dissociation and strand displacement rates displayed by short oligonucleotides (5nt–12nt) that exhibit no expected secondary structure as simple functions of oligonucleotide length, CG content, ΔG of hybridization and ΔG of toehold binding. We then show that the resultant calculated kinetic parameters are consistent with the experimentally observed time dependent changes in concentrations of the different species present in mixtures of multiple competing RNA strands. We show that by changing the mixture composition, it is possible to create and tune kinetic traps that extend by orders of magnitude the typical sub-second hybridization timescale of two complementary oligonucleotides. We suggest that the slow equilibration of complex oligonucleotide mixtures may have facilitated the nonenzymatic replication of RNA during the origin of life.

Rolling Circles as a Means of Encoding Genes in the RNA World

The rolling circle mechanism found in viroids and some RNA viruses is a likely way that replication could have begun in the RNA World. Here, we consider simulations of populations of protocells, each containing multiple copies of rolling circle RNAs that can replicate non-enzymatically. The mechanism requires the presence of short self-cleaving ribozymes such as hammerheads, which can cleave and re-circularize RNA strands. A rolling circle must encode a hammerhead and the complement of a hammerhead, so that both plus and minus strands can cleave. Thus, the minimal functional length is twice the length of the hammerhead sequence. Selection for speed of replication will tend to reduce circles to this minimum length. However, if sequence errors occur when copying the hammerhead sequence, this prevents cleavage at one point, but still allows cleavage on the next passage around the rolling circle. Thus, there is a natural doubling mechanism that creates strands that are multiple times the length of the minimal sequence. This can provide space for the origin of new genes with beneficial functions. We show that if a beneficial gene appears in this new space, the longer sequence with the beneficial function can be selected, even though it replicates more slowly. This provides a route for the evolution of longer circles encoding multiple genes.

Computer simulations of Template-Directed RNA Synthesis driven by temperature cycling in diverse sequence mixtures

We present simulations of non-enzymatic template-directed RNA synthesis that incorporate primer extension, ligation, melting, and reannealing. Strand growth occurs over multiple heating/cooling cycles, producing strands of several hundred nucleotides in length, starting with random oligomers of 4 to 10 nucleotides. A strand typically grows by only 1 or 2 nucleotides in each cycle. Therefore, a strand is copied from many different templates, not from one specific complementary strand. A diverse sequence mixture is produced, and there is no exact copying of sequences, even if single base additions are fully accurate (no mutational errors). It has been proposed that RNA systems may contain a virtual circular genome, in which sequences partially overlap in a way that is mutually catalytic. We show that virtual circles do not emerge naturally in our simulations, and that a system initiated with a virtual circle can only maintain itself if there are no mutational errors and there is no input of new sequences formed by random polymerization. Furthermore, if a virtual sequence and its complement contain repeated short words, new sequences can be produced that were not on the original virtual circle. Therefore the virtual circle sequence cannot maintain itself. Functional sequences with secondary structures contain complementary words on opposite sides of stem regions. Both these words are repeated in the complementary sequence; hence, functional sequences cannot be encoded on a virtual circle. Additionally, we consider sequence replication in populations of protocells. We suppose that functional ribozymes benefit the cell which contains them. Nevertheless, scrambling of sequences occurs, and the functional sequence is not maintained, even when under positive selection.

The earliest form of RNA replication may have been non-enzymatic, without requiring polymerase ribozymes. Non-enzymatic replication forms double strands that are unlikely to separate unless melting is driven by temperature cycling. However, re-annealing of existing strands occurs rapidly on cooling, and this prevents subsequent cycles of copying if there are multiple copies of similar sequences. In contrast, if there is a diverse mixture of sequences, partially matching sequences can reanneal in configurations that allow continued strand growth. We show that this allows continued synthesis of populations of random sequences that are quite long. We test the idea that a virtual circular genome could exist in such a mixture. We show that a virtual genome does not arise spontaneously and that it cannot be maintained except in unrealistic ideal cases. We conclude that functional sequence information cannot be encoded on the fragments of a virtual circle.

The Effect of Environment on the Evolution and Proliferation of Protocells of Increasing Complexity

The formation, growth, division and proliferation of protocells containing RNA strands is an important step in ensuring the viability of a mixed RNA-lipid world. Experiments and computer simulations indicate that RNA encapsulated inside protocells can favor the protocell, promoting its growth while protecting the system from being over-run by selfish RNA sequences. Recent work has also shown that the rolling-circle replication mechanism can be harnessed to ensure the rapid growth of RNA strands and the probabilistic emergence and proliferation of protocells with functionally diverse ribozymes. Despite these advances in our understanding of a primordial RNA-lipid world, key questions remain about the ideal environment for the formation of protocells and its role in regulating the proliferation of functionally complex protocells. The hot spring hypothesis suggests that mineral-rich regions near hot springs, subject to dry-wet cycles, provide an ideal environment for the origin of primitive protocells. We develop a computational model to study protocellular evolution in such environments that are distinguished by the occurrence of three distinct phases, a wet phase, followed by a gel phase, and subsequently by a dry phase. We determine the conditions under which protocells containing multiple types of ribozymes can evolve and proliferate in such regions. We find that diffusion in the gel phase can inhibit the proliferation of complex protocells with the extent of inhibition being most significant when a small fraction of protocells is eliminated during environmental cycling. Our work clarifies how the environment can shape the evolution and proliferation of complex protocells.

Template-Free Assembly of Functional RNAs by Loop-Closing Ligation

The first ribozymes are thought to have emerged at a time when RNA replication proceeded via nonenzymatic template copying processes. However, functional RNAs have stable folded structures, and such structures are much more difficult to copy than short unstructured RNAs. How can these conflicting requirements be reconciled? Also, how can the inhibition of ribozyme function by complementary template strands be avoided or minimized? Here, we show that short RNA duplexes with single-stranded overhangs can be converted into RNA stem loops by nonenzymatic cross-strand ligation. We then show that loop-closing ligation reactions enable the assembly of full-length functional ribozymes without any external template. Thus, one can envisage a potential pathway whereby structurally complex functional RNAs could have formed at an early stage of evolution when protocell genomes might have consisted only of collections of short replicating oligonucleotides.

Visualizing RNA polymers produced by hot wet-dry cycling

It is possible that the transition from abiotic systems to life relied on RNA polymers that served as ribozyme-like catalysts and for storing genetic information. The source of such polymers is uncertain, but previous investigations reported that wet-dry cycles simulating prebiotic hot springs provide sufficient energy to drive condensation reactions of mononucleotides to form oligomers and polymers. The aim of the study reported here was to verify this claim and visualize the products prepared from solutions composed of single mononucleotides and 1:1 mixture of two mononucleotides. Therefore, we designed experiments that allowed comparisons of all such mixtures representing six combinations of the four mononucleotides of RNA. We observed irregular stringy patches and crystal strands when wet-dry cycling was performed at room temperature (20 °C). However, when the same solutions were exposed to wet-dry cycles at 80 °C, we observed what appeared to be true polymers. Their thickness was consistent with RNA-like products composed of covalently bonded monomers, while irregular strings and crystal segments of mononucleotides dried or cycled at room temperature were consistent with structures assembled and stabilized by weak hydrogen bonds. In a few instances we observed rings with short polymer attachments. These observations are consistent with previous claims of polymerization during wet-dry cycling. We conclude that RNA-like polymers and rings could have been synthesized non-enzymatically in freshwater hot springs on the prebiotic Earth with sizes sufficient to fold into ribozymes and genetic molecules required for life to begin.

Thermodynamic and Kinetic Sequence Selection in Enzyme-Free Polymer Self-Assembly Inside a Non-Equilibrium RNA Reactor

The RNA world is one of the principal hypotheses to explain the emergence of living systems on the prebiotic Earth. It posits that RNA oligonucleotides acted as both carriers of information as well as catalytic molecules, promoting their own replication. However, it does not explain the origin of the catalytic RNA molecules. How could the transition from a pre-RNA to an RNA world occur? A starting point to answer this question is to analyze the dynamics in sequence space on the lowest level, where mononucleotide and short oligonucleotides come together and collectively evolve into larger molecules. To this end, we study the sequence-dependent self-assembly of polymers from a random initial pool of short building blocks via templated ligation. Templated ligation requires two strands that are hybridized adjacently on a third strand. The thermodynamic stability of such a configuration crucially depends on the sequence context and, therefore, significantly influences the ligation probability. However, the sequence context also has a kinetic effect, since non-complementary nucleotide pairs in the vicinity of the ligation site stall the ligation reaction. These sequence-dependent thermodynamic and kinetic effects are explicitly included in our stochastic model. Using this model, we investigate the system-level dynamics inside a non-equilibrium ‘RNA reactor’ enabling a fast chemical activation of the termini of interacting oligomers. Moreover, the RNA reactor subjects the oligomer pool to periodic temperature changes inducing the reshuffling of the system. The binding stability of strands typically grows with the number of complementary nucleotides forming the hybridization site. While shorter strands unbind spontaneously during the cold phase, larger complexes only disassemble during the temperature peaks. Inside the RNA reactor, strand growth is balanced by cleavage via hydrolysis, such that the oligomer pool eventually reaches a non-equilibrium stationary state characterized by its length and sequence distribution. How do motif-dependent energy and stalling parameters affect the sequence composition of the pool of long strands? As a critical factor for self-enhancing sequence selection, we identify kinetic stalling due to non-complementary base pairs at the ligation site. Kinetic stalling enables cascades of self-amplification that result in a strong reduction of occupied states in sequence space. Moreover, we discuss the significance of the symmetry breaking for the transition from a pre-RNA to an RNA world.

Kinetic explanations for the sequence biases observed in the nonenzymatic copying of RNA templates

The identification of nonenzymatic pathways for nucleic acid replication is a key challenge in understanding the origin of life. We have previously shown that nonenzymatic RNA primer extension using 2-aminoimidazole (2AI) activated nucleotides occurs primarily through an imidazolium-bridged dinucleotide intermediate. The reactive nature and preorganized structure of the intermediate increase the efficiency of primer extension but remain insufficient to drive extensive copying of RNA templates containing all four canonical nucleotides. To understand the factors that limit RNA copying, we synthesized all ten 2AI-bridged dinucleotide intermediates and measured the kinetics of primer extension in a model system. The affinities of the ten dinucleotides for the primer/template/helper complexes vary by over 7,000-fold, consistent with nearest neighbor energetic predictions. Surprisingly, the reaction rates at saturating intermediate concentrations still vary by over 15-fold, with the most weakly binding dinucleotides exhibiting a lower maximal reaction rate. Certain noncanonical nucleotides can decrease sequence dependent differences in affinity and primer extension rate, while monomers bridged to short oligonucleotides exhibit enhanced binding and reaction rates. We suggest that more uniform binding and reactivity of imidazolium-bridged intermediates may lead to the ability to copy arbitrary template sequences under prebiotically plausible conditions.

Rolling-circle and strand-displacement mechanisms for non-enzymatic RNA replication at the time of the origin of life

It is likely that RNA replication began non-enzymatically, and that polymerases were later selected to speed up the process. We consider replication mechanisms in modern viruses and ask which of these is possible non-enzymatically, using mathematical models and experimental data found in the literature to estimate rates of RNA synthesis and replication. Replication via alternating plus and minus strands is found in some single-stranded RNA viruses. However, if this occurred non-enzymatically it would lead to double-stranded RNA that would not separate. With some form of environmental cycling, such as temperature, salinity, or pH cycling, double-stranded RNA can be melted to form single-stranded RNA, although re-annealing of existing strands would then occur much faster than synthesis of new strands. We show that re-annealing blocks this form of replication at a very low concentration of strands. Other kinds of viruses synthesize linear double strands from single strands and then make new single strands from double strands via strand-displacement. This does not require environmental cycling and is not blocked by re-annealing. However, under non-enzymatic conditions, if strand-displacement occurs from a linear template, we expect the incomplete new strand to be almost always displaced by the tail end of the old strand through toehold-mediated displacement. A third kind of replication in viruses and viroids is rolling-circle replication which occurs via strand-displacement on a circular template. Rolling-circle replication does not require environmental cycling and is not prevented by toehold-mediated displacement. Rolling-circle replication is therefore expected to occur non-enzymatically and is a likely starting point for the evolution of polymerase-catalysed replication.

From RNA World to SARS-CoV-2: The Edited Story of RNA Viral Evolution

The current SARS-CoV-2 pandemic underscores the importance of understanding the evolution of RNA genomes. While RNA is subject to the formation of similar lesions as DNA, the evolutionary and physiological impacts RNA lesions have on viral genomes are yet to be characterized. Lesions that may drive the evolution of RNA genomes can induce breaks that are repaired by recombination or can cause base substitution mutagenesis, also known as base editing. Over the past decade or so, base editing mutagenesis of DNA genomes has been subject to many studies, revealing that exposure of ssDNA is subject to hypermutation that is involved in the etiology of cancer. However, base editing of RNA genomes has not been studied to the same extent. Recently hypermutation of single-stranded RNA viral genomes have also been documented though its role in evolution and population dynamics. Here, we will summarize the current knowledge of key mechanisms and causes of RNA genome instability covering areas from the RNA world theory to the SARS-CoV-2 pandemic of today. We will also highlight the key questions that remain as it pertains to RNA genome instability, mutations accumulation, and experimental strategies for addressing these questions.

Competition between bridged dinucleotides and activated mononucleotides determines the error frequency of nonenzymatic RNA primer extension

Nonenzymatic copying of RNA templates with activated nucleotides is a useful model for studying the emergence of heredity at the origin of life. Previous experiments with defined-sequence templates have pointed to the poor fidelity of primer extension as a major problem. Here we examine the origin of mismatches during primer extension on random templates in the simultaneous presence of all four 2-aminoimidazole-activated nucleotides. Using a deep sequencing approach that reports on millions of individual template-product pairs, we are able to examine correct and incorrect polymerization as a function of sequence context. We have previously shown that the predominant pathway for primer extension involves reaction with imidazolium-bridged dinucleotides, which form spontaneously by the reaction of two mononucleotides with each other. We now show that the sequences of correctly paired products reveal patterns that are expected from the bridged dinucleotide mechanism, whereas those associated with mismatches are consistent with direct reaction of the primer with activated mononucleotides. Increasing the ratio of bridged dinucleotides to activated mononucleotides, either by using purified components or by using isocyanide-based activation chemistry, reduces the error frequency. Our results point to testable strategies for the accurate nonenzymatic copying of arbitrary RNA sequences.