A class I ligase ribozyme with reduced Mg2+ dependence: Selection, sequence analysis, and identification of functional tertiary interactions

Structural basis for activation of fluorogenic dyes by an RNA aptamer lacking a G-quadruplex motif

The DIR2s RNA aptamer, a second-generation, in-vitro selected binder to dimethylindole red (DIR), activates the fluorescence of cyanine dyes, DIR and oxazole thiazole blue (OTB), allowing detection of two well-resolved emission colors. Using Fab BL3-6 and its cognate hairpin as a crystallization module, we solved the crystal structures of both the apo and OTB-SO₃ bound forms of DIR2s at 2.0 Å and 1.8 Å resolution, respectively. DIR2s adopts a compact, tuning fork-like architecture comprised of a helix and two short stem-loops oriented in parallel to create the ligand binding site through tertiary interactions. The OTB-SO₃ fluorophore binds in a planar conformation to a claw-like structure formed by a purine base-triple, which provides a stacking platform for OTB-SO_3, and an unpaired nucleotide, which partially caps the binding site from the top. The absence of a G-quartet or base tetrad makes the DIR2s aptamer unique among fluorogenic RNAs with known 3D structure.

The DIR2s RNA aptamer activates the fluorescence of cyanine dyes allowing detection of two well-resolved emission colors. Here authors solve the crystal structures of the apo and OTB-SO3 fluorophore-bound DIR2s and show how the fluorophore ligand is bound.

Handmade microfluidic device for biochemical applications in emulsion

A simple, inexpensive flow-focusing device has been developed to make uniform droplets for biochemical reactions, such as in vitro transcription and cell-free protein synthesis. The device was fabricated from commercially available components without special equipment. Using the emulsion droplets formed by the device, a class I ligase ribozyme, bcI 23, was successfully synthesized from DNA attached to magnetic microbeads by T7 RNA polymerase. It was also ligated with an RNA substrate on the same microbeads, and detected using flow cytometry with a fluorescent probe. In addition, a single-chain derivative of the lambda Cro protein was expressed using an Escherichia coli cell-free protein synthesis system in emulsion, which was prepared using the flow-focusing device. In both emulsified reactions, usage of the flow-focusing device was able to greatly reduce the coefficient of variation for the amount of RNA or protein displayed on the microbeads, demonstrating the device is advantageous for quantitative analysis in high-throughput screening.

Automated band annotation for RNA structure probing experiments with numerous capillary electrophoresis profiles

MOTIVATION

Capillary electrophoresis (CE) is a powerful approach for structural analysis of nucleic acids, with recent high-throughput variants enabling three-dimensional RNA modeling and the discovery of new rules for RNA structure design. Among the steps composing CE analysis, the process of finding each band in an electrophoretic trace and mapping it to a position in the nucleic acid sequence has required significant manual inspection and remains the most time-consuming and error-prone step. The few available tools seeking to automate this band annotation have achieved limited accuracy and have not taken advantage of information across dozens of profiles routinely acquired in high-throughput measurements.

RESULTS

We present a dynamic-programming-based approach to automate band annotation for high-throughput capillary electrophoresis. The approach is uniquely able to define and optimize a robust target function that takes into account multiple CE profiles (sequencing ladders, different chemical probes, different mutants) collected for the RNA. Over a large benchmark of multi-profile datasets for biological RNAs and designed RNAs from the EteRNA project, the method outperforms prior tools (QuSHAPE and FAST) significantly in terms of accuracy compared with gold-standard manual annotations. The amount of computation required is reasonable at a few seconds per dataset. We also introduce an 'E-score' metric to automatically assess the reliability of the band annotation and show it to be practically useful in flagging uncertainties in band annotation for further inspection.

AVAILABILITY AND IMPLEMENTATION

The implementation of the proposed algorithm is included in the HiTRACE software, freely available as an online server and for download at http://hitrace.stanford.edu.

CONTACT

sryoon@snu.ac.kr or rhiju@stanford.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Enhanced group II intron retrohoming in magnesium-deficient Escherichia coli via selection of mutations in the ribozyme core

Mobile group II introns are bacterial retrotransposons thought to be evolutionary ancestors of spliceosomal introns and retroelements in eukaryotes. They consist of a catalytically active intron RNA ("ribozyme") and an intron-encoded reverse transcriptase, which function together to promote RNA splicing and intron mobility via reverse splicing of the intron RNA into new DNA sites ("retrohoming"). Although group II introns are active in bacteria, their natural hosts, they function inefficiently in eukaryotes, where lower free Mg(2+) concentrations decrease their ribozyme activity and constitute a natural barrier to group II intron proliferation within nuclear genomes. Here, we show that retrohoming of the Ll.LtrB group II intron is strongly inhibited in an Escherichia coli mutant lacking the Mg(2+) transporter MgtA, and we use this system to select mutations in catalytic core domain V (DV) that partially rescue retrohoming at low Mg(2+) concentrations. We thus identified mutations in the distal stem of DV that increase retrohoming efficiency in the MgtA mutant up to 22-fold. Biochemical assays of splicing and reverse splicing indicate that the mutations increase the fraction of intron RNA that folds into an active conformation at low Mg(2+) concentrations, and terbium-cleavage assays suggest that this increase is due to enhanced Mg(2+) binding to the distal stem of DV. Our findings indicate that DV is involved in a critical Mg(2+)-dependent RNA folding step in group II introns and demonstrate the feasibility of selecting intron variants that function more efficiently at low Mg(2+) concentrations, with implications for evolution and potential applications in gene targeting.

