SignificanceThe columnar liquid crystal phases reported here are physical associations of the smallest molecular species to self-assemble into the duplex base-paired stacked columnar double-helical structures of nucleic acids. These assemblies of monomers can provide starting states capable of partitioning appropriate molecules from solution with a high degree of selectivity, acting as pathways for the prebiotic appearance of molecular selection, self-assembly, and, ultimately, of the sequence-directed assembly of polymers.

Self-synthesizing materials, in which supramolecular structuring enhances the formation of new molecules that participate to the process, represent an intriguing notion to account for the first appearance of biomolecules in an abiotic Earth. We present here a study of the abiotic formation of inter-chain phosphodiester bonds in solutions of short RNA oligomers in various states of supramolecular arrangement and their reaction kinetics. We found a spectrum of conditions in which RNA oligomers self-assemble and phase separate into highly concentrated ordered fluid liquid crystal (LC) microdomains. We show that such supramolecular state provides a template guiding their ligation into hundred-bases long chains. The quantitative analysis presented here demonstrates that nucleic acid LC boosts the rate of end-to-end ligation and suppresses the formation of the otherwise dominant cyclic oligomers. These results strengthen the concept of supramolecular ordering as an efficient pathway toward the emergence of the RNA World in the primordial Earth.

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.