General acid-base catalysis in the mechanism of a hepatitis delta virus ribozyme

Contribution of the closing base pair to exceptional stability in RNA tetraloops: roles for molecular mimicry and electrostatic factors

Hairpins are common RNA secondary structures that play multiple roles in nature. Tetraloops are the most frequent RNA hairpin loops and are often phylogenetically conserved. For both the UNCG and GNRA families, CG closing base pairs (cbps) confer exceptional thermodynamic stability but the molecular basis for this has remained unclear. We propose that, despite having very different overall folds, these two tetraloop families achieve stability by presenting the same functionalities to the major groove edge of the CG cbp. Thermodynamic contributions of this molecular mimicry were investigated using substitutions at the nucleobase and functional group levels. By either interrupting or deleting loop-cbp electrostatic interactions, which were identified by solving the nonlinear Poisson-Boltzmann (NLPB) equation, stability changed in a manner consistent with molecular mimicry. We also observed a linear relationship between DeltaG(o)(37) and log[Na(+)] for both families, and loops with a CG cbp had a decreased dependence of stability on salt. NLPB calculations revealed that, for both UUCG and GAAA tetraloops, the GC cbp form has a higher surface charge density, although it arises from changes in loop compaction for UUCG and changes in loop configuration for GAAA. Higher surface charge density leads to stronger interactions of GC cbp loops with solvent and salt, which explains the correlation between experimental and calculated trends of free energy with salt. Molecular mimicry as evidenced in these two stable but otherwise unrelated tetraloops may underlie common functional roles in other RNA and DNA motifs.

Identification of catalytic metal ion ligands in ribozymes

Site-bound metal ions participate in the catalytic mechanisms of many ribozymes. Understanding these mechanisms therefore requires knowledge of the specific ligands on both substrate and ribozyme that coordinate these catalytic metal ions. A number of different structural and biochemical strategies have been developed and refined for identifying metal ion binding sites within ribozymes, and for assessing the catalytic contributions of the metal ions bound at those sites. We review these approaches and provide examples of their application, focusing in particular on metal ion rescue experiments and their roles in the construction of the transition state models for the Tetrahymena group I and RNase P ribozymes.

Comparative enzymology and structural biology of RNA self-cleavage

Self-cleaving hammerhead, hairpin, hepatitis delta virus, and glmS ribozymes comprise a family of small catalytic RNA motifs that catalyze the same reversible phosphodiester cleavage reaction, but each motif adopts a unique structure and displays a unique array of biochemical properties. Recent structural, biochemical, and biophysical studies of these self-cleaving RNAs have begun to reveal how active site nucleotides exploit general acid-base catalysis, electrostatic stabilization, substrate destabilization, and positioning and orientation to reduce the free energy barrier to catalysis. Insights into the variety of catalytic strategies available to these model RNA enzymes are likely to have important implications for understanding more complex RNA-catalyzed reactions fundamental to RNA processing and protein synthesis.

Direct measurement of the ionization state of an essential guanine in the hairpin ribozyme

Active site guanines are critical for self-cleavage reactions of several ribozymes, but their precise functions in catalysis are unclear. To learn whether protonated or deprotonated forms of guanine predominate in the active site, microscopic pKa values were determined for ionization of 8-azaguanosine substituted for G8 in the active site of a fully functional hairpin ribozyme in order to determine microscopic pKa values for 8-azaguanine deprotonation from the pH dependence of fluorescence. Microscopic pKa values above 9 for deprotonation of 8-azaguanine in the active site were about 3 units higher than apparent pKa values determined from the pH dependence of self-cleavage kinetics. Thus, the increase in activity with increasing pH does not correlate with deprotonation of G8, and most of G8 is protonated at neutral pH. These results do not exclude a role in proton transfer, but a simple interpretation is that G8 functions in the protonated form, perhaps by donating hydrogen bonds.

Hepatitis D: thirty years after

The key to the discovery of the Hepatitis D Virus (HDV) was the description in Turin, Italy in the mid-1970s of the delta antigen and antibody in carriers of the hepatitis B surface antigen. The new antigen was first thought to be a marker of the Hepatitis B Virus (HBV) and in view of its intricate true nature, it would have possibly died away as another odd antigenic subtype of HBV, like many that were described in the 1970s. Fortunately, instead, a collaboration started in 1978 between the Turin group, and the National Institute of Health and Georgetown University in the US. With American facilities and expertise this collaboration led just a year later, in 1979, to the unfolding of an unexpected and amazing chapter in virology. Experiments in chimpanzees demonstrated that the delta antigen was not a component of the HBV but of a separate defective virus requiring HBV for its infection; it was named the hepatitis D virus to conform to the nomenclature of hepatitis viruses and classified within the genus Deltavirus. The animal experiments were also seminal in proposing to future clinical interpretation, the paradigm of a pathogenic infection (hepatitis D), that could develop only in HBV-infected patients, was mainly transmitted by superinfection of HDV on chronic HBV carriers and had the ability to strongly inhibit the helper HBV. The discovery of the HDV has driven three directions of further research: (1) The understanding of the replicative and infectious mechanisms of the HDV. (2) The assessment of its epidemiological and medical impact. (3) The search for a therapy for chronic hepatitis D (CHD). This review summarizes the progress achieved in each field of research in the thirty years that have passed since the discovery of HDV.

