Natural and artificial RNAs occupy the same restricted region of sequence space

Replication elongates short DNA, reduces sequence bias and develops trimer structure

The origin of molecular evolution required the replication of short oligonucleotides to form longer polymers. Prebiotically plausible oligonucleotide pools tend to contain more of some nucleobases than others. It has been unclear whether this initial bias persists and how it affects replication. To investigate this, we examined the evolution of 12-mer biased short DNA pools using an enzymatic model system. This allowed us to study the long timescales involved in evolution, since it is not yet possible with currently investigated prebiotic replication chemistries. Our analysis using next-generation sequencing from different time points revealed that the initial nucleotide bias of the pool disappeared in the elongated pool after isothermal replication. In contrast, the nucleotide composition at each position in the elongated sequences remained biased and varied with both position and initial bias. Furthermore, we observed the emergence of highly periodic dimer and trimer motifs in the rapidly elongated sequences. This shift in nucleotide composition and the emergence of structure through templated replication could help explain how biased prebiotic pools could undergo molecular evolution and lead to complex functional nucleic acids.

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.

The virtual circular genome model for primordial RNA replication

We propose a model for the replication of primordial protocell genomes that builds upon recent advances in the nonenzymatic copying of RNA. We suggest that the original genomes consisted of collections of oligonucleotides beginning and ending at all possible positions on both strands of one or more virtual circular sequences. Replication is driven by feeding with activated monomers and by the activation of monomers and oligonucleotides in situ. A fraction of the annealed configurations of the protocellular oligonucleotides would allow for template-directed oligonucleotide growth by primer extension or ligation. Rearrangements of these annealed configurations, driven either by environmental fluctuations or occurring spontaneously, would allow for continued oligonucleotide elongation. Assuming that shorter oligonucleotides were more abundant than longer ones, replication of the entire genome could occur by the growth of all oligonucleotides by as little as one nucleotide on average. We consider possible scenarios that could have given rise to such protocell genomes, as well as potential routes to the emergence of catalytically active ribozymes and thus the more complex cells of the RNA World.

On the emergence of structural complexity in RNA replicators

The RNA world hypothesis relies on the ability of ribonucleic acids to spontaneously acquire complex structures capable of supporting essential biological functions. Multiple sophisticated evolutionary models have been proposed for their emergence, but they often assume specific conditions. In this work, we explore a simple and parsimonious scenario describing the emergence of complex molecular structures at the early stages of life. We show that at specific GC content regimes, an undirected replication model is sufficient to explain the apparition of multibranched RNA secondary structures-a structural signature of many essential ribozymes. We ran a large-scale computational study to map energetically stable structures on complete mutational networks of 50-nt-long RNA sequences. Our results reveal that the sequence landscape with stable structures is enriched with multibranched structures at a length scale coinciding with the appearance of complex structures in RNA databases. A random replication mechanism preserving a 50% GC content may suffice to explain a natural enrichment of stable complex structures in populations of functional RNAs. In contrast, an evolutionary mechanism eliciting the most stable folds at each generation appears to help reaching multibranched structures at highest GC content.

Type-II tRNAs and Evolution of Translation Systems and the Genetic Code

Because tRNA is the core biological intellectual property that was necessary to evolve translation systems, tRNAomes, ribosomes, aminoacyl-tRNA synthetases, and the genetic code, the evolution of tRNA is the core story in evolution of life on earth. We have previously described the evolution of type-I tRNAs. Here, we use the same model to describe the evolution of type-II tRNAs, with expanded V loops. The models are strongly supported by inspection of typical tRNA diagrams, measuring lengths of V loop expansions, and analyzing the homology of V loop sequences to tRNA acceptor stems. Models for tRNA evolution provide a pathway for the inanimate-to-animate transition and for the evolution of translation systems, the genetic code, and cellular life.

RNA regulators responding to ribosomal protein S15 are frequent in sequence space

There are several natural examples of distinct RNA structures that interact with the same ligand to regulate the expression of homologous genes in different organisms. One essential question regarding this phenomenon is whether such RNA regulators are the result of convergent or divergent evolution. Are the RNAs derived from some common ancestor and diverged to the point where we cannot identify the similarity, or have multiple solutions to the same biological problem arisen independently? A key variable in assessing these alternatives is how frequently such regulators arise within sequence space. Ribosomal protein S15 is autogenously regulated via an RNA regulator in many bacterial species; four apparently distinct regulators have been functionally validated in different bacterial phyla. Here, we explore how frequently such regulators arise within a partially randomized sequence population. We find many RNAs that interact specifically with ribosomal protein S15 from Geobacillus kaustophilus with biologically relevant dissociation constants. Furthermore, of the six sequences we characterize, four show regulatory activity in an Escherichia coli reporter assay. Subsequent footprinting and mutagenesis analysis indicates that protein binding proximal to regulatory features such as the Shine–Dalgarno sequence is sufficient to enable regulation, suggesting that regulation in response to S15 is relatively easily acquired.

