Genotypes with phenotypes: adventures in an RNA toy world

Evolution has created the complexity of the animate world and deciphering the language of evolution is the key towards understanding nature. The dynamics of evolution is simplified by considering it as a superposition of three less sophisticated processes: population dynamics, population support dynamics, and genotype-phenotype mapping. Evolution of molecules in laboratory assays provides a sufficiently simple system for the quantitative analysis of the three phenomena. Coarse-grained notions of structures like RNA secondary structures are used as model phenotypes. They provide an excellent tool for a comprehensive analysis of the entire complex of molecular evolution. The mapping from RNA genotypes into secondary structures is highly redundant. In order to find at least one sequence for every common structures one need only search a (relatively) small part of sequence space. The existence of selectively neutral phenotypes plays an important role for the the success and the efficiency of evolutionary optimization. Molecular evolution found a highly promising technological application in the design of biomolecules with predefined properties.     

The Role of Complex Formation and Deleterious Mutations for the Stability of RNA-Like Replicator Systems

In the RNA world hypothesis, RNA(-like) self-replicators are suggested as the central player of prebiotic evolution. However, there is a serious problem in the evolution of complexity in such replicators, i.e., the problem of parasites. Parasites, which are replicated by catalytic replicators (catalysts), but do not replicate the others, can destroy a whole replicator system by exploitation. Recently, a theoretical study underlined complex formation between replicators--an often neglected but realistic process--as a stabilizing factor in a replicator system by demonstrating that complex formation can shift the viable range of diffusion intensity to higher values. In the current study, we extend the previous study of complex formation. Firstly, by investigating a well-mixed replicator system, we establish that complex formation gives parasites an implicit advantage over catalysts, which makes the system significantly more vulnerable to parasites. Secondly, by investigating a spatially extended replicator system, we show that the formation of traveling wave patterns plays a crucial role in the stability of the system against parasites, and that because of this the effect of complex formation is not straightforward; i.e., whether complex formation stabilizes or destabilizes the spatial system is a complex function of other parameters. We give a detailed analysis of the spatial system by considering the pattern dynamics of waves. Furthermore, we investigate the effect of deleterious mutations. Surprisingly, high mutation rates can weaken the exploitation of the catalyst by the parasite.     

Predicted functional RNAs within coding regions constrain evolutionary rates of yeast proteins

Functional RNAs (fRNAs) are being recognized as an important regulatory component in biological processes. Interestingly, recent computational studies suggest that the number and biological significance of functional RNAs within coding regions (coding fRNAs) may have been underestimated. We hypothesized that such coding fRNAs will impose additional constraint on sequence evolution because the DNA primary sequence has to simultaneously code for functional RNA secondary structures on the messenger RNA in addition to the amino acid codons for the protein sequence. To test this prediction, we first utilized computational methods to predict conserved fRNA secondary structures within multiple species alignments of Saccharomyces sensu strico genomes. We predict that as much as 5% of the genes in the yeast genome contain at least one functional RNA secondary structure within their protein-coding region. We then analyzed the impact of coding fRNAs on the evolutionary rate of protein-coding genes because a decrease in evolutionary rate implies constraint due to biological functionality. We found that our predicted coding fRNAs have a significant influence on evolutionary rates (especially at synonymous sites), independent of other functional measures. Thus, coding fRNA may play a role on sequence evolution. Given that coding regions of humans and flies contain many more predicted coding fRNAs than yeast, the impact of coding fRNAs on sequence evolution may be substantial in genomes of higher eukaryotes.     

A Fission-Fusion Origin for Life

To develop a comprehensive cells-first approach to the origin of life, we propose that protocells form spontaneously and that the fission and fusion of these protocells drives the dynamics of their evolution. The fitness criterion for this evolution is taken to be the the stability (conservation) of domains in the protocellular membrane as determined by non-covalent molecular associations between the amphiphiles of the membrane and a subset of the macromolecules in the protocell. In the presence of a source of free energy the macromolecular content of the protocell (co-)evolves as the result of (domain-dependent) membrane-catalysed polymerisation of the prebiotic constituents delivered to the protocell by fusion. The metabolism of the cell therefore (co-)evolves on a rugged fitness landscape. We indicate how domain evolution with the same fitness criterion can potentially give rise to coding. Membrane domains may therefore provide the link between protocells and the RNA/DNA-world.     

