Molecular modelling of the GIR1 branching ribozyme gives new insight into evolution of structurally related ribozymes

Progress toward SHAPE Constrained Computational Prediction of Tertiary Interactions in RNA Structure

As more sequencing data accumulate and novel puzzling genetic regulations are discovered, the need for accurate automated modeling of RNA structure increases. RNA structure modeling from chemical probing experiments has made tremendous progress, however accurately predicting large RNA structures is still challenging for several reasons: RNA are inherently flexible and often adopt many energetically similar structures, which are not reliably distinguished by the available, incomplete thermodynamic model. Moreover, computationally, the problem is aggravated by the relevance of pseudoknots and non-canonical base pairs, which are hardly predicted efficiently. To identify nucleotides involved in pseudoknots and non-canonical interactions, we scrutinized the SHAPE reactivity of each nucleotide of the 188 nt long lariat-capping ribozyme under multiple conditions. Reactivities analyzed in the light of the X-ray structure were shown to report accurately the nucleotide status. Those that seemed paradoxical were rationalized by the nucleotide behavior along molecular dynamic simulations. We show that valuable information on intricate interactions can be deduced from probing with different reagents, and in the presence or absence of Mg2+. Furthermore, probing at increasing temperature was remarkably efficient at pointing to non-canonical interactions and pseudoknot pairings. The possibilities of following such strategies to inform structure modeling software are discussed.

Use of a Lariat Capping Ribozyme to Study Cap Function In Vivo

A lariat cap is a naturally occurring substitute of a conventional mRNA cap and is found in a particular genomic setting in a few eukaryotic microorganisms. It is installed by the lariat capping ribozyme acting in cis. In principle, any RNA molecule in any organism can be equipped with a lariat cap in vivo when expressed downstream of a lariat capping ribozyme. Lariat capping is thus a versatile tool for studying the importance of the 5' end structure of RNA molecules. In this chapter, we present protocols to validate the presence of the lariat cap and measure the efficiency of in vivo cleavage by the lariat capping ribozyme.

RNA structure through multidimensional chemical mapping

The discoveries of myriad non-coding RNA molecules, each transiting through multiple flexible states in cells or virions, present major challenges for structure determination. Advances in high-throughput chemical mapping give new routes for characterizing entire transcriptomes in vivo, but the resulting one-dimensional data generally remain too information-poor to allow accurate de novo structure determination. Multidimensional chemical mapping (MCM) methods seek to address this challenge. Mutate-and-map (M2), RNA interaction groups by mutational profiling (RING-MaP and MaP-2D analysis) and multiplexed •OH cleavage analysis (MOHCA) measure how the chemical reactivities of every nucleotide in an RNA molecule change in response to modifications at every other nucleotide. A growing body of in vitro blind tests and compensatory mutation/rescue experiments indicate that MCM methods give consistently accurate secondary structures and global tertiary structures for ribozymes, ribosomal domains and ligand-bound riboswitch aptamers up to 200 nucleotides in length. Importantly, MCM analyses provide detailed information on structurally heterogeneous RNA states, such as ligand-free riboswitches that are functionally important but difficult to resolve with other approaches. The sequencing requirements of currently available MCM protocols scale at least quadratically with RNA length, precluding general application to transcriptomes or viral genomes at present. We propose a modify-cross-link-map (MXM) expansion to overcome this and other current limitations to resolving the in vivo 'RNA structurome'.

RNA-Puzzles Round II: assessment of RNA structure prediction programs applied to three large RNA structures

This paper is a report of a second round of RNA-Puzzles, a collective and blind experiment in three-dimensional (3D) RNA structure prediction. Three puzzles, Puzzles 5, 6, and 10, represented sequences of three large RNA structures with limited or no homology with previously solved RNA molecules. A lariat-capping ribozyme, as well as riboswitches complexed to adenosylcobalamin and tRNA, were predicted by seven groups using RNAComposer, ModeRNA/SimRNA, Vfold, Rosetta, DMD, MC-Fold, 3dRNA, and AMBER refinement. Some groups derived models using data from state-of-the-art chemical-mapping methods (SHAPE, DMS, CMCT, and mutate-and-map). The comparisons between the predictions and the three subsequently released crystallographic structures, solved at diffraction resolutions of 2.5-3.2 Å, were carried out automatically using various sets of quality indicators. The comparisons clearly demonstrate the state of present-day de novo prediction abilities as well as the limitations of these state-of-the-art methods. All of the best prediction models have similar topologies to the native structures, which suggests that computational methods for RNA structure prediction can already provide useful structural information for biological problems. However, the prediction accuracy for non-Watson-Crick interactions, key to proper folding of RNAs, is low and some predicted models had high Clash Scores. These two difficulties point to some of the continuing bottlenecks in RNA structure prediction. All submitted models are available for download at http://ahsoka.u-strasbg.fr/rnapuzzles/.

