Structurally complex and highly active RNA ligases derived from random RNA sequences

DNA nanodevices

The molecular recognition properties of DNA molecules combined with the distinct mechanical properties of single and double strands of DNA can be utilized for the construction of nanodevices, which can perform ever more tasks with possible applications ranging from nanoconstruction to intelligent drug delivery. With the help of DNA it is possible to construct machinelike devices that are capable of rotational motion, pulling and stretching, or even unidirectional motion. It is possible to devise autonomous nanodevices, to grab and release molecules, and also to perform simple information-processing tasks.

Rational modification of a selection strategy leads to deoxyribozymes that create native 3'-5' RNA linkages

We previously used in vitro selection to identify several classes of deoxyribozymes that mediate RNA ligation by attack of a hydroxyl group at a 5'-triphosphate. In these reactions, the nucleophilic hydroxyl group is located at an internal 2'-position of an RNA substrate, leading to 2',5'-branched RNA. To obtain deoxyribozymes that instead create linear 3'-5'-linked (native) RNA, here we strategically modified the selection approach by embedding the nascent ligation junction within an RNA:DNA duplex region. This approach should favor formation of linear rather than branched RNA because the two RNA termini are spatially constrained by Watson-Crick base pairing during the ligation reaction. Furthermore, because native 3'-5' linkages are more stable in a duplex than isomeric non-native 2'-5' linkages, this strategy is predicted to favor the formation of 3'-5' linkages. All of the new deoxyribozymes indeed create only linear 3'-5' RNA, confirming the effectiveness of the rational design. The new deoxyribozymes ligate RNA with k(obs) values up to 0.5 h(-1) at 37 degrees C and 40 mM Mg2+, pH 9.0, with up to 41% yield at 3 h incubation. They require several specific RNA nucleotides on either side of the ligation junction, which may limit their practical generality. These RNA ligase deoxyribozymes are the first that create native 3'-5' RNA linkages, which to date have been highly elusive via other selection approaches.

High-throughput ribozyme-based assays for detection of viral nucleic acids

Many reports have suggested that target-activated ribozymes hold potential value as detection reagents. We show that a "half"-ribozyme ligase is activated similarly by three unstructured oligoribonucleotides representing the major sequence variants of a hepatitis C virus 5'-untranslated region (5'-UTR) target and by a structured RNA corresponding to the entire 5'-UTR. Half-ribozyme ligation product was detected both in an ELISA-like assay and in an optical immunoassay through the use of hapten-carrying substrate RNAs. Both assay formats afford a limit of detection of approximately 1 x 10(6) HCV molecules (1.6 attomol, 330 fM), a sensitivity which compares favorably to that provided by standard immunoassays. These data suggest that target-activated ribozyme systems are a viable approach for the sensitive detection of viral nucleic acids using high-throughput platforms.

In vitro evolution and characterization of a ligase ribozyme adapted to acidic conditions: Effect of further rounds of evolution

A ligase ribozyme that accelerates the ligation reaction with an oligonucleotide under low pH conditions was identified by in vitro adaptation in a previous study. We examined the effects of further rounds of evolution to isolate a more active ribozyme. The ribozyme, which was obtained after four rounds of evolution, was randomly mutated, and the resultant RNA library was subjected to in vitro selection at low pH. One ribozyme isolated from the pool was found to react 8,000 times faster than the original b1 ribozyme at pH 4. The reaction rate of the isolated ribozyme was enhanced at various pH values, and its pH dependence was less than that of the original ribozyme or the ribozyme selected with four rounds of evolution. The reaction rate of the isolated ribozyme was reduced in the presence of 3' primer, the sequence of which is complementary to the 3' primer-binding site of the ligase ribozyme. This inhibition induced by the primer oligonucleotide binding to the ribozyme 3' region implies that the 3' region plays a role in the ligation reaction of the ribozyme.

Directed evolution of nucleic acid enzymes

Just as Darwinian evolution in nature has led to the development of many sophisticated enzymes, Darwinian evolution in vitro has proven to be a powerful approach for obtaining similar results in the laboratory. This review focuses on the development of nucleic acid enzymes starting from a population of random-sequence RNA or DNA molecules. In order to illustrate the principles and practice of in vitro evolution, two especially well-studied categories of catalytic nucleic acid are considered: RNA enzymes that catalyze the template-directed ligation of RNA and DNA enzymes that catalyze the cleavage of RNA. The former reaction, which involves attack of a 2'- or 3'-hydroxyl on the alpha-phosphate of a 5'-triphosphate, is more difficult. It requires a comparatively larger catalytic motif, containing more nucleotides than can be sampled exhaustively within a starting population of random-sequence RNAs. The latter reaction involves deprotonation of the 2'-hydroxyl adjacent to the cleavage site, resulting in cleaved products that bear a 2',3'-cyclic phosphate and 5'-hydroxyl. The difficulty of this reaction, and therefore the complexity of the corresponding DNA enzyme, depends on whether a catalytic cofactor, such as a divalent metal cation or small molecule, is present in the reaction mixture.

