Systematic substitution of individual bases in two important single-stranded regions of the HDV ribozyme for evaluation of the role of specific bases

New tools provide a second look at HDV ribozyme structure, dynamics and cleavage

The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA enzyme essential for processing viral transcripts during rolling circle viral replication. The first crystal structure of the cleaved ribozyme was solved in 1998, followed by structures of uncleaved, mutant-inhibited and ion-complexed forms. Recently, methods have been developed that make the task of modeling RNA structure and dynamics significantly easier and more reliable. We have used ERRASER and PHENIX to rebuild and re-refine the cleaved and cis-acting C75U-inhibited structures of the HDV ribozyme. The results correct local conformations and identify alternates for RNA residues, many in functionally important regions, leading to improved R values and model validation statistics for both structures. We compare the rebuilt structures to a higher resolution, trans-acting deoxy-inhibited structure of the ribozyme, and conclude that although both inhibited structures are consistent with the currently accepted hammerhead-like mechanism of cleavage, they do not add direct structural evidence to the biochemical and modeling data. However, the rebuilt structures (PDBs: 4PR6, 4PRF) provide a more robust starting point for research on the dynamics and catalytic mechanism of the HDV ribozyme and demonstrate the power of new techniques to make significant improvements in RNA structures that impact biologically relevant conclusions.

Calculation of pKas in RNA: on the structural origins and functional roles of protonated nucleotides

pK(a) calculations based on the Poisson-Boltzmann equation have been widely used to study proteins and, more recently, DNA. However, much less attention has been paid to the calculation of pK(a) shifts in RNA. There is accumulating evidence that protonated nucleotides can stabilize RNA structure and participate in enzyme catalysis within ribozymes. Here, we calculate the pK(a) shifts of nucleotides in RNA structures using numerical solutions to the Poisson-Boltzmann equation. We find that significant shifts are predicted for several nucleotides in two catalytic RNAs, the hairpin ribozyme and the hepatitis delta virus ribozyme, and that the shifts are likely to be related to their functions. We explore how different structural environments shift the pK(a)s of nucleotides from their solution values. RNA structures appear to use two basic strategies to shift pK(a)s: (a) the formation of compact structural motifs with structurally-conserved, electrostatic interactions; and (b) the arrangement of the phosphodiester backbone to focus negative electrostatic potential in specific regions.

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.

Inhibition of HCV NS3 protease by RNA aptamers in cells

Non-structural protein 3 (NS3) of hepatitis C virus (HCV) has two distinct activities, protease and helicase, which are essential for HCV proliferation. In previous work, we obtained RNA aptamers (G9-I, II and III) which specifically bound the NS3 protease domain (DeltaNS3), efficiently inhibiting protease activity in vitro. To utilize these aptamers in vivo, we constructed a G9 aptamer expression system in cultured cells, using the cytomegarovirus enhancer + chicken beta-actin globin (CAG) promoter. By conjugating the cis-acting genomic human hepatitis delta virus (HDV) ribozyme and G9-II aptamer, a chimeric HDV ribozyme-G9-II aptamer (HA) was constructed, which was used to produce stable RNA in vivo and to create tandem repeats of the functional unit. To target the transcribed RNA aptamers to the cytoplasm, the minimal mutant of constitutive transport element (CTE), derived from type D retroviruses, was conjugated at the 3' end of HA (HAC). Transcript RNAs from (HA)(n) and (HAC)(n) were processed into the G9-II aptamer unit by the cis-acting HDV ribozyme, both in vitro and in vivo. Efficient protease inhibition activity of HDV ribozyme-G9-II aptamer expression plasmid was demonstrated in HeLa cells. Protease inhibition activity level of tandem chimeric aptamers, (HA)(n) and (HAC)(n), rose with the increase of n from 1 to 4.

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.

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.

Requirement for canonical base pairing in the short pseudoknot structure of genomic hepatitis delta virus ribozyme

The tertiary structure of the 3'-cleaved product of the genomic hepatitis delta virus (HDV) ribozyme was solved by X-ray crystallographic analysis. In this structure, three single-stranded regions (SSrA, -B and -C) interact intricately with one another via hydrogen bonds between nucleotide bases, phosphate oxygens and 2'-OHs to form a nested double pseudoknot structure. Among these interactions, two Watson-Crick (W-C) base pairs, 726G-710C and 727G-709C, that form between SSrA and SSrC (P1.1) seem to be especially important for compact folding. To characterize the importance of these base pairs, ribozymes were subjected to in vitro selection from a pool of RNA molecules randomly substituted at positions 709, 710, 726 and 727. The results establish the importance of the two W-C base pairs for activity, although some mutants are active with one G-C base pair. In addition, the kinetic parameters were analyzed in all 16 combinations with two canonical base pairs. Comparison of variant ribozymes with the wild-type ribozyme reveals that the difference in reaction rates for these variants (DeltaDelta G (double dagger)) is not simply accounted for by the differences in the stability of P1.1 (DeltaDelta G (0)(37)). The role played by Mg(2+)ions in formation of the P1.1 structure is also discussed.

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.

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.

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.

