Evolutionary origins and directed evolution of RNA

In Vitro Selection of Varkud Satellite Ribozyme Variants that Cleave a Modified Stem-Loop Substrate

In vitro selection is an established approach to create artificial ribozymes with defined activities or to modify the properties of naturally occurring ribozymes. For the Varkud satellite ribozyme of Neurospora, an in vitro selection protocol based on its phosphodiester bond cleavage activity has not been previously reported. Here, we describe a simple protocol for cleavage-based in vitro selection that we recently used to identify variants of the Varkud satellite ribozyme able to target and cleave a non-natural stem-loop substrate derived from the HIV-1 TAR RNA. It allows quick selection of active ribozyme variants from the transcription reaction based on the size of the self-cleavage product without the need for RNA labeling. This results in a streamlined procedure that is easily adaptable to engineer ribozymes with new activities.

Promiscuous Ribozymes and Their Proposed Role in Prebiotic Evolution

The ability of enzymes, including ribozymes, to catalyze side reactions is believed to be essential to the evolution of novel biochemical activities. It has been speculated that the earliest ribozymes, whose emergence marked the origin of life, were low in activity but high in promiscuity, and that these early ribozymes gave rise to specialized descendants with higher activity and specificity. Here, we review the concepts related to promiscuity and examine several cases of highly promiscuous ribozymes. We consider the evidence bearing on the question of whether de novo ribozymes would be quantitatively more promiscuous than later evolved ribozymes or protein enzymes. We suggest that while de novo ribozymes appear to be promiscuous in general, they are not obviously more promiscuous than more highly evolved or active sequences. Promiscuity is a trait whose value would depend on selective pressures, even during prebiotic evolution.

The Alanine World Model for the Development of the Amino Acid Repertoire in Protein Biosynthesis

A central question in the evolution of the modern translation machinery is the origin and chemical ethology of the amino acids prescribed by the genetic code. The RNA World hypothesis postulates that templated protein synthesis has emerged in the transition from RNA to the Protein World. The sequence of these events and principles behind the acquisition of amino acids to this process remain elusive. Here we describe a model for this process by following the scheme previously proposed by Hartman and Smith, which suggests gradual expansion of the coding space as GC–GCA–GCAU genetic code. We point out a correlation of this scheme with the hierarchy of the protein folding. The model follows the sequence of steps in the process of the amino acid recruitment and fits well with the co-evolution and coenzyme handle theories. While the starting set (GC-phase) was responsible for the nucleotide biosynthesis processes, in the second phase alanine-based amino acids (GCA-phase) were recruited from the core metabolism, thereby providing a standard secondary structure, the α-helix. In the final phase (GCAU-phase), the amino acids were appended to the already existing architecture, enabling tertiary fold and membrane interactions. The whole scheme indicates strongly that the choice for the alanine core was done at the GCA-phase, while glycine and proline remained rudiments from the GC-phase. We suggest that the Protein World should rather be considered the Alanine World, as it predominantly relies on the alanine as the core chemical scaffold.

Virus-Inspired Function in Engineered Protein Cages

The structural and functional diversity of proteins combined with their genetic programmability has made them indispensable modern materials. Well-defined, hollow protein capsules have proven to be particularly useful due to their ability to compartmentalize macromolecules and chemical processes. To this end, viral capsids are common scaffolds and have been successfully repurposed to produce a suite of practical protein-based nanotechnologies. Recently, the recapitulation of viromimetic function in protein cages of nonviral origin has emerged as a strategy to both complement physical studies of natural viruses and produce useful scaffolds for diverse applications. In this perspective, we review recent progress toward generation of virus-like behavior in nonviral protein cages through rational engineering and directed evolution. These artificial systems can aid our understanding of the emergence of viruses from existing cellular components, as well as provide alternative approaches to tackle current problems, and open up new opportunities, in medicine and biotechnology.

Altered stoichiometry of an evolved RNA aptamer

Inhibitory aptamers against a protein are promising as antagonistic reagents and repressive genetic components. Typically, improvement of such aptamers is achieved by acquiring higher binding affinity. Here, we report an alternative mechanism for the improvement of aptamer activity. Recently, we reported a transcriptional activator based on an inhibitory RNA aptamer against lambda cI repressor. We improved the aptamer through in vitro selection (SELEX) from a randomly mutagenized aptamer pool, followed by in vivo screening and truncation. Biochemical analyses indicated that the activity improvement was achieved by alteration of the complex formation stoichiometry, rather than by higher affinity or expression. Our results suggest an alternative strategy for improving aptamer activity.

An inhibitory RNA aptamer against the lambda cI repressor shows transcriptional activator activity in vivo

RNA aptamers are one of the promising components for constructing artificial genetic circuits. In this study, we developed a transcriptional activator based on an RNA aptamer against one of the most frequently applied repressor proteins, lambda phage cI. In vitro selection (Systematic Evolution of Ligands by EXponential enrichment) and following in vivo screening identified an RNA aptamer with the intended transcriptional activator activity from an RNA pool containing a 40-nucleotide long random region. Quantitative analysis showed a 35-fold elevation of reporter expression upon aptamer expression. These results suggest that the diversity of artificial transcriptional activators can be extended by employing RNA aptamers against repressor proteins to broaden the parts for constructing genetic circuits.