A self-cleaving DNA enzyme modified with amines, guanidines and imidazoles operates independently of divalent metal cations (M2+)

The selection of modified DNAzymes represents an important endeavor in expanding the chemical and catalytic properties of catalytic nucleic acids. Few examples of such exist and to date, there is no example where three different modified bases have been simultaneously incorporated for catalytic activity. Herein, dCTP, dATP and dUTP bearing, respectively, a cationic amine, an imidazole and a cationic guanidine, were enzymatically polymerized on a DNA template for the selection of a highly functionalized DNAzyme, called DNAzyme 9-86, that catalyzed (M(2+))-independent self-cleavage under physiological conditions at a single ribo(cytosine)phosphodiester linkage with a rate constant of (0.134 +/- 0.026) min(-1). A pH rate profile analysis revealed pK(a)'s of 7.4 and 8.1, consistent with both general acid and base catalysis. The presence of guanidinium cations permits cleavage at significantly higher temperatures than previously observed for DNAzymes with only amines and imidazoles. Qualitatively, DNAzyme 9-86 presents an unprecedented ensemble of synthetic functionalities while quantitatively it expresses one of the highest reported values for any self-cleaving nucleic acid when investigated under M(2+)-free conditions at 37 degrees C.

Synthesis of novel C4-linked C2-imidazole ribonucleoside phosphoramidite and its application to probing the catalytic mechanism of a ribozyme

The synthesis of a novel C4-linked C2-imidazole ribonucleoside phosphoramidite (ICN-C2-PA 1) with a two-carbon linker between imidazole and ribose moieties is described. In the phosphoramidite, POM and 2-cyanoethyl groups were selected to protect the endocyclic amine function of imidazole and the 2'-hydroxyl function of D-ribose, respectively. The C2-imidazole nucleoside, a flexible structural mimic of a purine nucleobase, was successfully incorporated using ICN-C2-PA 1 into position 638 of the VS ribozyme through 2'-TBDMS chemistry to study the role of G638 in general acid-base catalysis. The modified VS ribozyme (G638C2Imz) exhibited significantly greater catalytic activity than observed with the C0-imidazole that has no carbon atoms linking the ribose and the C4-imidazole. Imidazole nucleoside analogues with variable spacer lengths could provide a valuable general methodology for exploring the catalytic mechanisms of ribozymes.

A catalytic metal ion interacts with the cleavage Site G.U wobble in the HDV ribozyme

The HDV ribozyme self-cleaves by a chemical mechanism involving general acid-base catalysis to generate 2',3'-cyclic phosphate and 5'-hydroxyl termini. Biochemical studies from several laboratories have implicated C75 as the general acid and hydrated magnesium as the general base. We have previously shown that C75 has a pK(a) shifted >2 pH units toward neutrality [Gong, B., Chen, J. H., Chase, E., Chadalavada, D. M., Yajima, R., Golden, B. L., Bevilacqua, P. C., and Carey, P. R. (2007) J. Am. Chem. Soc. 129, 13335-13342], while in crystal structures, it is well-positioned for proton transfer. However, no evidence for a hydrated magnesium poised to serve as a general base in the reaction has been observed in high-resolution crystal structures of various reaction states and mutants. Herein, we use solution kinetic experiments and parallel Raman crystallographic studies to examine the effects of pH on the rate and Mg(2+) binding properties of wild-type and 7-deazaguanosine mutants of the HDV ribozyme. These data suggest that a previously unobserved hydrated magnesium ion interacts with N7 of the cleavage site G.U wobble base pair. Integrating this metal ion binding site with the available crystal structures provides a new three-dimensional model for the active site of the ribozyme that accommodates all available biochemical data and appears competent for catalysis. The position of this metal is consistent with a role of a magnesium-bound hydroxide as a general base as dictated by biochemical data.

The Small Ribozymes: Common and Diverse Features Observed through the FRET Lens

The hammerhead, hairpin, HDV, VS and glmS ribozymes are the five known, naturally occurring catalytic RNAs classified as the "small ribozymes". They share common reaction chemistry in cleaving their own backbone by phosphodiester transfer, but are diverse in their secondary and tertiary structures, indicating that Nature has found at least five independent solutions to a common chemical task. Fluorescence resonance energy transfer (FRET) has been extensively used to detect conformational changes in these ribozymes and dissect their reaction pathways. Common and diverse features are beginning to emerge that, by extension, highlight general biophysical properties of non-protein coding RNAs.

Fluorine substituted adenosines as probes of nucleobase protonation in functional RNAs

