Search

Homologous recombination forms and resolves an entangled DNA Holliday Junction (HJ) crucial for achieving genetic reshuffling and genome repair. To maintain genomic integrity, specialized resolvase enzymes cleave the entangled DNA into two discrete DNA molecules. However, it is unclear how two similar stacking isomers are distinguished, and how a cognate sequence is found and recognized to achieve accurate recombination. We here use single-molecule fluorescence observation and cluster analysis to examine how prototypic bacterial resolvase RuvC singles out two of the four HJ strands and achieves sequence-specific cleavage. We find that RuvC first exploits, then constrains the dynamics of intrinsic HJ isomer exchange at a sampled branch position to direct cleavage toward the catalytically competent HJ conformation and sequence, thus controlling recombination output at minimal energetic cost. Our model of rapid DNA scanning followed by ‘snap-locking’ of a cognate sequence is strikingly consistent with the conformational proofreading of other DNA-modifying enzymes.

The sensitive and accurate quantification of protein biomarkers plays important roles in clinical diagnostics and biomedical research. Sandwich ELISA and its variants accomplish the capture and detection of a target protein via two antibodies that tightly bind at least two distinct epitopes of the same antigen and have been the gold standard for sensitive protein quantitation for decades. However, existing antibody-based assays cannot distinguish between signal arising from specific binding to the protein of interest and nonspecific binding to assay surfaces or matrix components, resulting in significant background signal even in the absence of the analyte. As a result, they generally do not achieve single-molecule sensitivity, and they require two high-affinity antibodies as well as stringent washing to maximize sensitivity and reproducibility. Here, we show that surface capture with a high-affinity antibody combined with kinetic fingerprinting using a dynamically binding, low-affinity fluorescent antibody fragment differentiates between specific and nonspecific binding at the single-molecule level, permitting the direct, digital counting of single protein molecules with femtomolar-to-attomolar limits of detection (LODs). We apply this approach to four exemplary antigens spiked into serum, demonstrating LODs 55- to 383-fold lower than commercially available ELISA. As a real-world application, we establish that endogenous interleukin-6 (IL-6) can be quantified in 2-µL serum samples from chimeric antigen receptor T cell (CAR-T cell) therapy patients without washing away excess serum or detection probes, as is required in ELISA-based approaches. This kinetic fingerprinting thus exhibits great potential for the ultrasensitive, rapid, and streamlined detection of many clinically relevant proteins.

The ability to control transcription in a time-dependent manner in vitro promises numerous applications in molecular biology and nanotechnology. Here we demonstrate an approach that enables precise, independent control over the production of multiple RNA transcripts in vitro using single guide RNA (sgRNA)-directed transcription blockades by catalytically dead Streptococcus pyogenes CRISPR-Cas9 enzyme (dCas9). We show that when bound to a DNA template, the dCas9:sgRNA complex forms a robust blockade to transcription by RNA polymerases (RNAPs) from bacteriophages SP6, T3, and T7 (>99.5% efficiency), and a partial blockade to transcription by Escherichia coli RNAP (∼70% efficiency). We find that all three bacteriophage RNAPs dissociate from the DNA template upon encountering the dCas9 blockade, while E. coli RNAP stays bound for at least the 90-min duration of our experiments. The blockade maintains >95% efficiency when four mismatches are introduced into the 5′ end of the sgRNA target sequence. Notably, when using such a mismatched blockade, production of specific RNA species can be activated on demand by addition of a double-stranded competitor DNA perfectly matching the sgRNA. This strategy enables the independent production of multiple RNA species in a temporally controlled fashion from the same DNA template, demonstrating a new approach for transcription control.

DEAD-box proteins are nonprocessive RNA helicases that can function as RNA chaperones by coupling ATP binding and hydrolysis to structural reorganization of RNA. Here, Jarmoskaite et al. quantify the ATP utilization of an RNA chaperone during refolding of a misfolded ribozyme substrate. Strikingly, 100 ATP hydrolysis events are needed per successfully refolded ribozyme, suggesting that each round of unfolding requires ten ATP molecules, since 90% of substrate unfolding cycles only lead back to the kinetically favored misfolded state. This near-Sisyphean effort reveals a potentially conserved model for RNA reorganization by RNA chaperones.

RNA plays a fundamental, ubiquitous role as either substrate or functional component of many large cellular complexes – ‘molecular machines' – used to maintain and control the readout of most genetic information, a functional landscape that we are only beginning to understand. The cellular mechanisms for the spatiotemporal organization of the plethora of RNAs involved in gene expression are particularly poorly understood. Intracellular single-molecule fluorescence microscopy provides a powerful emerging tool for probing the pertinent mechanistic parameters that govern cellular RNA functions, including those of protein coding messenger RNAs (mRNAs). Progress has been hampered, however, by the scarcity of efficient high-yield methods to fluorescently label RNA molecules without the need to drastically increase their molecular weight through artificial appendages that may result in altered behavior. Herein, we employ T7 RNA polymerase to body label an RNA with a cyanine dye, as well as yeast poly(A) polymerase to strategically place multiple 2'-azido-modifications for subsequent fluorophore labeling either between the body and tail or randomly throughout the tail. Using a combination of biochemical and single-molecule fluorescence microscopy approaches, we demonstrate that both yeast poly(A) polymerase labeling strategies result in fully functional mRNA, whereas protein coding is severely diminished in the case of body labeling. This article is protected by copyright. All rights reserved.

The detection and quantification of biomarkers have numerous applications in biological research and medicine. The most widely used methods to detect nucleic acids require amplification via the polymerase chain reaction (PCR). However, errors arising from the imperfect copying fidelity of DNA polymerases, limited specificity of primers, and heat-induced damage reduce the specificity of PCR-based methods, particularly for single-nucleotide variants. Furthermore, not all analytes can be amplified efficiently. While amplification-free methods avoid these pitfalls, the specificity of most such methods is strictly constrained by probe binding thermodynamics, which for example hampers detection of rare somatic mutations. In contrast, single-molecule recognition through equilibrium Poisson sampling (SiMREPS) provides ultraspecific detection with single-molecule and single-nucleotide sensitivity by monitoring the repetitive interactions of a fluorescent probe with surface-immobilized targets. In this review, we discuss SiMREPS in comparison with other analytical approaches, and describe its utility in quantifying a range of nucleic acids and other analytes.

