A toggle duplex in hepatitis delta virus self-cleaving RNA that stabilizes an inactive and a salt-dependent pro-active ribozyme conformation

Mapping the conformational landscape of the neutral network of RNA sequences that connect two functional distinctly different ribozymes

During evolution of an RNA world, the development of enzymatic function was essential. Such enzymatic function was linked to RNA sequences capable of adopting specific RNA folds that possess catalytic pockets to promote catalysis. Within this primordial RNA world, initially evolved self-replicating ribozymes presumably mutated to ribozymes with new functions. In 2000, Schultes and Bartel investigated such conversion from one ribozyme to a new ribozyme with distinctly different catalytic functions. Within a neutral network that linked these two prototype ribozymes, a single RNA chain could be identified that exhibited both enzymatic functions. This system serves as a paradigm for an evolutionary system that allows neutral drifts by stepwise mutation from one ribozyme into a different ribozyme without loss of intermittent function. Here, we investigated this complex functional diversification by analyzing several RNA sequences within this neutral network between two ribozymes with class III ligase activity and with self-cleavage reactivity. We utilized rapid RNA sample preparation for NMR spectroscopic studies together with SHAPE analysis and in-line probing to characterize secondary structure changes within the neutral network. Our investigations allowed delineation of the 2 nd structure space and by comparison with the previously determined catalytic function allowed correlation of the structure-function relation of ribozyme function in this neutral network.

RNA folding: structure prediction, folding kinetics and ion electrostatics

Beyond the "traditional" functions such as gene storage, transport and protein synthesis, recent discoveries reveal that RNAs have important "new" biological functions including the RNA silence and gene regulation of riboswitch. Such functions of noncoding RNAs are strongly coupled to the RNA structures and proper structure change, which naturally leads to the RNA folding problem including structure prediction and folding kinetics. Due to the polyanionic nature of RNAs, RNA folding structure, stability and kinetics are strongly coupled to the ion condition of solution. The main focus of this chapter is to review the recent progress in the three major aspects in RNA folding problem: structure prediction, folding kinetics and ion electrostatics. This chapter will introduce both the recent experimental and theoretical progress, while emphasize the theoretical modelling on the three aspects in RNA folding.

HDV family of self-cleaving ribozymes

The hepatitis delta virus (HDV) ribozymes are catalytic RNAs capable of cleaving their own sugar-phosphate backbone. The HDV virus possesses the ribozymes in both sense and antisense genomic transcripts, where they are essential for processing during replication. These ribozymes have been the subject of intense biochemical scrutiny and have yielded a wealth of mechanistic insights. In recent years, many HDV-like ribozymes have been identified in nearly all branches of life. The ribozymes are implicated in a variety of biological events, including episodic memory in mammals and retrotransposition in many eukaryotes. Detailed analysis of additional HDV-like ribozyme isolates will likely reveal many more biological functions and provide information about the evolution of this unique RNA.

Basin Hopping Graph: a computational framework to characterize RNA folding landscapes

MOTIVATION

RNA folding is a complicated kinetic process. The minimum free energy structure provides only a static view of the most stable conformational state of the system. It is insufficient to give detailed insights into the dynamic behavior of RNAs. A sufficiently sophisticated analysis of the folding free energy landscape, however, can provide the relevant information.

RESULTS

We introduce the Basin Hopping Graph (BHG) as a novel coarse-grained model of folding landscapes. Each vertex of the BHG is a local minimum, which represents the corresponding basin in the landscape. Its edges connect basins when the direct transitions between them are 'energetically favorable'. Edge weights endcode the corresponding saddle heights and thus measure the difficulties of these favorable transitions. BHGs can be approximated accurately and efficiently for RNA molecules well beyond the length range accessible to enumerative algorithms.

AVAILABILITY AND IMPLEMENTATION

The algorithms described here are implemented in C++ as standalone programs. Its source code and supplemental material can be freely downloaded from http://www.tbi.univie.ac.at/bhg.html.

Significantly improved rescue of rabies virus from cDNA plasmids