Ionized nucleobases are required for folding, conformational switching, or catalysis in a number of functional RNAs. A common strategy to study these sites employs nucleoside analogues with perturbed pKa, but the interpretation of these studies is often complicated by the chemical modification introduced, in particular modifications that add, remove, or translocate hydrogen bonding groups in addition to perturbing pKa values. In the present study we present a series of fluorine substituted adenosine analogues that produce large changes in N1 pKa values with minimal structural perturbation. These analogues include fluorine for hydrogen substitutions in the adenine ring of adenosine and 7-deaza-adenosine with resulting N1 pKa values spanning more than 4 pKa units. To demonstrate the utility of these analogues we have conducted a nucleotide analogue interference mapping (NAIM) study on a self-ligating construct of the Varkud Satellite (VS) ribozyme. We find that each of the analogues is readily incorporated by T7 RNA polymerase and produces fully active transcripts when substituted at the majority of sites. Strong interferences are observed for three sites known to be critical for VS ribozyme function, most notably A756. Substitutions at A756 lead to slight enhancements in activity for elevated pKa analogues and dramatic interferences in activity for reduced pKa analogues, supporting the proposed catalytic role for this base. The structural similarity of these analogues, combined with their even incorporation and selective interference, provides an improved method for identifying sites of adenosine protonation in a variety of systems.

Comparison of the acid-base properties of ribose and 2'-deoxyribose nucleotides

The extent to which the replacement of a ribose unit by a 2'-deoxyribose unit influences the acid-base properties of nucleotides has not hitherto been determined in detail. In this study, by potentiometric pH titrations in aqueous solution, we have measured the acidity constants of the 5'-di- and 5'-triphosphates of 2'-deoxyguanosine [i.e., of H(2)(dGDP)(-) and H(2)(dGTP)(2-)] as well as of the 5'-mono-, 5'-di-, and 5'-triphosphates of 2'-deoxyadenosine [i.e., of H(2)(dAMP)(+/-), H(2)(dADP)(-), and H(2)(dATP)(2-)]. These 12 acidity constants (of the 56 that are listed) are compared with those of the corresponding ribose derivatives (published data) measured under the same experimental conditions. The results show that all protonation sites in the 2'-deoxynucleotides are more basic than those in their ribose counterparts. The influence of the 2'-OH group is dependent on the number of 5'-phosphate groups as well as on the nature of the purine nucleobase. The basicity of N7 in guanine nucleotides is most significantly enhanced (by about 0.2 pK units), while the effect on the phosphate groups and the N1H or N1H(+) sites is less pronounced but clearly present. In addition, (1)H NMR chemical shift change studies in dependence on pD in D(2)O have been carried out for the dAMP, dADP, and dATP systems, which confirmed the results from the potentiometric pH titrations and showed the nucleotides to be in their anti conformations. Overall, our results are not only of relevance for metal ion binding to nucleotides or nucleic acids, but also constitute an exact basis for the calculation, determination, and understanding of perturbed pK(a) values in DNAzymes and ribozymes, as needed for the delineation of acid-base mechanisms in catalysis.

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.

The mechanism of acidic hydrolysis of esters explains the HDV ribozyme activity

The hepatitis delta virus (HDV) ribozyme is an RNA enzyme that catalyzes the site-specific trans-esterification reaction. Using high hydrostatic pressure (HHP) technique we showed that HDV ribozyme catalyzes the reaction of RNA cleavage in the absence of magnesium ions according to mechanism of acidic hydrolysis of esters. HHP induces changes of water structure, lowering pH and effect ribozyme catalytic site structure formation without magnesium. HHP, similarly to magnesium ion at ambient pressure stabilizes the higher order RNA structure of HDV, but Mg(2+) is not involved in the catalysis. Our results clearly support the new mechanism of HDV hydrolysis and show advantages of using HHP in analysis of macromolecules interaction.

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.

Mechanistic characterization of the HDV genomic ribozyme: the cleavage site base pair plays a structural role in facilitating catalysis

The hepatitis delta virus (HDV) ribozyme occurs in the genomic and antigenomic strands of the HDV RNA and within mammalian transcriptomes. Previous kinetic studies suggested that a wobble pair (G*U or A(+)*C) is preferred at the cleavage site; however, the reasons for this are unclear. We conducted sequence comparisons, which indicated that while G*U is the most prevalent combination at the cleavage site, G-C occurs to a significant extent in genomic HDV isolates, and G*U, G-C, and A-U pairs are present in mammalian ribozymes. We analyzed the folding of genomic HDV ribozymes by free energy minimization and found that variants with purine-pyrimidine combinations at the cleavage site are predicted to form native structures while pyrimidine-purine combinations misfold, consistent with earlier kinetic data and sequence comparisons. To test whether the cleavage site base pair contributes to catalysis, we characterized the pH and Mg(2+)-dependence of reaction kinetics of fast-folding genomic HDV ribozymes with cleavage site base pair purine-pyrimidine combinations: G*U, A-U, G-C, and A(+)*C. Rates for these native-folding ribozymes displayed highly similar pH and Mg(2+) concentration dependencies, with the exception of the A(+)*C ribozyme, which deviated at high pH. None of the four ribozymes underwent miscleavage. These observations support the A(+)*C ribozyme as being more active with a wobble pair at the cleavage site than with no base pair at all. Overall, the data support a model in which the cleavage site base pair provides a structural role in catalysis and does not need to be a wobble pair.