Frequency of RNA-RNA interaction in a model of the RNA World

The RNA World model for prebiotic evolution posits the selection of catalytic/template RNAs from random populations. The mechanisms by which these random populations could be generated de novo are unclear. Non-enzymatic and RNA-catalyzed nucleic acid polymerizations are poorly processive, which means that the resulting short-chain RNA population could contain only limited diversity. Nonreciprocal recombination of smaller RNAs provides an alternative mechanism for the assembly of larger species with concomitantly greater structural diversity; however, the frequency of any specific recombination event in a random RNA population is limited by the low probability of an encounter between any two given molecules. This low probability could be overcome if the molecules capable of productive recombination were redundant, with many nonhomologous but functionally equivalent RNAs being present in a random population. Here we report fluctuation experiments to estimate the redundancy of the set of RNAs in a population of random sequences that are capable of non-Watson-Crick interaction with another RNA. Parallel SELEX experiments showed that at least one in 10(6) random 20-mers binds to the P5.1 stem-loop of Bacillus subtilis RNase P RNA with affinities equal to that of its naturally occurring partner. This high frequency predicts that a single RNA in an RNA World would encounter multiple interacting RNAs within its lifetime, supporting recombination as a plausible mechanism for prebiotic RNA evolution. The large number of equivalent species implies that the selection of any single interacting species in the RNA World would be a contingent event, i.e., one resulting from historical accident.     

On the structural repertoire of pools of short, random RNA sequences

A detailed knowledge of the mapping between sequence and structure spaces in populations of RNA molecules is essential to better understand their present-day functional properties, to envisage a plausible early evolution of RNA in a prebiotic chemical environment and to improve the design of in vitro evolution experiments, among others. Analysis of natural RNAs, as well as in vitro and computational studies, show that certain RNA structural motifs are much more abundant than others, pointing out a complex relation between sequence and structure. Within this framework, we have investigated computationally the structural properties of a large pool (10⁸ molecules) of single-stranded, 35 nt-long, random RNA sequences. The secondary structures obtained are ranked and classified into structure families. The number of structures in main families is analytically calculated and compared with the numerical results. This permits a quantification of the fraction of structure space covered by a large pool of sequences. We further show that the number of structural motifs and their frequency is highly unbalanced with respect to the nucleotide composition: simple structures such as stem-loops and hairpins arise from sequences depleted in G, while more complex structures require an enrichment of G. In general, we observe a strong correlation between subfamilies—characterized by a fixed number of paired nucleotides—and nucleotide composition. Our results are compared to the structural repertoire obtained in a second pool where isolated base pairs are prohibited.     

On the origin of DNA genomes: evolution of the division of labor between template and catalyst in model replicator systems

The division of labor between template and catalyst is a fundamental property of all living systems: DNA stores genetic information whereas proteins function as catalysts. The RNA world hypothesis, however, posits that, at the earlier stages of evolution, RNA acted as both template and catalyst. Why would such division of labor evolve in the RNA world? We investigated the evolution of DNA-like molecules, i.e. molecules that can function only as template, in minimal computational models of RNA replicator systems. In the models, RNA can function as both template-directed polymerase and template, whereas DNA can function only as template. Two classes of models were explored. In the surface models, replicators are attached to surfaces with finite diffusion. In the compartment models, replicators are compartmentalized by vesicle-like boundaries. Both models displayed the evolution of DNA and the ensuing division of labor between templates and catalysts. In the surface model, DNA provides the advantage of greater resistance against parasitic templates. However, this advantage is at least partially offset by the disadvantage of slower multiplication due to the increased complexity of the replication cycle. In the compartment model, DNA can significantly delay the intra-compartment evolution of RNA towards catalytic deterioration. These results are explained in terms of the trade-off between template and catalyst that is inherent in RNA-only replication cycles: DNA releases RNA from this trade-off by making it unnecessary for RNA to serve as template and so rendering the system more resistant against evolving parasitism. Our analysis of these simple models suggests that the lack of catalytic activity in DNA by itself can generate a sufficient selective advantage for RNA replicator systems to produce DNA. Given the widespread notion that DNA evolved owing to its superior chemical properties as a template, this study offers a novel insight into the evolutionary origin of DNA.     