The rescue of recombinant rabies virus (RV) from cloned cDNA is an inefficient process because it relies on the de novo formation within cells of functional ribonucleoprotein (RNP) complexes from plasmid-expressed viral-like antigenome RNAs and three helper proteins. In the standard RV reverse genetics systems, bacteriophage T7 RNA polymerase drives the transcription of virus antigenome-like RNAs containing three nonviral G residues at the 5'-end and a correct 3'-end generated by the autocatalytic activity of an 85 nucleotides long hepatitis delta virus antigenomic "core" ribozyme (HDVagrz). Here, we show that employing optimized ribozyme sequences significantly improves RV rescue. Substitution of the "core" HDVagrz by a ribozyme with an enhanced cleavage activity resulted in an approximately 10-fold higher number of rescue events and faster initiation of an infectious cycle. The alternative use of a hammerhead ribozyme for the generation of an exact 5'-end similarly enhanced rescue efficiency. Notably, RV cDNA clones containing the combination of optimized 3'- and 5'-ribozymes were rescued at an at least 100-fold increase. In addition to virus rescue, reporter gene expression from transfected minigenome cDNAs was significantly enhanced by the novel ribozymes. The improved RV reverse genetics system greatly facilitates recovery of strongly attenuated viruses and vectors for biomedical applications.

BarMap: RNA folding on dynamic energy landscapes

Dynamical changes of RNA secondary structures play an important role in the function of many regulatory RNAs. Such kinetic effects, especially in time-variable and externally triggered systems, are usually investigated by means of extensive and expensive simulations of large sets of individual folding trajectories. Here we describe the theoretical foundations of a generic approach that not only allows the direct computation of approximate population densities but also reduces the efforts required to analyze the folding energy landscapes to a one-time preprocessing step. The basic idea is to consider the kinetics on individual landscapes and to model external triggers and environmental changes as small but discrete changes in the landscapes. A "barmap" links macrostates of temporally adjacent landscapes and defines the transfer of population densities from one "snapshot" to the next. Implemented in the BarMap software, this approach makes it feasible to study folding processes at the level of basins, saddle points, and barriers for many nonstationary scenarios, including temperature changes, cotranscriptional folding, refolding in consequence to degradation, and mechanically constrained kinetics, as in the case of the translocation of a polymer through a pore.

Predicting secondary structural folding kinetics for nucleic acids

We report a new computational approach to the prediction of RNA secondary structure folding kinetics. In this approach, each elementary kinetic step is represented as the transformation between two secondary structures that differ by a helix. Based on the free energy landscape analysis, we identify three types of dominant pathways and the rate constants for the kinetic steps: 1), formation; 2), disruption of a helix stem; and 3), helix formation with concomitant partial melting of a competing (incompatible) helix. The third pathway, termed the tunneling pathway, is the low-barrier dominant pathway for the conversion between two incompatible helices. Comparisons with experimental data indicate that this new method is quite reliable in predicting the kinetics for RNA secondary structural folding and structural rearrangements. The approach presented here may provide a robust first step for further systematic development of a predictive theory for the folding kinetics for large RNAs.

Folding kinetics of large RNAs

We introduce here a heuristic approach to kinetic RNA folding that constructs secondary structures by stepwise combination of building blocks. These blocks correspond to subsequences and their thermodynamically optimal structures. These are determined by the standard dynamic programming approach to RNA folding. Folding trajectories are modeled at base-pair resolution using the Morgan-Higgs heuristic and a barrier tree-based heuristic to connect combinations of the local building blocks. Implemented in the program Kinwalker, the algorithm allows co-transcriptional folding and can be used to fold sequences of up to about 1500 nucleotides in length. A detailed comparison with several well-studied examples from the literature, including the delayed folding of bacteriophage cloverleaf structures, the adenine sensing riboswitch, and the hok RNA, shows an excellent agreement of predicted trajectories and experimental evidence. The software is available as part of the ViennaRNA Package.

Long-range impact of peripheral joining elements on structure and function of the hepatitis delta virus ribozyme