The ionic environment determines ribozyme cleavage rate by modulation of nucleobase pK a

Several small ribozymes employ general acid-base catalysis as a mechanism to enhance site-specific RNA cleavage, even though the functional groups on the ribonucleoside building blocks of RNA have pK (a) values far removed from physiological pH. The rate of the cleavage reaction is strongly affected by the identity of the metal cation present in the reaction solution; however, the mechanism(s) by which different cations contribute to rate enhancement has not been determined. Using the Neurospora VS ribozyme, we provide evidence that different cations confer particular shifts in the apparent pK (a) values of the catalytic nucleobases, which in turn determines the fraction of RNA in the protonation state competent for general acid-base catalysis at a given pH, which determines the observed rate of the cleavage reaction. Despite large differences in observed rates of cleavage in different cations, mathematical models of general acid-base catalysis indicate that k (1), the intrinsic rate of the bond-breaking step, is essentially constant irrespective of the identity of the cation(s) in the reaction solution. Thus, in contrast to models that invoke unique roles for metal ions in ribozyme chemical mechanisms, we find that most, and possibly all, of the ion-specific rate enhancement in the VS ribozyme can be explained solely by the effect of the ions on nucleobase pK (a). The inference that k (1) is essentially constant suggests a resolution of the problem of kinetic ambiguity in favor of a model in which the lower pK (a) is that of the general acid and the higher pK (a) is that of the general base.

Structural effects of nucleobase variations at key active site residue Ade38 in the hairpin ribozyme

The hairpin ribozyme requires functional groups from Ade38 to achieve efficient bond cleavage or ligation. To identify molecular features that contribute to catalysis, structures of position 38 base variants 2,6-diaminopurine (DAP), 2-aminopurine (AP), cytosine (Cyt), and guanine (Gua) were determined between 2.2 and 2.8 A resolution. For each variant, two substrate modifications were compared: (1) a 2'-O-methyl-substituent at Ade-1 was used in lieu of the nucleophile to mimic the precatalytic state, and (2) a 3'-deoxy-2',5'-phosphodiester linkage between Ade-1 and Gua+1 was used to mimic a reaction-intermediate conformation. While the global fold of each variant remained intact, the results revealed the importance of Ade38 N1 and N6 groups. Absence of N6 resulting from AP38 coincided with failure to localize the precatalytic scissile phosphate. Cyt38 severely impaired catalysis in a prior study, and its structures here indicated an anti base conformation that sequesters the imino moiety from the scissile bond. Gua38 was shown to be even more deleterious to activity. Although the precatalytic structure was nominally affected, the reaction-intermediate conformation indicated a severe electrostatic clash between the Gua38 keto oxygen and the pro-Rp oxygen of the scissile bond. Overall, position 38 modifications solved in the presence of 2'-OMe Ade-1 deviated from in-line geometry, whereas variants with a 2',5' linkage exhibited S-turn destabilization, as well as base conformational changes from syn to anti. These findings demonstrate the importance of the Ade38 Watson-Crick face in attaining a reaction-intermediate state and the sensitivity of the RNA fold to restructuring when electrostatic and shape features fail to complement.

Catalytic strategies of self-cleaving ribozymes

[Structure: see text]. Five naturally occurring nucleolytic ribozymes have been identified: the hammerhead, hairpin, glmS, hepatitis delta virus (HDV), and Varkud satellite (VS) ribozymes. All of these RNA enzymes catalyze self-scission of the RNA backbone using a chemical mechanism equivalent to that of RNase A. RNase A uses four basic strategies to promote this reaction: geometric constraints, activation of the nucleophile, transition-state stabilization, and leaving group protonation. In this Account, we discuss the current thinking on how nucleolytic ribozymes harness RNase A's four sources of catalytic power. The geometry of the phosphodiester cleavage reaction constrains the nucleotides flanking the scissile phosphate so that they are unstacked from a canonical A-form helix and thus require alternative stabilization. Crystal structures and mutational analysis reveal that cross-strand base pairing, along with unconventional stacking and tertiary hydrogen-bonding interactions, work to stabilize the splayed conformation in nucleolytic ribozymes. Deprotonation of the 2'-OH nucleophile greatly increases its nucleophilicity in the strand scission reaction. Crystal structures of the hammerhead, hairpin, and glmS ribozymes reveal the N1 of a G residue within hydrogen-bonding distance of the 2'-OH. In each case, this residue has also been shown to be important for catalysis. In the HDV ribozyme, a hydrated magnesium has been implicated as the general base. Catalysis by the VS ribozyme requires both an A and a G, but the precise role of either has not been elucidated. Enzymes can lower the energy of a chemical reaction by binding more tightly to the transition state than to the ground states. Comparison of the hairpin ground- and transition-state mimic structures reveal greater hydrogen bonding to the transition-state mimic structure, suggesting transition-state stabilization as a possible catalytic strategy. However, the hydrogen-bonding pattern in the glmS ribozyme transition-state mimic structure and the ground-state structures are equivalent. Protonation of the 5'-O leaving group by a variety of functional groups can promote the cleavage reaction. In the HDV ribozyme, the general acid is a conserved C residue. In the hairpin ribozyme, a G residue has been implicated in protonation of the leaving group. An A in the hammerhead ribozyme probably plays a similar role. In the glmS ribozyme, an exogenous cofactor may provide the general acid. This diversity is in contrast to the relatively small number of functional groups that serve as a general base, where at least three of the nucleolytic ribozymes may use the N1 of a G.

