Terbium-mediated Footprinting Probes a Catalytic Conformational Switch in the Antigenomic Hepatitis Delta Virus Ribozyme

Inhibition of Cpeb3 ribozyme elevates CPEB3 protein expression and polyadenylation of its target mRNAs and enhances object location memory

A self-cleaving ribozyme that maps to an intron of the cytoplasmic polyadenylation element-binding protein 3 (Cpeb3) gene is thought to play a role in human episodic memory, but the underlying mechanisms mediating this effect are not known. We tested the activity of the murine sequence and found that the ribozyme's self-scission half-life matches the time it takes an RNA polymerase to reach the immediate downstream exon, suggesting that the ribozyme-dependent intron cleavage is tuned to co-transcriptional splicing of the Cpeb3 mRNA. Our studies also reveal that the murine ribozyme modulates maturation of its harboring mRNA in both cultured cortical neurons and the hippocampus: inhibition of the ribozyme using an antisense oligonucleotide leads to increased CPEB3 protein expression, which enhances polyadenylation and translation of localized plasticity-related target mRNAs, and subsequently strengthens hippocampal-dependent long-term memory. These findings reveal a previously unknown role for self-cleaving ribozyme activity in regulating experience-induced co-transcriptional and local translational processes required for learning and memory.

Inhibition of CPEB3 ribozyme elevates CPEB3 protein expression and polyadenylation of its target mRNAs, and enhances object location memory

A self-cleaving ribozyme that maps to an intron of the cytoplasmic polyadenylation element binding protein 3 (CPEB3) gene is thought to play a role in human episodic memory, but the underlying mechanisms mediating this effect are not known. We tested the activity of the murine sequence and found that the ribozyme's self-scission half-life matches the time it takes an RNA polymerase to reach the immediate downstream exon, suggesting that the ribozyme-dependent intron cleavage is tuned to co-transcriptional splicing of the CPEB3 mRNA. Our studies also reveal that the murine ribozyme modulates maturation of its harboring mRNA in both cultured cortical neurons and the hippocampus: inhibition of the ribozyme using an antisense oligonucleotide leads to increased CPEB3 protein expression, which enhances polyadenylation and translation of localized plasticity-related target mRNAs, and subsequently strengthens hippocampal-dependent long-term memory. These findings reveal a previously unknown role for self-cleaving ribozyme activity in regulating experience-induced co-transcriptional and local translational processes required for learning and memory.

Systematic detection of tertiary structural modules in large RNAs and RNP interfaces by Tb-seq

Compact RNA structural motifs control many aspects of gene expression, but we lack methods for finding these structures in the vast expanse of multi-kilobase RNAs. To adopt specific 3-D shapes, many RNA modules must compress their RNA backbones together, bringing negatively charged phosphates into close proximity. This is often accomplished by recruiting multivalent cations (usually Mg2+), which stabilize these sites and neutralize regions of local negative charge. Coordinated lanthanide ions, such as terbium (III) (Tb3+), can also be recruited to these sites, where they induce efficient RNA cleavage, thereby revealing compact RNA 3-D modules. Until now, Tb3+ cleavage sites were monitored via low-throughput biochemical methods only applicable to small RNAs. Here we present Tb-seq, a high-throughput sequencing method for detecting compact tertiary structures in large RNAs. Tb-seq detects sharp backbone turns found in RNA tertiary structures and RNP interfaces, providing a way to scan transcriptomes for stable structural modules and potential riboregulatory motifs.

The structural features of the ligand-free moaA riboswitch and its ion-dependent folding