The HDV ribozyme is an RNA enzyme from the human pathogenic hepatitis delta virus (HDV) that has recently also been identified in the human genome. It folds into a compact, nested double-pseudoknot. We examined here the functional relevance of the capping loop L4 and the helical crossover J1/2, which tightly interlace the two helical stacks of the ribozyme. Peripheral structural elements such as these are present in cis-acting, but not trans-acting ribozymes, which may explain the order-of-magnitude decrease in cleavage activity observed in trans-acting ribozymes with promise in gene therapy applications. Comparison of a systematic set of cis- and trans-acting HDV ribozymes shows that the absence of either L4 or J1/2 significantly and independently impacts catalytic activity. Using terbium(III) footprinting and affinity studies, as well as distance measurements based on time-resolved fluorescence resonance energy transfer, we find that J1/2 is most important for conferring structural properties similar to those of the cis-acting ribozyme. Our results are consistent with a model in which removal of either a helical crossover or surprisingly a capping loop induces greater dynamics and expansion of the catalytic core at long range, impacting local and global folding, as well as catalytic function.

HDV Ribozyme Activity in Monovalent Cations†

Activity of the two ribozymes from hepatitis delta virus in monovalent salts was examined and compared to activity in Mg2+. Both ribozymes self-cleaved in high concentrations of monovalent cations, and an active site cytosine was required for cleavage activity under those conditions. Cleavage rates were 30-50-fold higher for reactions in LiCl than for reactions in NaCl or NH4Cl, and a thio effect indicated that chemistry was rate-determining for cleavage of the HDV genomic ribozyme in LiCl. Still, in LiCl, there was a more than 100-fold increase in the rate when MgCl2 was included in the reaction. However, the pH-rate profiles for the reactions in LiCl with and without MgCl2 were both bell-shaped with the pH optima in the neutral range. These findings support the idea that monovalent cations can partially substitute for divalent metal ions in the HDV ribozymes, although a divalent metal ion is more effective in supporting catalysis. The absence of a dramatic change in the general shape of pH-rate profiles in LiCl, relative to the profile for reactions including Mg2+, is in contrast to earlier data for the reactions in NaCl and limits our interpretation of the specific role played by the divalent metal ion in the catalytic mechanism.

Molecular beacons as probes of RNA unfolding under native conditions

Hybridization of fluorescent molecular beacons provides real-time detection of RNA secondary structure with high specificity. We used molecular beacons to measure folding and unfolding rates of the Tetrahymena group I ribozyme under native conditions. A molecular beacon targeted against 15 nt in the 5' strand of the P3 helix specifically hybridized with misfolded forms of the ribozyme, without invading the native tertiary structure. The beacon associated with the misfolded ribozyme 300 times more slowly than with an unstructured oligonucleotide containing the same target sequence, suggesting that the misfolded ribozyme core remains structured in the absence of Mg2+. The rate of beacon hybridization under native conditions revealed a linear relationship between the free energy of unfolding and Mg2+ concentration. A small fraction of the RNA population unfolded very rapidly, suggesting parallel unfolding in one step or through misfolded intermediates.

Design of a highly reactive HDV ribozyme sequence uncovers facilitation of RNA folding by alternative pairings and physiological ionic strength

The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA that resides in the HDV genome and regulates its replication. The native fold of the ribozyme is complex, having two pseudoknots. Earlier work implicated four non-native pairings in slowing pseudoknot formation: Alt 1, Alt 2, Alt 3, and Alt P1. The goal of the present work was design of a kinetically simplified and maximally reactive construct for in vitro mechanistic and structural studies. The initial approach chosen was site-directed mutagenesis in which known alternative pairings were destabilized while leaving the catalytic core intact. Based on prior studies, the G11C/U27Delta double mutant was prepared. However, biphasic kinetics and antisense oligonucleotide response trends opposite those of the well-studied G11C mutant were observed suggesting that new alternative pairings with multiple registers, termed Alt X and Alt Y, had been created. Enzymatic structure mapping of oligonucleotide models supported this notion. This led to a model wherein Alt 2 and the phylogenetically conserved Alt 3 act as "folding guides", facilitating folding of the major population of the RNA molecules by hindering formation of the Alt X and Alt Y registers. Attempts to eliminate the strongest of the Alt X pairings by rational design of a quadruple mutant only resulted in more complex kinetic behavior. In an effort to simultaneously destabilize multiple alternative pairings, studies were carried out on G11C/U27Delta in the presence of urea or increased monovalent ion concentration. Inclusion of physiological ionic strength allowed the goal of monophasic, fast-folding (kobs approximately 60 min(-1)) kinetics to be realized. To account for this, a model is developed wherein Na+, which destabilizes secondary and tertiary structures in the presence of Mg2+, facilitates native folding by destabilizing the multiple alternative secondary structures with a higher-order dependence.

