In vitro selections with RNAs of variable length converge on a robust catalytic core.

In vitro selection is a powerful tool that can be used to understand basic principles of molecular evolution. We used in vitro selection to understand how changes in length and the accumulation of point mutations enable the evolution of functional RNAs. Using RNA populations of various lengths, we performed a series of in vitro experiments to select for ribozymes with RNA ligase activity. We identified a core ribozyme structure that was robust to changes in RNA length, high levels of mutagenesis, and increased selection pressure. Elaboration on this core structure resulted in improved activity which we show is consistent with a larger trend among functional RNAs in which increasing motif size can lead to an exponential improvement in fitness. We conclude that elaboration on conserved core structures is a preferred mechanism in RNA evolution. This conclusion, drawn from selections of RNAs from random sequences, is consistent with proposed evolutionary histories of specific biological RNAs. More generally, our results indicate that modern RNA structures can be used to infer ancestral structures. Our observations also suggest a mechanism by which structural outcomes of early RNA evolution would be largely reproducible even though RNA fitness landscapes consist of disconnected clusters of functional sequences.

Big on Change, Small on Innovation: Evolutionary Consequences of RNA Sequence Duplication.

The potential for biopolymers to evolve new structures has important consequences for their ability to optimize function and our attempts to reconstruct their evolutionary histories. Prior work with in vitro systems suggests that structural remodeling of RNA is difficult to achieve through the accumulation of point mutations or through recombination events. Sequence duplication may represent an alternative mechanism that can more readily lead to the evolution of new structures. Structural and sequence elements in many RNAs and proteins appear to be the products of duplication events, indicating that this mechanism plays a major role in molecular evolution. Despite the potential significance of this mechanism, little experimental data is available concerning the structural and evolutionary consequences of duplicating biopolymer sequences. To assess the structural consequences of sequence duplication on the evolution of RNA, we mutagenized an RNA sequence containing two copies of an ATP aptamer and subjected the resulting population to multiple in vitro evolution experiments. We identified multiple routes by which duplication, followed by the accumulation of functional point mutations, allowed our populations to sample two entirely different secondary structures. The two structures have no base pairs in common, but both structures contain two copies of the same ATP-binding motif. We do not observe the emergence of any other functional secondary structures beyond these two. Although this result suggests a limited capacity for duplication to support short-term functional innovation, major changes in secondary structure, like the one observed here, should be given careful consideration as they are likely to frustrate attempts to infer deep evolutionary histories of functional RNAs.

Evolution of ribozymes in the presence of a mineral surface.

Mineral surfaces are often proposed as the sites of critical processes in the emergence of life. Clay minerals in particular are thought to play significant roles in the origin of life including polymerizing, concentrating, organizing, and protecting biopolymers. In these scenarios, the impact of minerals on biopolymer folding is expected to influence evolutionary processes. These processes include both the initial emergence of functional structures in the presence of the mineral and the subsequent transition away from the mineral-associated niche. The initial evolution of function depends upon the number and distribution of sequences capable of functioning in the presence of the mineral, and the transition to new environments depends upon the overlap between sequences that evolve on the mineral surface and sequences that can perform the same functions in the mineral's absence. To examine these processes, we evolved self-cleaving ribozymes in vitro in the presence or absence of Na-saturated montmorillonite clay mineral particles. Starting from a shared population of random sequences, RNA populations were evolved in parallel, along separate evolutionary trajectories. Comparative sequence analysis and activity assays show that the impact of this clay mineral on functional structure selection was minimal; it neither prevented common structures from emerging, nor did it promote the emergence of new structures. This suggests that montmorillonite does not improve RNA's ability to evolve functional structures; however, it also suggests that RNAs that do evolve in contact with montmorillonite retain the same structures in mineral-free environments, potentially facilitating an evolutionary transition away from a mineral-associated niche.

In vitro evolution of distinct self-cleaving ribozymes in diverse environments.

In vitro evolution experiments have long been used to evaluate the roles of RNA in both modern and ancient biology, and as a tool for biotechnology applications. The conditions under which these experiments have been conducted, however, do not reflect the range of cellular environments in modern biology or our understanding of chemical environments on the early earth, when the atmosphere and oceans were largely anoxic and soluble Fe(2+) was abundant. To test the impact of environmental factors relevant to RNA's potential role in the earliest forms of life, we evolved populations of self-cleaving ribozymes in an anoxic atmosphere with varying pH in the presence of either Fe(2+) or Mg(2+). Populations evolved under these different conditions are dominated by different sequences and secondary structures, demonstrating global differences in the underlying fitness landscapes. Comparisons between evolutionary outcomes and catalytic activities also indicate that Mg(2+) can readily take the place of Fe(2+) in supporting the catalysis of RNA cleavage at neutral pH, but not at lower pH. These results highlight the importance of considering the specific environments in which functional biopolymers evolve when evaluating their potential roles in the origin of life, extant biology, or biotechnology.