Living cells contain diverse biopolymers, creating a heterogeneous crowding environment, the impact of which on RNA folding is poorly understood. Here, we have used single-molecule fluorescence resonance energy transfer to monitor tertiary structure formation of the hairpin ribozyme as a model to probe the effects of polyethylene glycol and yeast cell extract as crowding agents. As expected, polyethylene glycol stabilizes the docked, catalytically active state of the ribozyme, in part through excluded volume effects; unexpectedly, we found evidence that it additionally displays soft, non-specific interactions with the ribozyme. Yeast extract has a profound effect on folding at protein concentrations 1000-fold lower than found intracellularly, suggesting the dominance of specific interactions over volume exclusion. Gel shift assays and affinity pull-down followed by mass spectrometry identified numerous non-canonical RNA-binding proteins that stabilize ribozyme folding; the apparent chaperoning activity of these ubiquitous proteins significantly compensates for the low-counterion environment of the cell.

In numerous Gram-positive bacteria, the glmS ribozyme or catalytic riboswitch regulates the expression of glucosamine-6-phosphate (GlcN6P) synthase via site-specific cleavage of its sugar-phosphate backbone in response to GlcN6P ligand binding. Biochemical data have suggested a crucial catalytic role for an active site guanine (G40 in Thermoanaerobacter tengcongensis, G33 in Bacillus anthracis). We used hybrid quantum chemical/molecular mechanical (QM/MM) calculations to probe the mechanism where G40 is deprotonated and acts as a general base. The calculations suggest that the deprotonated guanine G40(-) is sufficiently reactive to overcome the thermodynamic penalty arising from its rare protonation state, and thus is able to activate the A-1(2-OH) group toward nucleophilic attack on the adjacent backbone. Furthermore, deprotonation of A-1(2-OH) and nucleophilic attack are predicted to occur as separate steps, where activation of A-1(2-OH) precedes nucleophilic attack. Conversely, the transition state associated with the rate-determining step corresponds to concurrent nucleophilic attack and protonation of the G1(O5) leaving group by the ammonium moiety of the GlcN6P cofactor. Overall, our calculations help to explain the crucial roles of G40 (as a general base) and GlcN6P (as a general acid) during glmS ribozyme self-cleavage. In addition, we show that the QM/MM description of the glmS ribozyme self-cleavage reaction is significantly more sensitive to the size of the QM region and the quality of the QM-MM coupling than that of other small ribozymes. (c) 2015 Wiley Periodicals, Inc.

The DNA damage response (DDR) preserves genomic integrity. Small non-coding RNAs termed DDRNAs are generated at DNA double-strand breaks (DSBs) and are critical for DDR activation. Here we show that active DDRNAs specifically localize to their damaged homologous genomic sites in a transcription-dependent manner. Following DNA damage, RNA polymerase II (RNAPII) binds to the MRE11-RAD50-NBS1 complex, is recruited to DSBs and synthesizes damage-induced long non-coding RNAs (dilncRNAs) from and towards DNA ends. DilncRNAs act both as DDRNA precursors and by recruiting DDRNAs through RNA-RNA pairing. Together, dilncRNAs and DDRNAs fuel DDR focus formation and associate with 53BP1. Accordingly, inhibition of RNAPII prevents DDRNA recruitment, DDR activation and DNA repair. Antisense oligonucleotides matching dilncRNAs and DDRNAs impair site-specific DDR focus formation and DNA repair. We propose that DDR signalling sites, in addition to sharing a common pool of proteins, individually host a unique set of site-specific RNAs necessary for DDR activation.

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.

Structured RNAs such as ribozymes must fold into specific 3D structures to carry out their biological functions. While it is well-known that architectural features such as flexible junctions between helices help guide RNA tertiary folding, the mechanisms through which junctions influence folding remain poorly understood. We combine computational modeling with single molecule Förster resonance energy transfer (smFRET) and catalytic activity measurements to investigate the influence of junction design on the folding and function of the hairpin ribozyme. Coarse-grained simulations of a wide range of junction topologies indicate that differences in sterics and connectivity, independent of stacking, significantly affect tertiary folding and appear to largely explain previously observed variations in hairpin ribozyme stability. We further use our simulations to identify stabilizing modifications of non-optimal junction topologies, and experimentally validate that a three-way junction variant of the hairpin ribozyme can be stabilized by specific insertion of a short single-stranded linker. Combined, our multi-disciplinary study further reinforces that junction sterics and connectivity are important determinants of RNA folding, and demonstrates the potential of coarse-grained simulations as a tool for rationally tuning and optimizing RNA folding and function.

