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

Structure and sequence at an RNA template 5' end influence insertion of transgenes by an R2 retrotransposon protein

R2 non-long terminal repeat retrotransposons insert site-specifically into ribosomal RNA genes (rDNA) in a broad range of multicellular eukaryotes. R2-encoded proteins can be leveraged to mediate transgene insertion at 28S rDNA loci in cultured human cells. This strategy, precise RNA-mediated insertion of transgenes (PRINT), relies on the codelivery of an mRNA encoding R2 protein and an RNA template encoding a transgene cassette of choice. Here, we demonstrate that the PRINT RNA template 5' module, which as a complementary DNA 3' end will generate the transgene 5' junction with rDNA, influences the efficiency and mechanism of gene insertion. Iterative design and testing identified optimal 5' modules consisting of a hepatitis delta virus-like ribozyme fold with high thermodynamic stability, suggesting that RNA template degradation from its 5' end may limit transgene insertion efficiency. We also demonstrate that transgene 5' junction formation can be either precise, formed by annealing the 3' end of first-strand complementary DNA with the upstream target site, or imprecise, by end-joining, but this difference in junction formation mechanism is not a major determinant of insertion efficiency. Sequence characterization of imprecise end-joining events indicates surprisingly minimal reliance on microhomology. Our findings expand the current understanding of the role of R2 retrotransposon transcript sequence and structure, and especially the 5' ribozyme fold, for retrotransposon mobility and RNA-templated gene synthesis in cells.

Experimental Resurrection of Ancestral Mammalian CPEB3 Ribozymes Reveals Deep Functional Conservation

Self-cleaving ribozymes are genetic elements found in all domains of life, but their evolution remains poorly understood. A ribozyme located in the second intron of the cytoplasmic polyadenylation binding protein 3 gene (CPEB3) shows high sequence conservation in mammals, but little is known about the functional conservation of self-cleaving ribozyme activity across the mammalian tree of life or during the course of mammalian evolution. Here, we use a phylogenetic approach to design a mutational library and a deep sequencing assay to evaluate the in vitro self-cleavage activity of numerous extant and resurrected CPEB3 ribozymes that span over 100 My of mammalian evolution. We found that the predicted sequence at the divergence of placentals and marsupials is highly active, and this activity has been conserved in most lineages. A reduction in ribozyme activity appears to have occurred multiple different times throughout the mammalian tree of life. The in vitro activity data allow an evaluation of the predicted mutational pathways leading to extant ribozyme as well as the mutational landscape surrounding these ribozymes. The results demonstrate that in addition to sequence conservation, the self-cleavage activity of the CPEB3 ribozyme has persisted over millions of years of mammalian evolution.

Topological constraints of structural elements in regulation of catalytic activity in HDV-like self-cleaving ribozymes

Self-cleaving ribozymes fold into intricate structures, which orient active site groups into catalytically competent conformations. Most ribozyme families have distinct catalytic cores stabilized by tertiary interactions between domains peripheral to those cores. We show that large hepatitis delta virus (HDV)-like ribozymes are activated by peripheral domains that bring two helical segments, P1 and P2, into proximity – a “pinch” that results in rate acceleration by almost three orders of magnitude. Kinetic analysis of ribozymes with systematically altered length and stability of the peripheral domain revealed that about one third of its free energy of formation is used to lower an activation energy barrier, likely related to a rate-limiting conformational change leading to the pre-catalytic state. These findings provide a quantitative view of enzyme regulation by peripheral domains and may shed light on the energetics of allosteric regulation.

A Two-Metal-Ion-Mediated Conformational Switching Pathway for HDV Ribozyme Activation

