Existence of efficient divalent metal ion-catalyzed and inefficient divalent metal ion-independent channels in reactions catalyzed by a hammerhead ribozyme

The hammerhead ribozyme is generally accepted as a well characterized metalloenzyme. However, the precise nature of the interactions of the RNA with metal ions remains to be fully defined. Examination of metal ion-catalyzed hammerhead reactions at limited concentrations of metal ions is useful for evaluation of the role of metal ions, as demonstrated in this study. At concentrations of Mn²⁺ ions from 0.3 to 3 mM, addition of the ribozyme to the reaction mixture under single-turnover conditions enhances the reaction with the product reaching a fixed maximum level. Further addition of the ribozyme inhibits the reaction, demonstrating that a certain number of divalent metal ions is required for proper folding and also for catalysis. At extremely high concentrations, monovalent ions, such as Na⁺ ions, can also serve as cofactors in hammerhead ribozyme-catalyzed reactions. However, the catalytic efficiency of monovalent ions is extremely low and, thus, high concentrations are required. Furthermore, addition of monovalent ions to divalent metal ion-catalyzed hammerhead reactions inhibits the divalent metal ion-catalyzed reactions, suggesting that the more desirable divalent metal ion–ribozyme complexes are converted to less desirable monovalent metal ion–ribozyme complexes via removal of divalent metal ions, which serve as a structural support in the ribozyme complex. Even though two channels appear to exist, namely an efficient divalent metal ion-catalyzed channel and an inefficient monovalent metal ion-catalyzed channel, it is clear that, under physiological conditions, hammerhead ribozymes are metalloenzymes that act via the significantly more efficient divalent metal ion-dependent channel. Moreover, the observed kinetic data are consistent with Lilley’s and DeRose’s two-phase folding model that was based on ground state structure analyses.     

Photocrosslinking detects a compact, active structure of the hammerhead ribozyme

The hammerhead ribozyme has been intensively studied for approximately 15 years, but its cleavage mechanism is not yet understood. Crystal structures reveal a Y-shaped molecule in which the cleavage site is not ideally aligned for an S(N)2 reaction and no RNA functional groups are positioned appropriately to perform the roles of acid and base or other functions in the catalysis. If the ribozyme folds to a more compact structure in the transition state, it probably does so only transiently. We have used photocrosslinking as a tool to trap hammerhead ribozyme-substrate complexes in various stages of folding. Results suggest that the two substrate residues flanking the cleavage site approach and stack upon two guanosines (G8 and G12) in domain 2, moving 10-15 A closer to domain 2 than they appear in the crystal structure. Most crosslinks obtained with the nucleotide analogues positioned in the ribozyme core are catalytically inactive; however, one cobalt(III) hexaammine-dependent crosslink of an unmodified ribozyme retains catalytic activity and confirms the close stacking of cleavage site residue C17 with nucleotide G8 in domain 2. These findings suggest that residues involved in the chemistry of hammerhead catalysis are likely located in that region containing G8 and G12.     

Disparate HDV ribozyme crystal structures represent intermediates on a rugged free-energy landscape

The hepatitis delta virus (HDV) ribozyme is a member of the class of small, self-cleaving catalytic RNAs found in a wide range of genomes from HDV to human. Both pre- and post-catalysis (precursor and product) crystal structures of the cis-acting genomic HDV ribozyme have been determined. These structures, together with extensive solution probing, have suggested that a significant conformational change accompanies catalysis. A recent crystal structure of a trans-acting precursor, obtained at low pH and by molecular replacement from the previous product conformation, conforms to the product, raising the possibility that it represents an activated conformer past the conformational change. Here, using fluorescence resonance energy transfer (FRET), we discovered that cleavage of this ribozyme at physiological pH is accompanied by a structural lengthening in magnitude comparable to previous trans-acting HDV ribozymes. Conformational heterogeneity observed by FRET in solution appears to have been removed upon crystallization. Analysis of a total of 1.8 µsec of molecular dynamics (MD) simulations showed that the crystallographically unresolved cleavage site conformation is likely correctly modeled after the hammerhead ribozyme, but that crystal contacts and the removal of several 2′-oxygens near the scissile phosphate compromise catalytic in-line fitness. A cis-acting version of the ribozyme exhibits a more dynamic active site, while a G-1 residue upstream of the scissile phosphate favors poor fitness, allowing us to rationalize corresponding changes in catalytic activity. Based on these data, we propose that the available crystal structures of the HDV ribozyme represent intermediates on an overall rugged RNA folding free-energy landscape.     