Mapping L1 ligase ribozyme conformational switch

L1 ligase (L1L) molecular switch is an in vitro optimized synthetic allosteric ribozyme that catalyzes the regioselective formation of a 5'-to-3' phosphodiester bond, a reaction for which there is no known naturally occurring RNA catalyst. L1L serves as a proof of principle that RNA can catalyze a critical reaction for prebiotic RNA self-replication according to the RNA world hypothesis. L1L crystal structure captures two distinct conformations that differ by a reorientation of one of the stems by around 80Å and are presumed to correspond to the active and inactive state, respectively. It is of great interest to understand the nature of these two states in solution and the pathway for their interconversion. In this study, we use explicit solvent molecular simulation together with a novel enhanced sampling method that utilizes concepts from network theory to map out the conformational transition between active and inactive states of L1L. We find that the overall switching mechanism can be described as a three-state/two-step process. The first step involves a large-amplitude swing that reorients stem C. The second step involves the allosteric activation of the catalytic site through distant contacts with stem C. Using a conformational space network representation of the L1L switch transition, it is shown that the connection between the three states follows different topographical patterns: the stem C swing step passes through a narrow region of the conformational space network, whereas the allosteric activation step covers a much wider region and a more diverse set of pathways through the network.

The structural role of Mg2+ ions in a class I RNA polymerase ribozyme: a molecular simulation study

According to the RNA world hypothesis, self-replicating ribozymes, storing the genetic information and being able to perform catalysis, were the constituents of the first living organisms. In particular, RNA polymerase ribozymes, similar to current proteinaceous enzymatic polymerases, may have been able to promote the synthesis of RNA strands in a primitive world. Polymerase catalysis is usually assisted by Mg(2+) ions, but it is not always trivial to find out experimentally the number of Mg(2+) ions placed in the active site as well as the identity and the number of their coordination ligands. Here, we addressed this issue in an artificial class I ligase ribozyme. On the basis of a recently solved crystal structure, we constructed computational models of reactant and product states of this ribozyme, considering monometallic and bimetallic species. Our models were relaxed by force field based molecular dynamics (MD) simulations and mixed quantum-classical (QM/MM) MD. The structural and dynamical properties of these models were consistent with experimental data and were validated by a comparison with the catalytic sites of proteinaceous DNA and RNA polymerases. Consistently with enzymatic polymerases, our results suggest that class I RNA ligases most probably contain two magnesium ions in the active site and they may, therefore, catalyze the junction of two RNA strands via "a two Mg(2+) ions" mechanism.

Crystal structure of an RNA polymerase ribozyme in complex with an antibody fragment

All models of the RNA world era invoke the presence of ribozymes that can catalyse RNA polymerization. The class I ligase ribozyme selected in vitro 15 years ago from a pool of random RNA sequences catalyses formation of a 3',5'-phosphodiester linkage analogous to a single step of RNA polymerization. Recently, the three-dimensional structure of the ligase was solved in complex with U1A RNA-binding protein and independently in complex with an antibody fragment. The RNA adopts a tripod arrangement and appears to use a two-metal ion mechanism similar to protein polymerases. Here, we discuss structural implications for engineering a true polymerase ribozyme and describe the use of the antibody framework both as a portable chaperone for crystallization of other RNAs and as a platform for exploring steps in evolution from the RNA world to the RNA-protein world.

The chemical origins of life and its early evolution: an introduction

Can we look at contemporary biology and couple this with chemical insight to propose some plausible mechanisms for the origin of life on the planet? In what follows, we examine some promising chemical reactions by which the building blocks for nucleic acids might have been created about a billion years after the Earth formed. This could have led to self-assembling systems that were based on an all-RNA metabolism, where RNA is both catalytic and informational. We consider the breadth of RNA enzymes presently existing in biology, and to what extent these might have covered a wider range of chemistry in the RNA world. Ultimately, the RNA world would probably have given way to protein-based life quite quickly, and the origins of peptidyl transferase activity are discussed below.

The structural basis of RNA-catalyzed RNA polymerization