The GC content of LSU rRNA evolves across topological and functional regions of the ribosome in all three domains of life

Large-subunit rRNA is the ribozyme that catalyzes protein synthesis by translation, and many of its features vary along a deep-to-superficial gradient. By measuring the G+C proportions in this rRNA in all three domains of life (60 bacteria, 379 eukaryote, and 23 archaean sequences), we tested whether the proportion of GC nucleotides varies along this in-out gradient. The rRNA regions used were several zones identified by Bokov and Steinberg (2009) as being arranged from deep to superficial within the LSU. To the Bokov-Steinberg zones, we added the most superficial zone of all, the divergent domains (expansion segments), which are greatly enlarged in eukaryotes. Regression lines constructed from the hundreds of species of organisms revealed the expected in-out gradient, showing that species with high %GC (or high %AT) in their rRNA distribute more of these abundant nucleotides into the peripheral zones. This could be explained by the evolutionary rates of replacement of all nucleotides (A, C, G, T), because these latter rates are fastest at the periphery and slowest near the conserved core. As an overall explanation, we propose that when extrinsic factors (whole-genome nucleotide composition, or environmental temperature) demand the percentage of GC in the rRNA of a species be high or low, then the deep-lying zones are buffered against GC variation because they are the slowest to evolve. The deep, conserved zones are also the most involved in translation, hinting that stabilizing selection there prevents a high GC variability that would diminish LSU rRNA's core functions. We found only a few domain-specific trends in rRNA-GC distribution, which relate to many Archaea living at high temperatures or to the highly complex genes and adaptations of Eukaryota. Use of rRNA sequences in molecular phylogenetic studies, for reconstructing the relationships of organisms across the tree of life, requires accurate models of how rRNA evolves. The demonstration that GC distributes in regular patterns across rRNA regions can improve these tree-reconstruction models in the future and should yield phylogenies of greater accuracy.

The paradox of dual roles in the RNA world: resolving the conflict between stable folding and templating ability

The hypothesized dual roles of RNA as both information carrier and biocatalyst during the earliest stages of life require a combination of features: good templating ability (for replication) and stable folding (for ribozymes). However, this poses the following paradox: well-folded sequences are poor templates for copying, but poorly folded sequences are unlikely to be good ribozymes. Here, we describe a strategy to overcome this dilemma through G:U wobble pairing in RNA. Unlike Watson-Crick base pairs, wobble pairs contribute highly to the energetic stability of the folded structure of their sequence, but only slightly, if at all, to the stability of the folded reverse complement. Sequences in the RNA World might thereby combine stable folding of the ribozyme with an unstructured, reverse-complementary genome, resulting in a "division of labor" between the strands. We demonstrate this strategy using computational simulations of RNA folding and an experimental model of early replication, nonenzymatic template-directed RNA primer extension. Additional study is needed to solve other problems associated with a complete replication cycle, including separation of strands after copying. Interestingly, viroid RNA sequences, which have been suggested to be relics of an RNA World (Diener, Proc Natl Acad Sci USA 86:9370-9374, 1989), also show significant asymmetry in folding energy between the infectious (+) and template (-) strands due to G:U pairing, suggesting that this strategy may even be used by replicators in the present day.

Prebiotically plausible mechanisms increase compositional diversity of nucleic acid sequences

During the origin of life, the biological information of nucleic acid polymers must have increased to encode functional molecules (the RNA world). Ribozymes tend to be compositionally unbiased, as is the vast majority of possible sequence space. However, ribonucleotides vary greatly in synthetic yield, reactivity and degradation rate, and their non-enzymatic polymerization results in compositionally biased sequences. While natural selection could lead to complex sequences, molecules with some activity are required to begin this process. Was the emergence of compositionally diverse sequences a matter of chance, or could prebiotically plausible reactions counter chemical biases to increase the probability of finding a ribozyme? Our in silico simulations using a two-letter alphabet show that template-directed ligation and high concatenation rates counter compositional bias and shift the pool toward longer sequences, permitting greater exploration of sequence space and stable folding. We verified experimentally that unbiased DNA sequences are more efficient templates for ligation, thus increasing the compositional diversity of the pool. Our work suggests that prebiotically plausible chemical mechanisms of nucleic acid polymerization and ligation could predispose toward a diverse pool of longer, potentially structured molecules. Such mechanisms could have set the stage for the appearance of functional activity very early in the emergence of life.

The prebiotic evolutionary advantage of transferring genetic information from RNA to DNA