Thymine DNA glycosylase exhibits negligible affinity for nucleobases that it removes from DNA

Thymine DNA Glycosylase (TDG) performs essential functions in maintaining genetic integrity and epigenetic regulation. Initiating base excision repair, TDG removes thymine from mutagenic G·T mispairs caused by 5-methylcytosine (mC) deamination and other lesions including uracil (U) and 5-hydroxymethyluracil (hmU). In DNA demethylation, TDG excises 5-formylcytosine (fC) and 5-carboxylcytosine (caC), which are generated from mC by Tet (ten–eleven translocation) enzymes. Using improved crystallization conditions, we solved high-resolution (up to 1.45 Å) structures of TDG enzyme–product complexes generated from substrates including G·U, G·T, G·hmU, G·fC and G·caC. The structures reveal many new features, including key water-mediated enzyme–substrate interactions. Together with nuclear magnetic resonance experiments, the structures demonstrate that TDG releases the excised base from its tight product complex with abasic DNA, contrary to previous reports. Moreover, DNA-free TDG exhibits no significant binding to free nucleobases (U, T, hmU), indicating a K_d >> 10 mM. The structures reveal a solvent-filled channel to the active site, which might facilitate dissociation of the excised base and enable caC excision, which involves solvent-mediated acid catalysis. Dissociation of the excised base allows TDG to bind the beta rather than the alpha anomer of the abasic sugar, which might stabilize the enzyme–product complex.     

Structure—function studies of the hammerhead ribozyme

Elucidation of the catalytic mechanism and structure—function relationship studies of the hammerhead ribozyme continue to be an area of intensive research. A combination of diverse approaches, such as X ray crystallography, spectral studies, chemical modifications, sequence variations and kinetic analyses, have provided valuable insight into the cleavage mechanism of this ribozyme. The hammerhead ribozyme crystal structures have provided valuable insight into conformational deformations needed to attain the catalytically active structure. Similarly, determination of ribozyme solution structure by spectroscopic analyses and the effect of divalent metal ions on RNA folding has further aided in the construction of a model for hammerhmead catalysis.     

Three conserved guanosines approach the reaction site in native and minimal hammerhead ribozymes

Native hammerhead ribozymes contain RNA domains that enable high catalytic activity under physiological conditions, where minimal hammerheads show little activity. However, little is known about potential differences in native versus minimal ribozyme folding. Here, we present results of photocross-linking analysis of native and minimal hammerheads containing photoreactive nucleobases 6-thioguanosine, 2,6-diaminopurine, 4-thiouridine, and pyrrolocytidine, introduced at specific sites within the catalytic core. Under conditions where catalytic activity is observed, the two substrate nucleobases spanning the cleavage site approach and stack upon G8 and G12 of the native hammerhead, two conserved nucleobases that show similar behavior in minimal constructs, have been implicated in general acid-base catalysis, and are >15 A from the cleavage site in the crystal structures. Pyrrolocytidine at cleavage site position 17 forms an efficient crosslink to G12, and the crosslinked RNA retains catalytic activity. Multiple cross-linked species point to a structural rearrangement within the U-turn, positioning residue G5 in the vicinity of cleavage site position 1.1. Intriguing crosslinks were triggered by nucleotide analogues at positions distal to the crosslinked residues; for example, 6-thioguanosine at position 5 induced a crosslink between G12 and C17, suggesting an intimate functional communication among these three nucleobases. Together, these results support a model in which the native hammerhead folds to an active structure similar to that of the minimal ribozyme, and significantly different from the crystallographic structures.     