The hairpin ribozyme is a prominent member of small ribozymes since it does not require metal ions to achieve catalysis. Guanine 8 (G8) and adenine 38 (A38) have been identified as key participants in self-cleavage and -ligation. We have carried out hybrid quantum-mechanical/molecular mechanical (QM/MM) calculations to evaluate the energy along several putative reaction pathways. The error of our DFT description of the QM region was tested and shown to be ~1 kcal/mol. We find that self-cleavage of the hairpin ribozyme may follow several competing microscopic reaction mechanisms, all with calculated activation barriers in good agreement with those from experiment (20-21 kcal/mol). The initial nucleophilic attack of the A-1(2'-OH) group on the scissile phosphate is predicted to be rate-limiting in all these mechanisms. An unprotonated G8(-) (together with A38H(+)) yields a feasible activation barrier (20.4 kcal/mol). Proton transfer to a nonbridging phosphate oxygen also leads to feasible reaction pathways. Finally, our calculations consider thio-substitutions of one or both nonbridging oxygens of the scissile phosphate and predict that they have only a negligible effect on the reaction barrier, as observed experimentally.

RNA interference is a conserved gene regulatory mechanism employed by most eukaryotes as a key component of their innate immune response to viruses and retrotransposons. During viral infection, the RNase-III-type endonuclease Dicer cleaves viral double-stranded RNA into small interfering RNAs (siRNAs) 21-24 nucleotides in length and helps load them into the RNA-induced silencing complex (RISC) to guide the cleavage of complementary viral RNA. As a countermeasure, many viruses have evolved viral RNA silencing suppressors (RSS) that tightly, and presumably quantitatively, bind siRNAs to thwart RNA-interference-mediated degradation. Viral RSS proteins also act across kingdoms as potential immunosuppressors in gene therapeutic applications. Here we report fluorescence quenching and electrophoretic mobility shift assays that probe siRNA binding by the dimeric RSS p19 from Carnation Italian Ringspot Virus, as well as by human Dicer and RISC assembly complexes. We find that the siRNA:p19 interaction is readily reversible, characterized by rapid binding [(1.69 ± 0.07) × 10(8) M(-)(1) s(-1)] and marked dissociation (k(off)=0.062 ± 0.002 s(-1)). We also observe that p19 efficiently competes with recombinant Dicer and inhibits the formation of RISC-related assembly complexes found in human cell extract. Computational modeling based on these results provides evidence for the transient formation of a ternary complex between siRNA, human Dicer, and p19. An expanded model of RNA silencing indicates that multiple turnover by reversible binding of siRNAs potentiates the efficiency of the suppressor protein. Our predictive model is expected to be applicable to the dosing of p19 as a silencing suppressor in viral gene therapy.

The glmS catalytic riboswitch is part of the 5'-untranslated region of mRNAs encoding glucosamine-6-phosphate (GlcN6P) synthetase (glmS) in numerous gram-positive bacteria. Binding of the cofactor GlcN6P induces site-specific self-cleavage of the RNA. However, the detailed reaction mechanism as well as the protonation state of the glmS reactive form still remains elusive. To probe the dominant protonation states of key active site residues, we carried out explicit solvent molecular dynamic simulations involving various protonation states of three crucial active site moieties observed in the available crystal structures: (i) guanine G40 (following the Thermoanaerobacter tengcongensis numbering), (ii) the GlcN6P amino/ammonium group, and (iii) the GlcN6P phosphate moiety. We found that a deprotonated G40(-) seems incompatible with the observed glmS active site architecture. Our data suggest that the canonical form of G40 plays a structural role by stabilizing an in-line attack conformation of the cleavage site A-1(2'-OH) nucleophile, rather than a more direct chemical role. In addition, we observe weakened cofactor binding upon protonation of the GlcN6P phosphate moiety, which explains the experimentally observed increase in K(m) with decreasing pH. Finally, we discuss a possible role of cofactor binding and its interaction with the G65 and G1 purines in structural stabilization of the A-1(2'-OH) in-line attack conformation. On the basis of the identified dominant protonation state of the reaction precursor, we propose a hypothesis of the self-cleavage mechanism in which A-1(2'-OH) is activated as a nucleophile by the G1(pro-R(p)) nonbridging oxygen of the scissile phosphate, whereas the ammonium group of GlcN6P acts as the general acid protonating the G1(O5') leaving group.

The hairpin ribozyme is a prominent member of the group of small catalytic RNAs (RNA enzymes or ribozymes) because it does not require metal ions to achieve catalysis. Biochemical and structural data have implicated guanine 8 (G8) and adenine 38 (A38) as catalytic participants in cleavage and ligation catalyzed by the hairpin ribozyme, yet their exact role in catalysis remains disputed. To gain insight into dynamics in the active site of a minimal self-cleaving hairpin ribozyme, we have performed extensive classical, explicit-solvent molecular dynamics (MD) simulations on time scales of 50-150 ns. Starting from the available X-ray crystal structures, we investigated the structural impact of the protonation states of G8 and A38, and the inactivating A-1(2'-methoxy) substitution employed in crystallography. Our simulations reveal that a canonical G8 agrees well with the crystal structures while a deprotonated G8 profoundly distorts the active site. Thus MD simulations do not support a straightforward participation of the deprotonated G8 in catalysis. By comparison, the G8 enol tautomer is structurally well tolerated, causing only local rearrangements in the active site. Furthermore, a protonated A38H(+) is more consistent with the crystallography data than a canonical A38. The simulations thus support the notion that A38H(+) is the dominant form in the crystals, grown at pH 6. In most simulations, the canonical A38 departs from the scissile phosphate and substantially perturbs the structures of the active site and S-turn. Yet, we occasionally also observe formation of a stable A-1(2'-OH)...A38(N1) hydrogen bond, which documents the ability of the ribozyme to form this hydrogen bond, consistent with a potential role of A38 as general base catalyst. The presence of this hydrogen bond is, however, incompatible with the expected in-line attack angle necessary for self-cleavage, requiring a rapid transition of the deprotonated 2'-oxyanion to a position more favorable for in-line attack after proton transfer from A-1(2'-OH) to A38(N1). The simulations revealed a potential force field artifact, occasional but irreversible formation of "ladder-like", underwound A-RNA structure in one of the external helices. Although it does not affect the catalytic center of the hairpin ribozyme, further studies are under way to better assess possible influence of such force field behavior on long RNA simulations.

Despite numerous structural and biochemical investigations, the catalytic mechanism of hairpin ribozyme self-cleavage remains elusive. To gain insight into the coupling of active site dynamics with activity of this small catalytic RNA, we analyzed a total of ∼300 ns of molecular dynamics (MD) simulations. Our simulations predict improved global stability for an in vitro selected “gain of function” mutation, which is validated by native gel electrophoretic mobility shift assay. We observe that active site nucleobases and water molecules stabilize a geometry favorable to catalysis through a dynamic hydrogen bonding network. Simulations in which A38 is unprotonated show its N1 move into close proximity of the active site 2′-OH, indicating that A38 may act as a general base during cleavage, a role that has generally been discounted due to the longer distances observed in crystal structures involving inactivating substrate analogs. By contrast, simulations in which N1 of A38 is protonated place N1 in close proximity to the 5′-oxygen leaving group, which supports the proposal that A38 serves as a general acid. In analogy to protein enzymes, we discuss a plausible mechanism in which A38 acts bifunctionally and shuttles a proton directly from the 2′-OH to the 5′-oxygen. Furthermore, our simulations suggest an important role for protonation of N1 of A38 in promoting a favorable geometry similar to that observed in transition-state analog crystal structures, and support previously proposed roles of A38, G8, and long residency water molecules in transition-state stabilization.