Architecture and folding mechanism of the Azoarcus Group I Pre-tRNA

Self-splicing RNAs must evolve to function in their specific exon context. The conformation of a group I pre-tRNA(ile) from the bacterium Azoarcus was probed by ribonuclease T(1) and hydroxyl radical cleavage, and by native gel electrophoresis. Biochemical data and three-dimensional models of the pre-tRNA showed that the tRNA is folded, and that the tRNA and intron sequences form separate tertiary domains. Models of the active site before steps 1 and 2 of the splicing reaction predict that exchange of the external G-cofactor and the 3'-terminal G is accomplished by a slight conformational change in P9.0 of the Azoarcus group I intron. Kinetic assays showed that the pre-tRNA folds in minutes, much more slowly than the intron alone. The dependence of the folding kinetics on Mg(2+) and the concentration of urea, and RNase T(1) experiments showed that formation of native pre-tRNA is delayed by misfolding of P3-P9, including mispairing between residues in P9 and the tRNA. Thus, although the intron and tRNA sequences form separate domains in the native pre-tRNA, their folding is coupled via metastable non-native base-pairs. This could help prevent premature processing of the 5' and 3' ends of unspliced pre-tRNA.

Catalytic strategies of the hepatitis delta virus ribozymes

The hepatitis delta virus (HDV) ribozymes are self-cleaving RNA sequences critical to the replication of a small RNA genome. A recently determined crystal structure together with biochemical and biophysical studies provides new insight into the possible catalytic mechanism of these ribozymes. The HDV ribozymes are examples of naturally occurring small ribozymes that catalyze cleavage of the RNA backbone with a rate enhancement of 10(6)- to 10(7)-fold over the uncatalyzed rate. To achieve this level of rate enhancement, the HDV ribozymes have been proposed to employ several catalytic strategies that include the use of metal ions, intrinsic binding energy, and a novel example of general acid-base catalysis with a cytosine side chain acting as a proton donor or acceptor.

Catalytic cleavage of cis- and trans-acting antigenomic delta ribozymes in the presence of various divalent metal ions

Catalytic activity of four structural variants of the antigenomic delta ribozyme, two cis- and two trans-acting, has been compared in the presence of selected divalent metal ions that effectively support catalysis. The ribozymes differ in regions that are not directly involved in formation of the ribozyme active site: the region immediately preceding the catalytic cleavage site, the P4 stem and a stretch of the viral RNA sequence extending the minimal ribozyme sequence at its 3'-terminus. The variants show high cleavage activity in the presence of Mg(2+), Ca(2+) and Mn(2+), lower with Co(2+) and Sr(2+) and some variants are also active with Cd(2+) and Zn(2+) ions. In the presence of a particular metal ion the ribozymes cleave, however with different initial rates, according to pseudo-first or higher order kinetics and to different final cleavage extents. On the other hand, relatively small differences are observed in the reactions induced by various metal ions. The cleavage of trans-acting ribozymes induced by Mg(2+) is partially inhibited in the presence of Na(+), spermidine and some other divalent metal ions. The inert Co(NH(3))(6)(3+) complex is unable to support catalysis, as reported earlier for the genomic ribozyme. The results are discussed in terms of the influence of structural elements peripheral to the ribozyme active site on its cleavage rate and efficiency as well as the role of metal ions in the cleavage mechanism. Some implications concerning further studies and possible applications of delta ribozymes are also considered.

Energetic contribution of non-essential 5' sequence to catalysis in a hepatitis delta virus ribozyme

Hepatitis delta virus (HDV) ribozymes employ multiple catalytic strategies to achieve overall rate enhancement of RNA cleavage. These strategies include general acid-base catalysis by a cytosine side chain and involvement of divalent metal ions. Here we used a trans-acting form of the antigenomic ribozyme to examine the contribution of the 5' sequence in the substrate to HDV ribozyme catalysis. The cleavage rate constants increased for substrates with 5' sequence alterations that reduced ground-state binding to the ribozyme. Quantitatively, a plot of activation free energy of chemical conversion versus Gibb's free energy of substrate binding revealed a linear relationship with a slope of -1. This relationship is consistent with a model in which components of the substrate immediately 5' to the cleavage site in the HDV ribozyme-substrate complex destabilize ground-state binding. The intrinsic binding energy derived from the ground-state destabilization could contribute up to 2 kcal/mol toward the total 8.5 kcal/mol reduction in activation free energy for RNA cleavage catalyzed by the HDV ribozyme.

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.