Molecular characterization of a new member of the lariat capping twin-ribozyme introns

BACKGROUND

Twin-ribozyme introns represent a complex class of mobile group I introns that harbour a lariat capping (LC) ribozyme and a homing endonuclease gene embedded in a conventional self-splicing group I ribozyme (GIR2). Twin-ribozyme introns have so far been confined to nucleolar DNA in Naegleria amoeboflagellates and the myxomycete Didymium iridis.

RESULTS

We characterize structural organization, catalytic properties and molecular evolution of a new twin-ribozyme intron in Allovahlkampfia (Heterolobosea). The intron contains two ribozyme domains with different functions in ribosomal RNA splicing and homing endonuclease mRNA maturation. We found Allovahlkampfia GIR2 to be a typical group IC1 splicing ribozyme responsible for addition of the exogenous guanosine cofactor (exoG), exon ligation and circularization of intron RNA. The Allovahlkampfia LC ribozyme, by contrast, represents an efficient self-cleaving ribozyme that generates a small 2',5' lariat cap at the 5' end of the homing endonuclease mRNA, and thus contributes to intron mobility.

CONCLUSIONS

The discovery of a twin-ribozyme intron in a member of Heterolobosea expands the distribution pattern of LC ribozymes. We identify a putative regulatory RNA element (AP2.1) in the Allovahlkampfia LC ribozyme that involves homing endonuclease mRNA coding sequences as an important structural component.

Speciation of a group I intron into a lariat capping ribozyme

The lariat-capping (LC) ribozyme is a natural ribozyme isolated from eukaryotic microorganisms. Despite apparent structural similarity to group I introns, the LC ribozyme catalyzes cleavage by a 2',5' branching reaction, leaving the 3' product with a 3-nt lariat cap that functionally substitutes for a conventional mRNA cap in the downstream pre-mRNA encoding a homing endonuclease. We describe the crystal structures of the precleavage and postcleavage LC ribozymes, which suggest that structural features inherited from group I ribozymes have undergone speciation due to profound changes in molecular selection pressure, ultimately giving rise to an original branching ribozyme family. The structures elucidate the role of key elements that regulate the activity of the LC ribozyme by conformational switching and suggest a mechanism by which the signal for branching is transmitted to the catalytic core. The structures also show how conserved interactions twist residues, forming the lariat to join chemical groups involved in branching.

Nuclear group I introns in self-splicing and beyond

Group I introns are a distinct class of RNA self-splicing introns with an ancient origin. All known group I introns present in eukaryote nuclei interrupt functional ribosomal RNA genes located in ribosomal DNA loci. The discovery of the Tetrahymena intron more than 30 years ago has been essential to our understanding of group I intron catalysis, higher-order RNA structure, and RNA folding, but other intron models have provided information about the biological role. Nuclear group I introns appear widespread among eukaryotic microorganisms, and the plasmodial slime molds (myxomycetes) contain an abundance of self-splicing introns. Here, we summarize the main conclusions from previous work on the Tetrahymena intron on RNA self-splicing catalysis as well as more recent work on myxomycete intron biology. Group I introns in myxomycetes that represent different evolutionary stages, biological roles, and functional settings are discussed.

Group I intron ribozymes

Group I intron ribozymes constitute one of the main classes of ribozymes and have been a particularly important model in the discovery of key concepts in RNA biology as well as in the development of new methods. Compared to other ribozyme classes, group I intron ribozymes display considerable variation both in their structure and the reactions they catalyze. The best described pathway is the splicing pathway that results in a spliced out intron and ligated exons. This is paralleled by the circularization pathway that leads to full-length circular intron and un-ligated exons. In addition, the intronic products of these pathways have the potential to integrate into targets and to form various types of circular RNA molecules. Thus, group I intron ribozymes and associated elements found within group I introns is a rich source of biological phenomena. This chapter provides a strategy and protocols for initial characterization of new group I intron ribozymes.

The role of robustness in phenotypic adaptation and innovation