Riboswitches are structural elements of mRNA involved in the regulation of gene expression by responding to specific cellular metabolites. To fulfil their regulatory function, riboswitches prefold into an active state, the so-called binding competent form, that guarantees metabolite binding and allows a consecutive refolding of the RNA. Here, we describe the folding pathway to the binding competent form as well as the ligand free structure of the moaA riboswitch of E. coli. This RNA proposedly responds to the molybdenum cofactor (Moco), a highly oxygen-sensitive metabolite, essential in the carbon and sulfur cycles of eukaryotes. K⁺- and Mg²⁺-dependent footprinting assays and spectroscopic investigations show a high degree of structure formation of this RNA already at very low ion-concentrations. Mg²⁺ facilitates additionally a general compaction of the riboswitch towards its proposed active structure. We show that this fold agrees with the earlier suggested secondary structure which included also a long-range tetraloop/tetraloop-receptor like interaction. Metal ion cleavage assays revealed specific Mg²⁺-binding pockets within the moaA riboswitch. These Mg²⁺ binding pockets are good indicators for the potential Moco binding site, since in riboswitches, Mg²⁺ was shown to be necessary to bind phosphate-carrying metabolites. The importance of the phosphate and of other functional groups of Moco is highlighted by binding assays with tetrahydrobiopterin, the reduced and oxygen-sensitive core moiety of Moco. We demonstrate that the general molecular shape of pterin by its own is insufficient for the recognition by the riboswitch.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Tb3+-Cleavage Assays Reveal Specific Mg2+ Binding Sites Necessary to Pre-fold the btuB Riboswitch for AdoCbl Binding

Riboswitches are RNA elements that bind specific metabolites in order to regulate the gene expression involved in controlling the cellular concentration of the respective molecule or ion. Ligand recognition is mostly facilitated by Mg²⁺ mediated pre-organization of the riboswitch to an active tertiary fold. To predict these specific Mg²⁺ induced tertiary interactions of the btuB riboswitch from E. coli, we here report Mg²⁺ binding pockets in its aptameric part in both, the ligand-free and the ligand-bound form. An ensemble of weak and strong metal ion binding sites distributed over the entire aptamer was detected by terbium(III) cleavage assays, Tb³⁺ being an established Mg²⁺ mimic. Interestingly many of the Mⁿ⁺ (n = 2 or 3) binding sites involve conserved bases within the class of coenzyme B₁₂-binding riboswitches. Comparison with the published crystal structure of the coenzyme B₁₂ riboswitch of S. thermophilum aided in identifying a common set of Mⁿ⁺ binding sites that might be crucial for tertiary interactions involved in the organization of the aptamer. Our results suggest that Mⁿ⁺ binding at strategic locations of the btuB riboswitch indeed facilitates the assembly of the binding pocket needed for ligand recognition. Binding of the specific ligand, coenzyme B₁₂ (AdoCbl), to the btuB aptamer does however not lead to drastic alterations of these Mⁿ⁺ binding cores, indicating the lack of a major rearrangement within the three-dimensional structure of the RNA. This finding is strengthened by Tb³⁺ mediated footprints of the riboswitch's structure in its ligand-free and ligand-bound state indicating that AdoCbl indeed induces local changes rather than a global structural rearrangement.

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.

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.

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.

HDV family of self-cleaving ribozymes

The hepatitis delta virus (HDV) ribozymes are catalytic RNAs capable of cleaving their own sugar-phosphate backbone. The HDV virus possesses the ribozymes in both sense and antisense genomic transcripts, where they are essential for processing during replication. These ribozymes have been the subject of intense biochemical scrutiny and have yielded a wealth of mechanistic insights. In recent years, many HDV-like ribozymes have been identified in nearly all branches of life. The ribozymes are implicated in a variety of biological events, including episodic memory in mammals and retrotransposition in many eukaryotes. Detailed analysis of additional HDV-like ribozyme isolates will likely reveal many more biological functions and provide information about the evolution of this unique RNA.

Monitoring global structural changes and specific metal-ion-binding sites in RNA by in-line probing and Tb(III) cleavage

In this chapter we describe the use of two methods, in-line probing as well as terbium(III) cleavage. Both methods can be applied to RNAs of any size, structure, and function. Aside from revealing directly metal ion-binding sites these techniques also provide structural information for longer RNA sequences that are out of range to be analyzed with other techniques such as NMR. The cleavage pattern derived from in-line probing experiments reflects local and overall conformational changes in RNA upon the addition of metal ions, metal complexes, or other ligands. On the other side, terbium(III) cleavage experiments are applied to locate specific metal ion-binding sites in RNA molecules.

Enhanced group II intron retrohoming in magnesium-deficient Escherichia coli via selection of mutations in the ribozyme core