An important role of G638 in the cis-cleavage reaction of the Neurospora VS ribozyme revealed by a novel nucleotide analog incorporation method

We describe a chemical coupling procedure that allows joining of two RNAs, one of which contains a site-specific base analog substitution, in the absence of divalent ions. This method allows incorporation of nucleotide analogs at specific positions even into large, cis-cleaving ribozymes. Using this method we have studied the effects of substitution of G638 in the cleavage site loop of the VS ribozyme with a variety of purine analogs having different functional groups and pK(a) values. Cleavage rate versus pH profiles combined with kinetic solvent isotope experiments indicate an important role for G638 in proton transfer during the rate-limiting step of the cis-cleavage reaction.

The Cech Symposium: a celebration of 25 years of ribozymes, 10 years of TERT, and 60 years of Tom

The Cech Symposium was held in Boulder, Colorado, on July 12-13, 2007, to celebrate a triple anniversary: 25 years since the first publication reporting RNA self-splicing, 10 years since the identification of reverse transcriptase motifs in the catalytic subunit of telomerase, and 60 years since the birth of Thomas R. Cech. Past and present members of the Cech laboratory presented on their current research, which branched into many categories of study including RNA-mediated catalysis, telomerase and telomeres, new frontiers in nucleic acids, alternative splicing, as well as scientific research with direct medical applications.

Ribozyme catalysis revisited: is water involved?

Enzymatic catalysis by RNA was discovered 25 years ago, yet mechanistic insights are emerging only slowly. Thought to be metalloenzymes at first, some ribozymes proved more versatile than anticipated when shown to utilize their own functional groups for catalysis. Recent evidence suggests that some may also judiciously place structural water molecules to shuttle protons in acid-base catalyzed reactions.

Wild-type is the optimal sequence of the HDV ribozyme under cotranscriptional conditions

RNA viruses are responsible for a variety of human diseases, and the pathogenicity of RNA viruses is often attributed to a high rate of mutation. Self-cleavage activity of the wild-type hepatitis delta virus (HDV) ribozyme as measured in standard divalent ion renaturation assays is biphasic and mostly slow and can be improved by multiple rational changes to ribozyme sequence or by addition of chemical denaturants. This is unusual in the sense that wild type is the most catalytically active sequence for the majority of protein enzymes, and RNA viruses are highly mutable. To see whether the ribozyme takes advantage of fast-reacting sequence changes in vivo, we performed alignment of 76 genomic and 269 antigenomic HDV isolates. Paradoxically, the sequence for the ribozyme was found to be essentially invariant in nature. We therefore tested whether three ribozyme sequence changes that improve self-cleavage under standard divalent ion renaturation assays also improve self-cleavage during transcription. Remarkably, wild type was as fast, or faster, than these mutants under cotranscriptional conditions. Slowing the rate of transcription or adding the hepatitis delta antigen protein only further stimulated cotranscriptional self-cleavage activity. Thus, the relative activity of HDV ribozyme mutants depends critically on whether the reaction is assayed under in vivo-like conditions. A model is presented for how wild-type ribozyme sequence and flanking sequence work in concert to promote efficient self-cleavage during transcription. Wild type being the optimal ribozyme sequence under in vivo-like conditions parallels the behavior of most protein enzymes.

RNA catalysis: ribozymes, ribosomes, and riboswitches

The catalytic mechanisms employed by RNA are chemically more diverse than initially suspected. Divalent metal ions, nucleobases, ribosyl hydroxyl groups, and even functional groups on metabolic cofactors all contribute to the various strategies employed by RNA enzymes. This catalytic breadth raises intriguing evolutionary questions about how RNA lost its biological role in some cases, but not in others, and what catalytic roles RNA might still be playing in biology.

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.

Impact of an extruded nucleotide on cleavage activity and dynamic catalytic core conformation of the hepatitis delta virus ribozyme

The self-cleaving hepatitis delta virus (HDV) ribozyme is essential for the replication of HDV, a liver disease causing pathogen in humans. The catalytically critical nucleotide C75 of the ribozyme is buttressed by a trefoil turn pivoting around an extruded G76. In all available crystal structures, the conformation of G76 is restricted by stacking with G76 of a neighboring molecule. To test whether this crystal contact introduces a structural perturbation into the catalytic core, we have analyzed approximately 200 ns of molecular dynamics (MD) simulations. In the absence of crystal packing, the simulated G76 fluctuates between several conformations, including one wherein G76 establishes a perpendicular base quadruplet in the major groove of the adjacent P1 stem. Second-site mutagenesis experiments suggest that the identity of the nucleotide in position 76 (N76) indeed contributes to the catalytic activity of a trans-acting HDV ribozyme through its capacity for hydrogen bonding with P1. By contrast, in the cis-cleaving genomic ribozyme the functional relevance of N76 is less pronounced and not correlated with the P1 sequence. Terbium(III) footprinting and additional MD show that the activity differences between N76 mutants of this ribozyme are related instead to changes in average conformation and modified cross-correlations in the trefoil turn.