Involvement of a cytosine side chain in proton transfer in the rate-determining step of ribozyme self-cleavage

Ribozymes of hepatitis delta virus have been proposed to use an active-site cytosine as an acid-base catalyst in the self-cleavage reaction. In this study, we have examined the role of cytosine in more detail with the antigenomic ribozyme. Evidence that proton transfer in the rate-determining step involved cytosine 76 (C76) was obtained from examining cleavage activity of the wild-type and imidazole buffer-rescued C76-deleted (C76 Delta) ribozymes in D(2)O and H(2)O. In both reactions, a similar kinetic isotope effect and shift in the apparent pKa indicate that the buffer is functionally substituting for the side chain in proton transfer. Proton inventory of the wild-type reaction supported a mechanism of a single proton transfer at the transition state. This proton transfer step was further characterized by exogenous base rescue of a C76 Delta mutant with cytosine and imidazole analogues. For the imidazole analogues that rescued activity, the apparent pKa of the rescue reaction, measured under k(cat)/K(M) conditions, correlated with the pKa of the base. From these data a Brønsted coefficient (beta) of 0.51 was determined for the base-rescued reaction of C76 Delta. This value is consistent with that expected for proton transfer in the transition state. Together, these data provide strong support for a mechanism where an RNA side chain participates directly in general acid or general base catalysis of the wild-type ribozyme to facilitate RNA cleavage.

A pH-sensitive RNA tertiary interaction affects self-cleavage activity of the HDV ribozymes in the absence of added divalent metal ion

Self-cleavage of the genomic and antigenomic ribozymes from hepatitis delta virus (HDV) requires divalent cation for optimal activity. Recently, the HDV genomic ribozyme has been shown to be active in NaCl in the absence of added divalent metal ion at low pH (apparent pKa 5.7). However, we find that the antigenomic ribozyme is 100 to 1000-fold less active under similar conditions. With deletion of a three-nucleotide sequence (C41-A42-A43) unique to the genomic ribozyme, the rate constant for cleavage decreased substantially, while activity of the antigenomic ribozyme was enhanced by introducing a CAA sequence. From the crystal structure, it has been proposed that C41 in this sequence is protonated. To investigate a possible connection between activity at low pH and protonation of C41, mutations were made that were predicted to either eliminate protonation or alter the nature of the tertiary interaction upon protonation. In the absence of added Mg2+, these mutations reduced activity and eliminated the observed pH-rate dependence. Thermal denaturation studies revealed a pH-sensitive structural feature in the genomic ribozyme, while unfolding of the mutant ribozymes was pH-independent. We propose that, in the absence of added Mg2+, protonation of C41 contributes to enhanced activity of the HDV genomic ribozyme at low pH.

Kinetic scheme for intermolecular RNA cleavage by a ribozyme derived from hepatitis delta virus RNA

A minimal kinetic mechanism for a trans-acting ribozyme derived from the HDV antigenomic RNA self-cleaving element was established from steady-state, pre-steady-state, single-turnover, and binding kinetics. Rate constants for individual steps, including substrate binding and dissociation, cleavage, and product release and binding, were measured at 37 degrees C at pH 8.0 in 10 mM Mg(2+) using oligonucleotides as either substrates, noncleavable analogues or 3' product mimics. A substrate containing a normal 3',5'-linkage was cleaved with a first-order rate constant (k(2)) of 0.91 min(-)(1). The association rate constant for the substrate to the ribozyme (2.1 x 10(7) M(-)(1) min(-)(1)) was at the lower range of the expected value for RNA duplex formation, and the substrate dissociated with a rate constant (1.4 min(-)(1)) slightly faster than that for cleavage. Thus the binary complex was not at equilibrium with free enzyme and substrate prior to the cleavage step. Following cleavage, product release was kinetically ordered in that the 5' product was released rapidly (>12 min(-)(1)) relative to the 3' product (6.0 x 10(-)(3) min(-)(1)). Rapid 5' product release and lack of a demonstrable binding site for the 5' product could contribute to the difficulty in establishing the ribozyme-catalyzed reverse reaction (ligation). Slow release of the 3' product was consistent with the extremely low turnover under steady-state conditions as 3' product dissociation was rate-limiting. The equilibrium dissociation constant for the substrate was 24-fold higher than that of the 3' cleavage product. A substrate with a 2',5'-linkage at the cleavage site was cleaved with a rate constant (k(2)) of 1.1 x 10(-)(2) min(-)(1). Thus, whereas cleavage of a 3',5'-linkage followed a Briggs-Haldane mechanism, 2', 5' cleavage followed a Michaelis-Menten mechanism.

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.