Engineered “aptazymes” fuse in vitro selected aptamers with ribozymes to create allosteric enzymes as biosensing components and artificial gene regulatory switches through ligand-induced conformational rearrangement and activation. By contrast, activating ligand is employed as an enzymatic cofactor in the only known natural aptazyme, the glmS ribozyme, which is devoid of any detectable conformational rearrangements. To better understand this difference in biosensing strategy, we monitored by single molecule fluorescence resonance energy transfer (FRET) and 2-aminopurine (AP) fluorescence the global conformational dynamics and local base (un)stacking, respectively, of a prototypical drug-sensing aptazyme, built from a theophylline aptamer and the hammerhead ribozyme. Single molecule FRET reveals that a catalytically active state with distal Stems I and III of the hammerhead ribozyme is accessed both in the theophylline-bound and, if less frequently, in the ligand-free state. The resultant residual activity (leakage) in the absence of theophylline contributes to a limited dynamic range of the aptazyme. In addition, site-specific AP labeling shows that rapid local theophylline binding to the aptamer domain leads to only slow allosteric signal transduction into the ribozyme core. Our findings allow us to rationalize the suboptimal biosensing performance of the engineered compared to the natural aptazyme and to suggest improvement strategies. Our single molecule FRET approach also monitors in real time the previously elusive equilibrium docking dynamics of the hammerhead ribozyme between several inactive conformations and the active, long-lived, Y-shaped conformer.

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.

Many recently discovered RNA functions rely on highly complex multi-step conformational transitions that occur in response to an array of cellular signals. These dynamics accompany and guide, for example, RNA co-transcriptional folding, ligand sensing and signaling, site-specific catalysis in ribozymes, and the hierarchically ordered assembly of ribonucleoproteins. RNA dynamics are encoded by both the inherent properties of RNA structure, spanning many motional modes with a large range of amplitudes and timescales, and external trigger factors, ranging from proteins, nucleic acids, metal ions, metabolites, and vitamins to temperature and even directional RNA biosynthesis itself. Here, we review recent advances in our understanding of RNA dynamics as highlighted by biophysical tools.

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.

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 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 ∼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. © 2007 Wiley Periodicals, Inc. Biopolymers 85: 392–406, 2007. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com

The ability of RNA to catalyze chemical reactions was first demonstrated 25 years ago with the discovery that group I introns and RNase P function as RNA enzymes (ribozymes). Several additional ribozymes were subsequently identified, most notably the ribosome, followed by intense mechanistic studies. More recently, the introduction of single molecule tools has dissected the kinetic steps of several ribozymes in unprecedented detail and has revealed surprising heterogeneity not evident from ensemble approaches. Still, many fundamental questions of how RNA enzymes work at the molecular level remain unanswered. This review surveys the current status of our understanding of RNA catalysis at the single molecule level and discusses the existing challenges and opportunities in developing suitable assays. © 2007 Wiley Periodicals, Inc. Biopolymers 87: 302–316, 2007. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com

I have a dream. I have been hooked on it ever since 25 years ago when I prepared for my German diploma exam by reading a chapter in Lubert Stryer's Biochemistry textbook about catalytic RNAs, or ribozymes. Around this time the journal RNA was born as a manifestation of a growing RNA field and society. My PhD mentor Manfred Eigen organized regular winter seminars in Klosters, Switzerland; hearing Tom Cech talk at one of them about the beauty of RNA catalysis underscored just how relevant his and Sid Altman's discovery of non-protein catalysts was, earning them the Nobel Prize in Chemistry just a few years earlier in 1989. I decided to become part of the discoveries in the RNA World since they would reveal from where we came through deciphering the Origin of Life; as well as make us appreciate our roots and each other through the interconnectedness embodied in the Tree of Life. Yet throughout my graduate and postdoctoral studies of ribozymes I would feel that a piece of the puzzle was missing. It seemed rather strange to think of the primordial RNA World as a place of ubiquitous RNA functions that entirely would be lost in translation. I dreamed—like probably many of us—that somehow there had to be more to life than passive, repository nucleic acids and busy protein doers. And it turns out there is, in a big way. It is only a little over 60 years ago that Alfred Hershey and Martha Chase discovered that DNA represents the genetic material that stores all information necessary for life. A mere decade ago or so, the ∼3 billion base pairs of our genome were first sequenced in their entirety through the human genome project, at a cost of ∼$1 per base pair. The expectation for this breakthrough was that the identification of genetically encoded protein aberrations would realize a new era of personalized medical treatment for intractable human diseases such as cancer and Alzheimer's. Instead, it revealed the stunning truth that

X-ray crystallography can provide important insights into the structure of RNA enzymes (ribozymes). However, the details of a ribozyme's active site architecture are often altered by the inactivating chemical modifications necessary to inhibit self-cleavage. Molecular dynamics (MD) simulations are able to complement crystallographic data and model the conformation of the ribozyme's active site in its native form. However, the performance of MD simulations is driven by the quality of the force field used. Force fields are primarily parametrized and tested for a description of canonical structures and thus may be less accurate for noncanonical RNA elements, including ribozyme catalytic cores. Here, we show that our recent reparametrization of ε/ζ torsions significantly improves the description of the hairpin ribozyme's scissile phosphate conformational behavior. In addition, we find that an imbalance in the force field description of the nonbonded interactions of the ribose 2'-OH contributes to the conformational behavior observed for the scissile phosphate in the presence of a deprotonated G8(-). On the basis of the new force field, we obtain a reactive conformation for the hairpin ribozyme active site that is consistent with the most recent mechanistic and structural data.

Transforming growth factor-β (TGF-β) signaling regulates cell proliferation and differentiation, which contributes to development and disease. Upon binding TGF-β, the type I receptor (TGFBRI) binds TGFBRII, leading to the activation of the transcription factors SMAD2 and SMAD3. Using an RNA interference screen of the human kinome and a live-cell reporter for TGFBR activity, we identified the kinase BUB1 (budding uninhibited by benzimidazoles-1) as a key mediator of TGF-β signaling. BUB1 interacted with TGFBRI in the presence of TGF-β and promoted the heterodimerization of TGFBRI and TGFBRII. Additionally, BUB1 interacted with TGFBRII, suggesting the formation of a ternary complex. Knocking down BUB1 prevented the recruitment of SMAD3 to the receptor complex, the phosphorylation of SMAD2 and SMAD3 and their interaction with SMAD4, SMAD-dependent transcription, and TGF-β-mediated changes in cellular phenotype including epithelial-mesenchymal transition (EMT), migration, and invasion. Knockdown of BUB1 also impaired noncanonical TGF-β signaling mediated by the kinases AKT and p38 MAPK (mitogen-activated protein kinase). The ability of BUB1 to promote TGF-β signaling depended on the kinase activity of BUB1. A small-molecule inhibitor of the kinase activity of BUB1 (2OH-BNPP1) and a kinase-deficient mutant of BUB1 suppressed TGF-β signaling and formation of the ternary complex in various normal and cancer cell lines. 2OH-BNPP1 administration to mice bearing lung carcinoma xenografts reduced the amount of phosphorylated SMAD2 in tumor tissue. These findings indicated that BUB1 functions as a kinase in the TGF-β pathway in a role beyond its established function in cell cycle regulation and chromosome cohesion.