Phenotypes that vary in response to DNA mutations are essential for evolutionary adaptation and innovation. Therefore, it seems that robustness, a lack of phenotypic variability, must hinder adaptation. The main purpose of this review is to show why this is not necessarily correct. There are two reasons. The first is that robustness causes the existence of genotype networks--large connected sets of genotypes with the same phenotype. I discuss why genotype networks facilitate phenotypic variability. The second reason emerges from the evolutionary dynamics of evolving populations on genotype networks. I discuss how these dynamics can render highly robust phenotypes more variable, using examples from protein and RNA macromolecules. In addition, robustness can help avoid an important evolutionary conflict between the interests of individuals and populations-a conflict that can impede evolutionary adaptation.

RNase P: at last, the key finds its lock

Apart from the ribosome, the crystal structure of the bacterial RNase P in complex with a tRNA, reported by Reiter and colleagues recently, constitutes the first example of a multiple turnover RNA enzyme. Except in rare exceptions, RNase P is ubiquitous and, like the ribosome, is older than the initial branch point of the phylogenetic tree. Importantly, the structure shows how the RNA and the protein moieties cooperate to process the pre-tRNA substrates. The catalytic site comprises some critical RNA residues spread over the secondary structure but gathered in a compact volume next to the protein, which helps recognize and orient the substrate. The discussion here outlines some important aspects of that crystal structure, some of which could apply to RNA molecules in general.

Intermolecular interaction between a branching ribozyme and associated homing endonuclease mRNA

RNA tertiary interactions involving docking of GNRA (N; any base; R; purine) hairpin loops into helical stem structures on other regions of the same RNA are one of the most common RNA tertiary interactions. In this study, we investigated a tertiary association between a GAAA hairpin tetraloop in a small branching ribozyme (DiGIR1) and a receptor motif (HEG P1 motif) present in a hairpin structure on a separate mRNA molecule. DiGIR1 generates a 2', 5' lariat cap at the 5' end of its downstream homing endonuclease mRNA by catalysing a self-cleavage branching reaction at an internal processing site. Upon release, the 5' end of the mRNA forms a distinct hairpin structure termed HEG P1. Our biochemical data, in concert with molecular 3D modelling, provide experimental support for an intermolecular tetraloop receptor interaction between the L9 GAAA in DiGIR1 and a GNRA tetraloop receptor-like motif (UCUAAG-CAAGA) found within the HEG P1. The biological role of this interaction appears to be linked to the homing endonuclease expression by promoting post-cleavage release of the lariat capped mRNA. These findings add to our understanding of how protein-coding genes embedded in nuclear ribosomal DNA are expressed in eukaryotes and controlled by ribozymes.

Dual-channel single-molecule fluorescence resonance energy transfer to establish distance parameters for RNA nanoparticles

The increasing interest in RNA nanotechnology and the demonstrated feasibility of using RNA nanoparticles as therapeutics have prompted the need for imaging systems with nanometer-scale resolution for RNA studies. Phi29 dimeric pRNAs can serve as building blocks in assembly into the hexameric ring of the nanomotors, as modules of RNA nanoparciles, and as vehicles for specific delivery of therapeutics to cancers or viral infected cells. The understanding of the 3D structure of this novel RNA dimeric particle is fundamentally and practically important. Although a 3D model of pRNA dimer has been proposed based on biochemical analysis, no distance measurements or X-ray diffraction data have been reported. Here we evaluated the application of our customized single-molecule dual-viewing system for distance measurement within pRNA dimers using single-molecule Fluorescence Resonance Energy Transfer (smFRET). Ten pRNA monomers labeled with single donor or acceptor fluorophores at various locations were constructed and eight dimers were assembled. smFRET signals were detected for six dimers. The tethered arm sizes of the fluorophores were estimated empirically from dual-labeled RNA/DNA standards. The distances between donor and acceptor were calculated and used as distance parameters to assess and refine the previously reported 3D model of the pRNA dimer. Distances between nucleotides in pRNA dimers were found to be different from those of the dimers bound to procapsid, suggesting a conformational change of the pRNA dimer upon binding to the procapsid.

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.

RNA reprogramming and repair based on trans-splicing group I ribozymes