From Prebiotics to Probiotics: The Evolution and Functions of tRNA Modifications

All nucleic acids in cells are subject to post-transcriptional chemical modifications. These are catalyzed by a myriad of enzymes with exquisite specificity and that utilize an often-exotic array of chemical substrates. In no molecule are modifications more prevalent than in transfer RNAs. In the present document, we will attempt to take a chemical rollercoaster ride from prebiotic times to the present, with nucleoside modifications as key players and tRNA as the centerpiece that drove the evolution of biological systems to where we are today. These ideas will be put forth while touching on several examples of tRNA modification enzymes and their modus operandi in cells. In passing, we submit that the choice of tRNA is not a whimsical one but rather highlights its critical function as an essential invention for the evolution of protein enzymes.

DNA Catalysis: The Chemical Repertoire of DNAzymes

Deoxyribozymes or DNAzymes are single-stranded catalytic DNA molecules that are obtained by combinatorial in vitro selection methods. Initially conceived to function as gene silencing agents, the scope of DNAzymes has rapidly expanded into diverse fields, including biosensing, diagnostics, logic gate operations, and the development of novel synthetic and biological tools. In this review, an overview of all the different chemical reactions catalyzed by DNAzymes is given with an emphasis on RNA cleavage and the use of non-nucleosidic substrates. The use of modified nucleoside triphosphates (dN*TPs) to expand the chemical space to be explored in selection experiments and ultimately to generate DNAzymes with an expanded chemical repertoire is also highlighted.

The dynamics of the RNA world: insights and challenges

The RNA world hypothesis of the origin of life, in which RNA emerged as both enzyme and information carrier, is receiving solid experimental support. The prebiotic synthesis of biomolecules, the catalytic aid offered by mineral surfaces, and the vast enzymatic repertoire of ribozymes are only pieces of the origin of life puzzle; the full picture can only emerge if the pieces fit together by either following from one another or coexisting with each other. Here, we review the theory of the origin, maintenance, and enhancement of the RNA world as an evolving population of dynamical systems. The dynamical view of the origin of life allows us to pinpoint the missing and the not fitting pieces: (1) How can the first self-replicating ribozyme emerge in the absence of template-directed information replication? (2) How can nucleotide replicators avoid competitive exclusion despite utilizing the very same resources (nucleobases)? (3) How can the information catastrophe be avoided? (4) How can enough genes integrate into a cohesive system in order to transition to a cellular stage? (5) How can the way information is stored and metabolic complexity coevolve to pave to road leading out of the RNA world to the present protein-DNA world?

Using droplet-based microfluidics to improve the catalytic properties of RNA under multiple-turnover conditions

In vitro evolution methodologies are powerful approaches to identify RNA with new functionalities. While Systematic Evolution of Ligands by Exponential enrichment (SELEX) is an efficient approach to generate new RNA aptamers, it is less suited for the isolation of efficient ribozymes as it does not select directly for the catalysis. In vitro compartmentalization (IVC) in aqueous droplets in emulsions allows catalytic RNAs to be selected under multiple-turnover conditions but suffers severe limitations that can be overcome using the droplet-based microfluidics workflow described in this paper. Using microfluidics, millions of genes in a library can be individually compartmentalized in highly monodisperse aqueous droplets and serial operations performed on them. This allows the different steps of the evolution process (gene amplification, transcription, and phenotypic assay) to be uncoupled, making the method highly flexible, applicable to the selection and evolution of a variety of RNAs, and easily adaptable for evolution of DNA or proteins. To demonstrate the method, we performed cycles of random mutagenesis and selection to evolve the X-motif, a ribozyme which, like many ribozymes selected using SELEX, has limited multiple-turnover activity. This led to the selection of variants, likely to be the optimal ribozymes that can be generated using point mutagenesis alone, with a turnover number under multiple-turnover conditions, k(ss) cat, ∼ 28-fold higher than the original X-motif, primarily due to an increase in the rate of product release, the rate-limiting step in the multiple-turnover reaction.

Evolution of the division of labor between genes and enzymes in the RNA world

The RNA world is a very likely interim stage of the evolution after the first replicators and before the advent of the genetic code and translated proteins. Ribozymes are known to be able to catalyze many reaction types, including cofactor-aided metabolic transformations. In a metabolically complex RNA world, early division of labor between genes and enzymes could have evolved, where the ribozymes would have been transcribed from the genes more often than the other way round, benefiting the encapsulating cells through this dosage effect. Here we show, by computer simulations of protocells harboring unlinked RNA replicators, that the origin of replicational asymmetry producing more ribozymes from a gene template than gene strands from a ribozyme template is feasible and robust. Enzymatic activities of the two modeled ribozymes are in trade-off with their replication rates, and the relative replication rates compared to those of complementary strands are evolvable traits of the ribozymes. The degree of trade-off is shown to have the strongest effect in favor of the division of labor. Although some asymmetry between gene and enzymatic strands could have evolved even in earlier, surface-bound systems, the shown mechanism in protocells seems inevitable and under strong positive selection. This could have preadapted the genetic system for transcription after the subsequent origin of chromosomes and DNA.

