General base catalysis for cleavage by the active-site cytosine of the hepatitis delta virus ribozyme: QM/MM calculations establish chemical feasibility

The hepatitis delta virus (HDV) ribozyme is an RNA motif embedded in human pathogenic HDV RNA. Previous experimental studies have established that the active-site nucleotide C75 is essential for self-cleavage of the ribozyme, although its exact catalytic role in the process remains debated. Structural data from X-ray crystallography generally indicate that C75 acts as the general base that initiates catalysis by deprotonating the 2'-OH nucleophile at the cleavage site, while a hydrated magnesium ion likely protonates the 5'-oxygen leaving group. In contrast, some mechanistic studies support the role of C75 acting as general acid and thus being protonated before the reaction. We report combined quantum chemical/molecular mechanical calculations for the C75 general base pathway, utilizing the available structural data for the wild type HDV genomic ribozyme as a starting point. Several starting configurations differing in magnesium ion placement were considered and both one-dimensional and two-dimensional potential energy surface scans were used to explore plausible reaction paths. Our calculations show that C75 is readily capable of acting as the general base, in concert with the hydrated magnesium ion as the general acid. We identify a most likely position for the magnesium ion, which also suggests it acts as a Lewis acid. The calculated energy barrier of the proposed mechanism, approximately 20 kcal/mol, would lower the reaction barrier by approximately 15 kcal/mol compared with the uncatalyzed reaction and is in good agreement with experimental data.     

Quantum Mechanical/Molecular Mechanical Study of the HDV Ribozyme: Impact of the Catalytic Metal Ion on the Mechanism

A recent crystal structure of the precleaved HDV ribozyme along with biochemical data support a mechanism for phosphodiester bond self-cleavage in which C75 acts as a general acid and bound Mg(2+) ion acts as a Lewis acid. Herein this precleaved crystal structure is used as the basis for quantum mechanical/molecular mechanical calculations. These calculations indicate that the self-cleavage reaction is concerted with a phosphorane-like transition state when a divalent ion, Mg(2+) or Ca(2+), is bound at the catalytic site but is sequential with a phosphorane intermediate when a monovalent ion, such as Na(+), is at this site. Electrostatic potential calculations suggest that the divalent metal ion at the catalytic site lowers the pK(a) of C75, leading to the concerted mechanism in which the proton is partially transferred to the leaving group in the phosphorane-like transition state. These observations are consistent with experimental data, including pK(a) measurements, reaction kinetics, and proton inventories with divalent and monovalent ions.     

Mechanistic characterization of the HDV genomic ribozyme: solvent isotope effects and proton inventories in the absence of divalent metal ions support C75 as the general acid

The hepatitis delta virus (HDV) ribozyme uses the nucleobase C75 and a hydrated Mg(2+) ion as the general acid-base catalysts in phosphodiester bond cleavage at physiological salt. A mechanistic framework has been advanced that involves one Mg(2+)-independent and two Mg(2+)-dependent channels. The rate-pH profile for wild-type (WT) ribozyme in the Mg(2+)-free channel is inverted relative to the fully Mg(2+)-dependent channel, with each having a near-neutral pKa. Inversion of the rate-pH profile was used as the crux of a mechanistic argument that C75 serves as general acid both in the presence and absence of Mg(2+). However, subsequent studies on a double mutant (DM) ribozyme suggested that the pKa observed for WT in the absence of Mg(2+) arises from ionization of C41, a structural nucleobase. To investigate this further, we acquired rate-pH/pD profiles and proton inventories for WT and DM in the absence of Mg(2+). Corrections were made for effects of ionic strength on hydrogen ion activity and pH meter readings. Results are accommodated by a model wherein the Mg(2+)-free pKa observed for WT arises from ionization of C75, and DM reactivity is compromised by protonation of C41. The Br?nsted base appears to be water or hydroxide ion depending on pH. The observed pKa's are related to salt-dependent pH titrations of a model oligonucleotide, as well as electrostatic calculations, which support the local environment for C75 in the absence of Mg(2+) being similar to that in the presence of Mg(2+) and impervious to bulk ions. Accordingly, the catalytic role of C75 as the general acid does not appear to depend on divalent ions or the identity of the Br?nsted base.     

