Importance of short pseudoknot base pairs between two single-stranded regions of HDV ribozyme.

Human hepatitis delta virus (HDV) ribozyme can catalyze self-cleavage reaction in the presence of Mg2+ ions, yielding products with 2',3'-cyclic phosphate and 5'-OH termini as do hammerhead and hairpin ribozymes. Recently, the tertiary structure of 3'-cleaved product of genomic HDV ribozyme was solved by X-ray crystallographic analysis. In this structure three single-stranded regions (SSrA, -B and -C) interacts intricately with hydrogen bonds between bases, phosphate oxygens and 2'-OHs to form nested double pseudoknot structure. Especially two Watson-Crick base pairs, 726G-710C and 727G-709C, between SSrA and SSrC, seems to be important for compact folding. To characterize the necessity of the two base pairs, we performed in vitro selection of active ribozymes using random RNA pool which mutated at 709, 710, 726 and 727. The result indicates that basically one G-C base pair is necessary for the activity.     

Formation of the P1.1 pseudoknot is critical for both the cleavage activity and substrate specificity of an antigenomic trans-acting hepatitis delta ribozyme

Hepatitis delta virus RNAs possess self-cleavage activities that produce 2′,3′-cyclic phosphate and 5′-hydroxyl termini (i.e. cis-acting delta ribozyme). Trans-acting delta ribozymes have been engineered by removing a junction from the cis version, thereby producing one molecule possessing the substrate sequence and the other the catalytic domain. According to the pseudoknot model, the secondary structure of the delta ribozyme includes a pseudoknot (i.e. P1.1 stem) formed by two base pairs from residues of the L3 loop and J1/4 junction. A collection of 48 P1.1 stem mutants was synthesized in order to provide an original characterization of both the importance and the structure of this pseudoknot in a trans-acting version of the ribozyme. Several structural differences were noted compared to the results reported for cis-acting ribozymes. For example, a combination of two stable Watson–Crick base pairs composing the essential P1.1 stem was demonstrated to be crucial for a significant level of activity, while the cis version required only one base pair. In addition, we present the first physical evidences revealing that the composition of the P1.1 stem affects the substrate specificity for ribozyme cleavage. Depending on the residues forming the J1/4 junction, non-productive ribozyme–substrate complexes can be observed. This phenomenon is proposed to be important for further development of a gene-inactivation system based on delta ribozyme.     

Noncanonical G(syn)-G(anti) base pairs stabilized by sulphate anions in two X-ray structures of the (GUGGUCUGAUGAGGCC) RNA duplex

The structures of two crystal forms of the RNA 16-mer with the sequence GUGGUCUGAUGAGGCC, grown in the presence of a high concentration of sulphate ions, have been determined using synchrotron radiation at 1.4- and 2.0-Å resolution. RNA with this sequence is known as one of the two strands of the noncleavable form of the hammerhead ribozyme. In both crystal structures, two G(syn)–G(anti) noncanonical base pairs are observed in the middle of a 14 base-pair (bp) duplex having 5′-dangling GU residues. Both structures contain sulphate anions interacting with the G–G bp stabilizing G in its syn conformation and bridging the two RNA strands. In both cases the interactions take place in the major groove, although the anions are accommodated within different helix geometries, most pronounced in the changing width of the major groove. In one structure, where a single sulphate spans both G–G pairs, the major groove is closed around the anion, while in the other structure, where each of the two G–G pairs is associated with a separate sulphate, the groove is open. This work provides the first examples of a G–G pair in syn-anti conformation, which minimizes the purine–purine clash in the center of the duplex, while utilizing its residual hydrogen bonding potential in specific interactions with sulphate anions.     

The structure of metallo-DNA with consecutive thymine–HgII–thymine base pairs explains positive entropy for the metallo base pair formation

We have determined the three-dimensional (3D) structure of DNA duplex that includes tandem Hg^II-mediated T–T base pairs (thymine–Hg^II–thymine, T–Hg^II–T) with NMR spectroscopy in solution. This is the first 3D structure of metallo-DNA (covalently metallated DNA) composed exclusively of ‘NATURAL’ bases. The T–Hg^II–T base pairs whose chemical structure was determined with the ¹⁵N NMR spectroscopy were well accommodated in a B-form double helix, mimicking normal Watson–Crick base pairs. The Hg atoms aligned along DNA helical axis were shielded from the bulk water. The complete dehydration of Hg atoms inside DNA explained the positive reaction entropy (ΔS) for the T–Hg^II–T base pair formation. The positive ΔS value arises owing to the Hg^II dehydration, which was approved with the 3D structure. The 3D structure explained extraordinary affinity of thymine towards Hg^II and revealed arrangement of T–Hg^II–T base pairs in metallo-DNA.     