De novo synthesis and development of an RNA enzyme

Arbitrary manipulation of molecular recognition at the atomic level has many applications. However, systematic design and de novo synthesis of an artificial enzyme based on such manipulation has been a long-standing challenge in the field of chemistry and biotechnology. In this report, we developed an artificial RNA ligase by implementing a synthetic strategy that fuses a series of 3D molecular modelings based on naturally occurring RNA-RNA recognition motifs with a small-scale combinatorial synthesis of a modular catalytic unit. The resulting ligase produces a 3'-5' linkage in a template-directed manner for any combinations of two nucleotides at the reaction site. The reaction rate is 10(6)-fold over that of the uncatalyzed reaction with a yield higher than those of previously reported ligase ribozymes. The strategy may be applicable to the synthesis and development of a variety of nonnatural functional RNAs with defined 3D structures.

Origins of introns based on the definition of exon modules and their conserved interfaces

Central to the unraveling of the early evolution of the genome is the origin and role of introns. The evolution of the genome can be characterized by a continuous expansion of functional modules that occurs without the interruption of existing processes. The design-by-contract methodology of software development offers a modular approach to design that seeks to increase flexibility by focusing on the design of constant interfaces between functional modules. Here, it is shown that design-by-contract can offer a framework for genome evolution. The definition of an ancient exon module with identical splice sites leads to a relatively simple sequence of events that explains the role of introns, intron phase differences and the evolution of multi-exon proteins in an RNA world. An interaction of the experimentally defined six-nucleotide splicing consensus sequence together with a limited number of primitive ribozymes can account for a rapid creation of protein diversity.

Informational complexity and functional activity of RNA structures

Very little is known about the distribution of functional DNA, RNA, and protein molecules in sequence space. The question of how the number and complexity of distinct solutions to a particular biochemical problem varies with activity is an important aspect of this general problem. Here we present a comparison of the structures and activities of eleven distinct GTP-binding RNAs (aptamers). By experimentally measuring the amount of information required to specify each optimal binding structure, we show that defining a structure capable of 10-fold tighter binding requires approximately 10 additional bits of information. This increase in information content is equivalent to specifying the identity of five additional nucleotide positions and corresponds to an approximately 1000-fold decrease in abundance in a sample of random sequences. We observe a similar relationship between structural complexity and activity in a comparison of two catalytic RNAs (ribozyme ligases), raising the possibility of a general relationship between the complexity of RNA structures and their functional activity. Describing how information varies with activity in other heteropolymers, both biological and synthetic, may lead to an objective means of comparing their functional properties. This approach could be useful in predicting the functional utility of novel heteropolymers.

Expanded divalent metal-ion tolerance of evolved ligase ribozymes

Class I ligases are artificial ribozymes that catalyze the joining of two single-stranded RNAs. These ribozymes are between 120 and 160 nucleotides in length, making them intermediate in size for catalytic RNAs. Previous characterization of the b1-207 ribozyme suggests that it behaves similar to larger ribozymes in terms of divalent metal-ion dependence. This molecule displays a strong preference for magnesium for catalysis, and is inactive in any other metal except manganese, which actually inhibits its operation in magnesium. Here, we sought to examine the metal-ion usages of two ligases that were obtained through continuous evolution in vitro from the b1-207 sequence framework. We found an expanded catalytic range for the E(100)(#3) and B(16)(#19) ribozymes, as they are both catalytically active in calcium and strontium, and less inhibited by manganese. Though not selected for activity in these salts, the evolved ribozymes exhibit several adaptations to in vitro catalysis, and their ability to accommodate metals other than magnesium can be viewed as an example of a molecular exaptation.

Practical and general synthesis of 5'-adenylated RNA (5'-AppRNA)