Mobile group II introns are bacterial retrotransposons thought to be evolutionary ancestors of spliceosomal introns and retroelements in eukaryotes. They consist of a catalytically active intron RNA ("ribozyme") and an intron-encoded reverse transcriptase, which function together to promote RNA splicing and intron mobility via reverse splicing of the intron RNA into new DNA sites ("retrohoming"). Although group II introns are active in bacteria, their natural hosts, they function inefficiently in eukaryotes, where lower free Mg(2+) concentrations decrease their ribozyme activity and constitute a natural barrier to group II intron proliferation within nuclear genomes. Here, we show that retrohoming of the Ll.LtrB group II intron is strongly inhibited in an Escherichia coli mutant lacking the Mg(2+) transporter MgtA, and we use this system to select mutations in catalytic core domain V (DV) that partially rescue retrohoming at low Mg(2+) concentrations. We thus identified mutations in the distal stem of DV that increase retrohoming efficiency in the MgtA mutant up to 22-fold. Biochemical assays of splicing and reverse splicing indicate that the mutations increase the fraction of intron RNA that folds into an active conformation at low Mg(2+) concentrations, and terbium-cleavage assays suggest that this increase is due to enhanced Mg(2+) binding to the distal stem of DV. Our findings indicate that DV is involved in a critical Mg(2+)-dependent RNA folding step in group II introns and demonstrate the feasibility of selecting intron variants that function more efficiently at low Mg(2+) concentrations, with implications for evolution and potential applications in gene targeting.

Secondary structure of bacteriophage T4 gene 60 mRNA: implications for translational bypassing

Translational bypassing is a unique phenomenon of bacteriophage T4 gene 60 mRNA wherein the bacterial ribosome produces a single polypeptide chain from a discontinuous open reading frame (ORF). Upon reaching the 50-nucleotide untranslated region, or coding gap, the ribosome either dissociates or bypasses the interruption to continue translating the remainder of the ORF, generating a subunit of a type II DNA topoisomerase. Mutational and computational analyses have suggested that a compact structure, including a stable hairpin, forms in the coding gap to induce bypassing, yet direct evidence is lacking. Here we have probed the secondary structure of gene 60 mRNA with both Tb³⁺ ions and the selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) reagent 1M7 under conditions where bypassing is observed. The resulting experimentally informed secondary structure models strongly support the presence of the predicted coding gap hairpin and highlight the benefits of using Tb³⁺ as a second, complementary probing reagent. Contrary to several previously proposed models, however, the rest of the coding gap is highly reactive with both probing reagents, suggesting that it forms only a short stem-loop. Mutational analyses coupled with functional assays reveal that two possible base-pairings of the coding gap with other regions of the mRNA are not required for bypassing. Such structural autonomy of the coding gap is consistent with its recently discovered role as a mobile genetic element inserted into gene 60 mRNA to inhibit cleavage by homing endonuclease MobA.

siRNA-like double-stranded RNAs are specifically protected against degradation in human cell extract

RNA interference (RNAi) is a set of intracellular pathways in eukaryotes that controls both exogenous and endogenous gene expression. The power of RNAi to knock down (silence) any gene of interest by the introduction of synthetic small-interfering (si)RNAs has afforded powerful insight into biological function through reverse genetic approaches and has borne a new field of gene therapeutics. A number of questions are outstanding concerning the potency of siRNAs, necessitating an understanding of how short double-stranded RNAs are processed by the cell. Recent work suggests unmodified siRNAs are protected in the intracellular environment, although the mechanism of protection still remains unclear. We have developed a set of doubly-fluorophore labeled RNAs (more precisely, RNA/DNA chimeras) to probe in real-time the stability of siRNAs and related molecules by fluorescence resonance energy transfer (FRET). We find that these RNA probes are substrates for relevant cellular degradative processes, including the RNase H1 mediated degradation of an DNA/RNA hybrid and Dicer-mediated cleavage of a 24-nucleotide (per strand) double-stranded RNA. In addition, we find that 21- and 24-nucleotide double-stranded RNAs are relatively protected in human cytosolic cell extract, but less so in blood serum, whereas an 18-nucleotide double-stranded RNA is less protected in both fluids. These results suggest that RNAi effector RNAs are specifically protected in the cellular environment and may provide an explanation for recent results showing that unmodified siRNAs in cells persist intact for extended periods of time.

Metal ions: supporting actors in the playbook of small ribozymes