Nucleic acid binding properties of the nucleic acid chaperone domain of hepatitis delta antigen

The N terminal region of hepatitis delta antigen (HDAg), referred to here as NdAg, has a nucleic acid chaperone activity that modulates the ribozyme activity of hepatitis delta virus (HDV) RNA and stimulates hammerhead ribozyme catalysis. We characterized the nucleic acid binding properties of NdAg, identified the structural and sequence domains important for nucleic acid binding, and studied the correlation between the nucleic acid binding ability and the nucleic acid chaperone activity. NdAg does not recognize the catalytic core of HDV ribozyme specifically. Instead, NdAg interacts with a variety of nucleic acids and has higher affinities to longer nucleic acids. The studies with RNA homopolymers reveal that the binding site size of NdAg is around nine nucleotides long. The extreme N terminal portion of NdAg, the following coiled-coil domain and the basic amino acid clusters in these regions are important for nucleic acid binding. The nucleic acid–NdAg complex is stabilized largely by electrostatic interactions. The formation of RNA–protein complex appears to be a prerequisite for facilitating hammerhead ribozyme catalysis of NdAg and its derivatives. Mutations that reduce the RNA binding activity or high ionic strength that destabilizes the RNA–protein complex, reduce the nucleic acid chaperone activity of NdAg.     

Sequence requirements of the hammerhead RNA self-cleavage reaction.

A previously well-characterized hammerhead catalytic RNA consisting of a 24-nucleotide substrate and a 19-nucleotide ribozyme was used to perform an extensive mutagenesis study. The cleavage rates of 21 different substrate mutations and 24 different ribozyme mutations were determined. Only one of the three phylogenetically conserved base pairs but all nine of the conserved single-stranded residues in the central core are needed for self cleavage. In most cases the mutations did not alter the ability of the hammerhead to assemble into a bimolecular complex. In the few cases where mutant hammerheads did not assemble, it appeared to be the result of the mutation stabilizing an alternate substrate or ribozyme secondary structure. All combinations of mutant substrate and mutant ribozyme were less active than the corresponding single mutations, suggesting that the hammerhead contains few, if any, replaceable tertiary interactions as are found in tRNA. The refined consensus hammerhead resulting from this work was used to identify potential hammerheads present in a variety of Escherichia coli gene sequences.     

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.     

An ultraviolet crosslink in the hammerhead ribozyme dependent on 2-thiocytidine or 4-thiouridine substitution.

The hammerhead domain is one of the smallest known ribozymes. Like other ribozymes it catalyzes site-specific cleavage of a phosphodiester bond. The hammerhead ribozyme has been the subject of a vast number of biochemical and structural studies aimed at determining the structure and mechanism of cleavage. Recently crystallographic analysis has produced a structure for the hammerhead. As the hammerhead is capable of undergoing cleavage within the crystal, it would appear that the crystal structure is representative of the catalytically active solution structure. However, the crystal structure conflicts with much of the biochemical data and reveals a catalytic metal ion binding site expected to be of very low affinity. Clearly, additional studies are needed to reconcile the discrepancies and provide a clear understanding of the structure and mechanism of the hammerhead ribozyme. Here we demonstrate that a unique crosslink can be induced in the hammerhead with 2-thiocytidine or 4-thiouridine substitution at different locations within the conserved core. Generation of the same crosslink with different modifications at different positions suggests that the structure trapped by the crosslink may be relevant to the catalytically active solution structure of the hammerhead ribozyme. As this crosslink appears to be incompatible with the crystal structure, this provides yet another indication that the active solution and crystal structures may differ significantly.     

Molecular dynamics investigations of hammerhead ribozyme RNA