A simple strategy is reported for 5'-adenylation of nearly any RNA sequence of indefinite length. The 5'-adenylated product (5'-AppRNA) is an activated RNA that is structurally similar to 5'-triphosphorylated RNA, which is usually prepared by in vitro transcription using T7 RNA polymerase. In the new 5'-adenylation strategy, the RNA substrate is first 5'-monophosphorylated either by T4 polynucleotide kinase, by in vitro transcription in the presence of excess GMP, or by appropriate derivatization during solid-phase synthesis. The RNA is then 5'-adenylated using ATP and T4 RNA ligase, in an interrupted version of the natural adenylation-ligation mechanism by which T4 RNA ligase joins two RNA substrates. Here, the final ligation step of the mechanism is inhibited with complementary DNA blocking oligonucleotide(s) that permit adenylation to occur with good yield. The 5'-AppRNA products of this approach should be valuable as activated RNAs for in vitro selection experiments as an alternative to 5'-triphosphorylated RNAs, among other likely applications. The 5'-terminal nucleotide of an RNA substrate to be adenylated using the new method is not restricted to guanosine, in contrast to 5'-triphosphorylated RNA prepared by in vitro transcription. Therefore, using the new approach, essentially any RNA obtained from solid-phase synthesis or other means can be activated by 5'-adenylation in a practical manner.

The three-dimensional architecture of the class I ligase ribozyme

The class I ligase ribozyme catalyzes a Mg(++)-dependent RNA-ligation reaction that is chemically analogous to a single step of RNA polymerization. Indeed, this ribozyme constitutes the catalytic domain of an accurate and general RNA polymerase ribozyme. The ligation reaction is also very rapid in both single- and multiple-turnover contexts and thus is informative for the study of RNA catalysis as well as RNA self-replication. Here we report the initial characterization of the three-dimensional architecture of the ligase. When the ligase folds, several segments become protected from hydroxyl-radical cleavage, indicating that the RNA adopts a compact tertiary structure. Ribozyme folding was largely, though not completely, Mg(++) dependent, with a K(1/2[Mg]) < 1 mM, and was observed over a broad temperature range (20 degrees C -50 degrees C). The hydroxyl-radical mapping, together with comparative sequence analyses and analogy to a region within 23S ribosomal RNA, were used to generate a three-dimensional model of the ribozyme. The predictive value of the model was tested and supported by a photo-cross-linking experiment.

Aptamer database

The aptamer database is designed to contain comprehensive sequence information on aptamers and unnatural ribozymes that have been generated by in vitro selection methods. Such data are not normally collected in 'natural' sequence databases, such as GenBank. Besides serving as a storehouse of sequences that may have diagnostic or therapeutic utility, the database serves as a valuable resource for theoretical biologists who describe and explore fitness landscapes. The database is updated monthly and is publicly available at http://aptamer. icmb.utexas.edu/.

Zeptomole detection of a viral nucleic acid using a target-activated ribozyme

We describe a strategy for the ultra-sensitive detection of nucleic acids using "half" ribozymes that are devoid of catalytic activity unless completed by a trans-acting target nucleic acid. The half-ribozyme concept was initially demonstrated using a construct derived from a multiple turnover Class I ligase. Iterative RNA selection was carried out to evolve this half-ribozyme into one activated by a conserved sequence present in the hepatitis C virus (HCV) genome. Following sequence optimization of substrate RNAs, this HCV-activated half-ribozyme displayed a maximal turnover rate of 69 min(-1) (pH 8.3) and was induced in rate by approximately 2.6 x 10(9)-fold by the HCV target. It detected the HCV target oligonucleotide in the zeptomole range (6700 molecules), a sensitivity of detection roughly 2.6 x 10(6)-fold greater than that previously demonstrated by oligonucleotide-activated ribozymes, and one that is sufficient for molecular diagnostic applications.

Primordial genetics: phenotype of the ribocyte

The idea that the ancestors of modern cells were RNA cells (ribocytes) can be investigated by asking whether all essential cellular functions might be performed by RNAs. This requires isolating suitable molecules by selection-amplification when the predicted molecules are presently extinct. In fact, RNAs with many properties required during a period in which RNA was the major macromolecular agent in cells (an RNA world) have been selected in modern experiments. There is, accordingly, reason to inquire how such a ribocyte might appear, based on the properties of the RNAs that composed it. Combining the intrinsic qualities of RNA with the fundamental characteristics of selection from randomized sequence pools, one predicts ribocytes with a cell cycle measured (roughly) in weeks. Such cells likely had a rapidly varying genome, composed of many small genetic and catalytic elements made of tens of ribonucleotides. There are substantial arguments that, at the mid-RNA era, a subset of these nucleotides are reproducibly available and resemble the modern four. Such cells are predicted to evolve rapidly. Instead of modifying preexisting genes, ribocytes frequently draw new functions from an internal pool containing zeptomoles (<1 attomole) of predominantly inactive random sequences.

Deoxyribozymes with 2'-5' RNA ligase activity

In vitro selection was used to identify deoxyribozymes that ligate two RNA substrates. In the ligation reaction, a 2'-5' RNA phosphodiester linkage is created from a 2',3'-cyclic phosphate and a 5'-hydroxyl group. The new Mg(2+)-dependent deoxyribozymes provide 50-60% yield of ligated RNA in overnight incubations at pH 7.5 and 37 degrees C, and they afford 40-50% yield in 1 h at pH 9.0 and 37 degrees C. Various RNA substrate sequences may be joined by simple Watson-Crick covaration of the DNA binding arms that interact with the two RNA substrates. The current deoxyribozymes have some RNA substrate sequence requirements at the nucleotides immediately surrounding the ligation junction (either UAUA GGAA or UAUN GGAA, where the arrow denotes the ligation site and N equals any nucleotide). One of the new deoxyribozymes was used to prepare by ligation the Tetrahymena group I intron RNA P4-P6 domain, a representative structured RNA. Nondenaturing gel electrophoresis revealed that a 2'-5' linkage between nucleotides A233 and G234 of P4-P6 does not disrupt its Mg(2+)-dependent folding (DeltaDeltaG degrees ' < 0.2 kcal/mol). This demonstrates that a 2'-5' linkage does not necessarily interfere with structure in a folded RNA. Therefore, these non-native linkages may be acceptable in modified RNAs when structure/function relationships are investigated. Deoxyribozymes that ligate RNA should be particularly useful for preparing site-specifically modified RNAs for studies of RNA structure, folding, and catalysis.