Crystallization and structure determination of a hepatitis delta virus ribozyme: use of the RNA-binding protein U1A as a crystallization module

Well-ordered crystals of a genomic hepatitis delta virus (HDV) ribozyme, a large, globular RNA, were obtained employing a new crystallization method. A high-affinity binding site for the spliceosomal protein U1A was engineered into a segment of the catalytic RNA that is dispensable for catalysis. Because molecular surfaces of proteins are more chemically varied than those of RNA, the presence of the protein moiety was expected to facilitate crystallization and improve crystal order. The HDV ribozyme-U1A complex crystallized readily, and its structure was solved using standard techniques for heavy-atom derivatization of protein crystals. Over 1200 A(2) of the solvent-accessible surface area of the complex are involved in crystal contacts. As protein-protein interactions comprise 85% of this buried area, these crystals appear to be held together predominantly by the protein component of the complex. Our crystallization method should be useful for the structure determination of other biochemically important RNAs for which protein partners do not exist or are experimentally intractable. The refined model of the complex (R-free=27.9% for all reflections between 20.0 and 2.3 A) reveals an RNA with a deep active site cleft. Well-ordered metal ions are not observed crystallographically in this cavity. Biochemical results of previous workers had suggested an important role in catalysis for cytosine 75. The pyrimidine base of this residue is buried at the bottom of the active site in an environment that could raise its pK(a) value. We propose that this highly conserved cytosine may be the general base that catalyzes the transesterification.

Satellite cereal yellow dwarf virus-RPV (satRPV) RNA requires a douXble hammerhead for self-cleavage and an alternative structure for replication

The 110 nt hammerhead ribozyme in the satellite RNA of cereal yellow dwarf virus-RPV (satRPV RNA) folds into an alternative conformation that inhibits self-cleavage. This alternative structure comprises a pseudoknot with base-pairing between loop (L1) and a single-stranded bulge (L2a), which are located in hammerhead stems I and II, respectively. Mutations that disrupt this base-pairing, or otherwise cause the ribozyme to more closely resemble a canonical hammerhead, greatly increase self-cleavage. In a more natural multimeric sequence context containing the full-length satRPV RNA and two copies of the hammerhead, wild-type RNA cleaves much more efficiently than in the 110 nt context. Mutations in the upstream hammerhead, including a knock-out in the catalytic core, affect cleavage at the downstream cleavage site, indicating that multimers of satRPV RNA cleave via a double hammerhead. The double hammerhead includes base-pairing between two copies of the L1 sequence which extends stem I. Disruption of L1-L1 base-pairing slows cleavage of the multimer. L1-L2a base-pairing is required for efficient replication of satRPV RNA in oat protoplasts. Mutations that affect self-cleavage of the multimer do not correlate with replication efficiency, indicating that the ability to self-cleave is not a primary determinant of replication. We present a replication model in which multimeric satRPV RNA folds into alternative conformations that cannot form in the monomer. One potential metastable intermediate conformation involves L1-L2a base-pairing that may facilitate formation of the double hammerhead. However, we conclude that L1-L2a also performs some other essential function in the satRPV RNA replication cycle, because the L1-L2a base-pairing is more important than efficient self-cleavage for replication.

Imidazole rescue of a cytosine mutation in a self-cleaving ribozyme

