Developing three-dimensional models of putative-folding intermediates of the HDV ribozyme

Plasmid template design and in vitro transcription of short RNAs within a "structure cassette" for structure probing experiments

Chemical and enzymatic RNA structure probing methods are important tools for examining RNA secondary and tertiary structures and their interactions with proteins, small molecules, and ions. The recently developed "Selective 2'-Hydroxyl Acylation Analyzed by Primer Extension" (SHAPE) approach has proven especially useful for uncovering details of secondary structures, complex tertiary interactions, and RNA dynamics. Analysis of short RNAs using SHAPE or other probing methods that require reverse transcription to detect RNA modifications presents a technical hurdle in that intense bands corresponding to abortive transcription during primer extension and the full-length RT product may obscure information corresponding to the 3' and 5' ends of the molecule, respectively. This chapter describes the design and use of an RNA "structure cassette" that addresses these issues. First, we describe methods by which any RNA of interest may be cloned into a new plasmid preloaded with sequences that encode a structure cassette surrounding the short internal target RNA. Next, we outline key considerations and analyses of the RNAs produced that should be performed prior to SHAPE or other structure probing experiments.

Programming a highly structured ribozyme into complex allostery using RNA oligonucleotides

RNA possesses great potential for expanding the toolbox currently available to synthetic biologists. Here, the modulation of the Hepatitis Delta Virus ribozyme's activity with a series of rationally designed aptamers and effector RNA oligonucleotides is described. The ribozyme was initially fused with an 18-nucleotide hairpin structure that abolished its self-cleaving activity. The binding of a 14-mer oligonucleotide to the hairpin rescued the self-cleavage in a concentration-dependent manner. This modified ribozyme was inserted into the 5' UTR of a reporter gene, and the resulting construct was used to demonstrate that it is possible to modulate the ribozyme activity in cellulo with the oligonucleotide. Subsequently, ribozymes possessing specialized aptamers respecting other logic gates were also successfully designed and found to be functional in vitro. To our knowledge, this is the first example of HDV ribozyme regulation by oligonucleotides, as well as the first allosteric regulation of HDV ribozyme in mammalian cells.

Characterization of the trans Watson-Crick GU base pair located in the catalytic core of the antigenomic HDV ribozyme

The HDV ribozyme’s folding pathway is, by far, the most complex folding pathway elucidated to date for a small ribozyme. It includes 6 different steps that have been shown to occur before the chemical cleavage. It is likely that other steps remain to be discovered. One of the most critical of these unknown steps is the formation of the trans Watson-Crick GU base pair within loop III. The U₂₃ and G₂₈ nucleotides that form this base pair are perfectly conserved in all natural variants of the HDV ribozyme, and therefore are considered as being part of the signature of HDV-like ribozymes. Both the formation and the transformation of this base pair have been studied mainly by crystal structure and by molecular dynamic simulations. In order to obtain physical support for the formation of this base pair in solution, a set of experiments, including direct mutagenesis, the site-specific substitution of chemical groups, kinetic studies, chemical probing and magnesium-induced cleavage, were performed with the specific goal of characterizing this trans Watson-Crick GU base pair in an antigenomic HDV ribozyme. Both U₂₃ and G₂₈ can be substituted for nucleotides that likely preserve some of the H-bond interactions present before and after the cleavage step. The formation of the more stable trans Watson-Crick base pair is shown to be a post-cleavage event, while a possibly weaker trans Watson-Crick/Hoogsteen interaction seems to form before the cleavage step. The formation of this unusually stable post-cleavage base pair may act as a driving force on the chemical cleavage by favouring the formation of a more stable ground state of the product-ribozyme complex. To our knowledge, this represents the first demonstration of a potential stabilising role of a post-cleavage conformational switch event in a ribozyme-catalyzed reaction.

The mechanisms of RNA SHAPE chemistry

The biological functions of RNA are ultimately governed by the local environment at each nucleotide. Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistry is a powerful approach for measuring nucleotide structure and dynamics in diverse biological environments. SHAPE reagents acylate the 2'-hydroxyl group at flexible nucleotides because unconstrained nucleotides preferentially sample rare conformations that enhance the nucleophilicity of the 2'-hydroxyl. The critical corollary is that some constrained nucleotides must be poised for efficient reaction at the 2'-hydroxyl group. To identify such nucleotides, we performed SHAPE on intact crystals of the Escherichia coli ribosome, monitored the reactivity of 1490 nucleotides in 16S rRNA, and examined those nucleotides that were hyper-reactive toward SHAPE and had well-defined crystallographic conformations. Analysis of these conformations revealed that 2'-hydroxyl reactivity is broadly facilitated by general base catalysis involving multiple RNA functional groups and by two specific orientations of the bridging 3'-phosphate group. Nucleotide analog studies confirmed the contributions of these mechanisms to SHAPE reactivity. These results provide a strong mechanistic explanation for the relationship between SHAPE reactivity and local RNA dynamics and will facilitate interpretation of SHAPE information in the many technologies that make use of this chemistry.

Determining RNA three-dimensional structures using low-resolution data