Since the 1980s, several small RNA motifs capable of chemical catalysis have been discovered. These small ribozymes, composed of between approximately 40 and 200 nucleotides, have been found to play vital roles in the replication of subviral and viral pathogens, as well as in gene regulation in prokaryotes, and have recently been discovered in noncoding eukaryotic RNAs. All of the known natural small ribozymes - the hairpin, hammerhead, hepatitis delta virus, Varkud satellite, and glmS ribozymes--catalyze the same self-cleavage reaction as RNase A, resulting in two products, one bearing a 2'-3' cyclic phosphate and the other a 5'-hydroxyl group. Although originally thought to be obligate metalloenzymes like the group I and II self-splicing introns, the small ribozymes are now known to support catalysis in a wide variety of cations that appear to be only indirectly involved in catalysis. Nevertheless, under physiologic conditions, metal ions are essential for the proper folding and function of the small ribozymes, the most effective of these being magnesium. Metal ions contribute to catalysis in the small ribozymes primarily by stabilizing the catalytically active conformation, but in some cases also by activating RNA functional groups for catalysis, directly participating in catalytic acid-base chemistry, and perhaps by neutralizing the developing negative charge of the transition state. Although interactions between the small ribozymes and cations are relatively nonspecific, ribozyme activity is quite sensitive to the types and concentrations of metal ions present in solution, suggesting a close evolutionary relationship between cellular metal ion homeostasis and cation requirements of catalytic RNAs, and perhaps RNA in general.

Nondenaturing purification of co-transcriptionally folded RNA avoids common folding heterogeneity

Due to the energetic frustration of RNA folding, tertiary structured RNA is typically characterized by a rugged folding free energy landscape where deep kinetic barriers separate numerous misfolded states from one or more native states. While most in vitro studies of RNA rely on (re)folding chemically and/or enzymatically synthesized RNA in its entirety, which frequently leads into kinetic traps, nature reduces the complexity of the RNA folding problem by segmental, co-transcriptional folding starting from the 5' end. We here have developed a simplified, general, nondenaturing purification protocol for RNA to ask whether avoiding denaturation of a co-transcriptionally folded RNA can reduce commonly observed in vitro folding heterogeneity. Our protocol bypasses the need for large-scale auxiliary protein purification and expensive chromatographic equipment and involves rapid affinity capture with magnetic beads and removal of chemical heterogeneity by cleavage of the target RNA from the beads using the ligand-induced glmS ribozyme. For two disparate model systems, the Varkud satellite (VS) and hepatitis delta virus (HDV) ribozymes, we achieve >95% conformational purity within one hour of enzymatic transcription, without the need for any folding chaperones. We further demonstrate that in vitro refolding introduces severe conformational heterogeneity into the natively-purified VS ribozyme but not into the compact, double-nested pseudoknot fold of the HDV ribozyme. We conclude that conformational heterogeneity in complex RNAs can be avoided by co-transcriptional folding followed by nondenaturing purification, providing rapid access to chemically and conformationally pure RNA for biologically relevant biochemical and biophysical studies.

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.

Structure and dynamics of a DNA-based model system for the study of electron spin-spin interactions

We report on the structure and dynamics of a model system for measuring long-range distances in biological macromolecules by saturation-recovery EPR. Four DNA duplexes that incorporate a paramagnetic dysprosium ion (Dy(III)) and a nitroxide spin-label were examined by electron paramagnetic resonance (EPR), circular dichroism (CD), and ultra-violet absorbance (UV) spectroscopy. Dy(III) is chelated by the modified base deoxythymidine-EDTA, (dT-EDTA). Electron spin-spin interactions between the Dy(III) ion and the nitroxide radical are observed at distances as great as approximately 5.3 nm. A slight change in the conformation of those nucleotides lying between the EDTA(Dy(III)) complex and the nitroxide spin-label results in a "stiffening" of the DNA helix on the EPR time scale. Changes in conformation and helix dynamics are due to the binding of the EDTA(Dy(III)) complex to the phosphodiester backbone of the complementary strand. Molecular mechanics calculations indicate that binding occurs in the 5' direction on the complementary strand, at a position 3 or 4 phosphates distant from the dT-EDTA(Dy(III))*dA base pair.

A rugged free energy landscape separates multiple functional RNA folds throughout denaturation