Molecular crowding and early evolution

The environment of protocells might have been crowded with small molecules and functional and non-specific polymers. In addition to altering conformational equilibria, affecting reaction rates and changing the structure and activity of water, crowding might have enhanced the capabilities of protocells for evolutionary innovation through the creation of extended neutral networks in the fitness landscape.

Identification of RNA aptamers against recombinant proteins with a hexa-histidine tag

Artificial riboswitches that respond to the concentrations of intracellular proteins are a promising tool with a variety of applications. They can be designed and engineered using existing RNA aptamers that target proteins. Aptamers are generated via an iterative selection-amplification process, known as systematic evolution of ligands by exponential enrichment (SELEX). This chapter describes a SELEX procedure for the identification of RNA aptamers against hexa-histidine-tagged proteins. For the efficient enrichment of higher affinity aptamers, the selection stringency should be gradually increased. Undesired RNA species that bind to affinity resins can be eliminated from the pool by using a negative selection step and alternating different types of resins.

Experimental fitness landscapes to understand the molecular evolution of RNA-based life

In evolutionary biology, the relationship between genotype and Darwinian fitness is known as a fitness landscape. These landscapes underlie natural selection, so understanding them would greatly improve quantitative prediction of evolutionary outcomes, guiding the development of synthetic living systems. However, the structure of fitness landscapes is essentially unknown. Our ability to experimentally probe these landscapes is physically limited by the number of different sequences that can be identified. This number has increased dramatically in the last several years, leading to qualitatively new investigations. Several approaches to illuminate fitness landscapes are possible, ranging from tight focus on a single peak to random speckling or even comprehensive coverage of an entire landscape. We discuss recent experimental studies of fitness landscapes, with a special focus on functional RNA, an important system for both synthetic cells and the origin of life.

Continuous in vitro evolution of a ribozyme ligase: a model experiment for the evolution of a biomolecule

Evolution is a defining criterion of life and is central to understanding biological systems. However, the timescale of evolutionary shifts in phenotype limits most classroom evolution experiments to simple probability simulations. In vitro directed evolution (IVDE) frequently serves as a model system for the study of Darwinian evolution but produces noticeable phenotypic shifts in a matter of hours. An IVDE demonstration lab would serve to both directly demonstrate how Darwinian selection can act on a pool of variants and introduce students to an essential method of modern molecular biology. To produce an IVDE demonstration lab, continuous IVDE of a T500 ribozyme ligase population has been paired with a fluorescent strand displacement reporter system to visualize the selection of improved catalytic function. A ribozyme population is taken through rounds of isothermal amplification dependent on the self-ligation of a T7 promoter. As the population is selectively enriched with better ligase activity, the strand displacement system allows for the monitoring of the population's ligation rate. The strand displacement reporter system permits the detection of ligated ribozyme. Once ligated with the T7 promoter, the 5' end of the ribozyme displaces paired fluorophore-quencher oligonucleotides, in turn, generating visible signal upon UV light excitation. As the ligation rate of the population increases, due to the selection for faster ligating species, the fluorescent signal develops more rapidly. The pairing of the continuous isothermal system with the fluorescent reporting scheme allows any user, provided with minimal materials, to model the continuous directed evolution of a biomolecule.

In-ice evolution of RNA polymerase ribozyme activity

Mechanisms of molecular self-replication have the potential to shed light on the origins of life. In particular, self-replication through RNA-catalysed templated RNA synthesis is thought to have supported a primordial 'RNA world'. However, existing polymerase ribozymes lack the capacity to synthesize RNAs approaching their own size. Here, we report the in vitro evolution of such catalysts directly in the RNA-stabilizing medium of water ice, which yielded RNA polymerase ribozymes specifically adapted to sub-zero temperatures and able to synthesize RNA in ices at temperatures as low as -19 °C. The combination of cold-adaptive mutations with a previously described 5' extension operating at ambient temperatures enabled the design of a first polymerase ribozyme capable of catalysing the accurate synthesis of an RNA sequence longer than itself (adding up to 206 nucleotides), an important stepping stone towards RNA self-replication.

Fold or hold: experimental evolution in vitro

We introduce a system for experimental evolution consisting of populations of short oligonucleotides (Oli populations) evolving in a modified quantitative polymerase chain reaction (qPCR). It is tractable at the genetic, genomic, phenotypic and fitness levels. The Oli system uses DNA hairpins designed to form structures that self-prime under defined conditions. Selection acts on the phenotype of self-priming, after which differences in fitness are amplified and quantified using qPCR. We outline the methodological and bioinformatics tools for the Oli system here and demonstrate that it can be used as a conventional experimental evolution model system by test-driving it in an experiment investigating adaptive evolution under different rates of environmental change.

Comprehensive experimental fitness landscape and evolutionary network for small RNA