Ribozymes are catalytically competent examples of highly structured noncoding RNAs, which are ubiquitous in the processing and regulation of genetic information. Combining explicit-solvent molecular dynamics simulation and single molecule fluorescence spectroscopy approaches, we find that a ribozyme from a subviral plant pathogen exhibits a coupled hydrogen bonding network that communicates dynamic structural rearrangements throughout the catalytic core in response to site-specific chemical modification. Trapped long-residency water molecules are critical for this network and only occasionally exchange with bulk solvent as they pass through a breathing interdomain base stack. These highly structured water molecules line up in a string that may potentially also be involved in specific base catalysis. Our observations suggest important, still underappreciated roles for specifically bound water molecules in the structural dynamics and function of noncoding RNAs.

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.

The hepatitis delta virus (HDV) ribozyme is a catalytic RNA motif embedded in the human pathogenic HDV RNA. It catalyzes self-cleavage of its sugar-phosphate backbone with direct participation of the active site cytosine C75. Biochemical and structural data support a general acid role of C75. Here, we used hybrid quantum mechanical/molecular mechanical (QM/MM) calculations to probe the reaction mechanism and changes in Gibbs energy along the ribozyme's reaction pathway with an N3-protonated C75H(+) in the active site, which acts as the general acid, and a partially hydrated Mg(2+) ion with one deprotonated, inner-shell coordinated water molecule that acts as the general base. We followed eight reaction paths with a distinct position and coordination of the catalytically important active site Mg(2+) ion. For six of them, we observed feasible activation barriers ranging from 14.2 to 21.9 kcal mol(-1), indicating that the specific position of the Mg(2+) ion in the active site is predicted to strongly affect the kinetics of self-cleavage. The deprotonation of the U-1(2'-OH) nucleophile and the nucleophilic attack of the resulting U-1(2'-O(-)) on the scissile phosphodiester are found to be separate steps, as deprotonation precedes the nucleophilic attack. This sequential mechanism of the HDV ribozyme differs from the concerted nucleophilic activation and attack suggested for the hairpin ribozyme. We estimate the pKa of the U-1(2'-OH) group to range from 8.8 to 11.2, suggesting that it is lowered by several units from that of a free ribose, comparable to and most likely smaller than the pKa of the solvated active site Mg(2+) ion. Our results thus support the notion that the structure of the HDV ribozyme, and particularly the positioning of the active site Mg(2+) ion, facilitate deprotonation and activation of the 2'-OH nucleophile.

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.

The hairpin ribozyme is a minimalist paradigm for studying RNA folding and function. In this enzyme, two domains dock by induced fit to form a catalytic core that mediates a specific backbone cleavage reaction. Here, we have fully dissected its reversible reaction pathway, which comprises two structural transitions (docking/undocking) and a chemistry step (cleavage/ligation), by applying a combination of single-molecule fluorescence resonance energy transfer (FRET) assays, ensemble cleavage assays, and kinetic simulations. This has allowed us to quantify the effects that modifications of essential functional groups remote from the site of catalysis have on the individual rate constants. We find that all ribozyme variants show similar fractionations into effectively noninterchanging molecule subpopulations of distinct undocking rate constants. This leads to heterogeneous cleavage activity as commonly observed for RNA enzymes. A modification at the domain junction additionally leads to heterogeneous docking. Surprisingly, most modifications not only affect docking/undocking but also significantly impact the internal chemistry rate constants over a substantial distance from the site of catalysis. We propose that a network of coupled molecular motions connects distant parts of the RNA with its reaction site, which suggests a previously undescribed analogy between RNA and protein enzymes. Our findings also have broad implications for applications such as the action of drugs and ligands distal to the active site or the engineering of allostery into RNA.

The hepatitis delta virus (HDV), an infectious human pathogen affecting millions of people worldwide, leads to intensified disease symptoms, including progression to liver cirrhosis upon coinfection with its helper virus, HBV. Both the circular RNA genome of HDV and its complementary antigenome contain a common cis-cleaving catalytic RNA motif, the HDV ribozyme, which plays a crucial role in viral replication. Previously, the crystal structure of the product form of the cis-acting genomic HDV ribozyme has been determined, and the precursor form has been suggested to be structurally similar. In contrast, solution studies by fluorescence resonance energy transfer (FRET) on a trans-cleaving form of the ribozyme have shown significant global conformational changes upon catalysis, while 2-aminopurine (AP) fluorescence assays have detected concomitant local conformational changes in the catalytic core. Here, we augment these studies by using terbium(III) to probe the structure of the trans-acting HDV ribozyme at nucleotide resolution. We observe significant structural differences between the precursor and product forms, especially in the P1.1 helix and the trefoil turn in the single-stranded region connecting P4 and P2 (termed J4/2) of the catalytic core. We show, using terbium(III) footprinting and sensitized luminescence spectroscopy as well as steady-state, time-resolved, and gel-mobility FRET assays on a systematic set of substrates, that the substrate sequence immediately 5 ' to the cleavage site significantly modulates these local as well as resultant global structural differences. Our results suggest a structural basis for the previously observed impact of the 5' substrate sequence on catalytic activity.