The dynamic mechanisms by which RNAs acquire biologically functional structures are of increasing importance to the rapidly expanding fields of RNA therapeutics and biotechnology. Large energy barriers separating misfolded and functional states arising from alternate base pairing are a well-appreciated characteristic of RNA. In contrast, it is typically assumed that functionally folded RNA occupies a single native basin of attraction that is free of deeply dividing energy barriers (ergodic hypothesis). This assumption is widely used as an implicit basis to interpret experimental ensemble-averaged data. Here, we develop an experimental approach to isolate persistent sub-populations of a small RNA enzyme and show by single molecule fluorescence resonance energy transfer (smFRET), biochemical probing and high-resolution mass spectrometry that commitment to one of several catalytically active folds occurs unexpectedly high on the RNA folding energy landscape, resulting in partially irreversible folding. Our experiments reveal the retention of molecular heterogeneity following the complete loss of all native secondary and tertiary structure. Our results demonstrate a surprising longevity of molecular heterogeneity and advance our current understanding beyond that of non-functional misfolds of RNA kinetically trapped on a rugged folding-free energy landscape.

Assay for glucosamine 6-phosphate using a ligand-activated ribozyme with fluorescence resonance energy transfer or CE-laser-induced fluorescence detection

A naturally occurring aptazyme, the glmS ribozyme, is adapted to an assay for glucosamine 6-phosphate, an effector molecule for the aptazyme. In the assay, binding of analyte allosterically activates aptazyme to cleave a fluorescently labeled oligonucleotide substrate. The extent of reaction, and hence analyte concentration, is detected by either fluorescence resonance energy transfer (FRET) or capillary electrophoresis with laser-induced fluorescence (CE-LIF). With FRET, assay signal is the rate of increase in FRET in presence of analyte. With CE-LIF, the assay signal is the peak height of cleavage product formed after a fixed incubation time. The assay has a linear response up to 100 (CE-LIF) or 500 microM (FRET) and detection limit of approximately 500 nM for glucosamine 6-phosphate under single-turnover conditions. When substrate is present in excess of the aptazyme, it is possible to amplify the signal by multiple turnovers to achieve a 13-fold improvement in sensitivity and detection limit of 50 nM. Successful signal amplification requires a temperature cycle to alternately dissociate cleaved substrate and allow fresh substrate to bind aptazyme. The results show that aptazymes have potential utility as analytical reagents for quantification of effector molecules.

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.

Chemical modification resolves the asymmetry of siRNA strand degradation in human blood serum

Small interfering (si)RNAs have recently been used to therapeutically silence genes in vivo after intravenous systemic delivery. Further progress in the development of siRNA therapeutics will in part rely on tailoring site-specific chemical modifications of siRNAs to optimize their pharmacokinetic properties. Advances are particularly needed to improve the nucleolytic stability of these double-stranded RNA drugs in vivo and suppress adverse off-target gene silencing effects. Here we demonstrate that specific chemical 2'-O-methylation, which has already been shown to ameliorate the omnipresent off-target toxicity of siRNAs, selectively protects the particularly vulnerable 5'-end of the guide strand against exonucleolytic degradation in human blood serum. Specific chemical modification thus resolves the asymmetric degradation of the guide and passenger strands, which is inherent to the thermodynamic asymmetry of the siRNA termini as required for proper utilization of the guide strand in RNA interference.

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.

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.

Structural Features of Target RNA Molecules Greatly Modulate the Cleavage Efficiency of trans-Acting Delta Ribozymes

The aim of this work was to shed some more light on factors influencing the effectiveness of delta ribozyme cleavage of structured RNA molecules. An oligoribonucleotide that corresponds to the 3'-terminal region X of HCV RNA and yeast tRNAPhe were used as representative RNA targets. Only a few sites susceptible to ribozyme cleavage were identified in these targets using a combinatorial library of ribozyme variants, in which the region responsible for ribozyme-target interaction was randomized. On the other hand, the targets were fairly accessible for binding of complementary oligonucleotides, as was shown by 6-mer DNA libraries and RNase H approach. Moreover, the specifically acting ribozymes cleaved the targets precisely but with unexpectedly modest efficacy. To explain these observations, six model RNA molecules were designed, in which the same seven nucleotide long sequence recognized by the delta ribozyme was always single stranded but was embedded into different RNA structural context. These molecules were cleaved with differentiated rates, and the corresponding k2 values were in the range of 0.91-0.021 min-1; thus they differed almost 50-fold. This clearly shows that cleavage of structured RNAs might be much slower than cleavage of a short unstructured oligoribonucleotide, despite full accessibility of the targeted regions for hybridization. Restricted possibilities of conformational transitions, which are necessary to occur on the cleavage reaction trajectory, seem to be responsible for these differences. Their magnitude, which was evaluated in this work, should be taken into account while considering the use of delta ribozymes for practical applications.