High-throughput sequence analysis reveals structural diversity and improved potency among RNA inhibitors of HIV reverse transcriptase.

Systematic evolution of ligands through exponential enrichment (SELEX) is a well-established method for generating nucleic acid populations that are enriched for specified functions. High-throughput sequencing (HTS) enhances the power of comparative sequence analysis to reveal details of how RNAs within these populations recognize their targets. We used HTS analysis to evaluate RNA populations selected to bind type I human immunodeficiency virus reverse transcriptase (RT). The populations are enriched in RNAs of independent lineages that converge on shared motifs and in clusters of RNAs with nearly identical sequences that share common ancestry. Both of these features informed inferences of the secondary structures of enriched RNAs, their minimal structural requirements and their stabilities in RT-aptamer complexes. Monitoring population dynamics in response to increasing selection pressure revealed RNA inhibitors of RT that are more potent than the previously identified pseudoknots. Improved potency was observed for inhibition of both purified RT in enzymatic assays and viral replication in cell-based assays. Structural and functional details of converged motifs that are obscured by simple consensus descriptions are also revealed by the HTS analysis. The approach presented here can readily be generalized for the efficient and systematic post-SELEX development of aptamers for down-stream applications.

Broad-spectrum aptamer inhibitors of HIV reverse transcriptase closely mimic natural substrates.

A detailed understanding of how aptamers recognize biological binding partners is of considerable importance in the development of oligonucleotide therapeutics. For antiviral nucleic acid aptamers, current models predict a correlation between broad-spectrum inhibition of viral proteins and suppression of emerging viral resistance, but there is little understanding of how aptamer structures contribute to recognition specificity. We previously established that two independent single-stranded DNA aptamers, R1T and RT1t49(-5), are potent inhibitors of reverse transcriptases (RTs) from diverse branches of the primate lentiviral family, including HIV-1, HIV-2 and SIV(cpz). In contrast, class 1 RNA pseudoknots, such as aptamer T1.1, are specific for RTs from only a few viral clades. Here, we map the binding interfaces of complexes formed between RT and aptamers R1T, RT1t49(-5) and T1.1, using mass spectrometry-based protein footprinting of RT and hydroxyl radical footprinting of the aptamers. These complementary methods reveal that the broad-spectrum aptamers make contacts throughout the primer-template binding cleft of RT. The double-stranded stems of these aptamers closely mimic natural substrates near the RNase H domain, while their binding within the polymerase domain significantly differs from RT substrates. These results inform our perspective on how sustained, broad-spectrum inhibition of RT can be achieved by aptamers.

Molecular dynamics suggest multifunctionality of an adenine imino group in acid-base catalysis of the hairpin ribozyme.

Despite numerous structural and biochemical investigations, the catalytic mechanism of hairpin ribozyme self-cleavage remains elusive. To gain insight into the coupling of active site dynamics with activity of this small catalytic RNA, we analyzed a total of approximately 300 ns of molecular dynamics (MD) simulations. Our simulations predict improved global stability for an in vitro selected "gain of function" mutation, which is validated by native gel electrophoretic mobility shift assay. We observe that active site nucleobases and water molecules stabilize a geometry favorable to catalysis through a dynamic hydrogen bonding network. Simulations in which A38 is unprotonated show its N1 move into close proximity of the active site 2'-OH, indicating that A38 may act as a general base during cleavage, a role that has generally been discounted due to the longer distances observed in crystal structures involving inactivating substrate analogs. By contrast, simulations in which N1 of A38 is protonated place N1 in close proximity to the 5'-oxygen leaving group, which supports the proposal that A38 serves as a general acid. In analogy to protein enzymes, we discuss a plausible mechanism in which A38 acts bifunctionally and shuttles a proton directly from the 2'-OH to the 5'-oxygen. Furthermore, our simulations suggest an important role for protonation of N1 of A38 in promoting a favorable geometry similar to that observed in transition-state analog crystal structures, and support previously proposed roles of A38, G8, and long residency water molecules in transition-state stabilization.