In the early 'RNA world' stage of life, RNA stored genetic information and catalyzed chemical reactions. However, the RNA world eventually gave rise to the DNA-RNA-protein world, and this transition included the 'genetic takeover' of information storage by DNA. We investigated evolutionary advantages for using DNA as the genetic material. The error rate of replication imposes a fundamental limit on the amount of information that can be stored in the genome, as mutations degrade information. We compared misincorporation rates of RNA and DNA in experimental non-enzymatic polymerization and calculated the lowest possible error rates from a thermodynamic model. Both analyses found that RNA replication was intrinsically error-prone compared to DNA, suggesting that total genomic information could increase after the transition to DNA. Analysis of the transitional RNA/DNA hybrid duplexes showed that copying RNA into DNA had similar fidelity to RNA replication, so information could be maintained during the genetic takeover. However, copying DNA into RNA was very error-prone, suggesting that attempts to return to the RNA world would result in a considerable loss of information. Therefore, the genetic takeover may have been driven by a combination of increased chemical stability, increased genome size and irreversibility.

Protein folding absent selection

Biological proteins are known to fold into specific 3D conformations. However, the fundamental question has remained: Do they fold because they are biological, and evolution has selected sequences which fold? Or is folding a common trait, widespread throughout sequence space? To address this question arbitrary, unevolved, random-sequence proteins were examined for structural features found in folded, biological proteins. Libraries of long (71 residue), random-sequence polypeptides, with ensemble amino acid composition near the mean for natural globular proteins, were expressed as cleavable fusions with ubiquitin. The structural properties of both the purified pools and individual isolates were then probed using circular dichroism, fluorescence emission, and fluorescence quenching techniques. Despite this necessarily sparse "sampling" of sequence space, structural properties that define globular biological proteins, namely collapsed conformations, secondary structure, and cooperative unfolding, were found to be prevalent among unevolved sequences. Thus, for polypeptides the size of small proteins, natural selection is not necessary to account for the compact and cooperative folded states observed in nature.

Computational approaches to RNA structure prediction, analysis, and design

RNA molecules are important cellular components involved in many fundamental biological processes. Understanding the mechanisms behind their functions requires RNA tertiary structure knowledge. Although modeling approaches for the study of RNA structures and dynamics lag behind efforts in protein folding, much progress has been achieved in the past two years. Here, we review recent advances in RNA folding algorithms, RNA tertiary motif discovery, applications of graph theory approaches to RNA structure and function, and in silico generation of RNA sequence pools for aptamer design. Advances within each area can be combined to impact many problems in RNA structure and function.

Motif frequency and evolutionary search times in RNA populations

RNA molecules, through their dual identity as sequence and structure, are an appropriate experimental and theoretical model to study the genotype-phenotype map and evolutionary processes taking place in simple replicator populations. In this computational study, we relate properties of the sequence-structure map, in particular the abundance of a given secondary structure in a random pool, with the number of replicative events that an initially random population of sequences needs to find that structure through mutation and selection. For common structures, this search process turns out to be much faster than for rare structures. Furthermore, search and fixation processes are more efficient in a wider range of mutation rates for common structures, thus indicating that evolvability of RNA populations is not simply determined by abundance. We also find significant differences in the search and fixation processes for structures of same abundance, and relate them with the number of base pairs forming the structure. Moreover, the influence of the nucleotide content of the RNA sequences on the search process is studied. Our results advance in the understanding of the distribution and attainability of RNA secondary structures. They hint at the fact that, beyond sequence length and sequence-to-function redundancy, the mutation rate that permits localization and fixation of a given phenotype strongly depends on its relative abundance and global, in general non-uniform, distribution in sequence space.

Nucleotides that are essential but not conserved; a sufficient L-tryptophan site in RNA

Conservation is often used to define essential sequences within RNA sites. However, conservation finds only invariant sequence elements that are necessary for function, rather than finding a set of sequence elements sufficient for function. Biochemical studies in several systems-including the hammerhead ribozyme and the purine riboswitch-find additional elements, such as loop-loop interactions, required for function yet not phylogenetically conserved. Here we define a critical test of sufficiency: We embed a minimal, apparently sufficient motif for binding the amino acid tryptophan in a random-sequence background and ask whether we obtain functional molecules. After a negative result, we use a combination of three-dimensional structural modeling, selection, designed mutations, high-throughput sequencing, and bioinformatics to explore functional insufficiency. This reveals an essential unpaired G in a diverse structural context, varied sequence, and flexible distance from the invariant internal loop binding site identified previously. Addition of the new element yields a sufficient binding site by the insertion criterion, binding tryptophan in 22 out of 23 tries. Random insertion testing for site sufficiency seems likely to be broadly revealing.