Early life presumably required polymerase ribozymes capable of replicating RNA. Known polymerase ribozymes best approximating such replicases use as their catalytic engine an RNA-ligase ribozyme originally selected from random RNA sequences. Here, we report 3.15 Å crystal structures of this ligase trapped in catalytically viable pre-ligation states, with the 3′-hydroxyl nucleophile positioned for in-line attack on the 5′-triphosphate. Guided by metal and solvent-mediated interactions, the 5′-triphosphate hooks into the major groove of the adjoining RNA duplex in an unanticipated conformation. Two phosphates and the nucleophile jointly coordinate an active-site metal ion. Atomic mutagenesis experiments demonstrate that active-site nucleobase and hydroxyl groups also participate directly in catalysis, collectively playing a role that in proteinaceous polymerases is performed by a second metal ion. Thus artificial ribozymes can employ complex catalytic strategies that differ dramatically from those of analogous biological enzymes.

A portable RNA sequence whose recognition by a synthetic antibody facilitates structural determination

RNA crystallization and phasing represent major bottlenecks in RNA structure determination. Seeking to exploit antibody fragments as RNA crystallization chaperones, we have used an arginine-enriched synthetic Fab library displayed on phage to obtain Fabs against the class I ligase ribozyme. We solved the structure of a Fab:ligase complex at 3.1Å using molecular replacement with Fab coordinates, confirming the ribozyme architecture and revealing the chaperone’s role in RNA recognition and crystal contacts. The epitope resides in the GAAACAC sequence that caps the P5 helix and retains high-affinity Fab binding within the context of other structured RNAs. This portable epitope provides a new RNA crystallization chaperone system that easily can be screened in parallel to the U1A RNA-binding protein, with the advantages of the smaller size of the loop and high molecular weight, large surface area, and phasing power provided by Fabs.

Identification of dynamical hinge points of the L1 ligase molecular switch

The L1 ligase is an in vitro selected ribozyme that uses a noncanonically base-paired ligation site to catalyze regioselectively and regiospecifically the 5' to 3' phosphodiester bond ligation, a reaction relevant to origin of life hypotheses that invoke an RNA world scenario. The L1 ligase crystal structure revealed two different conformational states that were proposed to represent the active and inactive forms. It remains an open question as to what degree these two conformers persist as stable conformational intermediates in solution, and along what pathway are they able to interconvert. To explore these questions, we have performed a series of molecular dynamics simulations in explicit solvent of the inactive-active conformational switch in L1 ligase. Four simulations were performed departing from both conformers in both the reactant and product states, in addition to a simulation where local unfolding in the active state was induced. From these simulations, along with crystallographic data, a set of four virtual torsion angles that span two evolutionarily conserved and restricted regions were identified as dynamical hinge points in the conformational switch transition. The ligation site visits three distinct states characterized by hydrogen bond patterns that are correlated with the formation of specific contacts that may promote catalysis. The insights gained from these simulations contribute to a more detailed understanding of the coupled catalytic/conformational switch mechanism of L1 ligase that may facilitate the design and engineering of new catalytic riboswitches.

Quasispecies-like behavior observed in catalytic RNA populations evolving in a test tube

BACKGROUND

During the RNA World, molecular populations were probably very small and highly susceptible to the force of strong random drift. In conjunction with Muller's Ratchet, this would have imposed difficulties for the preservation of the genetic information and the survival of the populations. Mechanisms that allowed these nascent populations to overcome this problem must have been advantageous.

RESULTS

Using continuous in vitro evolution experimentation with an increased mutation rate imposed by MnCl2, it was found that clonal 100-molecule populations of ribozymes clearly exhibit certain characteristics of a quasispecies. This is the first time this has been seen with a catalytic RNA. Extensive genotypic sampling from two replicate lineages was gathered and phylogenetic networks were constructed to elucidate the structure of the evolving RNA populations. A common distribution was found in which a mutant sequence was present at high frequency, surrounded by a cloud of mutant with lower frequencies. This is a typical distribution of quasispecies. Most of the mutants in these clouds were connected by short Hamming distance values, indicating their close relatedness.

CONCLUSIONS

The quasispecies nature of mutant RNA clouds facilitates the recovery of genotypes under pressure of being removed from the population by random drift. The empirical populations therefore evolved a genotypic resiliency despite a high mutation rate by adopting the characteristics of quasispecies, implying that primordial RNA pools could have used this strategy to avoid extinction.

Exploring RNA structure by integrative molecular modelling

RNA molecular modelling is adequate to rapidly tackle the structure of RNA molecules. With new structured RNAs constituting a central class of cellular regulators discovered every year, the need for swift and reliable modelling methods is more crucial than ever. The pragmatic method based on interactive all-atom molecular modelling relies on the observation that specific structural motifs are recurrently found in RNA sequences. Once identified by a combination of comparative sequence analysis and biochemical data, the motifs composing the secondary structure of a given RNA can be extruded in three dimensions (3D) and used as building blocks assembled manually during a bioinformatic interactive process. Comparing the models to the corresponding crystal structures has validated the method as being powerful to predict the RNA topology and architecture while being less accurate regarding the prediction of base-base interactions. These aspects as well as the necessary steps towards automation will be discussed.