While many traditional gene therapy strategies attempt to deliver new copies of wild-type genes back to cells harboring the defective genes, RNA-directed strategies offer a range of novel therapeutic applications. Revision or reprogramming of mRNA is a form of gene therapy that modifies mRNA without directly changing the transcriptional regulation or the genomic gene sequence. Group I ribozymes can be engineered to act in trans by recognizing a separate RNA molecule in a sequence-specific manner, and to covalently link a new RNA sequence to this separate RNA molecule. Group I ribozymes have been shown to repair defective transcripts that cause human genetic or malignant diseases, as well as to replace transcript sequences by foreign RNA resulting in new cellular functions. This review provides an overview of current strategies using trans-splicing group I ribozymes in RNA repair and reprogramming.

A conformational switch in the DiGIR1 ribozyme involved in release and folding of the downstream I-DirI mRNA

DiGIR1 is a group I-like cleavage ribozyme found as a structural domain within a nuclear twin-ribozyme group I intron. DiGIR1 catalyzes cleavage by branching at an Internal Processing Site (IPS) leading to formation of a lariat cap at the 5'-end of the 3'-cleavage product. The 3'-cleavage product is subsequently processed into an mRNA encoding a homing endonuclease. By analysis of combinations of 5'- and 3'-deletions, we identify a hairpin in the 5'-UTR of the mRNA (HEG P1) that is formed by conformational switching following cleavage. The formation of HEG P1 inhibits the reversal of the branching reaction, thus giving it directionality. Furthermore, the release of the mRNA is a consequence of branching rather than hydrolytic cleavage. A model is put forward that explains the release of the I-DirI mRNA with a lariat cap and a structured 5'-UTR as a direct consequence of the DiGIR1 branching reaction. The role of HEG P1 in GIR1 branching is reminiscent of that of hairpin P-1 in splicing of the Tetrahymena rRNA group I intron and illustrates a general principle in RNA-directed RNA processing.

A natural ribozyme with 3',5' RNA ligase activity

Using electrophoresis, sequencing and enzymatic digestion, we show that the group I intron from the cyanobacterium Anabaena sp. PCC 7120 catalyzes phosphodiester bond formation using a triphosphate on the 5'-terminal nucleotide, much like protein polymerases and engineered ribozymes. In the process, this ribozyme forms a unique circular RNA that incorporates the exogenous guanosine cofactor added during self-splicing. This finding may have relevance to a prebiotic RNA world and to modern biology.

A structural basis for the recognition of 2'-deoxyguanosine by the purine riboswitch

Riboswitches are noncoding RNA elements that are commonly found in the 5'-untranslated region of bacterial mRNA. Binding of a small-molecule metabolite to the riboswitch aptamer domain guides the folding of the downstream sequence into one of two mutually exclusive secondary structures that directs gene expression. The purine riboswitch family, which regulates aspects of purine biosynthesis and transport, contains three distinct classes that specifically recognize guanine/hypoxanthine, adenine, or 2'-deoxyguanosine (dG). Structural analysis of the guanine and adenine classes revealed a binding pocket that almost completely buries the nucleobase within the core of the folded RNA. Thus, it is somewhat surprising that this family of RNA elements also recognizes dG. We have used a combination of structural and biochemical techniques to understand how the guanine riboswitch could be converted into a dG binder and the structural basis for dG recognition. These studies reveal that a limited number of sequence changes to a guanine-sensing RNA are required to cause a specificity switch from guanine to 2'-deoxyguanosine, and to impart an altered structure for accommodating the additional deoxyribose sugar moiety.

Toward predicting self-splicing and protein-facilitated splicing of group I introns

In the current era of massive discoveries of noncoding RNAs within genomes, being able to infer a function from a nucleotide sequence is of paramount interest. Although studies of individual group I introns have identified self-splicing and nonself-splicing examples, there is no overall understanding of the prevalence of self-splicing or the factors that determine it among the >2300 group I introns sequenced to date. Here, the self-splicing activities of 12 group I introns from various organisms were assayed under six reaction conditions that had been shown previously to promote RNA catalysis for different RNAs. Besides revealing that assessing self-splicing under only one condition can be misleading, this survey emphasizes that in vitro self-splicing efficiency is correlated with the GC content of the intron (>35% GC was generally conductive to self-splicing), and with the ability of the introns to form particular tertiary interactions. Addition of the Neurospora crassa CYT-18 protein activated splicing of two nonself-splicing introns, but inhibited the second step of self-splicing for two others. Together, correlations between sequence, predicted structure and splicing begin to establish rules that should facilitate our ability to predict the self-splicing activity of any group I intron from its sequence.