Molecular dynamics and quantum mechanics of RNA: conformational and chemical change we can believe in.

Structure and dynamics are both critical to RNA's vital functions in biology. Numerous techniques can elucidate the structural dynamics of RNA, but computational approaches based on experimental data arguably hold the promise of providing the most detail. In this Account, we highlight areas wherein molecular dynamics (MD) and quantum mechanical (QM) techniques are applied to RNA, particularly in relation to complementary experimental studies. We have expanded on atomic-resolution crystal structures of RNAs in functionally relevant states by applying explicit solvent MD simulations to explore their dynamics and conformational changes on the submicrosecond time scale. MD relies on simplified atomistic, pairwise additive interaction potentials (force fields). Because of limited sampling, due to the finite accessible simulation time scale and the approximated force field, high-quality starting structures are required. Despite their imperfection, we find that currently available force fields empower MD to provide meaningful and predictive information on RNA dynamics around a crystallographically defined energy minimum. The performance of force fields can be estimated by precise QM calculations on small model systems. Such calculations agree reasonably well with the Cornell et al. AMBER force field, particularly for stacking and hydrogen-bonding interactions. A final verification of any force field is accomplished by simulations of complex nucleic acid structures. The performance of the Cornell et al. AMBER force field generally corresponds well with and augments experimental data, but one notable exception could be the capping loops of double-helical stems. In addition, the performance of pairwise additive force fields is obviously unsatisfactory for inclusion of divalent cations, because their interactions lead to major polarization and charge-transfer effects neglected by the force field. Neglect of polarization also limits, albeit to a lesser extent, the description accuracy of other contributions, such as interactions with monovalent ions, conformational flexibility of the anionic sugar-phosphate backbone, hydrogen bonding, and solute polarization by solvent. Still, despite limitations, MD simulations are a valid tool for analyzing the structural dynamics of existing experimental structures. Careful analysis of MD simulations can identify problematic aspects of an experimental RNA structure, unveil structural characteristics masked by experimental constraints, reveal functionally significant stochastic fluctuations, evaluate the structural role of base ionization, and predict structurally and potentially functionally important details of the solvent behavior, including the presence of tightly bound water molecules. Moreover, combining classical MD simulations with QM calculations in hybrid QM/MM approaches helps in the assessment of the plausibility of chemical mechanisms of catalytic RNAs (ribozymes). In contrast, the reliable prediction of structure from sequence information is beyond the applicability of MD tools. The ultimate utility of computational studies in understanding RNA function thus requires that the results are neither blindly accepted nor flatly rejected, but rather considered in the context of all available experimental data, with great care given to assessing limitations through the available starting structures, force field approximations, and sampling limitations. The examples given in this Account showcase how the judicious use of basic MD simulations has already served as a powerful tool to help evaluate the role of structural dynamics in biological function of RNA.

A rugged free energy landscape separates multiple functional RNA folds throughout denaturation.

The dynamic mechanisms by which RNAs acquire biologically functional structures are of increasing importance to the rapidly expanding fields of RNA therapeutics and biotechnology. Large energy barriers separating misfolded and functional states arising from alternate base pairing are a well-appreciated characteristic of RNA. In contrast, it is typically assumed that functionally folded RNA occupies a single native basin of attraction that is free of deeply dividing energy barriers (ergodic hypothesis). This assumption is widely used as an implicit basis to interpret experimental ensemble-averaged data. Here, we develop an experimental approach to isolate persistent sub-populations of a small RNA enzyme and show by single molecule fluorescence resonance energy transfer (smFRET), biochemical probing and high-resolution mass spectrometry that commitment to one of several catalytically active folds occurs unexpectedly high on the RNA folding energy landscape, resulting in partially irreversible folding. Our experiments reveal the retention of molecular heterogeneity following the complete loss of all native secondary and tertiary structure. Our results demonstrate a surprising longevity of molecular heterogeneity and advance our current understanding beyond that of non-functional misfolds of RNA kinetically trapped on a rugged folding-free energy landscape.

Selection of 2'-Deoxy-2'-Fluoroarabino Nucleic Acid (FANA) Aptamers That Bind HIV-1 Integrase with Picomolar Affinity.