A chemo-genetic approach for the study of nucleobase participation in nucleolytic ribozymes

A novel chemo-genetic approach for the analysis of general acid-base catalysis by nucleobases in ribozymes is reviewed. This involves substitution of a C-nucleoside with imidazole in place of a natural nucleobase. The Varkud satellite ribozyme in which the nucleobase at the critical 756 position has been replaced by imidazole is active in both cleavage and ligation reactions. Similarly, a modified hairpin ribozyme with the nucleobase at position 8 substituted by imidazole is active in cleavage and ligation reactions. Although the rates are lower than those of the natural ribozymes, they are significantly greater than other variants at these positions. The dependence of the hairpin ribozyme reaction rates on pH has been studied. Both cleavage and ligation reactions display a bell-shaped pH dependence, consistent with general acid-base catalysis involving the nucleotide at position 8.

Reactions of phosphate and phosphorothiolate diesters with nucleophiles: comparison of transition state structures

A series of methyl aryl phosphorothiolate esters (SP) were synthesized and their reactions with pyridine derivatives were compared to those for methyl aryl phosphate esters (OP). Results show that SP esters react with pyridine nucleophiles via a concerted S(N)2(P) mechanism. Brønsted analysis suggests that reactions of both SP and OP esters proceed via transition states with dissociative character. The overall similarity of the transition state structures supports the use of phosphorothiolates as substrate analogues to probe mechanisms of enzyme-catalyzed phosphoryl transfer reactions.

Distinct reaction pathway promoted by non-divalent-metal cations in a tertiary stabilized hammerhead ribozyme

Divalent ion sensitivity of hammerhead ribozymes is significantly reduced when the RNA structure includes appropriate tertiary stabilization. Therefore, we investigated the activity of the tertiary stabilized "RzB" hammerhead ribozyme in several nondivalent ions. Ribozyme RzB is active in spermidine and Na(+) alone, although the cleavage rates are reduced by more than 1,000-fold relative to the rates observed in Mg(2+) and in transition metal ions. The trivalent cobalt hexammine (CoHex) ion is often used as an exchange-inert analog of hydrated magnesium ion. Trans-cleavage rates exceeded 8 min(-1) in 20 mM CoHex, which promoted cleavage through outersphere interactions. The stimulation of catalysis afforded by the tertiary structural interactions within RzB does not require Mg(2+), unlike other extended hammerhead ribozymes. Site-specific interaction with at least one Mg(2+) ion is suggested by CoHex competition experiments. In the presence of a constant, low concentration of Mg(2+), low concentrations of CoHex decreased the rate by two to three orders of magnitude relative to the rate in Mg(2+) alone. Cleavage rates increased as CoHex concentrations were raised further, but the final fraction cleaved was lower than what was observed in CoHex or Mg(2+) alone. These observations suggest that Mg(2+) and CoHex compete for binding and that they cause misfolded structures when they are together. The results of this study support the existence of an alternate catalytic mechanism used by nondivalent ions (especially CoHex) that is distinct from the one promoted by divalent metal ions, and they imply that divalent metals influence catalysis through a specific nonstructural role.

A guanine nucleobase important for catalysis by the VS ribozyme

A guanine (G638) within the substrate loop of the VS ribozyme plays a critical role in the cleavage reaction. Replacement by any other nucleotide results in severe impairment of cleavage, yet folding of the substrate is not perturbed, and the variant substrates bind the ribozyme with similar affinity, acting as competitive inhibitors. Functional group substitution shows that the imino proton on the N1 is critical, suggesting a possible role in general acid-base catalysis, and this in accord with the pH dependence of the reaction rate for the natural and modified substrates. We propose a chemical mechanism for the ribozyme that involves general acid-base catalysis by the combination of the nucleobases of guanine 638 and adenine 756. This is closely similar to the probable mechanism of the hairpin ribozyme, and the active site arrangements for the two ribozymes appear topologically equivalent. This has probably arisen by convergent evolution.

Asymmetric catalysis by chiral hydrogen-bond donors

Hydrogen bonding is responsible for the structure of much of the world around us. The unusual and complex properties of bulk water, the ability of proteins to fold into stable three-dimensional structures, the fidelity of DNA base pairing, and the binding of ligands to receptors are among the manifestations of this ubiquitous noncovalent interaction. In addition to its primacy as a structural determinant, hydrogen bonding plays a crucial functional role in catalysis. Hydrogen bonding to an electrophile serves to decrease the electron density of this species, activating it toward nucleophilic attack. This principle is employed frequently by Nature's catalysts, enzymes, for the acceleration of a wide range of chemical processes. Recently, organic chemists have begun to appreciate the tremendous potential offered by hydrogen bonding as a mechanism for electrophile activation in small-molecule, synthetic catalyst systems. In particular, chiral hydrogen-bond donors have emerged as a broadly applicable class of catalysts for enantioselective synthesis. This review documents these advances, emphasizing the structural and mechanistic features that contribute to high enantioselectivity in hydrogen-bond-mediated catalytic processes.