Structure–function analysis of the ribozymes of chrysanthemum chlorotic mottle viroid: a loop–loop interaction motif conserved in most natural hammerheads

Loop–loop tertiary interactions play a key role in the folding and catalytic activity of natural hammerhead ribozymes. Using a combination of NMR spectroscopy, site-directed mutagenesis and kinetic and infectivity analyses, we have examined the structure and function of loops 1 and 2 of the (+) and (–) hammerheads of chrysanthemum chlorotic mottle viroid RNA. In both hammerheads, loop 1 is a heptanucleotide hairpin loop containing an exposed U at its 5′ side and an extrahelical U at its 3′-side critical for the catalytic activity of the ribozyme in vitro and for viroid infectivity in vivo, whereas loop 2 has a key opened A at its 3′-side. These structural features promote a specific loop–loop interaction motif across the major groove. The essential features of this tertiary structure element, base pairing between the 5′ U of loop 1 and the 3′ A of loop 2, and interaction of the extrahelical pyrimidine of loop 1 with loop 2, are likely shared by a significant fraction of natural hammerheads.     

Hydration of short DNA, RNA and 2'-OMe oligonucleotides determined by osmotic stressing

Studies on hydration are important for better understanding of structure and function of nucleic acids. We compared the hydration of self-complementary DNA, RNA and 2′-O-methyl (2′-OMe) oligonucleotides GCGAAUUCGC, (UA)₆ and (CG)₃ using the osmotic stressing method. The number of water molecules released upon melting of oligonucleotide duplexes, Δn_W, was calculated from the dependence of melting temperature on water activity and the enthalpy, both measured with UV thermal melting experiments. The water activity was changed by addition of ethylene glycol, glycerol and acetamide as small organic co-solutes. The Δn_W was 3–4 for RNA duplexes and 2–3 for DNA and 2′-OMe duplexes. Thus, the RNA duplexes were hydrated more than the DNA and the 2′-OMe oligonucleotide duplexes by approximately one to two water molecules depending on the sequence. Consistent with previous studies, GC base pairs were hydrated more than AU pairs in RNA, whereas in DNA and 2′-OMe oligonucleotides the difference in hydration between these two base pairs was relatively small. Our data suggest that the better hydration of RNA contributes to the increased enthalpic stability of RNA duplexes compared with DNA duplexes.     

A ribozyme that triphosphorylates RNA 5′-hydroxyl groups

The RNA world hypothesis describes a stage in the early evolution of life in which RNA served as genome and as the only genome-encoded catalyst. To test whether RNA world organisms could have used cyclic trimetaphosphate as an energy source, we developed an in vitro selection strategy for isolating ribozymes that catalyze the triphosphorylation of RNA 5′-hydroxyl groups with trimetaphosphate. Several active sequences were isolated, and one ribozyme was analyzed in more detail. The ribozyme was truncated to 96 nt, while retaining full activity. It was converted to a trans-format and reacted with rates of 0.16 min⁻¹ under optimal conditions. The secondary structure appears to contain a four-helical junction motif. This study showed that ribozymes can use trimetaphosphate to triphosphorylate RNA 5′-hydroxyl groups and suggested that RNA world organisms could have used trimetaphosphate as their energy source.     

Structure-function relationships of two closely related group IC3 intron ribozymes from Azoarcus and Synechococcus pre-tRNA

The two group IC3 pre-tRNA introns from Azoarcus and Synechococcus share very analogous secondary structures. They are small group I ribozymes that possess only two peripheral domains, P2 and P9. However, the 3′-splice site hydrolysis activity of the Synechococcus ribozyme critically depends on P2 whereas that of Azoarcus does not, indicating that the structure–function relationships of the two ribozymes are strikingly different despite their structural resemblance. To identify the element(s) that determines the catalytic properties of these ribozymes, we undertook analyses of chimeric ribozymes prepared by swapping their structural elements. We found that the difference can be attributed to a small number of nucleotides within the conserved core region. Further analysis by employing in vitro selection revealed that a base triple interaction (P4bp3 × J6/7-2) is a critical element for determining activity and suggests the existence of a novel base quintuple involving the base triple P4bp5 × J8/7-5.     