Direct measurement of a pK(a) near neutrality for the catalytic cytosine in the genomic HDV ribozyme using Raman crystallography

The hepatitis delta virus (HDV) ribozyme uses a cytosine to facilitate general acid-base catalysis. Biochemical studies suggest that C75 has a pKa perturbed to near neutrality. To measure this pKa directly, Raman spectra were recorded on single ribozyme crystals using a Raman microscope. A spectral feature arising from a single neutral cytosine was identified at 1528 cm(-1). At low pH, this mode was replaced with a new spectral feature. Monitoring these features as a function of pH revealed pKa values for the cytosine that couple anticooperatively with Mg2+ binding, with values of 6.15 and 6.40 in the presence of 20 and 2 mM Mg2+, respectively. These pKa values agree well with those obtained from ribozyme activity experiments in solution. To correlate the observed pKa with a specific nucleotide, crystals of C75U, which is catalytically inactive, were examined. The Raman difference spectra show that this mutation does not affect the conformation of the ribozyme. However, crystals of C75U did not produce a signal from a protonatable cytosine, providing strong evidence that protonation of C75 is being monitored in the wild-type ribozyme. These studies provide the first direct physical measurement of a pKa near neutrality for a catalytic residue in a ribozyme and show that ribozymes, like their protein enzyme counterparts, can optimize the pKa of their side chains for proton transfer.     

Theoretical examination of two opposite mechanisms proposed for hepatitis delta virus ribozyme

Two mechanisms were previously proposed for the hepatitis delta virus (HDV) ribozyme where an active-site cytosine residue (C75) either functioned as a general base to deprotonate the 2'-OH at the rupture site or as a general acid to protonate the O5' leaving group. Here, we reported the first theoretical examination of the two mechanisms using a combination of the quantum mechanics (QM)/molecular mechanics (MM), molecular dynamics (MD), and near-attack-conformation (NAC) techniques. Our theoretical results supported the C75-acid mechanism, which was demonstrated to have an unfavorable starting geometry (in agreement with the crystallographic data) but a significantly lower energy barrier as compared to the C75-base mechanism. Therefore, the chemical details of the transition state in the HDV ribozyme may dramatically differ from those inferred from the structural studies.     

Role of Mg2+ in hammerhead ribozyme catalysis from molecular simulation

Molecular dynamics simulations have been performed to investigate the role of Mg2+ in the full-length hammerhead ribozyme cleavage reaction. In particular, the aim of this work is to characterize the binding mode and conformational events that give rise to catalytically active conformations and stabilization of the transition state. Toward this end, a series of eight 12 ns molecular dynamics simulations have been performed with different divalent metal binding occupations for the reactant, early and late transition state using recently developed force field parameters for metal ions and reactive intermediates in RNA catalysis. In addition, hybrid QM/MM calculations of the early and late transition state were performed to study the proton-transfer step in general acid catalysis that is facilitated by the catalytic Mg2+ ion. The simulations suggest that Mg2+ is profoundly involved in the hammerhead ribozyme mechanism both at structural and catalytic levels. Binding of Mg2+ in the active site plays a key structural role in the stabilization of stem I and II and to facilitate formation of near attack conformations and interactions between the nucleophile and G12, the implicated general base catalyst. In the transition state, Mg2+ binds in a bridging position where it stabilizes the accumulated charge of the leaving group while interacting with the 2'OH of G8, the implicated general acid catalyst. The QM/MM simulations provide support that, in the late transition state, the 2'OH of G8 can transfer a proton to the leaving group while directly coordinating the bridging Mg2+ ion. The present study provides evidence for the role of Mg2+ in hammerhead ribozyme catalysis. The proposed simulation model reconciles the interpretation of available experimental structural and biochemical data, and provides a starting point for more detailed investigation of the chemical reaction path with combined QM/MM methods.     

Probing general base catalysis in the hammerhead ribozyme