Evidence for proton transfer in the rate-limiting step of a fast-cleaving Varkud satellite ribozyme

A fast-cleaving version of the Varkud satellite ribozyme, called RG, shows an apparent cis-cleavage rate constant of 5 sec(-1), similar to the rates of protein enzymes that catalyze similar reactions. Here, we describe mutational, pH-rate, and kinetic solvent isotope experiments that investigate the identity and rate constant of the rate-limiting step in this reaction. Self-cleavage of RG exhibits a bell-shaped rate vs. pH profile with apparent pK(a)s of 5.8 and 8.3, consistent with the protonation state of two nucleotides being important for the rate of cleavage. Cleavage experiments in heavy water (D(2)O) revealed a kinetic solvent isotope effect consistent with proton transfer in the rate-limiting step. A mutant RNA that disrupts a peripheral loop-loop interaction involved in RNA folding exhibits pH- and D(2)O-independent cleavage approximately 10(3)-fold slower than wild type, suggesting that this mutant is limited by a different step than wild type. Substitution of adenosine 756 in the putative active-site loop with cytosine also decreases the cleavage rate approximately 10(3)-fold, but the A756C mutant retains pH- and D(2)O-sensitivity similar to wild type, consistent with this mutant and wild type being limited by the chemical step of the reaction. These results suggest that the RG ribozyme provides a good experimental system to investigate the nature of fast, rate-limiting steps in a ribozyme cleavage reaction.

The genomic HDV ribozyme utilizes a previously unnoticed U-turn motif to accomplish fast site-specific catalysis

The genome of the human hepatitis delta virus (HDV) harbors a self-cleaving catalytic RNA motif, the genomic HDV ribozyme, whose crystal structure shows the dangling nucleotides 5′ of the cleavage site projecting away from the catalytic core. This 5′-sequence contains a clinically conserved U − 1 that we find to be essential for fast cleavage, as the order of activity follows U − 1 > C − 1 > A − 1 > G − 1, with a >25-fold activity loss from U − 1 to G − 1. Terbium(III) footprinting detects conformations for the P1.1 stem, the cleavage site wobble pair and the A-minor motif of the catalytic trefoil turn that depend on the identity of the N − 1 base. The most tightly folded catalytic core, resembling that of the reaction product, is found in the U − 1 wild-type precursor. Molecular dynamics simulations demonstrate that a U − 1 forms the most robust kink around the scissile phosphate, exposing it to the catalytic C75 in a previously unnoticed U-turn motif found also, for example, in the hammerhead ribozyme and tRNAs. Strikingly, we find that the common structural U-turn motif serves distinct functions in the HDV and hammerhead ribozymes.

Structural Roles of Monovalent Cations in the HDV Ribozyme

The hepatitis delta virus (HDV) ribozyme catalyzes viral RNA self-cleavage through general acid-base chemistry in which an active-site cytidine and at least one metal ion are involved. Monovalent metal ions support slow catalysis and were proposed to substitute for structural, but not catalytic, divalent metal ions in the RNA. To investigate the role of monovalent cations in ribozyme structure and function, we determined the crystal structure of the precursor HDV ribozyme in the presence of thallium ions (Tl(+)). Two Tl(+) ions can occupy a previously observed divalent metal ion hexahydrate-binding site located near the scissile phosphate, but are easily competed away by cobalt hexammine, a magnesium hexahydrate mimic and potent reaction inhibitor. Intriguingly, a third Tl(+) ion forms direct inner-sphere contacts with the ribose 2'-OH nucleophile and the pro-S(p) scissile phosphate oxygen. We discuss possible structural and catalytic implications of monovalent cation binding for the HDV ribozyme mechanism.

Molecular dynamics simulations of RNA: an in silico single molecule approach

RNA molecules are now known to be involved in the processing of genetic information at all levels, taking on a wide variety of central roles in the cell. Understanding how RNA molecules carry out their biological functions will require an understanding of structure and dynamics at the atomistic level, which can be significantly improved by combining computational simulation with experiment. This review provides a critical survey of the state of molecular dynamics (MD) simulations of RNA, including a discussion of important current limitations of the technique and examples of its successful application. Several types of simulations are discussed in detail, including those of structured RNA molecules and their interactions with the surrounding solvent and ions, catalytic RNAs, and RNA-small molecule and RNA-protein complexes. Increased cooperation between theorists and experimentalists will allow expanded judicious use of MD simulations to complement conceptually related single molecule experiments. Such cooperation will open the door to a fundamental understanding of the structure-function relationships in diverse and complex RNA molecules. .

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.

Long-range impact of peripheral joining elements on structure and function of the hepatitis delta virus ribozyme