Trans-acting glmS catalytic riboswitch: locked and loaded

A recently discovered class of gene regulatory RNAs, coined riboswitches, are commonly found in noncoding segments of bacterial and some eukaryotic mRNAs. Gene up- or down-regulation is triggered by binding of a small organic metabolite, which typically induces an RNA conformational change. Unique among these noncoding RNAs is the glmS catalytic riboswitch, or ribozyme, found in the 5'-untranslated region of the glmS gene in Gram-positive bacteria. It is activated by glucosamine-6-phosphate (GlcN6P), leading to site-specific backbone cleavage of the mRNA and subsequent repression of the glmS gene, responsible for cellular GlcN6P production. Recent biochemical and structural evidence suggests that the GlcN6P ligand acts as a coenzyme and participates in the cleavage reaction without inducing a conformational change. To better understand the role of GlcN6P in solution structural dynamics and function, we have separated the glmS riboswitch core from Bacillus subtilis into a trans-cleaving ribozyme and an externally cleaved substrate. We find that trans cleavage is rapidly activated by nearly 5000-fold to a rate of 4.4 min(-1) upon addition of 10 mM GlcN6P, comparable to the cis-acting ribozyme. Fluorescence resonance energy transfer suggests that this ribozyme-substrate complex does not undergo a global conformational change upon ligand binding in solution. In addition, footprinting at nucleotide resolution using terbium(III) and RNase V1 indicates no significant changes in secondary and tertiary structure upon ligand binding. These findings suggest that the glmS ribozyme is fully folded in solution prior to binding its activating ligand, supporting recent observations in the crystalline state.

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.

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.

Catalytic core structure of the trans-acting HDV ribozyme is subtly influenced by sequence variation outside the core

The human pathogenic hepatitis delta virus (HDV) employs a unique self-cleaving catalytic RNA motif, the HDV ribozyme, during double-rolling circle replication. Fluorescence spectroscopy, circular dichroism, terbium(III) footprinting, and X-ray crystallography of precursor and product forms have revealed that a conformational change accompanies catalysis. In addition, fluorescence resonance energy transfer (FRET) has previously been used on a trans-acting HDV ribozyme to demonstrate surprisingly significant catalytic and global conformational effects of substrate analogues with varying 5' sequences, which reside as dangling overhangs outside the catalytic core. Here, we use the fluorescent guanine analogue 2-aminopurine (AP) in nucleotide position 76, immediately downstream of the catalytically involved C75, to monitor the relative structural effects of these substrate analogues on the ribozyme's trefoil turn of the catalytic core. Steady-state and time-resolved AP fluorescence spectroscopies show that the binding of each substrate analogue induces a unique local conformation with a specific AP76 stacking equilibrium. Binding of the 3' product results in a relative increase in AP fluorescence, suggesting that AP76 becomes more unstacked upon catalysis. These local conformational changes are kinetically concomitant with global conformational changes monitored by FRET. Finally, the rate constant of the local conformational change upon 3' product binding is fast and independent of 3' product concentration yet Mg2+ dependent. Our results demonstrate that the trefoil turn of the HDV ribozyme catalytic core is in a state of dynamic equilibrium not captured by static crystal structures and is highly sensitive to the identity of the 5' sequence and Mg2+ ions.

Cations and hydration in catalytic RNA: molecular dynamics of the hepatitis delta virus ribozyme

The hepatitis delta virus (HDV) ribozyme is an RNA enzyme from the human pathogenic HDV. Cations play a crucial role in self-cleavage of the HDV ribozyme, by promoting both folding and chemistry. Experimental studies have revealed limited but intriguing details on the location and structural and catalytic functions of metal ions. Here, we analyze a total of approximately 200 ns of explicit-solvent molecular dynamics simulations to provide a complementary atomistic view of the binding of monovalent and divalent cations as well as water molecules to reaction precursor and product forms of the HDV ribozyme. Our simulations find that an Mg2+ cation binds stably, by both inner- and outer-sphere contacts, to the electronegative catalytic pocket of the reaction precursor, in a position to potentially support chemistry. In contrast, protonation of the catalytically involved C75 in the precursor or artificial placement of this Mg2+ into the product structure result in its swift expulsion from the active site. These findings are consistent with a concerted reaction mechanism in which C75 and hydrated Mg2+ act as general base and acid, respectively. Monovalent cations bind to the active site and elsewhere assisted by structurally bridging long-residency water molecules, but are generally delocalized.