Constructing an efficient trans-acting genomic HDV ribozyme

We have engineered a genomic HDV ribozyme to construct several trans-acting ribozymes for use in trans to cleave target RNAs. Among the 10 different combinations attempted, only HDV88-Trans had cleavage activity on the 13-nucleotide substrate, R13, in vitro. To improve the cleavage efficiency, at least in vitro, of the HDV88-Trans ribozyme (kclv = 0.022 min(-1)), we have constructed several variants that differ in forming stem II (length) in the pseudoknot secondary structure model. When cleavage rate constants were analyzed and compared among variants of HDV88-Trans, HDV88-Trans-4 yielded kclv = 1.7 min(-1). HDV88-Trans-4 thus represents the highest active genomic HDV ribozyme that functions in trans thus far constructed, and has activity under physiological conditions (pH 7.1 at 37 degrees C with 1 mM of MgCl2).

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

In an effort to identify the functional structure as well as new active variants of the trans-acting genomic ribozyme of human hepatitis delta virus (HDV), we applied an in vitro selection procedure. A total of 14 rounds of selection and amplification was repeated and various mutant ribozymes in G10 and G14 pools analyzed. Active ribozymes which were isolated in the present study (from G10 and G14) all possessed conserved bases (that were identified earlier) in the cis-acting molecule. A dominant clone G10-68 variant was accumulated in generation 14. Interestingly, when base substitutions were analyzed in G10-68 variant, we found that this variant appears to be close to antigenome-like HDV ribozyme molecule. Further investigations of G10-68 confirmed that each mutated base was the most appropriate nucleotide at every position of the HDV ribozyme.

3-D models of the antigenomic ribozyme of the hepatitis delta agent with eight new contacts suggested by sequence analysis of 188 cDNA clones

We mapped 359 mutations at 25 positions in synthetic variants of the antigenomic ribozyme of the hepatitis delta agent by analyzing the sequences of 188 cDNA clones. These data were used to identify three features of the ribozyme: highly conserved nucleotides, positions with restricted nucleotide substitutions and three-dimensional relationships between nucleotides. The distribution of mutations at the 25 positions was as follows: G-11 (the eleventh nucleotide from the cleavage site) was mutated in 56 clones; G-12 in 36; U-15 in 33; C-13 in 26; G-28 in 23; C-27 in 21; C-29 in 19; U-26 in 17; C-18 in 14; A-14 in 13; C-16 in 13; C-19 in 12; U-17 in 11; A-20 in 10; G-42 in 9; G-40 in 7; G-41 in 7; C-24 in 6; U-32 in 6; U-23 in 5; C-25 in 4; C-21 in 3; G-30 in 3; G-31 in 3; C-22 in 1. All clones containing a mutation at C-25 had an A at this position, suggesting that the extra cyclic amino group present in adenine and cytosine may function during the cleavage event. Mutations at certain positions were common in simple clones (containing only one or two mutations), while mutations at other positions were over-represented in more complex clones. Both compensatory base changes and co-mutational frequencies were used to identify eight pairs of nucleotides which may interact with each other: G-11 and C-18, G-12 and C-27, C-13 and G-28, C-21 and U-23/C-24, C-21 and G-30, U-23 and G-31/U-32, C24 and G-30, C-27 and G-42. These pairs, which involve some of the most conserved positions in the molecule, suggest interactions among nucleotides previously depicted in open-loop structures. The newly proposed points of contact between pairs of nucleotides are compatible with both the axehead and pseudoknot secondary structural models and were combined with previously proposed Watson-Crick base paired helices to produce two three dimensional models. In both of these, C-25 and C-76 are placed near the cleavage site.

An RNA tertiary structure of the hepatitis delta agent contains UV-sensitive bases U-712 and U-865 and can form in a bimolecular complex