A ribozyme with michaelase activity: synthesis of the substrate precursors

The ability to generate RNA molecules that can catalyze complex organic transformations not only facilitates the reconstruction and plausibility of possible prebiotic reaction pathways but is also crucial for elucidating the potential of the application of RNA catalysts in organic syntheses. Iterative RNA selection previously identified a ribozyme that catalyzes the Michael addition of a cysteine thiol to an alpha,beta-unsaturated amide. This reaction is chemically similar to the rate limiting step of the thymidylate synthase reaction, which is the corresponding reaction of a cysteine thiol to the double-bond of the uracil nucleobase. Here we provide a detailed description of the synthesis of the ribozyme substrates and the substrate oligonucleotides used for its characterization and the investigation of the background reaction. We also describe the further characterization of the ribozyme with respect to substrate specificity. We show that the thiol group of the cysteine nucleophile is essential for the reaction to proceed. When substituted for a thiomethyl group, no reaction takes place.

In search of an RNA replicase ribozyme

The theory that an RNA world played a pivotal role in life's evolutionary past has prompted investigations into the scope of RNA catalysis. These efforts have attempted to demonstrate the plausibility of an RNA-based genetic system, which would require RNA molecules that catalyze their own replication. The mechanistic features of modern protein polymerases have been used to guide the laboratory evolution of catalytic RNAs (ribozymes) that exhibit polymerase-like activity. Ribozymes have been developed that recognize a primer-template complex in a general way and catalyze the template-directed polymerization of mononucleotides. These experiments demonstrate that RNA replicase behavior is likely within the catalytic repertoire of RNA, although many obstacles remain to be overcome in order to demonstrate that RNA can catalyze its own replication in a manner that could have sustained a genetic system on the early Earth.

RNA-cleaving DNA enzymes with altered regio- or enantioselectivity

In vitro evolution methods were used to obtain DNA enzymes that cleave either a 2',5'-phosphodiester following a D-ribonucleotide or a 3',5'-phosphodiester following an L-ribonucleotide. Both enzymes can operate in an intermolecular reaction format with multiple turnover. The DNA enzyme that cleaves a 2',5'-phosphodiester exhibits a k(cat) of approximately 0.01 min(-1) and catalytic efficiency, k(cat)/K(m), of approximately 10(8) M(-1) min(-1). The enzyme that cleaves an L-ribonucleotide is about 10-fold slower and has a catalytic efficiency of approximately 4 x 10(5) M(-1) min(-1). Both enzymes require a divalent metal cation for their activity and have optimal catalytic rate at pH 7-8 and 35-50 degrees C. In a comparison of each enzyme's activity with either its corresponding substrate that contains an unnatural ribonucleotide or a substrate that instead contains a standard ribonucleotide, the 2',5'-phosphodiester-cleaving DNA enzyme exhibited a regioselectivity of 6000-fold, while the L-ribonucleotide-cleaving DNA enzyme exhibited an enantioselectivity of 40-fold. These molecules demonstrate how in vitro evolution can be used to obtain regio- and enantioselective catalysts that exhibit specificities for nonnatural analogues of biological compounds.

Nucleic acid evolution and minimization by nonhomologous random recombination

We have developed a simple method for exploring nucleic acid sequence space by nonhomologous random recombination (NRR) that enables DNA fragments to randomly recombine in a length-controlled manner without the need for sequence homology. We compared the results of using NRR and error-prone PCR to evolve DNA aptamers that bind streptavidin. Starting with two parental sequences of modest avidin affinity, evolution using NRR resulted in aptamers with 15- to 20-fold higher affinity than the highest-affinity aptamers evolved using error-prone PCR, and 27- or 46-fold higher affinities than parental sequences derived using systematic evolution of ligands by exponential enrichment (SELEX). NRR also facilitates the identification of functional regions within evolved sequences. Inspection of a small number of NRR-evolved clones identified a 40-base DNA sequence, present in multiple copies in each clone, that binds streptavidin. Our findings suggest that NRR may enhance the effectiveness of nucleic acid evolution and the ease of identifying structure-activity relationships among evolved sequences.

