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.

Equilibrium and non-equilibrium furanose selection in the ribose isomerisation network

The exclusive presence of β-D-ribofuranose in nucleic acids is still a conundrum in prebiotic chemistry, given that pyranose species are substantially more stable at equilibrium. However, a precise characterisation of the relative furanose/pyranose fraction at temperatures higher than about 50 °C is still lacking. Here, we employ a combination of NMR measurements and statistical mechanics modelling to predict a population inversion between furanose and pyranose at equilibrium at high temperatures. More importantly, we show that a steady temperature gradient may steer an open isomerisation network into a non-equilibrium steady state where furanose is boosted beyond the limits set by equilibrium thermodynamics. Moreover, we demonstrate that nonequilibrium selection of furanose is maximum at optimal dissipation, as gauged by the temperature gradient and energy barriers for isomerisation. The predicted optimum is compatible with temperature drops found in hydrothermal vents associated with extremely fresh lava flows on the seafloor.

Anomalous thermal fluctuation distribution sustains proto-metabolic cycles and biomolecule synthesis

An environment far from equilibrium is thought to be a necessary condition for the origin and persistence of life. In this context we report open-flow simulations of a non-enzymic proto-metabolic system, in which hydrogen peroxide acts both as oxidant and driver of thermochemical cycling. We find that a Gaussian perturbed input produces a non-Boltzmann output fluctuation distribution around the mean oscillation maximum. Our main result is that net biosynthesis can occur under fluctuating cyclical but not steady drive. Consequently we may revise the necessary condition to "dynamically far from equilibrium".

Stability of 2',3' and 3',5' cyclic nucleotides in formamide and in water: a theoretical insight into the factors controlling the accumulation of nucleic acid building blocks in a prebiotic pool

Synthesis of the first RNAs represents one of the cornerstones of the emergence of life. Recent studies demonstrated powerful scenarios of prebiotic synthesis of cyclic nucleotides in aqueous and formamide environments. This raised a question about their thermodynamic stability, a decisive factor determining their accumulation in a prebiotic pool. Here we performed ab initio molecular dynamics simulations at various temperatures in formamide and water to study the relative stabilities of the 2',3' and 3',5' isomers of cyclic nucleotides. The computations show that in an aqueous environment 2',3' cyclic nucleotides are more stable than their 3',5' counterparts at all temperatures up to the boiling point. In contrast, in formamide higher temperatures favor the accumulation of the 3',5' cyclic form, whereas below about 400 K the 2',3' cyclic form becomes more stable. The latter observation is consistent with a formamide-based origin scenario, suggesting that 3',5' cyclic nucleotides accumulated at higher temperatures subsequently allowed oligomerization reactions after fast cooling to lower temperatures. A statistical analysis of the geometrical parameters of the solutes indicates that thermodynamics of cyclic nucleotides in aqueous and formamide environments are dictated by the floppiness of the molecules rather than by the ring strain of the cyclic phosphodiester linkages.