The origin of life is believed to have progressed through an RNA world, in which RNA acted as both genetic material and functional molecules. The structure of the evolutionary fitness landscape of RNA would determine natural selection for the first functional sequences. Fitness landscapes are the subject of much speculation, but their structure is essentially unknown. Here we describe a comprehensive map of a fitness landscape, exploring nearly all of sequence space, for short RNAs surviving selection in vitro. With the exception of a small evolutionary network, we find that fitness peaks are largely isolated from one another, highlighting the importance of historical contingency and indicating that natural selection would be constrained to local exploration in the RNA world.

Cascade of reduced speed and accuracy after errors in enzyme-free copying of nucleic acid sequences

Nonenzymatic, template-directed synthesis of nucleic acids is a paradigm for self-replicating systems. The evolutionary dynamics of such systems depend on several factors, including the mutation rates, relative replication rates, and sequence characteristics of mutant sequences. We measured the kinetics of correct and incorrect monomer insertion downstream of a primer-template mismatch (mutation), using a range of backbone structures (RNA, DNA, and LNA templates and RNA and DNA primers) and two types of 5'-activated nucleotides (oxyazabenzotriazolides and imidazolides, i.e., nucleoside 5'-phosphorimidazolides). Our study indicated that for all systems studied, an initial mismatch was likely to be followed by another error (54-75% of the time), and extension after a single mismatch was generally 10-100 times slower than extension without errors. If the mismatch was followed by a matched base pair, the extension rate recovered to nearly normal levels. On the basis of these data, we simulated nucleic acid replication in silico, which indicated that a primer suffering an initial error would lag behind properly extended counterparts due to a cascade of subsequent errors and kinetic stalling, with the typical mutational event consisting of several consecutive errors. Our study also included different sequence contexts, which suggest the presence of cooperativity among monomers affecting both absolute rate (by up to 2 orders of magnitude) and fidelity. The results suggest that molecular evolution in enzyme-free replication systems would be characterized by large "leaps" through sequence space rather than isolated point mutations, perhaps enabling rapid exploration of diverse sequences. The findings may also be useful for designing self-replicating systems combining high fidelity with evolvability.

Fixation and accumulation of thermotolerant catalytic competence of a pair of ligase ribozymes through complex formation and cross ligation

In the early stages of the hypothetical RNA world, some primitive RNA catalysts (ribozymes) may have emerged through self-assembly of short RNA oligomers. Although they may be unstable against temperature fluctuations and other environmental changes, ligase ribozymes (ribozymes with RNA strand-joining activity) may resolve structural instability of self-assembling RNAs by converting them to the corresponding unimolecular formats. To investigate this possibility, we constructed a model system using a cross-ligation system composed of a pair of self-assembling ligase ribozymes. Their abilities to act as catalysts, substrates, and a cross-ligation system were analyzed with or without thermal pretreatment before the reactions. A pair of self-assembling ligase ribozymes, each of which can form multiple conformations, demonstrated that thermotolerance was acquired and accumulated through complex-formation that stabilized the active forms of the bimolecular ribozymes and also cross-ligation that produced the unimolecular ribozymes.

Ocean-atmosphere interactions in the emergence of complexity in simple chemical systems

The prebiotic conversion of simple organic molecules into complex biopolymers necessary for life can only have emerged on a stage set by geophysics. The transition between "prebiotic soup," the diverse mixture of small molecules, and complex, self-replicating organisms requires passing through the bottleneck of fundamental chemistry. In this Account, we examine how water-air interfaces, namely, the surfaces of lakes, oceans, and atmospheric aerosols on ancient Earth, facilitated the emergence of complex structures necessary for life. Aerosols are liquid or solid suspensions in air with a broad, power law size distribution. Collectively, these globally distributed atmospheric particles have an enormous surface area. Organic films at the interface between water and air offer advantages for biomolecular synthesis compared with the bulk and can simultaneously participate in the folding of biopolymers into primitive enclosed structures. We survey the advantages of the water-air interface for prebiotic chemistry in a geophysical context from three points of view. We examine the formation of biopolymers from simple organic precursors and describe the necessity and availability of enclosures. In addition, we provide a statistical mechanical approach to natural selection and emergence of complexity that proposes a link between these molecular mechanisms and macroscopic scales. Very large aerosol populations were ubiquitous on ancient Earth, and the surfaces of lakes, oceans, and atmospheric aerosols would have provided an auspicious environment for the emergence of complex structures necessary for life. These prebiotic reactors would inevitably have incorporated the products of chemistry into their anhydrous, two-dimensional organic films in the three-dimensional fluids of the gaseous atmosphere and the liquid ocean. The untrammeled operation of natural selection on these aerosols provided the likely location where condensation reactions could form biopolymers by elimination of water. The fluctuating exposure of the large, recycling aerosol populations to radiation, pressure, temperature, and humidity over geological time allows complexity to emerge from simple molecular precursors. We propose an approach that connects chemical statistical thermodynamics and the macroscopic world of the planetary ocean and atmosphere.

Evolutionary dynamics of RNA-like replicator systems: A bioinformatic approach to the origin of life