Novel RNA structural features of an alternatively splicing group II intron from Clostridium tetani

Group II introns are ribozymes in bacterial and organellar genomes that function as self-splicing introns and as retroelements. Previously, we reported that the group II intron C.te.I1 of Clostridium tetani alternatively splices in vivo to produce five distinct coding mRNAs. Accurate fusion of upstream and downstream reading frames requires a shifted 5′ splice site located 8 nt upstream of the usual 5′ GUGYG motif. This site is specified by the ribozyme through an altered intron/exon-binding site 1 (IBS1–EBS1) pairing. Here we use mutagenesis and self-splicing assays to investigate in more detail the significance of the structural features of the C.te.I1 ribozyme. The shifted 5′ splice site is shown to be affected by structures in addition to IBS1–EBS1, and unlike other group II introns, C.te.I1 appears to require a spacer between IBS1 and the GUGYG motif. In addition, the mechanism of 3′ exon recognition is modified from the ancestral IIB mechanism to a IIA-like mechanism that appears to be longer than the typical single base-pair interaction and may extend up to 4 bp. The novel ribozyme properties that have evolved for C.te.I1 illustrate the plasticity of group II introns in adapting new structural and catalytic properties that can be utilized to affect gene expression.     

Relationship between Intracellular Phosphate, Proton Motive Force, and Rate of Nongrowth Energy Dissipation (Energy Spilling) in Streptococcus bovis JB1

When the rate of glucose addition to nongrowing Streptococcus bovis cell suspensions was increased, the fermentation was homolactic, fructose-1,6-diphosphate (FDP) increased, intracellular inorganic phosphate (P_i) declined, and the energy-spilling rate increased. ATP and ADP were not significantly affected by glucose consumption rate, but the decrease in P_i was sufficient to cause an increase in the free energy of ATP hydrolysis (ΔG′p). The increase in ΔG′p was correlated with an increase in proton motive force (Δp). S. bovis continuous cultures (dilution rate of 0.65 h⁻¹) that were provided with ammonia as the sole nitrogen source also had high rates of lactate production and energy spilling. When Trypticase was added as a source of amino acids, lactate production decreased; a greater fraction of the glucose was converted to acetate, formate, and ethanol; and the energy-spilling rate decreased. Trypticase also caused a decrease in FDP, an increase in P_i, and a decrease in Δp. The change in Δp could be explained by P_i-dependent changes in the ΔG′p. When P_i declined, ΔG′p and Δp increased. The ratio of ΔG′p to Δp (millivolt per millivolt) was always high (>4) at low rates of energy spilling but declined when the energy-spilling rate increased. Based on these results, it appears that Δp and the energy-spilling rate are responsive to fluctuations in the intracellular P_i concentration.     

New insights into Hoogsteen base pairs in DNA duplexes from a structure-based survey

Hoogsteen (HG) base pairs (bps) provide an alternative pairing geometry to Watson–Crick (WC) bps and can play unique functional roles in duplex DNA. Here, we use structural features unique to HG bps (syn purine base, HG hydrogen bonds and constricted C1′–C1′ distance across the bp) to search for HG bps in X-ray structures of DNA duplexes in the Protein Data Bank. The survey identifies 106 A•T and 34 G•C HG bps in DNA duplexes, many of which are undocumented in the literature. It also uncovers HG-like bps with syn purines lacking HG hydrogen bonds or constricted C1′–C1′ distances that are analogous to conformations that have been proposed to populate the WC-to-HG transition pathway. The survey reveals HG preferences similar to those observed for transient HG bps in solution by nuclear magnetic resonance, including stronger preferences for A•T versus G•C bps, TA versus GG steps, and also suggests enrichment at terminal ends with a preference for 5′-purine. HG bps induce small local perturbations in neighboring bps and, surprisingly, a small but significant degree of DNA bending (∼14°) directed toward the major groove. The survey provides insights into the preferences and structural consequences of HG bps in duplex DNA.     