Recent structural and computational studies have shed new light on the catalytic mechanism and active site structure of the RNA cleaving hammerhead ribozyme. Consequently, specific ribozyme functional groups have been hypothesized to be directly involved in general/acid base catalysis. In order to test this hypothesis, we have developed an affinity label to identify the functional general base in the S. mansoni hammerhead ribozyme. The ribozyme was reacted with a substrate analogue bearing a 2'-bromoacetamide group in place of the nucleophilic 2'-hydroxyl group which would normally be deprotonated by a general base. The electrophilic 2'-bromoacetamide group is poised to alkylate the general base, which is subsequently identified by footprinting analysis. Herein, we demonstrate alkylation of N1 of G12 in the hammerhead ribozyme in a pH and [Mg(2+)] dependent manner that is consistent with the native cleavage reaction. These results provide substantial evidence that deprotonated N1 of G12 functions directly as a general base in the hammerhead ribozyme; moreover, our experiments provide evidence that the pKa of G12 is perturbed downward in the context of the active site structure. We also observed other pH-independent alkylations, which do not appear to reflect the catalytic mechanism, but offer further insight into ribozyme conformation and structure.     

Acid-base and metal-ion binding properties of the RNA dinucleotide uridylyl-(5'-->3')-[5']uridylate (pUpU3-)

It is well known that Mg2+ and other divalent metal ions bind to the phosphate groups of nucleic acids. Subtle differences in the coordination properties of these metal ions to RNA, especially to ribozymes, determine whether they either promote or inhibit catalytic activity. The ability of metal ions to coordinate simultaneously with two neighboring phosphate groups is important for ribozyme structure and activity. However, such an interaction has not yet been quantified. Here, we have performed potentiometric pH titrations to determine the acidity constants of the protonated dinucleotide H2(pUpU)-, as well as the binding properties of pUpU3- towards Mg2+, Mn2+, Cd2+, Zn2+, and Pb2+. Whereas Mg2+, Mn2+, and Cd2+ only bind to the more basic 5'-terminal phosphate group, Pb2+, and to a certain extent also Zn2+, show a remarkably enhanced stability of the [M(pUpU)]- complex. This can be attributed to the formation of a macrochelate by bridging the two phosphate groups within this dinucleotide by these metal ions. Such a macrochelate is also possible in an oligonucleotide, because the basic structural units are the same, despite the difference in charge. The formation degrees of the macrochelated species of [Zn(pUpU)]- and [Pb(pUpU)]- amount to around 25 and 90 %, respectively. These findings are important in the context of ribozyme and DNAzyme catalysis, and explain, for example, why the leadzyme could be selected in the first place, and why this artificial ribozyme is inhibited by other divalent metal ions, such as Mg2+.     

Coordination environment of a site-bound metal ion in the hammerhead ribozyme determined by 15N and 2H ESEEM spectroscopy

Although site-bound Mg2+ ions have been proposed to influence RNA structure and function, establishing the molecular properties of such sites has been challenging due largely to the unique electrostatic properties of the RNA biopolymer. We have previously determined that, in solution, the hammerhead ribozyme (a self-cleaving RNA) has a high-affinity metal ion binding site characterized by a K(d,app) < 10 microM for Mn2+ in 1 M NaCl and speculated that this site has functional importance in the ribozyme cleavage reaction. Here we determine both the precise location and the hydration level of Mn2+ in this site using ESEEM (electron spin-echo envelope modulation) spectroscopy. Definitive assignment of the high-affinity site to the activity-sensitive A9/G10.1 region is achieved by site-specific labeling of G10.1 with 15N guanine. The coordinated metal ion retains four water ligands as measured by 2H ESEEM spectroscopy. The results presented here show that a functionally important, specific metal binding site is uniquely populated in the hammerhead ribozyme even in a background of high ionic strength. Although it has a relatively high thermodynamic affinity, this ion remains partially hydrated and is chelated to the RNA by just two ligands.     

Capreomycin and hygromycin B modulate the catalytic activity of the delta ribozyme in a manner that depends on the protonation and complexation with Cu2+ ions of these antibiotics

Catalytic RNA molecules (ribozymes) have often been used for the testing of interactions of antibiotics with ribonucleic acids. We showed that the impact of capreomycin and hygromycin B on delta ribozyme catalysis might change dramatically, from stimulation to inhibition, depending on conditions. In order to evaluate possible mechanisms of modulation of the ribozyme catalytic activity we used our earlier data on species distribution for protonated forms of capreomycin and hygromycin B and their complexes with Cu(2+) ions at different pH values. We proposed that, upon inhibition, the protonated amino group of capreomycin was located in the ribozyme catalytic cleft interfering with binding catalytic Mg(2+). Such a mechanism was also supported by the results of ribozyme inhibition with capreomycin complexed with Cu(2+). The effects of stimulation of the delta ribozyme activity by capreomycin and hygromycin B were less pronounced than inhibition. Possibly, the amino functions of these antibiotics might be involved in a general acid-base catalysis performed by the ribozyme, acting as proton acceptors/donors.     

Detection of innersphere interactions between magnesium hydrate and the phosphate backbone of the HDV ribozyme using Raman crystallography

A Raman microscope and Raman difference spectroscopy are used to detect the vibrational signature of RNA-bound magnesium hydrate in crystals of hepatitis delta virus (HDV) ribozyme and to follow the effects of magnesium hydrate binding to the nonbridging phosphate oxygens in the phosphodiester backbone. There is a correlation between the Raman intensity of the innersphere magnesium hydrate signature peak, near 322 cm-1, and the intensity of the PO2- symmetric stretch, near 1100 cm-1, perturbed by magnesium binding, demonstrating direct observation of -PO2-...Mg2+(H2O)x innersphere complexes. The complexes may be pentahydrates (x = 5) and tetrahydrates (x = 4). The assignment of the Raman feature near 322 cm-1 to a magnesium hydrate species is confirmed by isotope shifts observed in D2O and H218O that are semiquantitatively reproduced by calculations. The standardized intensity changes in the 1100 cm-1 PO2- feature seen upon magnesium hydrate binding indicates that there are approximately 5 innersphere Mg2+...-O2P contacts per HDV molecule when the crystal is exposed to a solution containing 20 mM magnesium.     

Fast cleavage kinetics of a natural hammerhead ribozyme

The hammerhead ribozyme is a small RNA motif that catalyzes the cleavage and ligation of RNA. The well-studied minimal hammerhead motif is inactive under physiological conditions and requires high Mg(2+) concentrations for efficient cleavage. In contrast, natural hammerheads are active under physiological conditions and contain motifs outside the catalytic core that lower the requirement for Mg(2+). Single-turnover kinetics were used here to characterize the Mg(2+) and pH dependence for cleavage of a trans-cleaving construct of the Schistosoma mansoni natural hammerhead ribozyme. Compared to the minimal hammerhead motif, the natural Schistosoma ribozyme requires 100-fold less Mg(2+) to achieve a cleavage rate of 1 min(-1). The improved catalysis results from tertiary interactions between loops in stems I and II and likely arises from increasing the population of the active conformation. Under optimum pH and Mg(2+) conditions this ribozyme cleaves at over 870 min(-1) at 25 degrees C, further demonstrating the impressive catalytic power of this ribozyme.     

Active site labeling of G8 in the hairpin ribozyme: implications for structure and mechanism

There is mounting evidence that suggests that general acid/base catalysis is operative in the hairpin ribozyme, with analogy to the protein enzyme RNaseA. Nevertheless, the extent of general base catalysis as well as the identity of the specific chemical groups responsible remains the subject of some controversy. An affinity label has previously been used to alkylate histidine 12 (His12), the active general base in RNaseA. To date, no such experiment has been applied to a ribozyme. We have synthesized the analogous affinity label for the hairpin ribozyme with an electrophilic 2'-bromoacetamide group in lieu of the 2'-hydroxyl (2'OH) at the substrate cleavage site and show that guanosine 8 (G8) of the hairpin ribozyme is specifically alkylated, most likely at the N1 position. This evidence strongly implicates N1 of G8 in active site chemistry. By direct analogy to RNase A, these findings could be consistent with the hypothesis that deprotonated G8 residue functions as a general base in the hairpin ribozyme. Other mechanistic possibilities for N1 of G8 such as indirect general base catalysis mediated by a water molecule or transition state stabilization could also be consistent with our findings.     

Probing the binding of Tb(III) and Eu(III) to the hammerhead ribozyme using luminescence spectroscopy

BACKGROUND: Divalent metal ions serve as structural as well as catalytic cofactors in the hammerhead ribozyme reaction. The natural cofactor in these reactions is Mg(II), but its spectroscopic silence makes it difficult to study. We previously showed that a single Tb(III) ion inhibits the hammerhead ribozyme by site-specific competition for a Mg(II) ion and therefore can be used as a spectroscopic probe for the Mg(II) it replaces. RESULTS: Lanthanide luminescence spectroscopy was used to study the coordination environment around Tb(III) and Eu(III) ions bound to the structurally well-characterized site on the hammerhead ribozyme. Sensitized emission and direct excitation experiments show that a single lanthanide ion binds to the ribozyme under these conditions and that three waters of hydration are displaced from the Tb(III) upon binding the RNA. Furthermore, we show that these techniques allow the comparison of binding affinities for a series of ions to this site. The binding affinities for ions at the G5 site correlates linearly with the function Z(2)/r of the aqua ion (where Z is the charge and r is the radius of the ion). CONCLUSIONS: This study compares the crystallographic nature of the G5 metal-binding site with solution measurements and gives a clearer picture of the coordination environment of this ion. These results provide one of the best characterized metal-binding sites from a ribozyme, so we use this information to compare the RNA site with that of typical metalloproteins.     

Metal-ion-coordinating properties of the dinucleotide 2'-deoxyguanylyl(5'-->3')-2'-deoxy-5'-guanylate (d(pGpG)3-): isomeric equilibria including macrochelated complexes relevant for nucleic acids

The interaction between divalent metal ions and nucleic acids is well known, yet knowledge about the strength of binding of labile metal ions at the various sites is very scarce. We have therefore studied the stabilities of complexes formed between the nucleic acid model d(pGpG) and the essential metal ions Mg2+ and Zn2+ as well as with the generally toxic ions Cd2+ and Pb2+ by potentiometric pH titrations; all four ions are of relevance in ribozyme chemistry. A comparison of the present results with earlier data obtained for M(pUpU)- complexes allows the conclusion that phosphate-bound Mg2+ and Cd2+ form macrochelates by interaction with N7, whereas the also phosphate-coordinated Pb2+ forms a 10-membered chelate with the neighboring phosphate diester bridge. Zn2+ forms both types of chelates with formation degrees of about 91% and 2.4% for Zn[d(pGpG)]cl/N7 and Zn[d(pGpG)]-cl/PO, respectively; the open form with Zn2+ bound only to the terminal phosphate group, Zn[d(pGpG)]-op, amounts to about 6.8 %. The various intramolecular equilibria have also been quantified for the other metal ions. Zn2+, Cu2+, and Cd2+ also form macrochelates in the monoprotonated M[H;d(pGpG)] species (the proton being at the terminal phosphate group), that is, the metal ion at N7 interacts to some extent with the P(O)2(OH)- group. Thus, this study demonstrates that the coordinating properties of the various metal ions toward a pGpG unit in a nucleic acid differ: Mg2+, Zn2+, and Cd2+ have a significant tendency to bridge the distance between N7 and the phosphate group of a (d)GMP unit, although to various extents, whereas Pb2+ (and possibly Ca2+) prefer a pure phosphate coordination.     

Determination of the intrinsic affinities of multiple site-specific Mg(2+) ions coordinated to domain 6 of a group II intron ribozyme

Group II introns are large metallo-ribozymes that use divalent metal ions in folding and catalysis. The 3'-terminal domain 6 (D6) contains a conserved adenosine whose 2'-OH group acts as the nucleophile in the first splicing step. In the hierarchy of folding, D6 binds last into the active site. In order to investigate and understand the folding process to the catalytically active intron structure, it is important to know the individual binding affinities of Mg2+ ions to D6. We recently studied the solution structure of a 27 nucleotide long D6 (D6-27) from the mitochondrial yeast group II intron Sc.ai5gamma, also identifying five Mg2+ binding sites including the one at the 5'-terminal phosphate residues. Mg2+ coordination to the 5'-terminal di- and triphosphate groups is strongest (e.g., log KA,TP = 4.55 +/- 0.10) and is evaluated here in detail for the first time. The other four binding sites within D6-27 are filled simultaneously (e.g., log KA,BR = 2.38 +/- 0.06) and thus compete for the free Mg2+ ions in solution, having a distinct influence on the individual affinities of the various sites. For the first time, we take this competition into account to obtain the intrinsic binding constants, describing a method that is generally applicable. Our data illustrates that any RNA molecule undergoing tertiary contacts to a second RNA molecule first needs to be loaded evenly and specifically with metal ions to compensate for the repulsion between the negatively charged RNA molecules.     

Significant kinetic solvent isotope effects in folding of the catalytic RNA from the hepatitis delta virus

The exchange of deuterium for hydrogen in water often produces solvent kinetic isotope effects (KSIEs) on the rate constants associated with enzyme reactions, including those catalyzed by RNA. Recently, KSIEs have been used to show that proton transfer occurs in the rate-limiting step of cleavage by the hepatitis delta virus (HDV) ribozyme and other catalytic RNAs. To test the underlying assumption that KSIEs are related to the chemistry step of ribozyme-mediated cleavage reactions, we developed fluorescence resonance energy transfer assays to measure KSIEs on the rate constants of conformational changes associated with substrate binding and dissociation by a trans-acting HDV ribozyme. We observe comparable KSIEs ( approximately 2-2.5-fold) of rate constants of conformational change and cleavage, while proton inventory experiments are consistent with a shift in the ensemble of transition states upon increase of D2O in the solvent. Taken together, these results challenge the common assumption that pL profiles of RNA-catalyzed reactions yielding a pKa and KSIE necessarily provide evidence for an ionization (chemistry) step to be rate-limiting. They also suggest that an unusual proton inventory may provide a signature for a conformational change contributing to the rate-limiting step.     

Complex study of nickel skeleton catalysts: V. Study of novel non-pyrophoric nickel skeleton catalysts by derivatograph

The structure and properties of non-pyrophoric skeleton catalysts prepared from NiSi, NiAlSi, NiMg and NiZn alloys have been studied—apart from other thermal methods—by means of the derivatograph.Our experimental results have contributed to the explanation of the non-pyrophoric behaviour of these catalysts. We demonstrated that the desorption of the hydrogen content in our catalysts is not accompanied by the oxidation of active nickel. This oxidation takes place only at higher temperatures, above 200°C, at a rate proportional to the amount of active nickel.Other constituents of the catalysts (adsorbed water, hydroxide content) were also determined from the experimental data. The outstandingly high Mg(OH)₂ content of the NiMg catalyst indicates that its structure is dissimilar: Mg(OH)₂ also acts as support for the catalyst.     

Separate metal requirements for loop interactions and catalysis in the extended hammerhead ribozyme

An extended hammerhead ribozyme derived from Schistosoma mansoni, including conserved loops in stems I and II, has been examined to directly monitor the relationship between docking of loops and its activity using site-directed spin labeling (SDSL) and EPR spectroscopy. Dynamics with EPR spectroscopy and fast-quench kinetics measurements have shown that the docking of stems I and II occurs at low Mg2+ concentrations ([Mg2+]1/2,dock = 0.7 mM, 0.1 M NaCl), but a much weaker Mg2+ interaction ([Mg2+]1/2,cat approximately 90 mM) increases activity to very high maximum rates of approximately 1 s-1 at 0.1 M Na+ and pH 7.0.