We review computational studies on prebiotic evolution, focusing on informatic processes in RNA-like replicator systems. In particular, we consider the following processes: the maintenance of information by replicators with and without interactions, the acquisition of information by replicators having a complex genotype-phenotype map, the generation of information by replicators having a complex genotype-phenotype-interaction map, and the storage of information by replicators serving as dedicated templates. Focusing on these informatic aspects, we review studies on quasi-species, error threshold, RNA-folding genotype-phenotype map, hypercycle, multilevel selection (including spatial self-organization, classical group selection, and compartmentalization), and the origin of DNA-like replicators. In conclusion, we pose a future question for theoretical studies on the origin of life.

Evolution of an inhibitory RNA aptamer against T7 RNA polymerase

Aptamers are promising gene components that can be used for the construction of synthetic gene circuits. In this study, we isolated an RNA aptamer that specifically inhibits transcription of T7 RNA polymerase (RNAP). The 38-nucleotide aptamer, which was a shortened variant of an initial SELEX isolate, showed moderate inhibitory activity. By stepwise doped-SELEX, we isolated evolved variants with strong inhibitory activity. A 29-nucleotide variant of a doped-SELEX isolate showed 50% inhibitory concentration at 11 nM under typical in vitro transcription conditions. Pull-down experiments revealed that the aptamer inhibited the association of T7 RNAP with T7 promoter DNA.

A computational approach to predict suitable target sites for trans-acting minimal hammerhead ribozymes

Trans-acting hammerhead ribozymes are challenging tools for diagnostic, therapeutic, and biosensoristic purposes, owing to their specificity, efficiency, and great flexibility of use. One of the main problems in their application is related to the difficulties in the design of active molecules and identification of suitable target sites.The aim of this chapter is to describe ALADDIN, "SeArch computing tooL for hAmmerheaD ribozyme DesIgN," an open-access tool able to automatically identify suitable cleavage sites and provide a set of hammerhead ribozymes putatively active against the selected target.ALADDIN is a fast, cheap, helpful, and accurate tool designed to overcome the problems in the design of trans-acting minimal hammerhead ribozymes.

Selectionist and evolutionary approaches to brain function: a critical appraisal

We consider approaches to brain dynamics and function that have been claimed to be Darwinian. These include Edelman’s theory of neuronal group selection, Changeux’s theory of synaptic selection and selective stabilization of pre-representations, Seung’s Darwinian synapse, Loewenstein’s synaptic melioration, Adam’s selfish synapse, and Calvin’s replicating activity patterns. Except for the last two, the proposed mechanisms are selectionist but not truly Darwinian, because no replicators with information transfer to copies and hereditary variation can be identified in them. All of them fit, however, a generalized selectionist framework conforming to the picture of Price’s covariance formulation, which deliberately was not specific even to selection in biology, and therefore does not imply an algorithmic picture of biological evolution. Bayesian models and reinforcement learning are formally in agreement with selection dynamics. A classification of search algorithms is shown to include Darwinian replicators (evolutionary units with multiplication, heredity, and variability) as the most powerful mechanism for search in a sparsely occupied search space. Examples are given of cases where parallel competitive search with information transfer among the units is more efficient than search without information transfer between units. Finally, we review our recent attempts to construct and analyze simple models of true Darwinian evolutionary units in the brain in terms of connectivity and activity copying of neuronal groups. Although none of the proposed neuronal replicators include miraculous mechanisms, their identification remains a challenge but also a great promise.

Inhibitory RNA aptamer against SP6 RNA polymerase

Aptamers are attractive tools for modulating function of a desired target. In this study, we isolated an RNA aptamer that specifically inhibits transcription of SP6 RNA polymerase. The dissociation constant and 50% inhibitory concentration of the aptamer were estimated 9.5 nM and 24.8 nM, respectively. Doped-SELEX and mutational analysis revealed that the aptamer adopts the structure including two stems, two loops, and 5' single-stranded region. Based on the results, the aptamer could be engineered to circular permutant and binary construct forms without decreasing the activity. The aptamer would be applicable for the construction of expression regulation systems.

An in vivo selection method to optimize trans-splicing ribozymes

Group I intron ribozymes can repair mutated mRNAs by replacing the 3'-terminal portion of the mRNA with their own 3'-exon. This trans-splicing reaction has the potential to treat genetic disorders and to selectively kill cancer cells or virus-infected cells. However, these ribozymes have not yet been used in therapy, partially due to a low in vivo trans-splicing efficiency. Previous strategies to improve the trans-splicing efficiencies focused on designing and testing individual ribozyme constructs. Here we describe a method that selects the most efficient ribozymes from millions of ribozyme variants. This method uses an in vivo rescue assay where the mRNA of an inactivated antibiotic resistance gene is repaired by trans-splicing group I intron ribozymes. Bacterial cells that express efficient trans-splicing ribozymes are able to grow on medium containing the antibiotic chloramphenicol. We randomized a 5'-terminal sequence of the Tetrahymena thermophila group I intron and screened a library with 9 × 10⁶ ribozyme variants for the best trans-splicing activity. The resulting ribozymes showed increased trans-splicing efficiency and help the design of efficient trans-splicing ribozymes for different sequence contexts. This in vivo selection method can now be used to optimize any sequence in trans-splicing ribozymes.