Systematic Evolution of Ligands by Exponential Enrichment (SELEX) is the iterative process by which nucleic acids that can bind with high affinity and specificity (termed aptamers) to specific protein targets are selected. Using a SELEX protocol adapted for Xeno-Nucleic Acid (XNA) as a suitable substrate for aptamer generation, 2'-fluoroarabinonucleic acid (FANA) was used to select several related aptamers to HIV-1 integrase (IN). IN bound FANA aptamers with equilibrium dissociation constants (KD,app) of ∼50-100 pM in a buffer with 200 mM NaCl and 6 mM MgCl2. Comparisons to published HIV-1 IN RNA and DNA aptamers as well as IN genomic binding partners indicated that FANA aptamers bound more than 2 orders of magnitude more tightly to IN. Using a combination of RNA folding algorithms and covariation analysis, all strong binding aptamers demonstrated a common four-way junction structure, despite significant sequence variation. IN aptamers were selected from the same starting library as FA1, a FANA aptamer that binds with pM affinity to HIV-1 Reverse Transcriptase (RT). It contains a 20-nucleotide 5' DNA sequence followed by 59 FANA nucleotides. IN-1.1 (one of the selected aptamers) potently inhibited IN activity and intasome formation in vitro. Replacing the FANA nucleotides of IN-1.1 with 2'-fluororibonucleic acid (F-RNA), which has the same chemical formula but with a ribose rather than arabinose sugar conformation, dramatically reduced binding, suggesting that FANA adopts unique structural conformations that promote binding to HIV-1 IN.

Potent Inhibition of HIV-1 Reverse Transcriptase and Replication by Nonpseudoknot, "UCAA-motif" RNA Aptamers.

RNA aptamers that bind the reverse transcriptase (RT) of human immunodeficiency virus (HIV) compete with nucleic acid primer/template for access to RT, inhibit RT enzymatic activity in vitro, and suppress viral replication when expressed in human cells. Numerous pseudoknot aptamers have been identified by sequence analysis, but relatively few have been confirmed experimentally. In this work, a screen of nearly 100 full-length and >60 truncated aptamer transcripts established the predictive value of the F1Pk and F2Pk pseudoknot signature motifs. The screen also identified a new, nonpseudoknot motif with a conserved unpaired UCAA element. High-throughput sequence (HTS) analysis identified 181 clusters capable of forming this novel element. Comparative sequence analysis, enzymatic probing and RT inhibition by aptamer variants established the essential requirements of the motif, which include two conserved base pairs (AC/GU) on the 5' side of the unpaired UCAA. Aptamers in this family inhibit RT in primer extension assays with IC(50) values in the low nmol/l range, and they suppress viral replication with a potency that is comparable with that of previously studied aptamers. All three known anti-RT aptamer families (pseudoknots, the UCAA element, and the recently described "(6/5)AL" motif) are therefore suitable for developing aptamer-based antiviral gene therapies.Molecular Therapy - Nucleic Acids (2013) 2, e71; doi:10.1038/mtna.2012.62; published online 5 February 2013.

Conformational dynamics of single pre-mRNA molecules during in vitro splicing.

The spliceosome is a complex small nuclear RNA (snRNA)-protein machine that removes introns from pre-mRNAs via two successive phosphoryl transfer reactions. The chemical steps are isoenergetic, yet splicing requires at least eight RNA-dependent ATPases responsible for substantial conformational rearrangements. To comprehensively monitor pre-mRNA conformational dynamics, we developed a strategy for single-molecule FRET (smFRET) that uses a small, efficiently spliced yeast pre-mRNA, Ubc4, in which donor and acceptor fluorophores are placed in the exons adjacent to the 5' and 3' splice sites. During splicing in vitro, we observed a multitude of generally reversible time- and ATP-dependent conformational transitions of individual pre-mRNAs. The conformational dynamics of branchpoint and 3'-splice site mutants differ from one another and from wild type. Because all transitions are reversible, spliceosome assembly appears to be occurring close to thermal equilibrium.

Focus on function: single molecule RNA enzymology.

The ability of RNA to catalyze chemical reactions was first demonstrated 25 years ago with the discovery that group I introns and RNase P function as RNA enzymes (ribozymes). Several additional ribozymes were subsequently identified, most notably the ribosome, followed by intense mechanistic studies. More recently, the introduction of single molecule tools has dissected the kinetic steps of several ribozymes in unprecedented detail and has revealed surprising heterogeneity not evident from ensemble approaches. Still, many fundamental questions of how RNA enzymes work at the molecular level remain unanswered. This review surveys the current status of our understanding of RNA catalysis at the single molecule level and discusses the existing challenges and opportunities in developing suitable assays.