The antiquity of RNA-based evolution

All life that is known to exist on Earth today and all life for which there is evidence in the geological record seems to be of the same form--one based on DNA genomes and protein enzymes. Yet there are strong reasons to conclude that DNA- and protein-based life was preceded by a simpler life form based primarily on RNA. This earlier era is referred to as the 'RNA world', during which the genetic information resided in the sequence of RNA molecules and the phenotype derived from the catalytic properties of RNA.

A ribozyme that ligates RNA to protein

We have used a combination of in vitro selection and rational design to generate ribozymes that form a stable phosphoamide bond between the 5' terminus of an RNA and a specific polypeptide. This reaction differs from that of previously identified ribozymes, although the product is analogous to the enzyme-nucleotidyl intermediates isolated during the reactions of certain proteinaceous enzymes, such as guanyltransferase, DNA ligase, and RNA ligase. Comparative sequence analysis of the isolated ribozymes revealed that they share a compact secondary structure containing six stems arranged in a four-helix junction and branched pseudoknot. An optimized version of the ribozyme reacts with substrate-fusion proteins, allowing it to be used to attach RNA tags to proteins both in vitro and within bacterial cells, suggesting a simple way to tag a specific protein with amplifiable information.

Continuous in vitro evolution of a ribozyme that catalyzes three successive nucleotidyl addition reactions

Variants of the class I ligase ribozyme, which catalyzes joining of the 3' end of a template bound oligonucleotide to its own 5' end, have been made to evolve in a continuous manner by a simple serial transfer procedure that can be carried out indefinitely. This process was expanded to allow the evolution of ribozymes that catalyze three successive nucleotidyl addition reactions, two template-directed mononucleotide additions followed by RNA ligation. During the development of this behavior, a population of ribozymes was maintained against an overall dilution of more than 10(406). The resulting ribozymes were capable of catalyzing the three-step reaction pathway, with nucleotide addition occurring in either a 5'-->3' or a 3'-->5' direction. This purely chemical system provides a functional model of a multi-step reaction pathway that is undergoing Darwinian evolution.

RNA-catalyzed RNA ligation on an external RNA template

Variants of the hc ligase ribozyme, which catalyzes ligation of the 3' end of an RNA substrate to the 5' end of the ribozyme, were utilized to evolve a ribozyme that catalyzes ligation reactions on an external RNA template. The evolved ribozyme catalyzes the joining of an oligonucleotide 3'-hydroxyl to the 5'-triphosphate of an RNA hairpin molecule. The ribozyme can also utilize various substrate sequences, demonstrating a largely sequence-independent mechanism for substrate recognition. The ribozyme also carries out the ligation of two oligonucleotides that are bound at adjacent positions on a complementary template. Finally, it catalyzes addition of mononucleoside 5'-triphosphates onto the 3' end of an oligonucleotide primer in a template-dependent manner. The development of ribozymes that catalyze polymerase-type reactions contributes to the notion that an RNA world could have existed during the early history of life on Earth.

In vitro adaptation of a ligase ribozyme for activity under a low-pH condition

A ligase ribozyme accelerating a ligation reaction with oligonucleotide under a low-pH condition was selected by in vitro adaptation. A ribozyme active at pH 7 was randomly mutated, and the resultant RNA library was subjected to in vitro adaptation under a low-pH reaction condition. At pH 4, the adapted RNAs reacted with the oligonucleotide substrates about 200 times faster than the original ribozyme. When the ribozyme was cloned and sequenced, 10 of the 30 clones sequenced had identical sequences. The differences in sequence from the original ribozyme were found at four positions in the middle region and at the 3' end. A few sequential differences dominated the activity of the ribozyme under the extreme condition. The adapted ribozyme had one repeating sequence that was critical for the activity.

Investigating and Engineering Enzymes by Genetic Selection

Natural enzymes have arisen over millions of years by the gradual process of Darwinian evolution. The fundamental steps of evolution-mutation, selection, and amplification-can also be exploited in the laboratory to create and characterize protein catalysts on a human timescale. In vivo genetic selection strategies enable the exhaustive analysis of protein libraries with 10(10) different members, and even larger ensembles can be studied with in vitro methods. Evolutionary approaches can consequently yield statistically meaningful insight into the complex and often subtle interactions that influence protein folding, structure, and catalytic mechanism. Such methods are also being used increasingly as an adjunct to design, thus providing access to novel proteins with tailored catalytic activities and selectivities.