Solution probing of metal ion binding by helix 27 from Escherichia coli 16S rRNA

Helix (H)27 from Escherichia coli 16S ribosomal (r)RNA is centrally located within the small (30S) ribosomal subunit, immediately adjacent to the decoding center. Bacterial 30S subunit crystal structures depicting Mg(2+) binding sites resolve two magnesium ions within the vicinity of H27: one in the major groove of the G886-U911 wobble pair, and one within the GCAA tetraloop. Binding of such metal cations is generally thought to be crucial for RNA folding and function. To ask how metal ion-RNA interactions in crystals compare with those in solution, we have characterized, using solution NMR spectroscopy, Tb(3+) footprinting and time-resolved fluorescence resonance energy transfer (tr-FRET), location, and modes of metal ion binding in an isolated H27. NMR and Tb(3+) footprinting data indicate that solution secondary structure and Mg(2+) binding are generally consistent with the ribosomal crystal structures. However, our analyses also suggest that H27 is dynamic in solution and that metal ions localize within the narrow major groove formed by the juxtaposition of the loop E motif with the tandem G894-U905 and G895-U904 wobble pairs. In addition, tr-FRET studies provide evidence that Mg(2+) uptake by the H27 construct results in a global lengthening of the helix. We propose that only a subset of H27-metal ion interactions has been captured in the crystal structures of the 30S ribosomal subunit, and that small-scale structural dynamics afforded by solution conditions may contribute to these differences. Our studies thus highlight an example for differences between RNA-metal ion interactions observed in solution and in crystals.

Structural dynamics of precursor and product of the RNA enzyme from the hepatitis delta virus as revealed by molecular dynamics simulations

The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA enzyme involved in the replication of a human pathogen, the hepatitis delta virus. Recent crystal structures of the precursor and product of self-cleavage, together with detailed kinetic analyses, have led to hypotheses on the catalytic strategies employed by the HDV ribozyme. We report molecular dynamics (MD) simulations (approximately 120 ns total simulation time) to test the plausibility that specific conformational rearrangements are involved in catalysis. Site-specific self-cleavage requires cytidine in position 75 (C75). A precursor simulation with unprotonated C75 reveals a rather weak dynamic binding of C75 in the catalytic pocket with spontaneous, transient formation of a H-bond between U-1(O2') and C75(N3). This H-bond would be required for C75 to act as the general base. Upon protonation in the precursor, C75H+ has a tendency to move towards its product location and establish a firm H-bonding network within the catalytic pocket. However, a C75H+(N3)-G1(O5') H-bond, which would be expected if C75 acted as a general acid catalyst, is not observed on the present simulation timescale. The adjacent loop L3 is relatively dynamic and may serve as a flexible structural element, possibly gated by the closing U20.G25 base-pair, to facilitate a conformational switch induced by a protonated C75H+. L3 also controls the electrostatic environment of the catalytic core, which in turn may modulate C75 base strength and metal ion binding. We find that a distant RNA tertiary interaction involving a protonated cytidine (C41) becomes unstable when left unprotonated, leading to disruptive conformational rearrangements adjacent to the catalytic core. A Na ion temporarily compensates for the loss of the protonated hydrogen bond, which is strikingly consistent with the experimentally observed synergy between low pH and high Na+ concentrations in mediating residual self-cleavage of the HDV ribozyme in the absence of divalents.

Metal ions and RNA folding: a highly charged topic with a dynamic future

Metal ions are required to stabilize RNA tertiary structure and to begin the folding process. How different metal ions enable RNAs to fold depends on the electrostatic potential of the RNA and correlated fluctuations in the positions of the ions themselves. Theoretical models, fluorescence spectroscopy, small angle scattering and structural biology reveal that metal ions alter the RNA dynamics and folding transition states. Specifically coordinated divalent metal ions mediate conformational rearrangements within ribozyme active sites.