The phylogenomic roots of modern biochemistry: origins of proteins, cofactors and protein biosynthesis

The complexity of modern biochemistry developed gradually on early Earth as new molecules and structures populated the emerging cellular systems. Here, we generate a historical account of the gradual discovery of primordial proteins, cofactors, and molecular functions using phylogenomic information in the sequence of 420 genomes. We focus on structural and functional annotations of the 54 most ancient protein domains. We show how primordial functions are linked to folded structures and how their interaction with cofactors expanded the functional repertoire. We also reveal protocell membranes played a crucial role in early protein evolution and show translation started with RNA and thioester cofactor-mediated aminoacylation. Our findings allow elaboration of an evolutionary model of early biochemistry that is firmly grounded in phylogenomic information and biochemical, biophysical, and structural knowledge. The model describes how primordial α-helical bundles stabilized membranes, how these were decorated by layered arrangements of β-sheets and α-helices, and how these arrangements became globular. Ancient forms of aminoacyl-tRNA synthetase (aaRS) catalytic domains and ancient non-ribosomal protein synthetase (NRPS) modules gave rise to primordial protein synthesis and the ability to generate a code for specificity in their active sites. These structures diversified producing cofactor-binding molecular switches and barrel structures. Accretion of domains and molecules gave rise to modern aaRSs, NRPS, and ribosomal ensembles, first organized around novel emerging cofactors (tRNA and carrier proteins) and then more complex cofactor structures (rRNA). The model explains how the generation of protein structures acted as scaffold for nucleic acids and resulted in crystallization of modern translation.

Downward causation by information control in micro-organisms

The concepts of functional equivalence classes and information control in living systems are useful to characterize downward (or top-down) causation by feedback information control in synthetic biology. Herein, we re-analyse published experiments of microbiology and synthetic biology that demonstrate the existence of several classes of functional equivalence in microbial organisms. Classes of functional equivalence from the bacterial operating system, which processes and controls the information encoded in the genome, can readily be interpreted as strong evidence, if not demonstration, of top-down causation (TDC) by information control. The proposed biological framework reveals how this type of causality is put in action in the cellular operating system. Considerations on TDC by information control and adaptive selection can be useful for synthetic biology by delineating the irreducible set of properties that characterizes living systems. Through a 'retro-synthetic' biology approach, these considerations could contribute to identifying the constraints behind the emergence of molecular complexity during the evolution of an ancient RNA/peptide world into a modern DNA/RNA/protein world. In conclusion, we propose TDCs by information control and adaptive selection as the two types of downward causality absolutely necessary for life.

Ribozymes: Flexible molecular devices at work

The discovery of ribozymes, RNAs with catalytic activity, revealed the extraordinary characteristic of this molecule, and corroborated the idea that RNA was the first informative polymer. The "RNA world" hypothesis asserts that the DNA/RNA/PROTEIN world arose from an earlier RNA world in which were present only RNA molecules able to perform both of the two functions performed separately by DNA and proteins in the present-day cells: the ability to transfer genetic information and to carry out catalytic activity. The catalytic properties of ribozymes are exclusively due to the capacity of RNA molecules to assume particular structures. Moreover, the structural versatility of RNA can allow to a single RNA sequence to fold in more than one structure, able to perform more than one function. In the first part of this work we will discuss the RNA plasticity, focusing on "bifunctional" ribozymes isolated by in vitro selection experiments, and on the consequences of this plasticity in the prospective of the emergence of new specific functions. The possibility that one sequence could have more than one structure/function, greatly increase the evolutionary potential of RNA, and the capacity of RNA to switch from a structure/function to another is probably one of the reasons of the evolutionary success also in modern-day cells. Naturally occurring ribozymes discovered in contemporary cells, demonstrate the crucial role that ribozymes still have in the modern protein world. In the second part of this paper we will discuss the capacity of natural ribozymes to modulate gene expression making use of their exclusive catalytic properties. Moreover, we will consider the possibility of their ancient origin.

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.

At the core of all biological systems lies the division of labor between the storage of genetic information and its phenotypic implementation, in other words, the functional differentiation between templates (DNA) and catalysts (proteins). This fundamental property of life is believed to have been absent at the earliest stages of evolution. The RNA world hypothesis, the most realistic current scenario for the origin of life, posits that, in primordial replicating systems, RNA functioned both as template and as catalyst. How would such division of labor emerge through Darwinian evolution? We investigated the evolution of DNA-like molecules in minimal computational models of RNA replicator systems. Two models were considered: one where molecules are adsorbed on surfaces and another one where molecules are compartmentalized by dividing cellular boundaries. Both models exhibit the evolution of DNA and the ensuing division of labor, revealing the simple governing principle of these processes: DNA releases RNA from the trade-off between template and catalyst that is inevitable in the RNA world and thereby enhances the system's resistance against parasitic templates. Hence, this study offers a novel insight into the evolutionary origin of the division of labor between templates and catalysts in the RNA world.

Protonation of base pairs in RNA: context analysis and quantum chemical investigations of their geometries and stabilities