Strong RNA splicing enhancers identified by a modified method of cycled selection interact with SR protein

A modified method of cycled selection was used to characterize splicing enhancers for exon inclusion from a pool of beta-globin-based three exon/two intron pre-mRNAs with a variable number of random nucleotides incorporated in the internal exon. The pre-mRNAs generated by this method contained random sequences ranging from 0 to 18 nucleotides in length. This method was used to isolate particular splicing enhancer motifs from a previously enriched pool of extremely diverse enhancers. After four cycles of selection for mRNA containing the internal exon, a distinct enhancer motif (GACGAC...CAGCAG) was highly enriched. This motif served as strong splicing enhancers in a heterogeneous exon. We have shown here that the selected enhancer motif promotes exon inclusion through specific interaction with SRp30. We have also shown that although present in many of our selected splicing enhancers conforming to this motif, a typical purine-rich enhancer sequence is dispensable for either enhancer activity or binding with SRp30.

Evolution in silico and in vitro: the RNA model

Theoretical concepts and experiments dealing with the evolution of molecules in vitro reached a state that allows for direct applications to the design of biomolecules with predefined properties. RNA evolution in vitro represents a basis for the development of a new and comprehensive model of evolution, focusing on the phenotype and its fitness relevant properties. Relations between genotypes and phenotypes are described by mappings from genotype space onto a space of phenotypes, which are many-to-one and thus give ample room for neutrality as expressed by the existence of extended neutral networks in genotype space. The RNA model reduces genotype-phenotype relations to mappings from sequences into secondary structures of minimal free energies and allows for derivation of otherwise inaccessible quantitative results. Continuity and discontinuity in evolution are defined through a new notion of accessibility in phenotype space that provides a basis for straight forward interpretation of computer simulations on RNA optimization; furthermore, it reveals the constructive role of random genomic drift in the search for phenotypes of higher fitness. The effects of population size on the course of evolutionary optimization can be predicted quantitatively by means of a simple stochastic model based on a birth-anddeath process with immigration.

RNA-catalyzed RNA polymerization: accurate and general RNA-templated primer extension

The RNA world hypothesis regarding the early evolution of life relies on the premise that some RNA sequences can catalyze RNA replication. In support of this conjecture, we describe here an RNA molecule that catalyzes the type of polymerization needed for RNA replication. The ribozyme uses nucleoside triphosphates and the coding information of an RNA template to extend an RNA primer by the successive addition of up to 14 nucleotides-more than a complete turn of an RNA helix. Its polymerization activity is general in terms of the sequence and the length of the primer and template RNAs, provided that the 3' terminus of the primer pairs with the template. Its polymerization is also quite accurate: when primers extended by 11 nucleotides were cloned and sequenced, 1088 of 1100 sequenced nucleotides matched the template.

Critical analysis of antibody catalysis

Antibody molecules elicited with rationally designed transition-state analogs catalyze numerous reactions, including many that cannot be achieved by standard chemical methods. Although relatively primitive when compared with natural enzymes, these catalysts are valuable tools for probing the origins and evolution of biological catalysis. Mechanistic and structural analyses of representative antibody catalysts, generated with a variety of strategies for several different reaction types, suggest that their modest efficiency is a consequence of imperfect hapten design and indirect selection. Development of improved transition-state analogs, refinements in immunization and screening protocols, and elaboration of general strategies for augmenting the efficiency of first-generation catalytic antibodies are identified as evident, but difficult, challenges for this field. Rising to these challenges and more successfully integrating programmable design with the selective forces of biology will enhance our understanding of enzymatic catalysis. Further, it should yield useful protein catalysts for an enhanced range of practical applications in chemistry and biology.

Natural 2',5'-phosphodiester bonds found at the ligation sites of peach latent mosaic viroid

Peach latent mosaic viroid (PLMVd) is a circular RNA pathogen that replicates in a DNA-independent fashion via a rolling circle mechanism. PLMVd has been shown to self-ligate in vitro primarily via the formation of 2',5'-phosphodiester bonds; however, in vivo the occurrence and necessity of this nonenzymatic mechanism are not evident. Here, we unequivocally report the presence of 2', 5'-phosphodiester bonds at the ligation site of circular PLMVd strands isolated from infected peach leaves. These bonds serve to close the linear conformers (i.e., intermediates), yielding circular ones. Furthermore, these bonds are shown to stabilize the replicational circular templates, resulting in a significant advantage in terms of viroid viability. Although the mechanism responsible for the formation of these 2',5'-phosphodiester bonds remains to be elucidated, a hypothesis describing in vivo nonenzymatic self-ligation is proposed. Most significantly, our results clearly show that 2',5'-phosphodiester bonds are still present in nature and that they are of biological importance.