Spatial Models of Prebiotic Evolution: Soup Before Pizza?

The problem of information integration andresistance to the invasion of parasitic mutants in prebiotic replicator systemsis a notorious issue of research on the origin of life.Almost all theoretical studies published so far havedemonstrated that some kind of spatial structure is indispensable forthe persistence and/or the parasite resistance of any feasible replicator system.Based on a detailed critical survey of spatial models on prebiotic informationintegration, we suggest a possible scenario for replicator system evolution leadingto the emergence of the first protocells capable of independent life.We show that even the spatial versions of the hypercycle model are vulnerable toselfish parasites in heterogeneous habitats. Contrary, the metabolic system remainspersistent and coexistent with its parasites both on heterogeneous surfaces andin chaotically mixing flowing media. Persistent metabolic parasites can beconverted to metabolic cooperators, or they can gradually obtain replicase activity.Our simulations show that, once replicase activity emerged, a gradual and simultaneousevolutionary improvement of replicase functionality (speed and fidelity) andtemplate efficiency is possible only on a surface that constrains the mobility ofmacromolecule replicators. Based on the results of the models reviewed, we suggestthat open chaotic flows (`soup') and surface dynamics (`pizza') both played keyroles in the sequence of evolutionary events ultimately concluding in theappearance of the first living cell on Earth.     

4 Towards non-enzymatic RNA replication

The direct path from prebiotic chemistry to the RNA World requires a plausible route for the synthesis of activated ribonucleotides and RNA templates, along with a means for the complete replication of potentially useful RNA sequences. However, many apparent roadblocks make non-enzymatic RNA replication look quite difficult, if not impossible. These problems include the slow rate, low accuracy and poor regioselectivity of non-enzymatic template copying, the hydrolysis of activated monomers and absence of good re-activation chemistry, the difficulty of strand separation after template copying, the rapidity of strand reannealing, the absence of primers in any realistic replication scenario, and the apparent incompatibly of RNA copying chemistry (which requires a high Mg²⁺ concentration) with fatty acid-based protocell membranes, which are destroyed by low Mg²⁺ concentrations. I will discuss recent progress from my laboratory on four of these issues. We have found that functional RNAs such as aptamers and ribozymes can tolerate moderate levels of 2′–5′ linkages without great loss of activity. It therefore appears that the presence of 10–25% of such linkages in the products of non-enzymatic copying would not prevent the evolution of functional RNAs. Furthermore, 2′–5′ linkages can be helpful, as they decrease the melting temperature of RNA duplexes enough to allow strand separation to occur under geophysically plausible conditions. Recently, we have found that small chemical changes to the nucleobases can greatly increase the fidelity of non-enzymatic template copying, and we have found conditions that render RNA copying chemistry compatible with vesicle integrity, thereby allowing RNA copying to occur inside fatty acid-based model protocell membranes. I will discuss potential approaches to solving the remaining issues that stand in the way of complete RNA replication. If all of the problems with RNA replication can be overcome, it should be possible to construct functioning protocells in the laboratory.     

The integration of the earliest genetic information

New levels of evolutionary units have emerged a number of times. One major pathway of emergence is the integration of information dispersed in originally competitive lower-level units. Such a transition is thought to have occurred during the origin of life, leading from RNA replicators to protocells. Recent alternative ideas, and some of their connections with population biology, are reviewed.     

Ro's role in RNA reconnaissance

The Ro 60 kDa autoantigen binds misfolded RNAs and likely functions in small RNA quality control. In this issue of Cell, Stein et al. (2005) present crystal structures of Ro alone and bound to both double- and single-stranded RNA, revealing two distinct RNA binding sites that suggest how Ro may distinguish between native and misfolded small RNAs.     

Use of linear regression model to compare RNA secondary structures

With more and more ribonucleic acid (RNA) secondary structures accumulated, the need for comparing different RNA secondary structures often arises in function prediction and evolutionary analysis. Numerous efficient algorithms were developed for comparing different RNA secondary structures, but challenges remain. In this paper, six new models based on the linear regression model were proposed for the comparison of RNA secondary structures. The proposed models were tested on a mixed data, containing six secondary structures from RNase P RNAs, three secondary structures from SSU rRNA and five secondary structures from 16S ribosomal RNAs. The results have shown the effectiveness of the proposed models. Moreover, the time complexity of our models is favorable by comparing with that of the existing methods which solve the similar problem.     

Multilevel Selection in Models of Prebiotic Evolution II: A Direct Comparison of Compartmentalization and Spatial Self-Organization

Multilevel selection has been indicated as an essential factor for the evolution of complexity in interacting RNA-like replicator systems. There are two types of multilevel selection mechanisms: implicit and explicit. For implicit multilevel selection, spatial self-organization of replicator populations has been suggested, which leads to higher level selection among emergent mesoscopic spatial patterns (traveling waves). For explicit multilevel selection, compartmentalization of replicators by vesicles has been suggested, which leads to higher level evolutionary dynamics among explicitly imposed mesoscopic entities (protocells). Historically, these mechanisms have been given separate consideration for the interests on its own. Here, we make a direct comparison between spatial self-organization and compartmentalization in simulated RNA-like replicator systems. Firstly, we show that both mechanisms achieve the macroscopic stability of a replicator system through the evolutionary dynamics on mesoscopic entities that counteract that of microscopic entities. Secondly, we show that a striking difference exists between the two mechanisms regarding their possible influence on the long-term evolutionary dynamics, which happens under an emergent trade-off situation arising from the multilevel selection. The difference is explained in terms of the difference in the stability between self-organized mesoscopic entities and externally imposed mesoscopic entities. Thirdly, we show that a sharp transition happens in the long-term evolutionary dynamics of the compartmentalized system as a function of replicator mutation rate. Fourthly, the results imply that spatial self-organization can allow the evolution of stable folding in parasitic replicators without any specific functionality in the folding itself. Finally, the results are discussed in relation to the experimental synthesis of chemical Darwinian systems and to the multilevel selection theory of evolutionary biology in general. To conclude, novel evolutionary directions can emerge through interactions between the evolutionary dynamics on multiple levels of organization. Different multilevel selection mechanisms can produce a difference in the long-term evolutionary trend of identical microscopic entities.     

A new Motzkin class for joint RNA secondary structures

In general RNA prediction problem includes genetic mapping, physical mapping and structure prediction. The ultimate goal of structure prediction is to obtain the three dimensional structure of bimolecules through computation. The key concept for solving the above mentioned problem is the appropriate representation of the biological structures. Even though, the problems that concern representations of certain biological structures like secondary structures either are characterized as NP-complete or with high complexity, few approximation algorithms and techniques had been constructed, mainly with polynomial complexity, concerning the prediction of RNA secondary structures. In this paper, a new class of Motzkin paths is introduced, the so-called semi-elevated inverse Motzkin peakless paths for the representation of two interacting RNA molecules. The basic combinatorial interpretations on single RNA secondary structures are extended via these new Motzkin paths on two RNA molecules and can be applied to the prediction methods of joint structures formed by interacting RNAs.     

Phenotypic effect of mutations in evolving populations of RNA molecules

Background  

The secondary structure of folded RNA sequences is a good model to map phenotype onto genotype, as represented by the RNA sequence. Computational studies of the evolution of ensembles of RNA molecules towards target secondary structures yield valuable clues to the mechanisms behind adaptation of complex populations. The relationship between the space of sequences and structures, the organization of RNA ensembles at mutation-selection equilibrium, the time of adaptation as a function of the population parameters, the presence of collective effects in quasispecies, or the optimal mutation rates to promote adaptation all are issues that can be explored within this framework.     

Modular evolution and increase of functional complexity in replicating RNA molecules

At early stages of biochemical evolution, the complexity of replicating molecules was limited by unavoidably high mutation rates. In an RNA world, prior to the appearance of cellular life, an increase in molecular length, and thus in functional complexity, could have been mediated by modular evolution. We describe here a scenario in which short, replicating RNA sequences are selected to perform a simple function. Molecular function is represented through the secondary structure corresponding to each sequence, and a given target secondary structure yields the optimal function in the environment where the population evolves. The combination of independently evolved populations may have facilitated the emergence of larger molecules able to perform more complex functions (including RNA replication) that could arise as a combination of simpler ones. We quantitatively show that modular evolution has relevant advantages with respect to the direct evolution of large functional molecules, among them the allowance of higher mutation rates, the shortening of evolutionary times, and the very possibility of finding complex structures that could not be otherwise directly selected.