RNA enzymes serve as a potentially powerful platform from which to design catalysts and engineer new biotechnology. A fundamental understanding of these systems provides insight to guide design. The hepatitis delta virus ribozyme (HDVr) is a small, self-cleaving RNA motif widely distributed in nature, that has served as a paradigm for understanding basic principles of RNA catalysis. Nevertheless, questions remain regarding the precise roles of divalent metal ions and key nucleotides in catalysis. In an effort to establish a reaction mechanism model consistent with available experimental data, we utilize molecular dynamics simulations to explore different conformations and metal ion binding modes along the HDVr reaction path. Building upon recent crystallographic data, our results provide a dynamic model of the HDVr reaction mechanism involving a conformational switch between multiple non-canonical G25:U20 base pair conformations in the active site. These local nucleobase dynamics play an important role in catalysis by modulating the metal binding environments of two Mg²⁺ ions that support catalysis at different steps of the reaction pathway. The first ion plays a structural role by inducing a base pair flip necessary to obtain the catalytic fold in which C75 moves towards to the scissile phosphate in the active site. Ejection of this ion then permits a second ion to bind elsewhere in the active site and facilitate nucleophile activation. The simulations collectively describe a mechanistic scenario that is consistent with currently available experimental data from crystallography, phosphorothioate substitutions, and chemical probing studies. Avenues for further experimental verification are suggested.

Wobble pairs of the HDV ribozyme play specific roles in stabilization of active site dynamics