Plasticity, evolvability, and modularity in RNA

RNA folding from sequences into secondary structures is a simple yet powerful, biophysically grounded model of a genotype-phenotype map in which concepts like plasticity, evolvability, epistasis, and modularity can not only be precisely defined and statistically measured but also reveal simultaneous and profoundly non-independent effects of natural selection. Molecular plasticity is viewed here as the capacity of an RNA sequence to assume a variety of energetically favorable shapes by equilibrating among them at constant temperature. Through simulations based on experimental designs, we study the dynamics of a population of RNA molecules that evolve toward a predefined target shape in a constant environment. Each shape in the plastic repertoire of a sequence contributes to the overall fitness of the sequence in proportion to the time the sequence spends in that shape. Plasticity is costly, since the more shapes a sequence can assume, the less time it spends in any one of them. Unsurprisingly, selection leads to a reduction of plasticity (environmental canalization). The most striking observation, however, is the simultaneous slow-down and eventual halting of the evolutionary process. The reduction of plasticity entails genetic canalization, that is, a dramatic loss of variability (and hence a loss of evolvability) to the point of lock-in. The causal bridge between environmental canalization and genetic canalization is provided by a correlation between the set of shapes in the plastic repertoire of a sequence and the set of dominant (minimum free energy) shapes in its genetic neighborhood. This statistical property of the RNA genotype-phenotype map, which we call plastogenetic congruence, traps populations in regions where most genetic variation is phenotypically neutral. We call this phenomenon neutral confinement. Analytical models of neutral confinement, made tractable by the assumption of perfect plastogenetic congruence, formally connect mutation rate, the topography of phenotype space, and evolvability. These models identify three mutational regimes: that corresponding to neutral confinement, an exploration threshold corresponding to a breakdown of neutral confinement with the simultaneous persistence of the dominant phenotype, and a classic error threshold corresponding to the loss of the dominant phenotype. In a final step, we analyze the structural properties of canalized phenotypes. The reduction of plasticity leads to extreme modularity, which we analyze from several perspectives: thermophysical (melting--the RNA version of a norm of reaction), kinetic (folding pathways--the RNA version of development), and genetic (transposability--the insensitivity to genetic context). The model thereby suggests a possible evolutionary origin of modularity as a side effect of environmental canalization.

A system of two polymerases--a model for the origin of life

What was the first living molecule--RNA or protein? This question embodies the major disagreement in studies on the origin of life. The fact that in contemporary cells RNA polymerase is a protein and peptidyl transferase consists of RNA suggests the existence of a mutual catalytic dependence between these two kinds of biopolymers. I suggest that this dependence is a 'frozen accident', a remnant from the first living system. This system is proposed to be a combination of an RNA molecule capable of catalyzing amino acid polymerization and the resulting protein functioning as an RNA-dependent RNA polymerase. The specificity of the protein synthesis is thought to be achieved by the composition of the surrounding medium and the specificity of the RNA synthesis--by Watson-Crick base pairing. Despite its apparent simplicity, the system possesses a great potential to evolve into a primitive ribosome and further to life, as it is seen today. This model provides a possible explanation for the origin of the interaction between nucleic acids and protein. Based on the suggested system, I propose a new definition of life as a system of nucleic acid and protein polymerases with a constant supply of monomers, energy and protection.

Molecular evolution of catalysis

In this paper, we consider the evolutionary dynamics of catalytically active species with a distinct genotype-phenotype relationship. Folding landscapes of RNA molecules serve as a paradigm for this relationship with essential neutral properties. The landscape itself is partitioned by phenotypes (realized as RNA secondary structures). To each genotype (represented as a sequence) a structure is assigned in a unique way. The set of all sequences which map into a particular structure is modeled as a random graph in sequence space (the so-called neutral network). A catalytic network is realized as a random digraph with maximal out-degree two and secondary structures as vertex sets. A population of catalytic RNA molecules shows significantly different behavior compared to a deterministic description: hypercycles are able to co-exist and out-compete a parasite with superior catalytic support. A "switching" between different dynamic organizations of the network can be observed, dynamical stability of hypercyclic organizations against errors and the existence of an error-threshold of catalysis can be reported.

RNA structure. Ribozyme evolution at the crossroads

Molecules that have similar sequences usually adopt the same structures and have the same functions. In his Perspective, Joyce explains that this is not always the case. In a remarkable study (Schultes and Bartel), an RNA sequence has been designed that can adopt two different structures, each with a different catalytic function. Joyce details how this study sheds light on the evolution of enzymes.

One sequence, two ribozymes: implications for the emergence of new ribozyme folds