Genomic RNA of the hepatitis delta agent has a highly conserved element of local tertiary structure. This element contains two nucleotides which become covalently crosslinked to each other upon irradiation with UV light. Using direct RNA analysis, we now identify the two nucleotides as U-712 and U-865 and show that the UV-induced crosslink can be broken by re-exposure to a 254 nm peak UV light source. In the rod-like secondary structural model of delta RNA, nucleotides U-712 and U-865 are off-set from each other by 5-6 bases, a distance too great to permit crosslinking. This model needs to be modified. Our data indicate that bases U-712 and U-865 closely approximate each other and suggest that the smooth helical contour proposed for delta RNA is interrupted by the UV-sensitive element. The nucleotide sequence shows that the UV-sensitive site does not have a particularly high density of conventional Watson-Crick base pairs compared to the rest of the genome. However, this element may have a number of non-Watson-Crick bonds which confer stability. Following UV-crosslinking and digestion with 1 mg/ml of RNase T1 at 37 degrees C for 45 min in 10 mM Tris-HCl, 1 mM EDTA (conditions expected to give complete digestion), this element can be isolated as part of a 54 nucleotide long partial digestion product containing at least 16 internal G residues. UV-crosslinking analysis shows that this unusual tertiary structural element can form in a bimolecular complex.

Identification of phosphate oxygens that are important for self-cleavage activity of the HDV ribozyme by phosphorothioate substitution interference analysis

A phosphorothioate substitution interference assay was used to investigate the role of the pro-Rp oxygens of phosphate groups in the self-cleavage reaction of the genomic human hepatitis delta virus (HDV) ribozyme. Incorporation of several different phosphorothioates (NTP alpha S) into the HDV ribozyme inhibited the self-cleavage activity. Incorporation of uridine 5' phosphorothioate or adenosine 5' phosphorothioate maintained 72% of the original self-cleavage activity whereas incorporation of guanosine 5' phosphorothioate or cytosine 5' phosphorothioate into the precursor reduced self-cleavage activity to about 20% in each case. Using partially substituted phosphorothioate-modified transcripts, we identified the pro-Rp oxygens that are important for the ribozyme activity, and they are located at positions 0, 1, 4, 5, 21, 24, 25, 27, 28, 30-34, 40, 43 and 75. In particular, the pro-Rp oxygens at positions 0, 1 and 21 are appear to be critical for the self-cleavage activity of the HDV ribozyme.

RNA pseudoknots

Many new RNA pseudoknot structures have been detected and proposed in the past year. Although we are still waiting for the first detailed structure of a pseudoknot, their role in processes such as translational autoregulation or ribosomal frameshifting has been extensively studied and is now well established. Pseudoknot structures appear to play a pivotal role in small subunit ribosomal RNA and in the noncoding regions of viral RNAs. There are also strong indications that RNA pseudoknots are highly suitable structural motifs for the recognition and binding of proteins.

Chemical probing studies of variants of the genomic hepatitis delta virus ribozyme by primer extension analysis

We have investigated in detail the higher order structure of the genomic hepatitis delta virus (HDV) ribozyme using various base-specific chemical probes under native, semi-denaturing, and denaturing conditions. The bases of the HDV ribozyme were probed by treatment with dimethyl sulfate [which reacts with A (at N1) and C (at N3)] and a carbodiimide [which reacts with U (at N3) and G (at N1)]. In addition, for probing G residues (at N7), RNA samples were treated with NaBH4 and aniline after modification by treatment with dimethyl sulfate. The sites of modified positions were identified by primer extension analysis with reverse transcriptase. In general, our results are consistent with the proposed pseudoknot model of secondary structure, a model that is based on data from ribonucleolytic cleavage experiments. Our results provide clues to the identification of interacting bases in the HDV ribozyme. Furthermore, using this method we identified local conformational changes in several stem variants.

Identification of important bases in a single-stranded region (SSrC) of the hepatitis delta (delta) virus ribozyme

Models for the secondary structure of genomic and antigenomic self-cleaving RNAs of human hepatitis delta (delta) virus (HDV) have been proposed by several groups. Our recent results support a pseudoknot structure and have allowed us to identify functionally important nucleotides in single-stranded regions [nucleotides 726-731 (SSrA) and nucleotides 762-766 (SSrB)]. For the identification of the important residues in the remaining single-stranded region, nucleotides 708-715 (SSrC), of the genomic HDV ribozyme, we made derivatives with a single-base substitution in the SSrC region. To screen inactive mutants rapidly, we use a simplified in-vitro selection method. Among the various base substitutions in mutants in the SSrC, U708A, C709(A/G/U) and G713C variants had less than 10% of the cleavage activity of the wild-type SSrC (HDV86). By analyzing the self-cleavage activities of various mutants, we determined the base requirements for SSrC as 5'-(U/C/G)-C-N-N-(C/A/G)-(G/A/U)-N-N-3'.