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

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.

RNA Complexes with Nicks and Gaps: Thermodynamic and Kinetic Effects of Coaxial Stacking and Dangling Ends

Multiple RNA strands can interact in solution and assume a large variety of configurations dictated by their potential for base pairing. Although duplex formation from two complementary oligonucleotides has been studied in detail, we still lack a systematic characterization of the behavior of higher order complexes. Here, we focus on the thermodynamic and kinetic effects of an upstream oligonucleotide on the binding of a downstream oligonucleotide to a common template, as we vary the sequence and structure of the contact interface. We show that coaxial stacking in RNA is well correlated with but much more stabilizing than helix propagation over an analogous intact double helix step (median ΔΔG°37 °C ≈ 1.7 kcal/mol). Consequently, approximating coaxial stacking in RNA with the helix propagation term leads to large discrepancies between predictions and our experimentally determined melting temperatures, with an offset of ≈10 °C. Our kinetic study reveals that the hybridization of the downstream probe oligonucleotide is impaired (lower kon) by the presence of the upstream oligonucleotide, with the thermodynamic stabilization coming entirely from an extended lifetime (lower koff) of the bound downstream oligonucleotide, which can increase from seconds to months. Surprisingly, we show that the effect of nicks is dependent on the length of the stacking oligonucleotides, and we discuss the binding of ultrashort (1-4 nt) oligonucleotides that are relevant in the context of the origin of life. The thermodynamic and kinetic data obtained in this work allow for the prediction of the formation and stability of higher-order multistranded complexes.

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.

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.