We describe a single RNA sequence that can assume either of two ribozyme folds and catalyze the two respective reactions. The two ribozyme folds share no evolutionary history and are completely different, with no base pairs (and probably no hydrogen bonds) in common. Minor variants of this sequence are highly active for one or the other reaction, and can be accessed from prototype ribozymes through a series of neutral mutations. Thus, in the course of evolution, new RNA folds could arise from preexisting folds, without the need to carry inactive intermediate sequences. This raises the possibility that biological RNAs having no structural or functional similarity might share a common ancestry. Furthermore, functional and structural divergence might, in some cases, precede rather than follow gene duplication.

In vitro selection of functional nucleic acids

In vitro selection allows rare functional RNA or DNA molecules to be isolated from pools of over 10(15) different sequences. This approach has been used to identify RNA and DNA ligands for numerous small molecules, and recent three-dimensional structure solutions have revealed the basis for ligand recognition in several cases. By selecting high-affinity and -specificity nucleic acid ligands for proteins, promising new therapeutic and diagnostic reagents have been identified. Selection experiments have also been carried out to identify ribozymes that catalyze a variety of chemical transformations, including RNA cleavage, ligation, and synthesis, as well as alkylation and acyl-transfer reactions and N-glycosidic and peptide bond formation. The existence of such RNA enzymes supports the notion that ribozymes could have directed a primitive metabolism before the evolution of protein synthesis. New in vitro protein selection techniques should allow for a direct comparison of the frequency of ligand binding and catalytic structures in pools of random sequence polynucleotides versus polypeptides.

Replicability and recurrence in the experimental evolution of a group I ribozyme

In order to explore the variety of possible responses available to a ribozyme population evolving a novel phenotype, five Tetrahymena thermophila group I intron ribozyme pools were evolved in parallel for cleavage of a DNA oligonucleotide. These ribozyme populations were propagated under identical conditions and characterized when they reached apparent phenotypic plateaus; the populations that reached the highest plateau showed a near 100-fold improvement in DNA cleavage activity. A detailed characterization of the evolved response in these populations reveals at least two distinct phenotypic trajectories emerging as a result of the imposed selection. Not only do these distinct solutions exhibit differential DNA cleavage activity, but they also exhibit a very different correlation with a related, but unselected, phenotype: RNA cleavage activity. In turn, each of these trajectories is underwritten by differing genotypic profiles. This study underscores the complex network of possible trajectories through sequence space available to an evolving population and uncovers the diversity of solutions that result when the process of experimental evolution is repeated multiple times in a simple, engineered system.

Physico-chemical constraints connected with the coding properties of the genetic system

New insights on the origin of the genetic code, based on the analysis of the physico-chemical properties of its molecular constituents (RNA and amino acids), are reported in this paper. We point out a symmetry in the genetic code table and show that it can be explained by the nature of the anticodon-codon interaction. The importance of the strength of this interaction is examined and a correlation is found between the free-energy change (DeltaG(0)) of anticodon-codon association and the volume of the corresponding amino acids. This correlation is investigated in conjunction with the well-known one linking the hydrophobicity of the anticodons with that of the amino acids. We show that they can be considerated separately and that the energy vs. volume correlation may be explained by the process implicating the peptide bond formation between two successive amino acids during translation. This interpretation is supported by a statistical pattern of bases (purines or pyrimidines), observed in present coding genes, and by considerations involving the availability of the different kinds of amino acids. Finally, we try to explain the hydrophobicity correlation when reconstructing the events at the time of the so-called "RNA World". The whole of our investigation shows that the genetic code might be sufficiently robust to exist without the participation of pre-existing proteins, and that this robustness is a consequence of the physico-chemical properties of the four bases of the genetic system.

A complex ligase ribozyme evolved in vitro from a group I ribozyme domain

Like most proteins, complex RNA molecules often are modular objects made up of distinct structural and functional domains. The component domains of a protein can associate in alternative combinations to form molecules with different functions. These observations raise the possibility that complex RNAs also can be assembled from preexisting structural and functional domains. To test this hypothesis, an in vitro evolution procedure was used to isolate a previously undescribed class of complex ligase ribozymes, starting from a pool of 10(16) different RNA molecules that contained a constant region derived from a large structural domain that occurs within self-splicing group I ribozymes. Attached to this constant region were three hypervariable regions, totaling 85 nucleotides, that gave rise to the catalytic motif within the evolved catalysts. The ligase ribozymes catalyze formation of a 3',5'-phosphodiester linkage between adjacent template-bound oligonucleotides, one bearing a 3' hydroxyl and the other a 5' triphosphate. Ligation occurs in the context of a Watson-Crick duplex, with a catalytic rate of 0.26 min(-1) under optimal conditions. The constant region is essential for catalytic activity and appears to retain the tertiary structure of the group I ribozyme. This work demonstrates that complex RNA molecules, like their protein counterparts, can share common structural domains while exhibiting distinct catalytic functions.