The hammerhead ribozyme, a small catalytic RNA molecule, cleaves, in the presence of magnesium ions, a specific phosphodiester bond within its own backbone, leading to 23-cyclic phosphate and 5-OH extremities. In order to study the dynamical flexibility of the hammerhead RNA, we performed molecular dynamics simulations of the solvated crystal structure of an active hammerhead ribozyme, obtained after flash-freezing crystals soaked with magnesium. Because of a careful equilibration protocol and the use of the Ewald summation in calculating the electrostatic interactions, the RNA structure remained close to the crystal structure, as attested by a root-mean-square deviation below 2.5 A after 750 ps of simulation. All Watson-Crick base pairs were intact at the end of the simulations. The tertiary interactions, such as the sheared G.A pairs and the U-turn, important for the stabilisation of the three-dimensional RNA fold, were also retained. The results demonstrate that molecular dynamics simulations can be successfully used to investigate the dynamical behaviour of a ribozyme, thus, opening a road to study the role of transient structural changes involved in ribozyme catalysis.     

Probing the hammerhead ribozyme structure with ribonucleases.

Susceptibility to RNase digestion has been used to probe the conformation of the hammerhead ribozyme structure prepared from chemically synthesised RNAs. Less than about 1.5% of the total sample was digested to obtain a profile of RNase digestion sites. The observed digestion profiles confirmed the predicted base-paired secondary structure for the hammerhead. Digestion profiles of both cis and trans hammerhead structures were nearly identical which indicated that the structural interactions leading to self-cleavage were similar for both systems. Furthermore, the presence or absence of Mg2+ did not affect the RNase digestion profiles, thus indicating that Mg2+ did not modify the hammerhead structure significantly to induce self-cleavage. The base-paired stems I and II in the hammerhead structure were stable whereas stem III, which was susceptible to digestion, appeared to be an unstable region. The single strand domains separating the stems were susceptible to digestion with the exception of sites adjacent to guanosines; GL2.1 in the stem II loop and G12 in the conserved GAAAC sequence, which separates stems II and III. The absence of digestion at GL2.1 in the stem II hairpin loop of the hammerhead complex was maintained in uncomplexed ribozyme and in short oligonucleotides containing only the stem II hairpin region. In contrast, the G12 site became susceptible when the ribozyme was not complexed with its substrate. Overall the results are consistent with the role of Mg2+ in the hammerhead self-cleavage reaction being catalytic and not structural.     

Elucidation of Enzymatic Mechanism of Phenazine Biosynthetic Protein PhzF Using QM/MM and MD Simulations

The phenazine biosynthetic pathway is of considerable importance for the pharmaceutical industry. The pathway produces two products: phenazine-1,6-dicarboxylic acid and phenazine-1-carboxylic acid. PhzF is an isomerase that catalyzes trans-2,3-dihydro-3-hydroxyanthranilic acid isomerization and plays an essential role in the phenazine biosynthetic pathway. Although the PhzF crystal structure has been determined recently, an understanding of the detailed catalytic mechanism and the roles of key catalytic residues are still lacking. In this study, a computational strategy using a combination of molecular modeling, molecular dynamics simulations, and quantum mechanics/molecular mechanics simulations was used to elucidate these important issues. The Apo enzyme, enzyme–substrate complexes with negatively charged Glu45, enzyme–transition state analog inhibitor complexes with neutral Glu45, and enzyme–product complexes with negatively charged Glu45 structures were optimized and modeled using a 200 ns molecular dynamics simulation. Residues such as Gly73, His74, Asp208, Gly212, Ser213, and water, which play important roles in ligand binding and the isomerization reaction, were comprehensively investigated. Our results suggest that the Glu45 residue at the active site of PhzF acts as a general base/acid catalyst during proton transfer. This study provides new insights into the detailed catalytic mechanism of PhzF and the results have important implications for PhzF modification.     

Evidence for a hydroxide ion bridging two magnesium ions at the active site of the hammerhead ribozyme.