Recently, Breaker and coworkers engineered hammerhead ribozymes that rearrange from a catalytically inactive to an active conformation upon allosteric binding of a specific ligand. To monitor cleavage activity in real time, we have coupled a donor-acceptor fluorophore pair to the termini of the substrate RNA of such a hammerhead ribozyme, modified to cleave in trans in the presence of the bronchodilator theophylline. In the intact substrate, the fluorophores interact by fluorescence resonance energy transfer (FRET). The specific FRET signal breaks down as the effector ligand binds, the substrate is cleaved, and the products dissociate, with a rate constant dependent on the concentration of the ligand. Our biosensor cleaves substrate at 0.46 min(-1) in 1 mM theophylline and 0.04 min(-1) without effector, and discriminates against caffeine, a structural relative of theophylline. We have measured the theophylline-dependence profile of this biosensor, showing that concentrations as low as 1 microM can be distinguished from background. To probe the mechanism of allosteric regulation, a single nucleotide in the communication domain between the catalytic and ligand-binding domains was mutated to destabilize the inactive conformation of the ribozyme. As predicted, this mutant shows the same activity (0.3 min(-1)) in the presence and absence of theophylline. Additionally, time-resolved FRET measurements on the biosensor ribozyme in complex with a noncleavable substrate analog reveal no significant changes in fluorophore distance distribution upon binding of effector.

RNA performs a multitude of essential cellular functions involving the maintenance, transfer, and processing of genetic information. The reason probably is twofold: (a) Life started as a prebiotic RNA World, in which RNA served as the genetic information carrier and catalyzed all chemical reactions required for its proliferation and (b) some of the RNA World functions were conserved throughout evolution because neither DNA nor protein is as adept in fulfilling them. A particular advantage of RNA is its high propensity to form alternative structures as required in subsequent steps of a reaction pathway. Here I describe fluorescence resonance energy transfer (FRET) as a method to monitor a crucial conformational transition on the reaction pathway of the hairpin ribozyme, a small catalytic RNA motif from a self-replicating plant virus satellite RNA and well-studied paradigm of RNA folding. Steady-state FRET measurements in solution allow one to measure the kinetics and requirements of docking of its two independently folding domains; time-resolved FRET reveals the relative thermodynamic stability of the undocked (extended, inactive) and docked (active) ribozyme conformations; while single-molecule FRET experiments will highlight the dynamics of RNA at the individual molecule level. Similar domain docking events are expected to be at the heart of many biological functions of RNA, and the described FRET techniques promise to be adaptable to most of the involved RNA systems.

Global structure and flexibility of three different hairpin ribozyme constructs have been analyzed by measuring their electric dichroism decay in various buffers at temperatures between 2 and 30 degrees C. The hairpin ribozyme is characterized by two independently folding domains A and B that are connected through a hinge and have to interact to enable catalysis. The analyzed constructs feature extended terminal helices 1 and 4 with 27 and 25 bp, respectively, to increase the sensitivity of the molecular rotational diffusion time constants with respect to the interdomain bending angle. Constructs HP1 and HP2 cannot cleave because of a G+1A change at the 3'-side of the cleavage site; in HP1 the helices 2 and 3 that flank the hinge form a continuous double helical segment; in HP2 and HP3, a six nucleotide bulge confers flexibility to the expected bending site; HP3 is a cleavable form of HP2 with a G+1-base. For comparison, a standard RNA double helix with 72 bp was included in our analysis. The dichroism decay curves of the hairpin constructs after pulses of low electric field strengths can be fitted to single exponentials taus, whereas the curves after pulses of high field strengths require two exponentials. In all cases, time constants increase with RNA concentration, indicating intermolecular interactions. Extrapolation of the tausvalues measured in standard buffer (50 mM Tris (pH 7.5) and 12 mM MgCl2) to zero RNA concentration provide values of 112, 93, and 73 ns for HP1, HP2 and HP3, respectively, at 30 degrees C, indicating increasingly compact structures. The 72 bp RNA reference under corresponding conditions did not show a dependence of its decay time constant on the RNA concentration nor on the field strength; its time constant is 175 ns (standard buffer, 30 degrees C). The observation of two relaxation processes for the hairpin constructs at high field strengths indicates stretching to a more elongated state; the fast process with a time constant of the order of 50 ns is assigned to reversion of stretching, the slow process to overall rotation. The overall rotational time of the stretched state at 20 degrees C is close to that for a completely stretched rigid state; at 30 degrees C the experimental values are around 70 % of that expected for a completely stretched rigid state, indicating flexibility and/or residual bending. Bead models were constructed to simulate dichroism decay curves. The time constants observed for the 72 bp RNA are as expected for a rigid rod with a rise of 2.8 A per base-pair. Based on this rise per base-pair for models of a V and a Y-shape, we estimate average bending angles of 80(+/-20) degrees and 105 (+/-25) degrees, respectively, for the catalytically active hairpin ribozyme HP3. The energy required for stretching is of the order of the thermal energy.

The equilibrium distributions of hairpin ribozyme conformational isomers have been examined by time-resolved fluorescence resonance energy transfer. Ribozymes partition between active (docked) and inactive (extended) conformers, characterized by unique interdomain distance distributions, which define differences in folding free energy. The active tertiary structure is stabilized both by specific interactions between the catalytic and the substrate-binding domains and by the structure of the intervening helical junction. Under physiological conditions, the docking equilibrium of the natural four-way junction dramatically favors the active conformer, while those of a three-way and the two-way junction used in gene therapy applications favor the inactive conformer.

The complex between the hairpin ribozyme and its substrate consists of two domains that must interact in order to form a catalytic complex, yet experimental evidence concerning the points of interaction between the two domains has been lacking. Here, we report the use of hydroxyl radical footprinting to define the interface between the two domains. Cations that support very efficient ribozyme catalysis (magnesium and cobalt(III) hexammine) lead to the formation of a docked complex that features several regions of protection, indicating a solvent-inaccessible core within the tertiary structure of the complex. Cations that are suboptimal in cleavage reactions do not produce complexes with regions of reduced solvent accessibility. Nucleotides encompassing the substrate cleavage site (c-2, a-1, g+1, and u+2) are strongly protected, suggesting their internalization into the catalytic core. Four distinct segments of the ribozyme are protected, including G11-A14, C25-C27, A38, and U42-A43. Protection of these sites is eliminated when g+1, an essential base at the cleavage site, is replaced by A. In addition, mutations which are known to decrease the fraction of docked complexes decrease or eliminate formation of a solvent-inaccessible core. Taken together, these observations demonstrate that we have identified the catalytic core of the active hairpin ribozyme-substrate complex.