Random mutations to evaluate the role of bases at two important single-stranded regions of genomic HDV ribozyme.

In elucidating function of two important single-stranded regions [SSrA (726-731 nt) and SSrB (762-766 nt)] derived mainly from three secondary structure models in genomic hepatitis delta virus (HDV) ribozyme possessing self-cleavage activity, we have constructed several random mutants at those two regions on the HDV88 molecule (683-770 nt) by oligonucleotide-directed mutagenesis. When self-cleavage activities were compared among mutants, at the region SSrA, G726 was found to play an important role during cleavage reaction since substitutions of the base to A (mutant A20) or C (mutant A16) or U (mutant A23), reduced the ribozyme activity to very low levels suggesting the importance of G726 position. C763 at SSrB region was found to play a more significant role during catalysis than G726 (at region SSrA) since any substitutions at C763 completely inactivated the ribozyme. Other bases located in these two regions could be substituted to other bases at the expense of some self-cleavage activity. The results presented here together with our previous deletion analysis indicate that these two regions may play an important role during cleavage process.     

Effect of G-1 on histidine tRNA microhelix conformation

Histidine tRNAs (tRNA^His) are unique in that they possess an extra 5′-base (G-1) not found in other tRNAs. Deletion of G-1 results in at least a 250-fold reduction in the rate of histidine charging in vitro. To better understand the role of the G-1 nucleotide in defining the structure of tRNA^His, and to correlate structure with cognate amino acid charging, NMR and molecular dynamics (MD) studies were performed on the wild-type and a ΔG-1 mutant Escherichia coli histidine tRNA acceptor stem microhelix. Using NMR-derived distance restraints, global structural characteristics are described and interpreted to rationalize experimental observations with respect to aminoacylation activity. The quality of the NMR-derived solution conformations of the wild-type and ΔG-1 histidine microhelices (micro helix^His) is assessed using a variety of MD-based computational protocols. Most of the duplex regions of the acceptor stem and the UUCG tetraloop are well defined and effectively superimposable for the wild-type and ΔG-1 mutant microhelix^His. Differences, however, are observed at the end of the helix and in the single-stranded CCCA-3′ tail. The wild-type microhelix^His structure is more well defined than the mutant and folds into a ‘stacked fold-back’ conformation. In contrast, we observe fraying of the first two base pairs and looping back of the single-stranded region in the ΔG-1 mutant resulting in a much less well defined conformation. Thus the role of the extra G-1 base of the unique G-1:C73 base pair in tRNA^His may be to prevent end-fraying and stabilize the stacked fold-back conformation of the CCCA-3′ region.     

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.     

Pronounced instability of tandem IU base pairs in RNA

Optical melting was used to determine the stabilities of three series of RNA oligomers containing tandem XU base pairs, GGCXUGCC (5′XU3′), GGCUXGCC (5′UX3′) and GGCXXGGC/CCGUUCCG (5′XX3′), where X is either A, G or I (inosine). The helices containing tandem AU base pairs were the most stable in the first two series (5′XU3′ and 5′UX3′), with an average melting temperature ~11°C higher than the helices with tandem 5′GU3′ base pairs and 25°C higher than the helices with tandem 5′IU3′ base pairs. For the third series (5′XX3′), the helix containing tandem GG is the most stable, with an average melting temperature ~2°C higher than the helix with tandem AA base pairs and ~24°C higher than the helix with tandem II base pairs. The thermodynamic stability of the oligomers with tandem IU base pairs was also investigated as a function of magnesium ion concentration. As with normal A–U or G–U tandem duplexes, the data could best be interpreted as non-specific binding of magnesium ions to the inosine-containing RNA oligonucleotides.     

Biochemical analysis of pistol self-cleaving ribozymes

Pistol RNAs are members of a distinct class of self-cleaving ribozymes that was recently discovered by using a bioinformatics search strategy. Several hundred pistol ribozymes share a consensus sequence including 10 highly conserved nucleotides and many other modestly conserved nucleotides associated with specific secondary structure features, including three base-paired stems and a pseudoknot. A representative pistol ribozyme from the bacterium Lysinibacillus sphaericus was found to promote RNA strand scission with a rate constant of ∼10 min⁻¹ under physiological Mg²⁺ and pH conditions. The reaction proceeds via the nucleophilic attack of a 2′-oxygen atom on the adjacent phosphorus center, and thus adheres to the same general catalytic mechanism of internal phosphoester transfer as found with all other classes of natural self-cleaving ribozymes discovered to date. Analyses of the kinetic characteristics and the metal ion requirements of the cleavage reaction reveal that members of this ribozyme class likely use several catalytic strategies to promote the rapid cleavage of RNA.     

