Catalytic strategies of the hepatitis delta virus ribozymes

New Insights into the Dependence of CPEB3 Ribozyme Cleavage on Mn2+ and Mg2

CPEB3 ribozyme is a self-cleaving RNA that occurs naturally in mammals and requires divalent metal ions for efficient activity. Ribozymes exhibit preferences for specific metal ions, but the exact differences in the catalytic mechanisms of various metal ions on the CPEB3 ribozyme remain unclear. Our findings reveal that Mn2+ functions as a more effective cofactor for CPEB3 ribozyme catalysis compared to Mg2+, as confirmed by its stronger binding affinity to CPEB3 by EPR. Cleavage assays of CPEB3 mutants and molecular docking analyses further showed that excessive Mn2+ ions can bind to a second binding site near the catalytic site, hindering CPEB3 catalytic efficiency and contributing to the Mn2+ bell-shaped curve. These results implicate a pivotal role for the local nucleobase-Mn2+ interactions in facilitating RNA folding and modulating the directed attack of nucleophilic reagents. Our study provides new insights and experimental evidence for exploring the divalent cation dependent cleavage mechanism of the CPEB3 ribozyme.

Novel ribozymes: discovery, catalytic mechanisms, and the quest to understand biological function

Small endonucleolytic ribozymes promote the self-cleavage of their own phosphodiester backbone at a specific linkage. The structures of and the reactions catalysed by members of individual families have been studied in great detail in the past decades. In recent years, bioinformatics studies have uncovered a considerable number of new examples of known catalytic RNA motifs. Importantly, entirely novel ribozyme classes were also discovered, for most of which both structural and biochemical information became rapidly available. However, for the majority of the new ribozymes, which are found in the genomes of a variety of species, a biological function remains elusive. Here, we concentrate on the different approaches to find catalytic RNA motifs in sequence databases. We summarize the emerging principles of RNA catalysis as observed for small endonucleolytic ribozymes. Finally, we address the biological functions of those ribozymes, where relevant information is available and common themes on their cellular activities are emerging. We conclude by speculating on the possibility that the identification and characterization of proteins that we hypothesize to be endogenously associated with catalytic RNA might help in answering the ever-present question of the biological function of the growing number of genomically encoded, small endonucleolytic ribozymes.

Generation of 2',3'-Cyclic Phosphate-Containing RNAs as a Hidden Layer of the Transcriptome

Cellular RNA molecules contain phosphate or hydroxyl ends. A 2′,3′-cyclic phosphate (cP) is one of the 3′-terminal forms of RNAs mainly generated from RNA cleavage by ribonucleases. Although transcriptome profiling using RNA-seq has become a ubiquitous tool in biological and medical research, cP-containing RNAs (cP-RNAs) form a hidden transcriptome layer, which is infrequently recognized and characterized, because standard RNA-seq is unable to capture them. Despite cP-RNAs’ invisibility in RNA-seq data, increasing evidence indicates that they are not accumulated simply as non-functional degradation products; rather, they have physiological roles in various biological processes, designating them as noteworthy functional molecules. This review summarizes our current knowledge of cP-RNA biogenesis pathways and their catalytic enzymatic activities, discusses how the cP-RNA generation affects biological processes, and explores future directions to further investigate cP-RNA biology.

Evidence That Nucleophile Deprotonation Exceeds Bond Formation in the HDV Ribozyme Transition State