The HDV ribozyme is an RNA enzyme from the human pathogenic hepatitis delta virus (HDV) that has recently also been identified in the human genome. It folds into a compact, nested double-pseudoknot. We examined here the functional relevance of the capping loop L4 and the helical crossover J1/2, which tightly interlace the two helical stacks of the ribozyme. Peripheral structural elements such as these are present in cis-acting, but not trans-acting ribozymes, which may explain the order-of-magnitude decrease in cleavage activity observed in trans-acting ribozymes with promise in gene therapy applications. Comparison of a systematic set of cis- and trans-acting HDV ribozymes shows that the absence of either L4 or J1/2 significantly and independently impacts catalytic activity. Using terbium(III) footprinting and affinity studies, as well as distance measurements based on time-resolved fluorescence resonance energy transfer, we find that J1/2 is most important for conferring structural properties similar to those of the cis-acting ribozyme. Our results are consistent with a model in which removal of either a helical crossover or surprisingly a capping loop induces greater dynamics and expansion of the catalytic core at long range, impacting local and global folding, as well as catalytic function.

RNA folding and the origins of catalytic activity in the hairpin ribozyme

The nucleolytic ribozymes catalyse site-specific phosphodiester cleavage and ligation transesterification reactions in RNA. The hairpin ribozyme folds to generate an intimate loop-loop interaction to create the local environment in which catalysis can proceed. We have studied the ion-induced folding using single-molecule FRET experiments, showing that the four-way helical junction accelerates the folding 500-fold by introducing a discrete intermediate that juxtaposes the loops. Using FRET we can observe individual hairpin ribozyme molecules as they undergo multiple cycles of cleavage and ligation, and measure the rates of the internal reactions, free of uncertainties in the contributions of docking and substrate dissociation processes. On average, the cleaved ribozyme undergoes several docking-undocking events before a ligation reaction occurs. On the basis of these experiments, we have explored the role of the nucleobases G8 and A38 in the catalysis. Both cleavage and ligation reactions are pH dependent, corresponding to the titration of a group with pKa=6.2. We have used a novel ribonucleoside in which these bases are replaced by imidazole to investigate the role of acid-base catalysis in this ribozyme. We observe significant rates of cleavage and ligation, and a bell-shaped pH dependence for both.

Electron Hole Flow Patterns through the RNA-Cleaving 8-17 Deoxyribozyme Yield Unusual Information about Its Structure and Folding

DNA double helices have been shown to conduct electron holes over significant distances. Here, we report on the hole flow patterns within a more intricately folded DNA complex, the 8-17 deoxyribozyme bound to a DNA pseudosubstrate, incorporating three helical elements and two catalytically relevant loops. The observed hole flow patterns within the complex permitted a quantitative assessment of the stacking preferences of the three constituent helices and provided evidence for significant transitions within the complex's global geometry. The patterns further suggested varying levels of solvent exposure of the complex's constituent parts, and revealed that a catalytically relevant cytosine within the folded complex exists in an unusual structural/electronic environment. Our data suggest that the study of charge flow may provide novel perspectives on the structure and folding of intricately folded DNAs and RNAs.

Density Functional Theory Investigation on the Mechanism of the Hepatitis Delta Virus Ribozyme

Density functional theory methods have been used to investigate the hepatitis delta virus (HDV) ribozyme and its catalyzed phosphodiester cleavage. In particular, the effects of the environment's polarity and/or specific hydrogen-bond interactions on the proton affinity of the active site cytosine's N3 ring center have been considered. In addition, the basicities of possible hydrated Mg2+ ion species were also examined. The mechanism previously proposed for the HDV ribozyme in which the active site cytosine (C75) is protonated and thus acts as an acid while the Mg2+ species acts as the complementary base was then investigated. The possible role of tautomerization of C75 is also discussed.

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.

Structural basis of glmS ribozyme activation by glucosamine-6-phosphate

The glmS ribozyme is the only natural catalytic RNA known to require a small-molecule activator for catalysis. This catalytic RNA functions as a riboswitch, with activator-dependent RNA cleavage regulating glmS messenger RNA expression. We report crystal structures of the glmS ribozyme in precleavage states that are unliganded or bound to the competitive inhibitor glucose-6-phosphate and in the postcleavage state. All structures superimpose closely, revealing a remarkably rigid RNA that contains a preformed active and coenzyme-binding site. Unlike other riboswitches, the glmS ribozyme binds its activator in an open, solvent-accessible pocket. Our structures suggest that the amine group of the glmS ribozyme-bound coenzyme performs general acid-base and electrostatic catalysis.

A genomewide search for ribozymes reveals an HDV-like sequence in the human CPEB3 gene

Ribozymes are thought to have played a pivotal role in the early evolution of life, but relatively few have been identified in modern organisms. We performed an in vitro selection aimed at isolating self-cleaving RNAs from the human genome. The selection yielded several ribozymes, one of which is a conserved mammalian sequence that resides in an intron of the CPEB3 gene, which belongs to a family of genes regulating messenger RNA polyadenylation. The CPEB3 ribozyme is structurally and biochemically related to the human hepatitis delta virus (HDV) ribozymes. The occurrence of this ribozyme exclusively in mammals suggests that it may have evolved as recently as 200 million years ago. We postulate that HDV arose from the human transcriptome.