Base pairs involving protonated nucleobases play important roles in mediating global macromolecular conformational changes and in facilitation of catalysis in a variety of functional RNA molecules. Here we present our attempts at understanding the role of such base pairs by detecting possible protonated base pairs in the available RNA crystal structures using BPFind software, in their specific structural contexts, and by the characterization of their geometries, interaction energies, and stabilities using advanced quantum chemical computations. We report occurrences of 18 distinct protonated base pair combinations from a representative data set of RNA crystal structures and propose a theoretical model for one putative base pair combination. Optimization of base pair geometries was carried out at the B3LYP/cc-pVTZ level, and the BSSE corrected interaction energies were calculated at the MP2/aug-cc-pVDZ level of theory. The geometries for each of the base pairs were characterized in terms of H-bonding patterns observed, rmsd values observed on optimization, and base pair geometrical parameters. In addition, the intermolecular interaction in these complexes was also analyzed using Morokuma energy decomposition. The gas phase interaction energies of the base pairs range from -24 to -49 kcal/mol and reveal the dominance of Hartree-Fock component of interaction energy constituting 73% to 98% of the total interaction energy values. On the basis of our combined bioinformatics and quantum chemical analysis of different protonated base pairs, we suggest resolution of structural ambiguities and correlate their geometric and energetic features with their structural and functional roles. In addition, we also examine the suitability of specific base pairs as key elements in molecular switches and as nucleators for higher order structures such as base triplets and quartets.

Proteome evolution and the metabolic origins of translation and cellular life

The origin of life has puzzled molecular scientists for over half a century. Yet fundamental questions remain unanswered, including which came first, the metabolic machinery or the encoding nucleic acids. In this study we take a protein-centric view and explore the ancestral origins of proteins. Protein domain structures in proteomes are highly conserved and embody molecular functions and interactions that are needed for cellular and organismal processes. Here we use domain structure to study the evolution of molecular function in the protein world. Timelines describing the age and function of protein domains at fold, fold superfamily, and fold family levels of structural complexity were derived from a structural phylogenomic census in hundreds of fully sequenced genomes. These timelines unfold congruent hourglass patterns in rates of appearance of domain structures and functions, functional diversity, and hierarchical complexity, and revealed a gradual build up of protein repertoires associated with metabolism, translation and DNA, in that order. The most ancient domain architectures were hydrolase enzymes and the first translation domains had catalytic functions for the aminoacylation and the molecular switch-driven transport of RNA. Remarkably, the most ancient domains had metabolic roles, did not interact with RNA, and preceded the gradual build-up of translation. In fact, the first translation domains had also a metabolic origin and were only later followed by specialized translation machinery. Our results explain how the generation of structure in the protein world and the concurrent crystallization of translation and diversified cellular life created further opportunities for proteomic diversification.

Ice as a protocellular medium for RNA replication

A crucial transition in the origin of life was the emergence of an informational polymer capable of self-replication and its compartmentalization within protocellular structures. We show that the physicochemical properties of ice, a simple medium widespread on a temperate early Earth, could have mediated this transition prior to the advent of membraneous protocells. Ice not only promotes the activity of an RNA polymerase ribozyme but also protects it from hydrolytic degradation, enabling the synthesis of exceptionally long replication products. Ice furthermore relieves the dependence of RNA replication on prebiotically implausible substrate concentrations, while providing quasicellular compartmentalization within the intricate microstructure of the eutectic phase. Eutectic ice phases had previously been shown to promote the de novo synthesis of nucleotide precursors, as well as the condensation of activated nucleotides into random RNA oligomers. Our results support a wider role for ice as a predisposed environment, promoting all the steps from prebiotic synthesis to the emergence of RNA self-replication and precellular Darwinian evolution.

Sequence variability and sequence evolution: an explanation of molecular polymorphisms and why many molecular structures can be preserved although they are not predominant

The existence of many processes that regulate RNA expression poses a challenge to the idea that the cell is the culmination of a highly efficient interplay of individual proteins, each with specific, highly specialized functions. It will be demonstrated here the extent to which the cell may undergo evolutionary processes that also occur in the macrocosmos, specifically with reference to the rules of mutation and preservation. These molecular evolutionary processes could facilitate a better understanding of the development of molecular structures and the functioning of the cell and could give an explanation of the molecular polymorphisms and also explain why many molecular structures can be preserved although they are not predominant.

Closing the circle: replicating RNA with RNA

How life emerged on this planet is one of the most important and fundamental questions of science. Although nearly all details concerning our origins have been lost in the depths of time, there is compelling evidence to suggest that the earliest life might have exploited the catalytic and self-recognition properties of RNA to survive. If an RNA based replicating system could be constructed in the laboratory, it would be much easier to understand the challenges associated with the very earliest steps in evolution and provide important insight into the establishment of the complex metabolic systems that now dominate this planet. Recent progress into the selection and characterization of ribozymes that promote nucleotide synthesis and RNA polymerization are discussed and outstanding problems in the field of RNA-mediated RNA replication are summarized.

Intercalation as a means to suppress cyclization and promote polymerization of base-pairing oligonucleotides in a prebiotic world