Predicting loop–helix tertiary structural contacts in RNA pseudoknots

Tertiary interactions between loops and helical stems play critical roles in the biological function of many RNA pseudoknots. However, quantitative predictions for RNA tertiary interactions remain elusive. Here we report a statistical mechanical model for the prediction of noncanonical loop–stem base-pairing interactions in RNA pseudoknots. Central to the model is the evaluation of the conformational entropy for the pseudoknotted folds with defined loop–stem tertiary structural contacts. We develop an RNA virtual bond-based conformational model (Vfold model), which permits a rigorous computation of the conformational entropy for a given fold that contains loop–stem tertiary contacts. With the entropy parameters predicted from the Vfold model and the energy parameters for the tertiary contacts as inserted parameters, we can then predict the RNA folding thermodynamics, from which we can extract the tertiary contact thermodynamic parameters from theory–experimental comparisons. These comparisons reveal a contact enthalpy (ΔH) of −14 kcal/mol and a contact entropy (ΔS) of −38 cal/mol/K for a protonated C⁺•(G–C) base triple at pH 7.0, and (ΔH = −7 kcal/mol, ΔS = −19 cal/mol/K) for an unprotonated base triple. Tests of the model for a series of pseudoknots show good theory–experiment agreement. Based on the extracted energy parameters for the tertiary structural contacts, the model enables predictions for the structure, stability, and folding pathways for RNA pseudoknots with known or postulated loop–stem tertiary contacts from the nucleotide sequence alone.     

Multiple-turnover thio-ATP hydrolase and phospho-enzyme intermediate formation activities catalyzed by an RNA enzyme

Ribozymes that phosphorylate internal 2′-OH positions mimic the first mechanistic step of P-type ATPase enzymes by forming a phospho-enzyme intermediate. We previously described 2′-autophosphorylation and autothiophosphorylation by the 2PTmin3.2 ribozyme. In the present work we demonstrate that the thiophosphorylated form of this ribozyme can de-thiophosphorylate in the absence of ATPγS. Identical ionic conditions yield a thiophosphorylated strand when ATPγS is included, thus effecting a net ATPγS hydrolysis. The de-thiophosphorylation step is nearly independent of pH over the range of 6.3–8.5 and does not require a specifically folded RNA structure, but this step is greatly stimulated by transition metal ions. By monitoring thiophosphate release, we observe 29–46 ATPγS hydrolyzed per ribozyme strand in 24 h, corresponding to a turnover rate of 1.2–2.0 h⁻¹. The existence of an ATP- (or thio-ATP-)powered catalytic cycle raises the possibility of using ribozymes to transduce chemical energy into mechanical work for nucleic acid nanodevices.     

Imino proton NMR analysis of HDV ribozymes: nested double pseudoknot structure and Mg2+ ion-binding site close to the catalytic core in solution

Minimized trans-acting HDV ribozyme systems consisting of three (Rz-3) and two (Rz-2) RNA strands were prepared and their folding conformations were analyzed by NMR spectroscopy. The guanosine residues in one of the enzyme components of Rz-3 were labeled with ¹³C and ¹⁵N. Imino proton signals were assigned by analysis of NOESY and HSQC spectra. The results are consistent with the nested double pseudoknot model, which contains novel base pairs (P1.1), as observed in the crystal structure of a genomic HDV ribozyme. The NOE connectivities suggest an additional G:G pair at the bottom of P1.1 and at the top of P4. The effects of temperature and Mg²⁺ ions on base pairs for Rz-3 were examined. The temperature variation experiment on Rz-3 showed that P3 is the most stable and that P1.1 is as stable as P1 and P2. The imino proton signals of the G:U pair at the bottom of P1 and the top of P1.1, which are close to the cleavage site, showed the largest changes upon Mg²⁺ titration of Rz-3. The results suggest that the catalytic Mg²⁺ ion binds to the pocket formed by P1 and L3.