The hepatitis delta virus (HDV) is the only known human pathogen whose genome contains a catalytic RNA motif (ribozyme). The overall architecture of the HDV ribozyme is that of a double-nested pseudoknot, with two GU pairs flanking the active site. Although extensive studies have shown that mutation of either wobble results in decreased catalytic activity, little work has focused on linking these mutations to specific structural effects on catalytic fitness. Here we use molecular dynamics simulations based on an activated structure to probe the active site dynamics as a result of wobble pair mutations. In both wild-type and mutant ribozymes, the in-line fitness of the active site (as a measure of catalytic proficiency) strongly depends on the presence of a C75(N3H3+)N1(O5') hydrogen bond, which positions C75 as the general acid for the reaction. Our mutational analyses show that each GU wobble supports catalytically fit conformations in distinct ways; the reverse G25U20 wobble promotes high in-line fitness, high occupancy of the C75(N3H3+)G1(O5') general-acid hydrogen bond and stabilization of the G1U37 wobble, while the G1U37 wobble acts more locally by stabilizing high in-line fitness and the C75(N3H3+)G1(O5') hydrogen bond. We also find that stable type I A-minor and P1.1 hydrogen bonding above and below the active site, respectively, prevent local structural disorder from spreading and disrupting global conformation. Taken together, our results define specific, often redundant architectural roles for several structural motifs of the HDV ribozyme active site, expanding the known roles of these motifs within all HDV-like ribozymes and other structured RNAs.

Two Divalent Metal Ions and Conformational Changes Play Roles in the Hammerhead Ribozyme Cleavage Reaction

The hammerhead ribozyme is a self-cleaving RNA broadly dispersed across all kingdoms of life. Although it was the first of the small, nucleolytic ribozymes discovered, the mechanism by which it catalyzes its reaction remains elusive. The nucleobase of G12 is well positioned to be a general base, but it is unclear if or how this guanine base becomes activated for proton transfer. Metal ions have been implicated in the chemical mechanism, but no interactions between divalent metal ions and the cleavage site have been observed crystallographically. To better understand how this ribozyme functions, we have solved crystal structures of wild-type and G12A mutant ribozymes. We observe a pH-dependent conformational change centered around G12, consistent with this nucleotide becoming deprotonated. Crystallographic and kinetic analysis of the G12A mutant reveals a Zn(2+) specificity switch suggesting a direct interaction between a divalent metal ion and the purine at position 12. The metal ion specificity switch and the pH-rate profile of the G12A mutant suggest that the minor imino tautomer of A12 serves as the general base in the mutant ribozyme. We propose a model in which the hammerhead ribozyme rearranges prior to the cleavage reaction, positioning two divalent metal ions in the process. The first metal ion, positioned near G12, becomes directly coordinated to the O6 keto oxygen, to lower the pKa of the general base and organize the active site. The second metal ion, positioned near G10.1, bridges the N7 of G10.1 and the scissile phosphate and may participate directly in the cleavage reaction.

Determination of hepatitis delta virus ribozyme N(-1) nucleobase and functional group specificity using internal competition kinetics

Biological catalysis involves interactions distant from the site of chemistry that can position the substrate for reaction. Catalysis of RNA 2'-O-transphosphorylation by the hepatitis delta virus (HDV) ribozyme is sensitive to the identity of the N(-1) nucleotide flanking the reactive phosphoryl group. However, the interactions that affect the conformation of this position, and in turn the 2'O nucleophile, are unclear. Here, we describe the application of multiple substrate internal competition kinetic analyses to understand how the N(-1) nucleobase contributes to HDV catalysis and test the utility of this approach for RNA structure-function studies. Internal competition reactions containing all four substrate sequence variants at the N(-1) position in reactions using ribozyme active site mutations at A77 and A78 were used to test a proposed base-pairing interaction. Mutants A78U, A78G, and A79G retain significant catalytic activity but do not alter the specificity for the N(-1) nucleobase. Effects of nucleobase analog substitutions at N(-1) indicate that U is preferred due to the ability to donate an H-bond in the Watson-Crick face and avoid minor groove steric clash. The results provide information essential for evaluating models of the HDV active site and illustrate multiple substrate kinetic analyses as a practical approach for characterizing structure-function relationships in RNA reactions.

Identification of minimal HDV-like ribozymes with unique divalent metal ion dependence in the human microbiome

HDV-like self-cleaving ribozymes have been found in a wide variety of organisms, implicated in diverse biological processes, and their activity typically shows a strong divalent metal dependence, but little metal specificity. Recent studies suggested that very short variants of these ribozymes exist in nature, but their distribution and biochemical properties have not been established. To map out the distribution of small HDV-like ribozymes, the drz-Spur-3 sequence was minimized to yield a core construct for structure-based bioinformatic searches. These searches revealed several microbial ribozymes, particularly in the human microbiome. Kinetic profile of the smallest ribozyme revealed two distinct metal binding sites, only one of which promotes fast catalysis. Furthermore, this ribozyme showed markedly reduced activity in Ca(2+), even in the presence of physiological Mg(2+) concentrations. Our study substantially expands the number of microbial HDV-like ribozymes and provides an example of cleavage regulation by divalent metals.

Thio effects and an unconventional metal ion rescue in the genomic hepatitis delta virus ribozyme

Metal ion and nucleobase catalysis are important for ribozyme mechanism, but the extent to which they cooperate is unclear. A crystal structure of the hepatitis delta virus (HDV) ribozyme suggested that the pro-RP oxygen at the scissile phosphate directly coordinates a catalytic Mg(2+) ion and is within hydrogen bonding distance of the amine of the general acid C75. Prior studies of the genomic HDV ribozyme, however, showed neither a thio effect nor metal ion rescue using Mn(2+). Here, we combine experiment and theory to explore phosphorothioate substitutions at the scissile phosphate. We report significant thio effects at the scissile phosphate and metal ion rescue with Cd(2+). Reaction profiles with an SP-phosphorothioate substitution are indistinguishable from those of the unmodified substrate in the presence of Mg(2+) or Cd(2+), supporting the idea that the pro-SP oxygen does not coordinate metal ions. The RP-phosphorothioate substitution, however, exhibits biphasic kinetics, with the fast-reacting phase displaying a thio effect of up to 5-fold and the slow-reacting phase displaying a thio effect of ~1000-fold. Moreover, the fast- and slow-reacting phases give metal ion rescues in Cd(2+) of up to 10- and 330-fold, respectively. The metal ion rescues are unconventional in that they arise from Cd(2+) inhibiting the oxo substrate but not the RP substrate. This metal ion rescue suggests a direct interaction of the catalytic metal ion with the pro-RP oxygen, in line with experiments with the antigenomic HDV ribozyme. Experiments without divalent ions, with a double mutant that interferes with Mg(2+) binding, or with C75 deleted suggest that the pro-RP oxygen plays at most a redundant role in positioning C75. Quantum mechanical/molecular mechanical (QM/MM) studies indicate that the metal ion contributes to catalysis by interacting with both the pro-RP oxygen and the nucleophilic 2'-hydroxyl, supporting the experimental findings.

A Role for the 2' OH of peptidyl-tRNA substrate in peptide release on the ribosome revealed through RF-mediated rescue

The 2' OH of the peptidyl-tRNA substrate is thought to be important for catalysis of both peptide bond formation and peptide release in the ribosomal active site. The release reaction also specifically depends on a release factor protein (RF) to hydrolyze the ester linkage of the peptidyl-tRNA upon recognition of stop codons in the A site. Here, we demonstrate that certain amino acid substitutions (in particular those containing hydroxyl or thiol groups) in the conserved GGQ glutamine of release factor RF1 can rescue defects in the release reaction associated with peptidyl-tRNA substrates lacking a 2' OH. We explored this rescue effect through biochemical and computational approaches that support a model where the 2' OH of the P-site substrate is critical for orienting the nucleophile in a hydrogen-bonding network productive for catalysis.

Two distinct catalytic strategies in the hepatitis δ virus ribozyme cleavage reaction

The hepatitis delta virus (HDV) ribozyme and related RNAs are widely dispersed in nature. This RNA is a small nucleolytic ribozyme that self-cleaves to generate products with a 2',3'-cyclic phosphate and a free 5'-hydroxyl. Although small ribozymes are dependent on divalent metal ions under biologically relevant buffer conditions, they function in the absence of divalent metal ions at high ionic strengths. This characteristic suggests that a functional group within the covalent structure of small ribozymes is facilitating catalysis. Structural and mechanistic analyses have demonstrated that the HDV ribozyme active site contains a cytosine with a perturbed pK(a) that serves as a general acid to protonate the leaving group. The reaction of the HDV ribozyme in monovalent cations alone never approaches the velocity of the Mg(2+)-dependent reaction, and there is significant biochemical evidence that a Mg(2+) ion participates directly in catalysis. A recent crystal structure of the HDV ribozyme revealed that there is a metal binding pocket in the HDV ribozyme active site. Modeling of the cleavage site into the structure suggested that this metal ion can interact directly with the scissile phosphate and the nucleophile. In this manner, the Mg(2+) ion can serve as a Lewis acid, facilitating deprotonation of the nucleophile and stabilizing the conformation of the cleavage site for in-line attack of the nucleophile at the scissile phosphate. This catalytic strategy had previously been observed only in much larger ribozymes. Thus, in contrast to most large and small ribozymes, the HDV ribozyme uses two distinct catalytic strategies in its cleavage reaction.

Mechanistic strategies in the HDV ribozyme: chelated and diffuse metal ion interactions and active site protonation

The crystal structure of the precleaved form of the hepatitis delta virus (HDV) ribozyme reveals two G•U wobbles near the active site: a rare reverse G•U wobble involving a syn G base, and a standard G•U wobble at the cleavage site. The catalytic mechanism for this ribozyme has been proposed to involve a Mg(2+) ion bound to the reverse G•U wobble, as well as a protonated C75 base. We carried out molecular dynamics simulations to analyze metal ion interaction with the reverse and standard G•U wobbles and to investigate the impact of C75 protonation on the structure and motions of the ribozyme. We identified two types of Mg(2+) ions associated with the ribozyme, chelated and diffuse, at the reverse and standard G•U wobbles, respectively, which appear to contribute to catalysis and stability, respectively. These two metal ion sites exhibit relatively independent behavior. Protonation of C75 was observed to locally organize the active site in a manner that facilitates the catalytic mechanism, in which C75(+) acts as a general acid and Mg(2+) as a Lewis acid. The simulations also indicated that the overall structure and thermal motions of the ribozyme are not significantly influenced by the catalytic Mg(2+) interaction or C75 protonation. This analysis suggests that the reaction pathway of the ribozyme is dominated by small local motions at the active site rather than large-scale global conformational changes. These results are consistent with a wealth of experimental data.

Isotope labeling strategies for NMR studies of RNA

The known biological functions of RNA have expanded in recent years and now include gene regulation, maintenance of sub-cellular structure, and catalysis, in addition to propagation of genetic information. As for proteins, RNA function is tightly correlated with structure. Unlike proteins, structural information for larger, biologically functional RNAs is relatively limited. NMR signal degeneracy, relaxation problems, and a paucity of long-range (1)H-(1)H dipolar contacts have limited the utility of traditional NMR approaches. Selective isotope labeling, including nucleotide-specific and segmental labeling strategies, may provide the best opportunities for obtaining structural information by NMR. Here we review methods that have been developed for preparing and purifying isotopically labeled RNAs, as well as NMR strategies that have been employed for signal assignment and structure determination.

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.

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.