The RNA world hypothesis proposes that nucleic acids were once responsible for both information storage and chemical catalysis, before the advent of coded protein synthesis. However, it is difficult to imagine how nucleic acid polymers first appeared, as the abiotic chemical formation of long nucleic acid polymers from mononucleotides or short oligonucleotides remains elusive, and barriers to achieving this goal are substantial. One specific obstacle to abiotic nucleic acid polymerization is strand cyclization. Chemically activated short oligonucleotides cyclize efficiently, which severely impairs polymer growth. We show that intercalation, which stabilizes and rigidifies nucleic acid duplexes, almost totally eliminates strand cyclization, allowing for chemical ligation of tetranucleotides into duplex polymers of up to 100 base pairs in length. In contrast, when these reactions are performed in the absence of intercalators, almost exclusively cyclic tetra- and octanucleotides are produced. Intercalator-free polymerization is not observed, even at tetranucleotide concentrations > 10,000-fold greater than those at which intercalators enable polymerization. We also demonstrate that intercalation-mediated polymerization is most favored if the size of the intercalator matches that of the base pair; intercalators that bind to Watson-Crick base pairs promote the polymerization of oligonucleotides that form these base pairs. Additionally, we demonstrate that intercalation-mediated polymerization is possible with an alternative, non-Watson-Crick-paired duplex that selectively binds a complementary intercalator. These results support the hypothesis that intercalators (acting as 'molecular midwives') could have facilitated the polymerization of the first nucleic acids and possibly helped select the first base pairs, even if only trace amounts of suitable oligomers were available.

Genotype-phenotype mapping and the end of the 'genes as blueprint' metaphor

In a now classic paper published in 1991, Alberch introduced the concept of genotype-phenotype (G-->P) mapping to provide a framework for a more sophisticated discussion of the integration between genetics and developmental biology that was then available. The advent of evo-devo first and of the genomic era later would seem to have superseded talk of transitions in phenotypic space and the like, central to Alberch's approach. On the contrary, this paper shows that recent empirical and theoretical advances have only sharpened the need for a different conceptual treatment of how phenotypes are produced. Old-fashioned metaphors like genetic blueprint and genetic programme are not only woefully inadequate but positively misleading about the nature of G-->P, and are being replaced by an algorithmic approach emerging from the study of a variety of actual G-->P maps. These include RNA folding, protein function and the study of evolvable software. Some generalities are emerging from these disparate fields of analysis, and I suggest that the concept of 'developmental encoding' (as opposed to the classical one of genetic encoding) provides a promising computational-theoretical underpinning to coherently integrate ideas on evolvability, modularity and robustness and foster a fruitful framing of the G-->P mapping problem.

The continuous evolution in vitro technique

In vitro experimentation techniques were developed in response to the necessity of exploring new molecular structures and functions and to better understand evolutionary phenomena that shape organismal and molecular populations. The advancement of these techniques has allowed further exploration of more complicated evolutionary dynamics. One such technique is the continuous evolution in vitro (CE) method, to which this unit is devoted. The CE method is characterized by continuous cycles of amplification of RNA molecules that occur without much participation of the researcher. This feature allows us to evolve lineages in which the evolutionary phenomena occurring at the molecular level more closely mimic what happens in organismal populations in the present, or what may have happened in RNA populations during the RNA world stage of life.

Evolutionary optimization of a modular ligase ribozyme: a small catalytic unit and a hairpin motif masking an element that could form an inactive structure

The YFL ribozyme is an artificial ligase ribozyme isolated by a ‘design and selection’ strategy, in which a modular catalytic unit was generated on a rationally designed modular scaffold RNA. This ligase ribozyme has a versatile catalytic unit that accepts not only β-nicotinamide mononucleotide (β-NMN) but also inorganic pyrophosphate as leaving groups for template-dependent RNA ligation. Although this property is interesting from an evolutionary viewpoint regarding primitive RNA ligation/polymerization systems in the RNA world, structural analysis of the YFL ribozyme has not been continued due to apparent structural nonuniformity of its folded state. To elucidate the active structure of the YFL ribozyme, we performed in vitro evolution experiments to improve its folding ability. Biochemical and phylogenetic analyses of evolved variants indicated that the catalytic unit of the YFL ribozyme is compact and the 3′ single-stranded region of the parent YFL-1 ribozyme contributes to mask an element that could form an inactive structure.

Strategies for the design of RNA-binding small molecules

Bacterial ribosomal RNA is the target of clinically important antibiotics, while biologically important RNAs in viral and eukaryotic genomes present a range of potential drug targets. The physicochemical properties of RNA present difficulties for medicinal chemistry, particularly when oral availability is needed. Peptidic ligands and analysis of their RNA-binding properties are providing insight into RNA recognition. RNA-binding ligands include far more chemical classes than just aminoglycosides. Chemical functionalities from known RNA-binding small molecules are being exploited in fragment- and ligand-based projects. While targeting of RNA for drug design is very challenging, continuing advances in our understanding of the principles of RNA–ligand interaction will be necessary to realize the full potential of this class of targets.