Ribozymes use a number of the same catalytic strategies as protein enzymes. However, general base catalysis by a ribozyme has not been demonstrated. In the hepatitis delta virus antigenomic ribozyme, imidazole buffer rescued activity of a mutant with a cytosine-76 (C76) to uracil substitution. In addition, a C76 to adenine substitution reduced the apparent pKa (where Ka is the acid constant) of the self-cleavage reaction by an amount consistent with differences in the pKa values of these two side chains. These results suggest that, in the wild-type ribozyme, C76 acts as a general base. This finding has implications for potential catalytic functions of conserved cytosines and adenines in other ribozymes and in ribonuclear proteins with enzymatic activity.

Facilitation of group I splicing in vivo: misfolding of the Tetrahymena IVS and the role of ribosomal RNA exons

Self-splicing of the group I IVS from Tetrahymena thermophila rDNA is limited by the time required for the RNA to reach its active conformation. In vitro, folding is slow because the pre-rRNA becomes kinetically trapped in inactive structures. In vivo, splicing is 50 times more rapid, implying that misfolding of the pre-rRNA is corrected. Exon mutations that inhibit self-splicing 100-fold in vitro were fully rescued when the pre-rRNA containing the IVS was expressed in E. coli. In contrast, IVS mutations that cause misfolding were only partially suppressed at 42 degrees C, and doubled the activation energy of splicing. These results suggest that intracellular folding of the pre-rRNA involves metastable intermediates similar to those observed in vitro. Precursors with natural rRNA exons were more active and less cold-sensitive than those with non-rRNA exons. This shows that the rRNA reduces misfolding of the IVS, thereby facilitating splicing of the pre-rRNA in vivo.

Sequential folding of the genomic ribozyme of the hepatitis delta virus: structural analysis of RNA transcription intermediates

The structures of the model oligoribonucleotides that mimic the consecutive stages in the transcription of genomic HDV ribozyme have been analyzed by the Pb(2+)-induced cleavage method, partial digestion with specific nucleases and chemical probing. In the transcription intermediates, the P1 and P4 helical segments are found to be present in the final folded forms in which they exist in the full-length transcript. However, the region corresponding to the central hairpin forms another thermodynamically stable hairpin structure. Its correct folding requires the presence of a ribozyme 3'-terminal sequence and the formation of helix P2. This confirms the ribozyme structure of the pseudoknot type and points to the crucial role of helix P2 in its overall folding. Moreover, we show that the J4/2 region can be specifically cleaved in the presence of selected divalent metal ions in the full-length transcript, but not in a shorter one lacking six 3'-terminal nucleotides, which cannot form the pseudoknotted structure. Thus, a particular RNA conformation around that cleavage site is required for specific hydrolysis, and the J4/2 region seems to be involved in the formation of a general metal ion binding site. Recently, it has been proposed that, in the antigenomic ribozyme, a four nucleotide sequence within the J1/2 region may contribute to the folding pathway, being part of a mechanism responsible for controlling ribozyme cleavage activity. Our study shows that in the genomic ribozyme the central hairpin region may contribute to a similar mechanism, providing a barrier to the formation of an active structure in the ribozyme folding pathway.

A conserved bulged adenosine in a peripheral duplex of the antigenomic HDV self-cleaving RNA reduceskinetic trapping of inactive conformations

In the ribozyme of hepatitis delta virus antigenomic RNA, two short duplexes, P2 and P2a, stabilize the active self-cleaving structure. However, P2a also promotes kinetic trapping of non-native structures. A bulged adenosine (A14) separates P2a and P2; this bulged A is conserved in clinical isolates of HDV but is unlikely to be physically close to the cleavage site phosphate in the ribozyme structure. Removing the bulge did not significantly slow the rate of cleavage but slowed the conversion of inactive to active conformations. In the absence of the bulged A, inactive conformations required higher urea concentrations or higher temperatures to be activated. Thus, the bulged-nucleotide in the P2-P2a duplex did not provide an essential kink or hinge between P2 and P2a that was required for cleavage activity but, rather, increased the rate of refolding from an inactive to an active ribozyme structure. These data also suggest a model in which P2 and P2a form a coaxial stacked helix of 9 bp, the most likely arrangement being one in which P2-P2a is roughly parallel to P1.