Binding of small RNAs by the RNA-dependent RNA polymerase of coliphage Qβ was studied utilizing a fluorometric assay. A DNA oligonucleotide probe of sequence 5′-d(TTTTTCC) was 5′-end-labeled with pyrene. In this construct, the proximal thymine residues efficiently quench the fluorophore emission in solution. Upon stoichiometric binding of one probe per polymerase molecule, the pyrene steady-state fluorescence increases by two orders of magnitude, the fluorescence anisotropy increases, and a long fluorescence lifetime component of 140 ns appears. With addition of replicable RNA, steady-state fluorescence decreases in a concentration dependent manner and the long lifetime component is lost. This observation most likely reflects displacement of the pyrene-labeled probe from the proposed nucleic acid binding site II of Qβ replicase. The effect was utilized to access binding affinities of different RNAs to this site in a reverse titration assay format. In 10 mM sodium phosphate (pH 7.0), 100 mM NaCl, at 16°C, equilibrium dissociation constants for different template midi- and minivariant RNAs were calculated to be in the nanomolar range. In general, the minus and plus strands, concomitantly synthesized by Qβ replicase during replication, exhibited discriminative affinities, while their hybrid bound less efficiently than either of the single strands. Different non-replicable tRNAs also bound to the polymerase with comparable dissociation constants. By titration with DNA homo-oligonucleotides it was shown that the probed site on Qβ replicase does not require a 2′ hydroxyl group for binding nucleic acids, but recognizes pyrimidine residues. Its interaction with thymine is lost in an A·T base-pair, while that with cytosine is retained after Watson-Crick base-pairing. These findings can explain the affinities of RNA-Qβ replicase interactions reported here and in earlier investigations. The sensitivity of the described fluorometric assay allows detection of RNA amplification by Qβ replicase in real-time.

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

Historically, all small catalytic RNAs have been shown to undergo global conformational changes upon phosphodiester cleavage. However, the most recent of numerous hepatitis delta virus (HDV) ribozyme crystal structures has challenged this trend, as this crystal structure suggests that the precursor structure is already product-like in conformation. To further investigate this unusual observation, we have extensively characterized the solution behavior of several three-stranded versions of the HDV ribozyme from the recent crystal structure. Fluorescence gel shift assays show that varying lengths of the 5'overhang sequence adjacent to the active site result in the same degree of cleavage, whereas noncleavable substrates exhibit significantly more heterogeneity. Complementary steady state and time-resolved FRET assays demonstrated that the length of the 5'overhang sequence adjacent to the cleavage site affects the rates of conformational change upon substrate binding and cleavage. Molecular dynamics (MD) simulations were also performed to gain insight into the atomic behavior and catalytic relevance of the HDV ribozyme from the Chen et al crystal structure. These simulations suggested that the dU-1dG1dG2 motif used in the crystal does not result in a catalytically fit ribozyme compared to an all-ribose construct. Furthermore, altered active site conditions also result in lowered catalytic fitness. These simulations suggested that catalytic fitness is greatly disrupted by deprotonation of C75, supporting the hypothesis of the role of C75 as a general acid. Simulations also showed that the magnesium ion resolved near the scissile phosphate results in favorable catalytic geometry compared to simulations neutralized with sodium. Our experimental results demonstrate that, despite previously published results, all forms of the HDV ribozyme undergo significant global conformational changes upon self-cleavage, and our simulations show that C75 is poised to act as a general acid during cleavage.

Evolution comprises cycles of amplification, mutagenesis and selection. To study evolutionary phenomena, isothermal strand displacement amplification (SDA) of double-stranded DNA as an in vitro model for rolling-circle replication of viruses has been coupled to a positive selection procedure. First, two subsequent amplification reactions utilizing exo- Klenow polymerase were performed under direct observation using the fluorescent dye thiazole orange. Under the chosen conditions, the mutation rate was 1.5 x 10(-3) and 0.4 x 10(-3) for base substitutions and deletions, respectively. Then, a 16mer oligodeoxynucleotide with an acridine moiety coupled to its 5' end was used to select for double strands that retained their ability to form a triple-helix with the oligodeoxynucleotide. Conditions for triple-helix formation were chosen such that only 10 to 40% of the SDA products were allowed to bind the third strand. Non-denaturing polyacrylamide gel electrophoresis was used to separate triple-helices from unmodified double strands, and only triplex strands were used to initiate a new round of error-prone amplification and selection. Nine such rounds with about 270 molecular generations were performed. The final mutant spectrum was characterized and compared with those of amplification reactions without additional selection pressure. While without selection pressure base substitutions and deletions throughout the initial wild-type rapidly produce a diverse mutant distribution, the consensus after nine selection rounds clearly shows two mutational hotspot positions. Using gel shift assays and a newly developed non-radioactive DNase I footprinting technique, it could be shown that both the initial wild-type and the final consensus do not differ significantly in their triplex formation ability. As opposed to this, they do show different amplification efficiencies. The final consensus sequence is amplified with the highest rate in the exponential reaction phase, while the most abundant clone, which is characterized by two additional point deletions, is the sequence with the highest amplification rate in the linear growth phase.

The hepatitis delta virus (HDV) is the only known human pathogen to contain a self-cleaving catalytic RNA motif (ribozyme) in its genome. The native structure of the HDV ribozyme consists of five helices (P1 - P4 and P1.1) that come together to form a double-nested pseudoknot. The active site with the catalytically involved C75 residue lies at the heart of this complex structure. The recent emergence of a product-like crystal structure of the precursor form of the HDV ribozyme1 has cast doubt on previous structural and biochemical evidence that a significant conformational change occurs along the reaction trajectory from precursor to product. We here are using MD simulations to investigate the impact that key structural motifs as well as chemical and structural modifications necessary to crystallize the ribozyme have on conformational dynamics of the active site in this recent crystal structure. Our simulations specifically probe the effects of varying base protonation, base identity in key structural motifs, and alternate conformations of the scissile phosphate on catalytic fitness2. Analysis of our MD simulations indicates that the precise position of C75, and thus catalytic fitness, is significantly affected by the modifications introduced into the recent crystal structure. We also investigate in detail the effects of force field choice on the observed ribozyme structural dynamics, by comparing several force field variants, including our recent reparametrization of the chi profile of the Cornell et al (AMBER) force field 3.