Knowing the 3-D structure of an RNA is fundamental to understand its biological function. Nowadays X-ray crystallography and NMR spectroscopy are systematically applied to newly discovered RNAs. However, the application of these high-resolution techniques is not always possible, and thus scientists must turn to lower resolution alternatives. Here, we introduce a pipeline to systematically generate atomic resolution 3-D structures that are consistent with low-resolution data sets. We compare and evaluate the discriminative power of a number of low-resolution experimental techniques to reproduce the structure of the Escherichia coli tRNA(VAL) and P4-P6 domain of the Tetrahymena thermophila group I intron. We test single and combinations of the most accessible low-resolution techniques, i.e. hydroxyl radical footprinting (OH), methidiumpropyl-EDTA (MPE), multiplexed hydroxyl radical cleavage (MOHCA), and small-angle X-ray scattering (SAXS). We show that OH-derived constraints are accurate to discriminate structures at the atomic level, whereas EDTA-based constraints apply to global shape determination. We provide a guide for choosing which experimental techniques or combination of thereof is best in which context. The pipeline represents an important step towards high-throughput low-resolution RNA structure determination.

Dynamic energy landscapes of riboswitches help interpret conformational rearrangements and function

Riboswitches are RNAs that modulate gene expression by ligand-induced conformational changes. However, the way in which sequence dictates alternative folding pathways of gene regulation remains unclear. In this study, we compute energy landscapes, which describe the accessible secondary structures for a range of sequence lengths, to analyze the transcriptional process as a given sequence elongates to full length. In line with experimental evidence, we find that most riboswitch landscapes can be characterized by three broad classes as a function of sequence length in terms of the distribution and barrier type of the conformational clusters: low-barrier landscape with an ensemble of different conformations in equilibrium before encountering a substrate; barrier-free landscape in which a direct, dominant “downhill” pathway to the minimum free energy structure is apparent; and a barrier-dominated landscape with two isolated conformational states, each associated with a different biological function. Sharing concepts with the “new view” of protein folding energy landscapes, we term the three sequence ranges above as the sensing, downhill folding, and functional windows, respectively. We find that these energy landscape patterns are conserved in various riboswitch classes, though the order of the windows may vary. In fact, the order of the three windows suggests either kinetic or thermodynamic control of ligand binding. These findings help understand riboswitch structure/function relationships and open new avenues to riboswitch design.

Riboswitches are RNAs that modulate gene expression by ligand-induced conformational changes. However, the way that sequence dictates alternative folding pathways of gene regulation remains unclear. In this study, we mimic transcription by computing energy landscapes which describe accessible secondary structures for a range of sequence lengths. Consistent with experimental evidence, we find that most riboswitch landscapes can be characterized by three broad classes as a function of sequence length in terms of the distribution and barrier type of the conformational clusters: Low-barrier landscape with an ensemble of conformations in equilibrium before encountering a substrate; barrier-free landscape with a dominant “downhill” pathway to the minimum free energy structure; and barrier-dominated landscape with two isolated conformational states with different functions. Sharing concepts with the “new view” of protein folding energy landscapes, we term the three sequence ranges above as the sensing, downhill folding, and functional windows, respectively. We find that these energy landscape patterns are conserved between riboswitch classes, though the order of the windows may vary. In fact, the order of the three windows suggests either kinetic or thermodynamic control of ligand binding. These findings help understand riboswitch structure/function relationships and open new avenues to riboswitch design.

Selection of the most potent specific on/off adaptor-hepatitis delta virus ribozymes for use in gene targeting

The Hepatitis Delta Virus (HDV) ribozyme, which is well adapted to the environment of the human cell, is an excellent candidate for the future development of gene-inactivation systems. On top of this, a new generation of HDV ribozymes now exists that benefits from the addition of a specific on/off adaptor (specifically the SOFA-HDV ribozymes) which greatly increases both the ribozyme's specificity and its cleavage activity. Unlike RNAi and hammerhead ribozymes, the designing of SOFA-HDV ribozymes to cleave, in trans, given RNA species has never been the object of a systematic optimization study, even with their recent use for the gene knockdown of various targets. This report aims at both improving and clarifying the design process of SOFA-HDV ribozymes. Both the ribozyme and the targeted RNA substrate were analyzed in order to provide new criteria that are useful in the selection of the most potent SOFA-HDV ribozymes. The crucial features present in both the ribozyme's biosensor and blocker, as well as at the target site, were identified and characterized. Simple rules were derived and tested using hepatitis C virus NS5B RNA as a model target. Overall, this method should promote the use of the SOFA-HDV ribozymes in a plethora of applications in both functional genomics and gene therapy.

Topological constraints: using RNA secondary structure to model 3D conformation, folding pathways, and dynamic adaptation

Accompanying recent advances in determining RNA secondary structure is the growing appreciation for the importance of relatively simple topological constraints, encoded at the secondary structure level, in defining the overall architecture, folding pathways, and dynamic adaptability of RNA. A new view is emerging in which tertiary interactions do not define RNA 3D structure, but rather, help select specific conformers from an already narrow, topologically pre-defined conformational distribution. Studies are providing fundamental insights into the nature of these topological constraints, how they are encoded by the RNA secondary structure, and how they interplay with other interactions, breathing new meaning to RNA secondary structure. New approaches have been developed that take advantage of topological constraints in determining RNA backbone conformation based on secondary structure, and a limited set of other, easily accessible constraints. Topological constraints are also providing a much-needed framework for rationalizing and describing RNA dynamics and structural adaptation. Finally, studies suggest that topological constraints may play important roles in steering RNA folding pathways. Here, we review recent advances in our understanding of topological constraints encoded by the RNA secondary structure.

Finding the path in an RNA folding landscape

In this issue of Structure, Reymond et al. (2010) combine molecular and computational biology approaches to provide structural details for intermediates in the folding pathway of the hepatitis delta virus ribozyme.