Steric constraints imposed by the active sites of protein and RNA enzymes pose major challenges to the investigation of structure-function relationships within these systems. As a strategy to circumvent such constraints in the HDV ribozyme, we have synthesized phosphoramidites from propanediol derivatives and incorporated them at the 5'-termini of RNA and DNA oligonucleotides to generate a series of novel substrates with nucleophiles perturbed electronically through geminal fluorination. In nonenzymatic, hydroxide-catalyzed intramolecular transphosphorylation of the DNA substrates, pH-rate profiles revealed that fluorine substitution reduces the maximal rate and the kinetic p K_a, consistent with the expected electron-withdrawing effect. In HDV ribozyme reactions, we observed that the RNA substrates undergo transphosphorylation relatively efficiently, suggesting that the conformational constraints imposed by a ribofuranose ring are not strictly required for ribozyme catalysis. In contrast to the nonenzymatic reactions, however, substrate fluorination modestly increases the ribozyme reaction rate, consistent with a mechanism in which (1) the 2'-hydroxyl nucleophile exists predominantly in its neutral, protonated form in the ground state and (2) the 2'-hydroxyl bears some negative charge in the rate-determining step, consistent with a transition state in which the extent of 2'-OH deprotonation exceeds the extent of P-O bond formation.

Analyses of a whole-genome inter-clade recombination map of hepatitis delta virus suggest a host polymerase-driven and viral RNA structure-promoted template-switching mechanism for viral RNA recombination

The genome of hepatitis delta virus (HDV) is a 1.7-kb single-stranded circular RNA that folds into an unbranched rod-like structure and has ribozyme activity. HDV redirects host RNA polymerase(s) (RNAP) to perform viral RNA-directed RNA transcription. RNA recombination is known to contribute to the genetic heterogeneity of HDV, but its molecular mechanism is poorly understood. Here, we established a whole-genome HDV-1/HDV-4 recombination map using two cloned sequences coexisting in cultured cells. Our functional analyses of the resulting chimeric delta antigens (the only viral-encoded protein) and recombinant genomes provide insights into how recombination promotes the genotypic and phenotypic diversity of HDV. Our examination of crossover distribution and subsequent mutagenesis analyses demonstrated that ribozyme activity on HDV genome, which is required for viral replication, also contributes to the generation of an inter-clade junction. These data provide circumstantial evidence supporting our contention that HDV RNA recombination occurs via a replication-dependent mechanism. Furthermore, we identify an intrinsic asymmetric bulge on the HDV genome, which appears to promote recombination events in the vicinity. We therefore propose a mammalian RNAP-driven and viral-RNA-structure-promoted template-switching mechanism for HDV genetic recombination. The present findings improve our understanding of the capacities of the host RNAP beyond typical DNA-directed transcription.

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.

Specific Reaction Parametrization of the AM1/d Hamiltonian for Phosphoryl Transfer Reactions: H, O, and P Atoms