In the presence of magnesium ions, cleavage by the hammerhead ribozyme RNA at a specific residue leads to 2'3'-cyclic phosphate and 5'-OH extremities. In the cleavage reaction an activated ribose 2'-hydroxyl group attacks its attached 3'-phosphate. Molecular dynamics simulations of the crystal structure of the hammerhead ribozyme, obtained after flash-freezing of crystals under conditions where the ribozyme is active, provide evidence that a mu-bridging OH-ion is located between two Mg2+ions close to the cleavable phosphate. Constrained simulations show further that a flip from the C3'- endo to the C2'- endo conformation of the ribose at the cleavable phosphate brings the 2'-hydroxyl in proximity to both the attacked phosphorous atom and the mu-bridging OH-ion. Thus, the simulations lead to a detailed new insight into the mechanism of hammerhead ribozyme cleavage where a mu-hydroxo bridged magnesium cluster, located on the deep groove side, provides an OH-ion that is able to activate the 2'-hydroxyl nucleophile after a minor and localized conformational change in the RNA.     

High hydrostatic pressure approach proves RNA catalytic activity without magnesium

High hydrostatic pressure (HHP) technique was used to evaluate a mechanism of RNA hydrolysis with RNA. We showed that hammerhead ribozyme specifically cleaves RNA substrate at HHP in the absence of Mg(2+). A deoxyribozyme "10-23" was active in the same conditions. These results pointed out that the hydrolytic activity of nucleic acid depends on proper tertiary structure of a complex with a substrate. They prove that magnesium ion is not directly involved in catalysis process. On that basis we show the mechanism of RNA hydrolysis catalyzed with nucleic acids at HHP.     

Oligonucleotide facilitators may inhibit or activate a hammerhead ribozyme.

Facilitators are oligonucleotides capable of affecting hammerhead ribozyme activity by interacting with the substrate at the termini of the ribozyme. Facilitator effects were determined in vitro using a system consisting of a ribozyme with 7 nucleotides in every stem sequence and two substrates with inverted facilitator binding sequences. The effects of 9mer and 12mer RNA as well as DNA facilitators which bind either adjacent to the 3'- or 5'-end of the ribozyme were investigated. A kinetic model was developed which allows determination of the apparent dissociation constant of the ribozyme-substrate complex from single turnover reactions. We observed a decreased dissociation constant of the ribozyme-substrate complex due to facilitator addition corresponding to an additional stabilization energy of delta delta G=-1.7 kcal/mol with 3'-end facilitators. The cleavage rate constant was increased by 3'-end facilitators and decreased by 5'-end facilitators. Values for Km were slightly lowered by all facilitators and kcat was increased by 3'-end facilitators and decreased by 5'-end facilitators in our system. Generally the facilitator effects increased with the length of the facilitators and RNA provided greater effects than DNA of the same sequence. Results suggest facilitator influences on several steps of the hammerhead reaction, substrate association, cleavage and dissociation of products. Moreover, these effects are dependent in different manners on ribozyme and substrate concentration. This leads to the conclusion that there is a concentration dependence whether activation or inhibition is caused by facilitators. Conclusions are drawn with regard to the design of hammerhead ribozyme facilitator systems.     

Recent advances in the elucidation of the mechanisms of action of ribozymes

The cleavage of RNA can be accelerated by a number of factors. These factors include an acidic group (Lewis acid) or a basic group that aids in the deprotonation of the attacking nucleophile, in effect enhancing the nucleophilicity of the nucleophile; an acidic group that can neutralize and stabilize the leaving group; and any environment that can stabilize the pentavalent species that is either a transition state or a short-lived intermediate. The catalytic properties of ribozymes are due to factors that are derived from the complicated and specific structure of the ribozyme–substrate complex. It was postulated initially that nature had adopted a rather narrowly defined mechanism for the cleavage of RNA. However, recent findings have clearly demonstrated the diversity of the mechanisms of ribozyme-catalyzed reactions. Such mechanisms include the metal-independent cleavage that occurs in reactions catalyzed by hairpin ribozymes and the general double-metal-ion mechanism of catalysis in reactions catalyzed by the Tetrahymena group I ribozyme. Furthermore, the architecture of the complex between the substrate and the hepatitis delta virus ribozyme allows perturbation of the pK_a of ring nitrogens of cytosine and adenine. The resultant perturbed ring nitrogens appear to be directly involved in acid/base catalysis. Moreover, while high concentrations of monovalent metal ions or polyamines can facilitate cleavage by hammerhead ribozymes, divalent metal ions are the most effective acid/base catalysts under physiological conditions.     

