Selection in vitro of trans-acting genomic human hepatitis delta virus (HDV) ribozymes

Wobble pairs of the HDV ribozyme play specific roles in stabilization of active site dynamics

The hepatitis delta virus (HDV) is the only known human pathogen whose genome contains a catalytic RNA motif (ribozyme). The overall architecture of the HDV ribozyme is that of a double-nested pseudoknot, with two GU pairs flanking the active site. Although extensive studies have shown that mutation of either wobble results in decreased catalytic activity, little work has focused on linking these mutations to specific structural effects on catalytic fitness. Here we use molecular dynamics simulations based on an activated structure to probe the active site dynamics as a result of wobble pair mutations. In both wild-type and mutant ribozymes, the in-line fitness of the active site (as a measure of catalytic proficiency) strongly depends on the presence of a C75(N3H3+)N1(O5') hydrogen bond, which positions C75 as the general acid for the reaction. Our mutational analyses show that each GU wobble supports catalytically fit conformations in distinct ways; the reverse G25U20 wobble promotes high in-line fitness, high occupancy of the C75(N3H3+)G1(O5') general-acid hydrogen bond and stabilization of the G1U37 wobble, while the G1U37 wobble acts more locally by stabilizing high in-line fitness and the C75(N3H3+)G1(O5') hydrogen bond. We also find that stable type I A-minor and P1.1 hydrogen bonding above and below the active site, respectively, prevent local structural disorder from spreading and disrupting global conformation. Taken together, our results define specific, often redundant architectural roles for several structural motifs of the HDV ribozyme active site, expanding the known roles of these motifs within all HDV-like ribozymes and other structured RNAs.

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.

Evolution of the R2 retrotransposon ribozyme and its self-cleavage site

R2 is a non-long terminal repeat retrotransposon that inserts site-specifically in the tandem 28S rRNA genes of many animals. Previously, R2 RNA from various species of Drosophila was shown to self-cleave from the 28S rRNA/R2 co-transcript by a hepatitis D virus (HDV)-like ribozyme encoded at its 5' end. RNA cleavage was at the precise 5' junction of the element with the 28S gene. Here we report that RNAs encompassing the 5' ends of R2 elements from throughout its species range fold into HDV-like ribozymes. In vitro assays of RNA self-cleavage conducted in many R2 lineages confirmed activity. For many R2s, RNA self-cleavage was not at the 5' end of the element but at 28S rRNA sequences up to 36 nucleotides upstream of the junction. The location of cleavage correlated well with the types of endogenous R2 5' junctions from different species. R2 5' junctions were uniform for most R2s in which RNA cleavage was upstream in the rRNA sequences. The 28S sequences remaining on the first DNA strand synthesized during retrotransposition are postulated to anneal to the target site and uniformly prime second strand DNA synthesis. In species where RNA cleavage occurred at the R2 5' end, the 5' junctions were variable. This junction variation is postulated to result from the priming of second strand DNA synthesis by chance microhomologies between the target site and the first DNA strand. Finally, features of R2 ribozyme evolution, especially changes in cleavage site and convergence on the same active site sequences, are discussed.

HDV-like self-cleaving ribozymes

HDV ribozymes catalyze their own scission from the transcript during rolling circle replication of the hepatitis delta virus. In vitro selection of self-cleaving ribozymes from a human genomic library revealed an HDV-like ribozyme in the second intron of the human CPEB3 gene and recent results suggest that this RNA affects episodic memory performance. Bioinformatic searches based on the secondary structure of the HDV/CPEB3 fold yielded numerous functional ribozymes in a wide variety of organisms. Genomic mapping of these RNAs suggested several biological roles, one of which is the 5' processing of non-LTR retrotransposons. The family of HDV-like ribozymes thus continues to grow in numbers and biological importance.

Antigenomic delta ribozyme variants with mutations in the catalytic core obtained by the in vitro selection method

We have used the in vitro selection method to search for catalytically active variants of the antigenomic delta ribozyme with mutations in the regions that constitute the ribozyme active site: L3, J1/4 and J4/2. In the initial combinatorial library 16 nt positions were randomized and the library contained a full representation of all possible sequences. Following ten cycles of selection-amplification several catalytically active ribozyme variants were identified. It turned out that one-third of the variants contained only single mutation G80U and their activity was similar to that of the wild-type ribozyme. Unexpectedly, in the next one-third of the variants the C76 residue, which was proposed to play a crucial role in the ribozyme cleavage mechanism, was mutated. In these variants, however, a cytosine residue was present in a neighboring position to the polynucleotide chain. It shows that the ribozyme catalytic core possesses substantial 'structural plasticity' and the capacity of functional adaptation. Four selected ribozyme variants were subjected to more detailed analysis. It turned out that the variants differed in their relative preferences towards Mg2+, Ca2+ and Mn2+ ions. Thus, the functional properties of the variants were dependent on both the structure of their catalytic sites and divalent metal ions performing catalysis.

Unbiased in vitro selection reveals the unique character of the self-cleaving antigenomic HDV RNA sequence

In order to revisit the architecture of the catalytic center of the antigenomic hepatitis delta virus (HDV) ribozyme we developed an unbiased in vitro selection procedure that efficiently selected novel variants from a relatively small set of sequences. Using this procedure we examined all possible variants from a pool of HDV ribozymes that had been randomized at 25 positions (4²⁵). The isolated set of sequences shows more variability than do the natural variants. Nucleotide variations were found at all randomized positions, even at positions when the general belief was that the specific base was absolutely required for catalytic activity. Covariation analysis supports the presence of several base pairs, although it failed to propose any new tertiary contacts. HDV ribozyme appears to possess a greater number of constraints, in terms of sequences capable of supporting the catalysed cleavage, than do other catalytic RNAs. This supports the idea that the appearance of this catalytic RNA structure has a low probability (i.e. is a rare event), which may explain why to date it has been found in nature only in the HDV. These contrasts with the hammerhead self-cleaving motif that is proposed to have multiple origins, and that is widespread among different organisms. Thus, just because a self-cleaving RNA motif is small does not imply that it occurs easily.

Site-specific modification of functional groups in genomic hepatitis delta virus (HDV) ribozyme

Human hepatitis delta (HDV) ribozyme is one of small ribozymes, such as hammerhead and hairpin ribozymes, etc. Its secondary structure shows pseudoknot structure composed of four stems (I to IV) and three single-stranded regions (SSrA, -B and -C). The 3D structure of 3'-cleaved product of genomic HDV ribozyme provided extensive information about tertiary hydrogen bonding interactions between nucleotide bases, phosphate oxygens and 2'OHs including new stem structure P1.1. To analyze the role of these hydrogen bond networks in the catalytic reaction, site-specific atomic-level modifications (such as deoxynucleotides, deoxyribosyl-2-aminopurine, deoxyribosylpurine, 7-deaza-ribonucleotide and inosine) were incorporated in the smallest trans-acting HDV ribozyme (47-mer). Kinetic analysis of these ribozyme variants demonstrated the importance of the two W-C base pairs of P1.1 for cleavage; in addition, the results suggest that all hydrogen bond interactions detected in the crystal structure involving 2'-OH and N7 atoms are present in the active ribozyme structure. In most of the variants, the relative reduction in kobs caused by substitution of the 2'-OH group correlated with the number of hydrogen bonds affected by the substitution. However G74 and C75 may have more than one hydrogen bond involving the 2'-OH in both the trans- and cis-acting HDV ribozyme. Moreover, in variants in which N7 was deleted, kobs was reduced 5- to 15-fold, it may suggest that N7 assists in coordinating Mg2+ ions or water molecules which bind with weak affinity in the active structure.

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.

In vitro selection of functional nucleic acids

In vitro selection allows rare functional RNA or DNA molecules to be isolated from pools of over 10(15) different sequences. This approach has been used to identify RNA and DNA ligands for numerous small molecules, and recent three-dimensional structure solutions have revealed the basis for ligand recognition in several cases. By selecting high-affinity and -specificity nucleic acid ligands for proteins, promising new therapeutic and diagnostic reagents have been identified. Selection experiments have also been carried out to identify ribozymes that catalyze a variety of chemical transformations, including RNA cleavage, ligation, and synthesis, as well as alkylation and acyl-transfer reactions and N-glycosidic and peptide bond formation. The existence of such RNA enzymes supports the notion that ribozymes could have directed a primitive metabolism before the evolution of protein synthesis. New in vitro protein selection techniques should allow for a direct comparison of the frequency of ligand binding and catalytic structures in pools of random sequence polynucleotides versus polypeptides.

Structural diversity of self-cleaving ribozymes

In vitro selection was used to isolate Mg(2+)-dependent self-cleaving ribozymes from random sequence. Characterization of representative clones revealed the emergence of at least 12 classes of ribozymes that adopt distinct secondary structure motifs. Only one class corresponds to a previously known structural motif, that of the naturally occurring hammerhead ribozyme. Each ribozyme promotes self-cleavage via an internal phosphoester transfer reaction involving the adjacent 2'-hydroxyl group with a chemical rate enhancement of between 10(3)- and 10(6)-fold greater than the corresponding uncatalyzed rate. These findings indicate that RNA can form a multitude of secondary and tertiary structures that promote cleavage by internal phosphoester transfer. Upon further in vitro selection, a class I ribozyme that adopts an "X motif" structure dominates over all other ribozymes in the population. Thus, self-cleaving RNAs isolated by in vitro selection from random-sequence populations can rival the catalytic efficiency of natural ribozymes.

Mimicry of the hepatitis delta virus replication cycle mediated by synthetic circular oligodeoxynucleotides

BACKGROUND

Hepatitis delta virus (HDV) is a circular single-stranded RNA pathogen whose monomeric form results from self-processing. Although studies have examined minimal HDV ribozyme activities, the mechanism for forming the circular virus remains unclear, and the trans catalytic properties of self-processed forms of HDV ribozymes have not been studied. In addition, HDV ribozymes have not previously been engineered to cleave a non-HDV sequence.

RESULTS

Long repeating RNAs have been produced from in vitro rolling-circle transcription of synthetic circular oligodeoxynucleotides encoding catalytically active subsets of the entire antigenomic RNA virus. Like full-length HDV, these multimeric RNAs undergo self-processing to monomer length; importantly, cyclization is found to occur efficiently, but only in the presence of the circular template. Linear and circular monomer ribozymes and engineered variants are shown to be active in cleaving HDV and HIV RNA targets in trans, despite having self-binding domains.

CONCLUSIONS

Mimicry of the rolling-circle replication pathway for HDV replication has led to a new proposal for cyclization of HDV RNA. Under these conditions, cyclization is mediated by the complementary circular template. In addition, it has been shown that self-processed HDV ribozymes can be catalytically active in trans despite the presence of antisense sequences built into their structure.

Detailed analysis of stem I and its 5' and 3' neighbor regions in the trans-acting HDV ribozyme

To determine the stem I structure of the human hepatitis delta virus (HDV) ribozyme, which is related to the substrate sequence in the trans -acting system, we kinetically studied stem I length and sequences. Stem I extension from 7 to 8 or 9 bp caused a loss of activity and a low amount of active complex with 9 bp in the trans -acting system. In a previous report, we presented cleavage in a 6 bp stem I. The observed reaction rates indicate that the original 7 bp stem I is in the most favorable location for catalytic reaction among the possible 6-8 bp stems. To test base specificity, we replaced the original GC-rich sequence in stem I with AU-rich sequences containing six AU or UA base pairs with the natural +1G.U wobble base pair at the cleavage site. The cis -acting AU-rich molecules demonstrated similar catalytic activity to that of the wild-type. In trans -acting molecules, due to stem I instability, reaction efficiency strongly depended on the concentration of the ribozyme-substrate complex and reaction temperature. Multiple turnover was observed at 37 degreesC, strongly suggesting that stem I has no base specificity and more efficient activity can be expected under multiple turnover conditions by substituting several UA or AU base pairs into stem I. We also studied the substrate damaging sequences linked to both ends of stem I for its development in therapeutic applications and confirmed the functions of the unique structure.

Core-associated non-duplex sequences distinguishing the genomic and antigenomic self-cleaving RNAs of hepatitis delta virus

The two ribozymes found in hepatitis delta virus RNA form related but non-identical secondary structures and display similar cleavage properties in vitro. Three of the non-duplex elements hypothesized to contribute nucleotides to the catalytic core vary slightly in length between the two ribozymes and the differences are conserved in clinical isolates. Possible functional relationships of the core sequence elements were tested by systematically exchanging sequences between the two ribozymes. It was found that switching two of the elements (L3 and J4/2) from one ribozyme to the other reduced cleavage activity in both. On the other hand, exchanging the third region (J1/4) resulted in enhanced activity for one ribozyme and a smaller increase in activity for the other. Combining exchanges did not reveal any compensatory interactions involving these particular elements nor did a pattern emerge that would suggest an optimal combination of core sequences for a generalized HDV ribozyme. Non-compensatory behavior reinforces the idea that the non-duplex sequences may form sequence-specific contacts with duplex portions of the ribozyme, but, in addition, these data suggest that there may be selective pressures on the ribozyme sequences in the virus that are not reflected in the in vitro self-cleavage assays.

Analysis of the cleavage reaction of a trans-acting human hepatitis delta virus ribozyme

The cleavage reaction catalyzed by the trans -acting genomic ribozyme of human hepatitis delta virus (HDV) was analyzed with a 13mer substrate (R13) and thio-substituted [SR13(Rp) and SR13(Sp)] substrates under single-turnover conditions. The cleavage of RNA by the trans -acting HDV ribozyme proceeded as a first order reaction. The logarithm of the rate of cleavage (kclv) increased linearly (with a slope of approximately 1) between pH 4.0 and 6.0, an indication that a single deprotonation reaction occurred. This result suggests that kclv reflects the rate of the chemical cleavage step, at least around pH 5. The amount of active complex with the SR13(Sp) substrate was almost as large as with R13 (60-80%), whereas the amount of the corresponding active complex formed with the SR13(Rp) substrate was, at most, 20% of this value (with 0.5-100 mM Mg2+ions) at pH 5.0. Nonetheless, the value of kclv for all substrates was almost the same (0.4-0.5 min-1). Neither a 'thio effect' nor a 'Mn2+rescue effect' were observed. These results suggest that Mg2+ions do not interact with pro-R oxygen directly but are essential to the formation of the active complex of the ribozyme and its substrate.

Self-cleaving ribozymes of hepatitis delta virus RNA

Hepatitis delta virus (HDV) is a small single-stranded RNA satellite of hepatitis B virus. Although it is a human pathogen, it shares a number of features with a subset of the small plant satellite RNA viruses, including self-cleaving sequences in the genomic and antigenomic sequences of the viral RNA. The self-cleaving sequence is critical to viral replication and is thought to function as a ribozyme in vivo to process the products of rolling-circle replication to unit-length molecules. A divalent cation is required for cleavage and while a structural role is implicated for metal ions, a more direct role for a metal ion in catalysis has not yet been proven. A minimal natural ribozyme sequence with proficient in vitro self-cleavage activity is about 85 nucleotides long and adopts a secondary structure with four paired regions (P1-P4). The two pairings that define the 5' and 3' boundaries of the ribozyme, P1 and P2, form an atypical pseudoknot arrangement. This secondary structure places a number of constraints on the possible tertiary folding of the sequence, which together with chemical probing, photo-cross-linking, mutagenesis and computer-assisted modeling provides clues to the three-dimensional structure. The data are consistent with a model in which the cleavage site, located at the 5' end of P1, is in close proximity to three single-stranded regions, consisting of a hairpin loop at the end of P3 and two sequences joining P1 to P4 and P4 to P2. While the natural forms of the HDV ribozymes appear to be prone to misfolding, biochemical and mutagenesis studies from a number of laboratories has allowed the production of trans-acting ribozymes and smaller more active cis-acting ribozymes, both of which will aid in further mechanistic and structural studies of this RNA.

The structure of the isolated, central hairpin of the HDV antigenomic ribozyme: novel structural features and similarity of the loop in the ribozyme and free in solution

The structure of an RNA hairpin containing a seven-nucleotide loop that is present in the self-cleaving sequence of hepatitis delta virus antigenomic RNA was determined by high resolution NMR spectroscopy. The loop, which is composed of only one purine and six pyrimidines, has a suprisingly stable structure, mainly supported by sugar hydroxyl hydrogen bonds and base-base and base-phosphate stacking interactions. Compared with the structurally well-determined, seven-membered anticodon loop in tRNA, the sharp turn which affects the required 180 degrees change in direction of the sugar-phosphate backbone in the loop is shifted one nucleotide in the 3' direction. This change in direction can be characterized as a reversed U-turn. It is expected that the reversed U-turn may be found frequently in other molecules as well. There is evidence for a new non-Watson-Crick UC base pair formed between the first and the last residue in the loop, while most of the other bases in the loop are pointing outwards making them accessible to solvent. From chemical modification, mutational and photocrosslinking studies, a similar picture develops for the structure of the hairpin in the active ribozyme indicating that the loop structure in the isolated hairpin and in the ribozyme is very similar.

Detailed analysis of base preferences at the cleavage site of a trans-acting HDV ribozyme: a mutation that changes cleavage site specificity

In our previous attempt at in vitro selection of a trans - acting human hepatitis delta virus (HDV) ribozyme, we found that one of the variants, G10-68-725G, cleaved a 13 nt substrate, HDVS1, at two sites [Nishikawa,F., Kawakami,J., Chiba,A., Shirai,M., Kumar,P.K.R. and Nishikawa,S. (1996) Eur. J. Biochem., 237, 712-718]. One site was the normal cleavage site and the other site was shifted 1 nt toward the 3'-end. To clarify the interactions between nucleotides around the cleavage site of the trans -acting HDV ribozyme, we analyzed the efficiency of the reaction for every possible base pair between the substrate and the ribozyme at positions -1 (-1N:726N) and +1 (+1N:725N) relative to the cleavage site using the genomic HDV ribozyme, TdS4(Xho), and derivatives of the most active variant, G10-68. These mutagenesis analyses revealed that the +1 base of the substrate affects the structure of the catalytic core in the complex with G10-68-725G, substrate and divalent metal ions, and it shifts the cleavage site. In a comparison with other variants of the trans -acting HDV ribozyme, we found that this cleavage site shift occurred only with G10-68-725G.