A semiempirical AM1/d Hamiltonian is developed to model phosphoryl transfer reactions catalyzed by enzymes and ribozymes for use in linear-scaling calculations and combined quantum mechanical/molecular mechanical simulations. The model, designated AM1/d-PhoT, is parametrized for H, O, and P atoms to reproduce high-level density-functional results from a recently constructed database of quantum calculations for RNA catalysis ( http://theory.chem.umn.edu/Database/QCRNA ), including geometries and relative energies of minima, transition states and reactive intermediates, dipole moments, proton affinities, and other relevant properties. The model is tested in the gas phase and in solution using a QM/MM potential. The results indicate that the method provides significantly higher accuracy than MNDO/d, AM1, and PM3 methods and, for the transphosphorylation reactions, is in close agreement with the density-functional calculations at the B3LYP/6-311++G(3df,2p) level with a reduction in computational cost of 3-4 orders of magnitude. The model is expected to have considerable impact on the application of semiempirical QM/MM methods to transphosphorylation reactions in solution, enzymes, and ribozymes and to ultimately facilitate the design of improved next-generation multiscale quantum models.

Force Field for Mg(2+), Mn(2+), Zn(2+), and Cd(2+) Ions That Have Balanced Interactions with Nucleic Acids

Divalent metal ions are of fundamental importance to the function and folding of nucleic acids. Divalent metal ion-nucleic acid interactions are complex in nature and include both territorial and site specific binding. Commonly employed nonbonded divalent ion models, however, are often parametrized against bulk ion properties and are subsequently utilized in biomolecular simulations without considering any data related to interactions at specific nucleic acid sites. Previously, we assessed the ability of 17 different nonbonded Mg(2+) ion models to reproduce different properties of Mg(2+) in aqueous solution including radial distribution functions, solvation free energies, water exchange rates, and translational diffusion coefficients. In the present work, we depart from the recently developed 12-6-4 potential models for divalent metal ions developed by Li and Merz and tune the pairwise parameters for Mg(2+), Mn(2+), Zn(2+), and Cd(2+) binding dimethyl phosphate, adenosine, and guanosine in order to reproduce experimental site specific binding free energies derived from potentiometric pH titration data. We further apply these parameters to investigate a metal ion migration previously proposed to occur during the catalytic reaction of the hammerhead ribozyme. The new parameters are shown to be accurate and balanced for nucleic acid binding in comparison with available experimental data and provide an important tool for molecular dynamics and free energy simulations of nucleic acids where these ions may exhibit different binding modes.

Therapeutic Strategies and New Intervention Points in Chronic Hepatitis Delta Virus Infection

Chronic hepatitis delta virus infection (CHD) is a condition arising from super-infection of hepatitis B virus (HBV)-infected patients, resulting in a more rapid advance in liver pathology and hepatocellular carcinoma than is observed for HBV mono-infection. Although hepatitis delta virus (HDV) is structurally simple, its life cycle involves the complex participation of host enzymes, HBV-derived surface antigen (HBsAg), and HDV-auto-ribozyme and hepatitis delta antigen (HDAg) activities. Unsatisfactory clinical trial results with interferon-based therapies are motivating researchers to adjust and redirect the approach to CHD drug development. This new effort will likely require additional structural and functional studies of the viral and cellular/host components involved in the HDV replication cycle. This review highlights recent work aimed at new drug interventions for CHD, with interpretation of key pre-clinical- and clinical trial outcomes and a discussion of promising new technological approaches to antiviral drug design.

Transition State Features in the Hepatitis Delta Virus Ribozyme Reaction Revealed by Atomic Perturbations

Endonucleolytic ribozymes constitute a class of non-coding RNAs that catalyze single-strand RNA scission. With crystal structures available for all of the known ribozymes, a major challenge involves relating functional data to the physically observed RNA architecture. In the case of the hepatitis delta virus (HDV) ribozyme, there are three high-resolution crystal structures, the product state of the reaction and two precursor variants, with distinct mechanistic implications. Here, we develop new strategies to probe the structure and catalytic mechanism of a ribozyme. First, we use double-mutant cycles to distinguish differences in functional group proximity implicated by the crystal structures. Second, we use a corrected form of the Brønsted equation to assess the functional significance of general acid catalysis in the system. Our results delineate the functional relevance of atomic interactions inferred from structure, and suggest that the HDV ribozyme transition state resembles the cleavage product in the degree of proton transfer to the leaving group.

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

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

Silencing of amyloid precursor protein expression using a new engineered delta ribozyme

Alzheimer's disease (AD) etiological studies suggest that an elevation in amyloid-β peptides (Aβ) level contributes to aggregations of the peptide and subsequent development of the disease. The major constituent of these amyloid peptides is the 1 to 40–42 residue peptide (Aβ₄₀₋₄₂) derived from amyloid protein precursor (APP). Most likely, reducing Aβ levels in the brain may block both its aggregation and neurotoxicity and would be beneficial for patients with AD. Among the several possible ways to lower Aβ accumulation in the cells, we have selectively chosen to target the primary step in the Aβ cascade, namely, to reduce APP gene expression. Toward this end, we engineered specific SOFA-HDV ribozymes, a new generation of catalytic RNA tools, to decrease APP mRNA levels. Additionally, we demonstrated that APP-ribozymes are effective at decreasing APP mRNA and protein levels as well as Aβ levels in neuronal cells. Our results could lay the groundwork for a new protective treatment for AD.

Characterization of the Structure and Dynamics of the HDV Ribozyme at Different Stages Along the Reaction Path

The structure and dynamics of the hepatitis delta virus ribozyme (HDVr) are studies using molecular dynamics simulations at several stages along its catalytic reaction path, including reactant, activated precursor, transition state mimic and product states, departing from an initial structure based on the C75U mutant crystal structure (PDB: 1VC7). Results of five 350 ns molecular dynamics simulations reveal a spontaneous rotation of U-1 that leads to an in-line conformation and support the role of protonated C75 as the general acid in the transition state. Our results provide rationale for the interpretation of several important experimental results, and make experimentally testable predictions regarding the roles of key active site residues that are not obvious from any available crystal structures.

Structure-function analysis from the outside in: long-range tertiary contacts in RNA exhibit distinct catalytic roles

The conserved catalytic core of the Tetrahymena group I ribozyme is encircled by peripheral elements. We have conducted a detailed structure-function study of the five long-range tertiary contacts that fasten these distal elements together. Mutational ablation of each of the tertiary contacts destabilizes the folded ribozyme, indicating a role of the peripheral elements in overall stability. Once folded, three of the five tertiary contact mutants exhibit defects in overall catalysis that range from 20- to 100-fold. These and the subsequent results indicate that the structural ring of peripheral elements does not act as a unitary element; rather, individual connections have distinct roles as further revealed by kinetic and thermodynamic dissection of the individual reaction steps. Ablation of P14 or the metal ion core/metal ion core receptor (MC/MCR) destabilizes docking of the substrate-containing P1 helix into tertiary interactions with the ribozyme's conserved core. In contrast, ablation of the L9/P5 contact weakens binding of the guanosine nucleophile by slowing its association, without affecting P1 docking. The P13 and tetraloop/tetraloop receptor (TL/TLR) mutations had little functional effect and small, local structural changes, as revealed by hydroxyl radical footprinting, whereas the P14, MC/MCR, and L9/P5 mutants show structural changes distal from the mutation site. These changes extended into regions of the catalytic core involved in docking or guanosine binding. Thus, distinct allosteric pathways couple the long-range tertiary contacts to functional sites within the conserved core. This modular functional specialization may represent a fundamental strategy in RNA structure-function interrelationships.

Ribozymes: Flexible molecular devices at work

The discovery of ribozymes, RNAs with catalytic activity, revealed the extraordinary characteristic of this molecule, and corroborated the idea that RNA was the first informative polymer. The "RNA world" hypothesis asserts that the DNA/RNA/PROTEIN world arose from an earlier RNA world in which were present only RNA molecules able to perform both of the two functions performed separately by DNA and proteins in the present-day cells: the ability to transfer genetic information and to carry out catalytic activity. The catalytic properties of ribozymes are exclusively due to the capacity of RNA molecules to assume particular structures. Moreover, the structural versatility of RNA can allow to a single RNA sequence to fold in more than one structure, able to perform more than one function. In the first part of this work we will discuss the RNA plasticity, focusing on "bifunctional" ribozymes isolated by in vitro selection experiments, and on the consequences of this plasticity in the prospective of the emergence of new specific functions. The possibility that one sequence could have more than one structure/function, greatly increase the evolutionary potential of RNA, and the capacity of RNA to switch from a structure/function to another is probably one of the reasons of the evolutionary success also in modern-day cells. Naturally occurring ribozymes discovered in contemporary cells, demonstrate the crucial role that ribozymes still have in the modern protein world. In the second part of this paper we will discuss the capacity of natural ribozymes to modulate gene expression making use of their exclusive catalytic properties. Moreover, we will consider the possibility of their ancient origin.

The pH dependence of hairpin ribozyme catalysis reflects ionization of an active site adenine

Understanding how self-cleaving ribozymes mediate catalysis is crucial in light of compelling evidence that human and bacterial gene expression can be regulated through RNA self-cleavage. The hairpin ribozyme catalyzes reversible phosphodiester bond cleavage through a mechanism that does not require divalent metal cations. Previous structural and biochemical evidence implicated the amidine group of an active site adenosine, A38, in a pH-dependent step in catalysis. We developed a way to determine microscopic pK(a) values in active ribozymes based on the pH-dependent fluorescence of 8-azaadenosine (8azaA). We compared the microscopic pK(a) for ionization of 8azaA at position 38 with the apparent pK(a) for the self-cleavage reaction in a fully functional hairpin ribozyme with a unique 8azaA at position 38. Microscopic and apparent pK(a) values were virtually the same, evidence that A38 protonation accounts for the decrease in catalytic activity with decreasing pH. These results implicate the neutral unprotonated form of A38 in a transition state that involves formation of the 5'-oxygen-phosphorus bond.

Evidence that the polymerase of respiratory syncytial virus initiates RNA replication in a nontemplated fashion

RNA virus polymerases must initiate replicative RNA synthesis with extremely high accuracy to maintain their genome termini and to avoid generating defective genomes. For the single-stranded negative-sense RNA viruses, it is not known how this accuracy is achieved. To investigate this question, mutations were introduced into the 3' terminal base of a respiratory syncytial virus (RSV) template, and the RNA products were examined to determine the impact of the mutation. To perform the assay, RNA replication was reconstituted using a modified minireplicon system in which replication was limited to a single step. Importantly, this system allowed analysis of RSV RNA generated intracellularly, but from a defined template that was not subject to selection by replication. Sequence analysis of RNA products generated from templates containing 1U-C and 1U-A substitutions showed that, in both cases, replication products were initiated with a nontemplated, WT A residue, rather than a templated G or U residue, indicating that the polymerase selects the terminal NTP independently of the template. Examination of a template in which the position 1 nucleotide was deleted supported these findings. This mutant directed efficient replication at approximately 60% of WT levels, and its product was found to be initiated at the WT position (-1 relative to the template) with a WT A residue. These findings show that the RSV replicase selects ATP and initiates at the correct position, independently of the first nucleotide of the template, suggesting a mechanism by which highly accurate replication initiation is achieved.

Monitoring of an RNA multistep folding pathway by isothermal titration calorimetry

Isothermal titration calorimetry was used to monitor the energetic landscape of a catalytic RNA, specifically that of the hepatitis delta virus ribozyme. Using mutants that isolated various tertiary interactions, the thermodynamic parameters of several ribozyme-substrate intermediates were determined. The results shed light on the impact of several tertiary interactions on the global structure of the ribozyme. In addition, the data indicate that the formation of the P1.1 pseudoknot is the limiting step of the molecular mechanism. Last, as illustrated here, isothermal titration calorimetry appears to be a method of choice for the elucidation of an RNA's folding pathway.

Competition between Co(NH(3)(6)3+ and inner sphere Mg2+ ions in the HDV ribozyme

Divalent cations play critical structural and functional roles in many RNAs. While the hepatitis delta virus (HDV) ribozyme can undergo self-cleavage in the presence of molar concentrations of monovalent cations, divalent cations such as Mg(2+) are required for efficient catalysis under physiological conditions. Moreover, the cleavage reaction can be inhibited with Co(NH(3))(6)(3+), an analogue of Mg(H(2)O)(6)(2+). Here, the binding of Mg(2+) and Co(NH(3))(6)(3+) to the HDV ribozyme is studied by Raman microscopic analysis of crystals. Raman difference spectra acquired at different metal ion conditions reveal changes in the ribozyme. When Mg(2+) alone is introduced to the ribozyme, inner sphere coordination of Mg(H(2)O)(x)(2+) (x </= 5) to nonbridging PO(2)(-) oxygen and changes in base stretches and phosphodiester group conformation are observed. In addition, binding of Mg(2+) induces deprotonation of a cytosine assigned to the general acid C75, consistent with solution studies. When Co(NH(3))(6)(3+) alone is introduced, deprotonation of C75 is again observed, as are distinctive changes in base vibrational ring modes and phosphodiester backbone conformation. In contrast to Mg(2+) binding, Co(NH(3))(6)(3+) binding does not perturb PO(2)(-) group vibrations, consistent with its ability to make only outer sphere contacts. Surprisingly, competitive binding studies reveal that Co(NH(3))(6)(3+) ions displace some inner sphere-coordinated magnesium species, including ions coordinated to PO(2)(-) groups or the N7 of a guanine, likely G1 at the active site. These observations contrast with the tenet that Co(NH(3))(6)(3+) ions displace only outer sphere magnesium ions. Overall, our data support two classes of inner sphere Mg(2+)-PO(2)(-) binding sites: sites that Co(NH(3))(6)(3+) can displace and others it cannot.

Modulating RNA structure and catalysis: lessons from small cleaving ribozymes

RNA is a key molecule in life, and comprehending its structure/function relationships is a crucial step towards a more complete understanding of molecular biology. Even though most of the information required for their correct folding is contained in their primary sequences, we are as yet unable to accurately predict both the folding pathways and active tertiary structures of RNA species. Ribozymes are interesting molecules to study when addressing these questions because any modifications in their structures are often reflected in their catalytic properties. The recent progress in the study of the structures, the folding pathways and the modulation of the small ribozymes derived from natural, self-cleaving, RNA motifs have significantly contributed to today's knowledge in the field.

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.

Natural and unnatural ribozymes: back to the primordial RNA world

We review natural and in vitro selected ribozymes, for which combined studies could provide us with both insight into the functions performed by ancient RNA molecules in a primitive RNA world and a hypothesis about evolutionary steps that led to the contemporary world.

The use of simple model systems to study spliceosomal catalysis

Since direct analysis of many aspects of spliceosomal function is greatly hindered by the daunting complexity of the spliceosome, the development of functionally validated simple model systems can be of great value. The critical role played by a base-paired complex of U6 and U2 snRNAs in splicing in vivo suggests that this complex could be a suitable starting point for the development of such a simple model system. However, several criteria must be satisfied before such a snRNA-based in vitro system can be considered a valid model for the spliceosomal catalytic core, including similarities at the level of reaction chemistry and cationic and sequence requirements. Previous functional analyses of in vitro assembled base-paired complexes of human U2 and U6 snRNAs have been promising, providing insight into catalysis. Furthermore, they strongly suggest that with further optimization, these RNAs might indeed be able to recapitulate the function of the spliceosomal catalytic core, thus opening the door to several lines of study not previously possible.

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.

How RNA unfolds and refolds

Understanding how RNA folds and what causes it to unfold has become more important as knowledge of the diverse functions of RNA has increased. Here we review the contributions of single-molecule experiments to providing answers to questions such as: How much energy is required to unfold a secondary or tertiary structure? How fast is the process? How do helicases unwind double helices? Are the unwinding activities of RNA-dependent RNA polymerases and of ribosomes different from other helicases? We discuss the use of optical tweezers to monitor the unfolding activities of helicases, polymerases, and ribosomes, and to apply force to unfold RNAs directly. We also review the applications of fluorescence and fluorescence resonance energy transfer to measure RNA dynamics.

Potassium ions modulate a G-quadruplex-ribozyme's activity

Hepatitis delta virus ribozyme folds into a tightly packed tertiary structure. However, unlike other ribozymes, it does not appear to be able to follow alternative folding pathways. Molecular engineering of the hepatitis delta virus ribozyme led to the development of a ribozyme possessing an endoribonuclease activity that is under the control of a G-quadruplex structure (i.e., a G-quartzyme). This latter species represents an entirely new class of ribozyme. Mutants of this ribozyme were then generated in order to shed light on the modulation of the cleavage activity caused by the presence of the G-quadruplex structure. Kinetic characterization of the G-quartzyme was performed under various single turnover conditions. It was found to be active only in the presence of potassium cations that act as counter ions in the positioning of the four coplanar guanines that form the building block of the G-quadruplex structure. The G-quartzyme behaves as an allosteric ribozyme, with the potassium cations acting as positive effectors with a Hill coefficient of 2.9 +/- 0.2. The conformation transition caused by the presence of the potassium ions is supported by enzymatic and chemical probing of both the inactive (off) and active (on) structures. This study shows that it is possible to interfere with the tight structure of the hepatitis delta virus ribozyme by adding an unusual, stable structure. To our knowledge, the G-quartzyme is the sole ribozyme that exhibits a monovalent cation-dependent activity.

Quantum mechanical/molecular mechanical simulation study of the mechanism of hairpin ribozyme catalysis

The molecular mechanism of hairpin ribozyme catalysis is studied with molecular dynamics simulations using a combined quantum mechanical and molecular mechanical (QM/MM) potential with a recently developed semiempirical AM1/d-PhoT model for phosphoryl transfer reactions. Simulations are used to derive one- and two-dimensional potentials of mean force to examine specific reaction paths and assess the feasibility of proposed general acid and base mechanisms. Density-functional calculations of truncated active site models provide complementary insight to the simulation results. Key factors utilized by the hairpin ribozyme to enhance the rate of transphosphorylation are presented, and the roles of A38 and G8 as general acid and base catalysts are discussed. The computational results are consistent with available experimental data, provide support for a general acid/base mechanism played by functional groups on the nucleobases, and offer important insight into the ability of RNA to act as a catalyst without explicit participation by divalent metal ions.

Structure-function relationships in RNA and RNP enzymes: recent advances

The structural biology of ribozymes and ribonucleoprotein (RNP) enzymes is now sufficiently advanced that a true dialogue between structural and functional studies is possible. In this review, we consider three important systems in which an integration of structural and biochemical data has recently led to major advances in mechanistic understanding. In the hammerhead ribozyme, application-driven biochemical studies led to the discovery of a key structural interaction that had been omitted from previously-studied constructs. A new crystal structure of the resulting, tertiary-stabilized hammerhead has resolved a remarkable number of longstanding paradoxes in the structure-function relationship of this ribozyme. In the Group I intron ribozyme, a flurry of high-resolution structures has largely confirmed, but in some cases refined or challenged, a detailed model of a metalloenzyme active site that had previously been derived by meticulous quantitative metal ion rescue experiments. Finally, for the peptidyl transferase center of the ribosome, recent biochemical and chemical results motivated by the pioneering crystal structures have suggested a picture of a catalytic mechanism dominated by proximity and orientation effects and substrate-assisted catalysis. These results refocus attention on catalysis as a property of the integrated RNP machinery as a whole, as opposed to a narrow concern with the RNA functional groups in immediate contact with the reactive center.

A novel structural rearrangement of hepatitis delta virus antigenomic ribozyme

A bioinformatic covariation analysis of a collection of 119 novel variants of the antigenomic, self-cleaving hepatitis delta virus (HDV) RNA motif supported the formation of all of the Watson–Crick base pairs (bp) of the catalytic centre except the C19–G81 pair located at the bottom of the P2 stem. In fact, a novel Watson–Crick bp between C19 and G80 is suggested by the data. Both chemical and enzymatic probing demonstrated that initially the C19–G81 pair is formed in the ribozyme (Rz), but upon substrate (S) binding and the formation of the P1.1 pseudoknot C19 switches its base-pairing partner from G81 to G80. As a result of this finding, the secondary structure of this ribozyme has been redrawn. The formation of the C19–G80 bp results in a J4/2 junction composed of four nucleotides, similar to that seen in the genomic counterpart, thereby increasing the similarities between these two catalytic RNAs. Additional mutagenesis, cleavage activity and probing experiments yield an original characterization of the structural features involving the residues of the J4/2 junction.

RNA recombination in hepatitis delta virus: Implications regarding the abilities of mammalian RNA polymerases

Hepatitis delta virus (HDV) requires the surface antigens of hepatitis B virus (HBV) for packaging and transmission, but replicates its RNA in an HBV-independent fashion. HDV contains a 1.7-kb circular RNA genome that is folded into an unbranched rod-like structure via intramolecular base-pairing, and possesses ribozyme activity. The HDV genome does not encode an RNA-dependent RNA polymerase (RdRp), but is instead replicated by host RNA polymerase(s) via a rolling-circle mechanism. As such, HDV is similar to the viroid plant pathogens. Recent findings suggest that HDV can also undergo template-switching recombination, a well-documented process that has been found in a large number of RdRp-encoding RNA viruses and is thought to affect viral evolution and pathogenesis. This mini-review examines HDV RNA recombination and how it may improve our understanding of the capacities of host RNA polymerases beyond typical DNA-directed transcription, and speculates on the role of host RNA polymerase-directed RNA template-switching in the origin of HDV.

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.

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. .

CHARMM force field parameters for simulation of reactive intermediates in native and thio-substituted ribozymes

Force field parameters specifically optimized for residues important in the study of RNA catalysis are derived from density-functional calculations, in a fashion consistent with the CHARMM27 all-atom empirical force field. Parameters are presented for residues that model reactive RNA intermediates and transition state analogs, thio-substituted phosphates and phosphoranes, and bound Mg(2+) and di-metal bridge complexes. Target data was generated via density-functional calculations at the B3LYP/6-311++G(3df,2p)// B3LYP/6-31++G(d,p) level. Partial atomic charges were initially derived from CHelpG electrostatic potential fitting and subsequently adjusted to be consistent with the CHARMM27 charges. Lennard-Jones parameters were determined to reproduce interaction energies with water molecules. Bond, angle, and torsion parameters were derived from the density-functional calculations and renormalized to maintain compatibility with the existing CHARMM27 parameters for standard residues. The extension of the CHARMM27 force field parameters for the nonstandard biological residues presented here will have considerable use in simulations of ribozymes, including the study of freeze-trapped catalytic intermediates, metal ion binding and occupation, and thio effects.

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.

Examination of the folding pathway of the antigenomic hepatitis delta virus ribozyme reveals key interactions of the L3 loop

With the goal of gaining insight into the tertiary structure of the hepatitis delta virus ribozyme, cross-linking experiments using 4-thiouridine residues introduced in either the 5'-end portion of the substrate, or at seven strategic positions within the ribozyme, were performed. Mapping of the newly formed covalent bonds in cross-linked species obtained under various conditions, as well as using several mutated ribozymes, permitted monitoring of the formation of the ribozyme-substrate complex as the ribozyme proceeded along the folding pathway. In order to aid visualization of the tertiary structure transformation, an in silico animation of the "on" folding pathway was developed. In combination with those of the cleavage assays of structured substrates, these data shed light on the key contribution of the L3 loop in the formation of an active tertiary complex.

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.

Inhibition of the catalytic activity of trans-acting antigenomic δ ribozyme by selected antibiotics and their Cu2+ complexes

The effect of selected 10 antibiotics and their complexes with Cu(2+) ions on the catalytic activity of the trans-acting antigenomic delta ribozyme was investigated. Sisomicin, vancomycin, and actinomycin D displayed weak inhibitory properties. However, much stronger effects were detected with complexes of these antibiotics with Cu(2+) ions. The strongest inhibition was observed with actinomycin D-Cu(2+) complex, for which the calculated K(i) value was reduced ca. 35-fold upon metal ion complexation. We postulate that the antibiotic-Cu(2+) complexes are guided to the ribozyme metal ion binding site(s) presumably displacing the catalytically important metal ion(s). Moreover, we assume that, once positioned in appropriate distances to RNA phosphate groups and bases, the coordinated Cu(2+) ions become positively charged factors that enhance the affinity of the antibiotics to the ribozyme. These observations indicate that coordination of metal ions to antibiotics substantially changes their properties which might also have a biological relevance inside the cell.