Ion-induced folding of the hammerhead ribozyme: a fluorescence resonance energy transfer study.

The ion-induced folding transitions of the hammerhead ribozyme have been analysed by fluorescence resonance energy transfer. The hammerhead ribozyme may be regarded as a special example of a three-way RNA junction, the global structure of which has been studied by comparing the distances (as energy transfer efficiencies) between the ends of pairs of labelled arms for the three possible end-to-end vectors as a function of magnesium ion concentration. The data support two sequential ion-dependent transitions, which can be interpreted in the light of the crystal structures of the hammerhead ribozyme. The first transition corresponds to the formation of a coaxial stacking between helices II and III; the data can be fully explained by a model in which the transition is induced by a single magnesium ion which binds with an apparent association constant of 8000-10 000 M-1. The second structural transition corresponds to the formation of the catalytic domain of the ribozyme, induced by a single magnesium ion with an apparent association constant of approximately 1100 M-1. The hammerhead ribozyme provides a well-defined example of ion-dependent folding in RNA.     

Design and specificity of hammerhead ribozymes against calretinin mRNA.

We obtained a partial sequence of mouse calretinin mRNA from cDNA clones, and designed hammerhead ribozymes to cleave positions within it. With a view to optimising hammerhead ribozymes for eliminating the mRNA in vivo, we varied the length and sequence of the three duplex 'arms' and measured the cleavage of long RNA substrates in vitro at 37 degrees C (as well as 50 degrees C). Precise cleavage occurred, but it could only go to completion with a large excess of ribozyme. The evidence suggests that the rate-limiting step with a large target is not the cleavage, but the formation of the active ribozyme: substrate complex. The efficiency varied unpredictably according to the target site, the length of the substrate RNA, and the length of the ribozyme; secondary structure in vitro may be responsible. We particularly investigated the degree of sequence-specificity. Some mismatches could be tolerated, but shortening of the total basepairing with the substrate to less than 14 bp drastically reduced activity, implying that interaction with weakly-matched RNAs is unlikely to be a serious problem in vivo. These results suggest that specific and complete cleavage of a mRNA in vivo should be possible, given high-level expression of a ribozyme against a favourable target site.     

Reaction pathway of the trans-acting hepatitis delta virus ribozyme: a conformational change accompanies catalysis.

The hepatitis delta virus (HDV), an infectious human pathogen and satellite of hepatitis B virus, leads to intensified disease symptoms, including progression to liver cirrhosis. Both the circular RNA genome of HDV and its complementary antigenome contain the same cis-cleaving catalytic RNA motif that plays a crucial role in virus replication. Previously, the high-resolution crystal structure of the product form of a cis-acting genomic HDV ribozyme has been determined, while a trans-acting version of the ribozyme was used to dissect the cleavage reaction pathway. Using fluorescence resonance energy transfer (FRET) on a synthetic trans-cleaving form of the ribozyme, we are able to directly observe substrate binding (at a rate constant k(on) of 7.8 x 10(6) M(-1) min(-1) at pH 7.5, 11 mM MgCl(2), and 25 degrees C) and dissociation (at 0.34 min(-1)). Steady-state and time-resolved FRET experiments in solution and in nondenaturing gels reveal that the substrate (precursor) complex is slightly more compact (by approximately 3 A) than the free ribozyme, yet becomes significantly extended (by approximately 15 A) upon cleavage and product complex formation. We also find that trans cleavage is characterized by a high transition-state entropy (-26 eu). We propose that the significant global conformational change that we observe between the precursor and product structures occurs on the reaction trajectory into a constrained product complex-like transition state. Our observations may present the structural basis of the recently described utilization of intrinsic substrate binding energy to the overall catalytic rate enhancement by the trans-acting HDV ribozyme.