Strand displacement amplification is an isothermal DNA amplification reaction based on a restriction endonuclease nicking its recognition site and a polymerase extending the nick at its 3' end, displacing the downstream strand. The reaction resembles rolling-circle replication of single-stranded phages and small plasmids. The displaced sense strand serves as target for an antisense reaction and vice versa, resulting in exponential growth and the autocatalytic nature of this in vitro reaction as long as the template is the limiting agent. We describe the optimization of strand displacement amplification for in vitro evolution experiments under serial transfer conditions. The reaction was followed and controlled by use of the fluorescent dye thiazole orange binding to the amplified DNA. We were able to maintain exponential growth conditions with a doubling time of 3.0 min throughout 100 transfers or approximately 350 molecular generations by using an automatic handling device. Homology of in vitro amplification with rolling-circle replication was mirrored by the occurring evolutionary processes. Deletion events most likely caused by a slipped mispairing mechanism as postulated for in vivo replication took place. Under our conditions, the mutation rate was high and a molecular quasi-species formed with a mutant lacking internal hairpin formation ability and thus outgrowing all other species under dGTP/dCTP deficiency.

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.

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

The hepatitis delta virus (HDV) is the only known human pathogen to contain a catalytic RNA motif (ribozyme) in its genome. The native structure of the HDV ribozyme consists of five helices (P1 - P4 and P1.1), which come together to form a double-nested pseudoknot. Within this complicated ribozyme are several landmarks of structural importance, and alterations of these landmarks often lead to decreased ribozymal activity. We have used extended explicit solvent molecular dynamics simulations to investigate structural dynamics of HDV ribozyme. We plan on including effects of base protonation and base substitutions, using the available X-ray structures and their modifications as the starting structures. Specific attention has been paid to nucleotides in the catalytic pocket and the flexible L3 region. We have in detail investigated the effect of the force field choice on the ribozyme structural dynamics, by comparing several simulation force field variants, including our recent reparametrization of the chi profile of the Cornell et al (AMBER) force field.

The two domains of the hairpin ribozyme-substrate complex, usually depicted as straight structural elements, must interact with one another in order to form an active conformation. Little is known about the internal geometry of the individual domains in an active docked complex. Using various crosslinking and structural approaches in conjunction with molecular modeling (constraint-satisfaction program MC-SYM), we have investigated the conformation of the substrate-binding domain in the context of the active docked ribozyme-substrate complex. The model generated by MC-SYM showed that the domain is not straight but adopts a bent conformation (D-shaped) in the docked state of the ribozyme, indicating that the two helices bounding the internal loop are closer than was previously assumed. This arrangement rationalizes the observed ability of hairpin ribozymes with a circularized substrate-binding strand to cleave a circular substrate, and provides essential information concerning the organization of the substrate in the active conformation. The internal geometry of the substrate-binding strand places G8 of the substrate-binding strand near the cleavage site, which has allowed us to predict the crucial role played by this nucleotide in the reaction chemistry., discussions, David Pecchia for the oligonucleotide synthesis. This work is supported by grants from the National Institutes of Health (AI 44186 to J.M.B. and GM52168 to G.D.G.) and from the Medical Research council of Canada (MT 14504 to F.M.). R.P. was supported by a postdoctoral fellowship from the Medical Research Council of Canada.

Metal ions interact with RNA to enhance folding, stabilize structure, and, in some cases, facilitate catalysis. Assigning functional roles to specifically bound metal ions presents a major challenge in analyzing the catalytic mechanisms of ribozymes. Bacillus subtilis ribonuclease P (RNase P), composed of a catalytically active RNA subunit (PRNA) and a small protein subunit (P protein), catalyzes the 5'-end maturation of precursor tRNAs (pre-tRNAs). Inner-sphere coordination of divalent metal ions to PRNA is essential for catalytic activity but not for the formation of the RNase P x pre-tRNA (enzyme-substrate, ES) complex. Previous studies have demonstrated that this ES complex undergoes an essential conformational change (to the ES* conformer) before the cleavage step. Here, we show that the ES* conformer is stabilized by a high-affinity divalent cation capable of inner-sphere coordination, such as Ca(II) or Mg(II). Additionally, a second, lower-affinity Mg(II) activates cleavage catalyzed by RNase P. Structural changes that occur upon binding Ca(II) to the ES complex were determined by time-resolved Forster resonance energy transfer measurements of the distances between donor-acceptor fluorophores introduced at specific locations on the P protein and pre-tRNA 5' leader. These data demonstrate that the 5' leader of pre-tRNA moves 4 to 6 A closer to the PRNA x P protein interface during the ES-to-ES* transition and suggest that the metal-dependent conformational change reorganizes the bound substrate in the active site to form a catalytically competent ES* complex.

Hybrid QM/MM methods combine the rigor of quantum mechanical (QM) calculations with the low computational cost of empirical molecular mechanical (MM) treatment allowing to capture dynamic properties to probe critical atomistic details of enzyme reactions. Catalysis by RNA enzymes (ribozymes) has only recently begun to be addressed with QM/MM approaches and is thus still a field under development. This review surveys methodology as well as recent advances in QM/MM applications to RNA mechanisms, including those of the HDV, hairpin, and hammerhead ribozymes, as well as the ribosome. We compare and correlate QM/MM results with those from QM and/or molecular dynamics (MD) simulations, and discuss scope and limitations with a critical eye on current shortcomings in available methodologies and computer resources. We thus hope to foster mutual appreciation and facilitate collaboration between experimentalists and theorists to jointly advance our understanding of RNA catalysis at an atomistic level.

Non-coding RNAs of complex tertiary structure are involved in numerous aspects of the replication and processing of genetic information in many organisms; however, an understanding of the complex relationship between their structural dynamics and function is only slowly emerging. The Neurospora Varkud Satellite (VS) ribozyme provides a model system to address this relationship. First, it adopts a tertiary structure assembled from common elements, a kissing loop and two three-way junctions. Second, catalytic activity of the ribozyme is essential for replication of VS RNA in vivo and can be readily assayed in vitro. Here we exploit single molecule FRET to show that the VS ribozyme exhibits previously unobserved dynamic and heterogeneous hierarchical folding into an active structure. Readily reversible kissing loop formation combined with slow cleavage of the upstream substrate helix suggests a model whereby the structural dynamics of the VS ribozyme favor cleavage of the substrate downstream of the ribozyme core instead. This preference is expected to facilitate processing of the multimeric RNA replication intermediate into circular VS RNA, which is the predominant form observed in vivo.