Conformational heterogeneity at position U37 of an all-RNA hairpin ribozyme with implications for metal binding and the catalytic structure of the S-turn

Two RNA Folds from One Sequence: A Ribozyme with Versatile Substrate Processing Abilities

We report the design of a single RNA sequence capable of adopting one of two ribozyme folds and catalyzing the cleavage and/or ligation of the respective substrates. The RNA is able to change its conformation in response to its environment, hence it is called chameleon ribozyme (CHR). Efficient RNA cleavage of two different substrates as well as RNA ligation by CHR is demonstrated in separate experiments and in a one pot reaction. Our study shows that sequence variants of the hairpin ribozyme intersect with the hammerhead ribozyme and that rather short RNA molecules can have comprehensive conformational flexibility, which is an important feature for the emergence of new functional folds in early evolution.

Full-Length NAD+-I Riboswitches Bind a Single Cofactor but Cannot Discriminate against Adenosine Triphosphate

Bacterial riboswitches are structured RNAs that bind small metabolites to control downstream gene expression. Two riboswitch classes have been reported to sense nicotinamide adenine dinucleotide (NAD+), which plays a key redox role in cellular metabolism. The NAD+-I (class I) riboswitch stands out because it comprises two homologous, tandemly arranged domains. However, previous studies examined the isolated domains rather than the full-length riboswitch. Crystallography and ligand binding analyses led to the hypothesis that each domain senses NAD+ but with disparate equilibrium binding constants (KD) of 127 μM (domain I) and 3.4 mM (domain II). Here, we analyzed individual domains and the full-length riboswitch by isothermal titration calorimetry to quantify the cofactor affinity and specificity. Domain I senses NAD+ with a KD of 24.6 ± 8.4 μM but with a reduced ligand-to-receptor stoichiometry, consistent with nonproductive domain self-association observed by gel-filtration chromatography; domain II revealed no detectable binding. By contrast, the full-length riboswitch binds a single NAD+ with a KD of 31.5 ± 1.5 μM; dinucleotides NADH and AP2-ribavirin also bind with one-to-one stoichiometry. Unexpectedly, the full-length riboswitch also binds a single ATP equivalent (KD = 11.0 ± 3.5 μM). The affinity trend of the full-length riboswitch is ADP = ATP > NAD+ = AP2-ribavirin > NADH. Although our results support riboswitch sensing of a single NAD+ at concentrations significantly below the intracellular levels of this cofactor, our findings do not support the level of specificity expected for a riboswitch that exclusively senses NAD+. Gene regulatory implications and future challenges are discussed.

Biological Catalysis and Information Storage Have Relied on N-Glycosyl Derivatives of β-D-Ribofuranose since the Origins of Life

Most naturally occurring nucleotides and nucleosides are N-glycosyl derivatives of β-d-ribose. These N-ribosides are involved in most metabolic processes that occur in cells. They are essential components of nucleic acids, forming the basis for genetic information storage and flow. Moreover, these compounds are involved in numerous catalytic processes, including chemical energy production and storage, in which they serve as cofactors or coribozymes. From a chemical point of view, the overall structure of nucleotides and nucleosides is very similar and simple. However, their unique chemical and structural features render these compounds versatile building blocks that are crucial for life processes in all known organisms. Notably, the universal function of these compounds in encoding genetic information and cellular catalysis strongly suggests their essential role in the origins of life. In this review, we summarize major issues related to the role of N-ribosides in biological systems, especially in the context of the origin of life and its further evolution, through the RNA-based World(s), toward the life we observe today. We also discuss possible reasons why life has arisen from derivatives of β-d-ribofuranose instead of compounds based on other sugar moieties.

Who stole the proton? Suspect general base guanine found with a smoking gun in the pistol ribozyme

The pistol ribozyme (Psr) is one among the most recently discovered classes of small nucleolytic ribozymes that catalyze site-specific RNA self-cleavage through 2'-O-transphosphorylation. The Psr contains a conserved guanine (G40) that in crystal structures is in a position suggesting it plays the role of the general base to abstract a proton from the nucleophile to activate the reaction. Although some functional data is consistent with this mechanistic role, a notable exception is 2-aminopurine (2AP) substitution which has no effect on the rate, unlike similar substitutions across other so-called "G + M" and "G + A" ribozyme classes. Herein we postulate that an alternate conserved guanine, G42, is the primary general base, and provide evidence from molecular simulations that the active site of Psr can undergo local refolding into a structure that is consistent with the common "L-platform/L-scaffold" architecture identified in G + M and G + A ribozyme classes with Psr currently the notable exception. We summarize the key currently available experimental data and present new classical and combined quantum mechanical/molecular mechanical simulation results that collectively suggest a new hypothesis. We hypothesize that there are two available catalytic pathways supported by different conformational states connected by a local refolding of the active site: (1) a primary pathway with an active site architecture aligned with the L-platform/L-scaffold framework where G42 acts as a general base, and (2) a secondary pathway with the crystallographic active site architecture where G40 acts as a general base. We go on to make several experimentally testable predictions, and suggest specific experiments that would ultimately bring closure to the mystery as to "who stole the proton in the pistol ribozyme?".

Investigation of the pKa of the Nucleophilic O2' of the Hairpin Ribozyme

Small ribozymes cleave their RNA phosphodiester backbone by catalyzing a transphosphorylation reaction wherein a specific O2' functions as the nucleophile. While deprotonation of this alcohol through its acidification would increase its nucleophilicity, little is known about the pK_a of this O2' in small ribozymes, in part because high pK_a's are not readily accessible experimentally. Herein, we turn to molecular dynamics to calculate the pK_a of the nucleophilic O2' in the hairpin ribozyme and to study interactions within the active site that may impact its value. We estimate the pK_a of the nucleophilic O2' in the wild-type hairpin ribozyme to be 18.5 ± 0.8, which is higher than the reference compound, and identify a correlation between proper positioning of the O2' for nucleophilic attack and elevation of its pK_a. We find that monovalent ions may play a role in depression of the O2' pK_a, while the exocyclic amine appears to be important for organizing the ribozyme active site. Overall, this study suggests that the pK_a of the O2' is raised in the ground state and lowers during the course of the reaction owing to positioning and metal ion interactions.

Context-sensitivity of isosteric substitutions of non-Watson-Crick basepairs in recurrent RNA 3D motifs

Sequence variation in a widespread, recurrent, structured RNA 3D motif, the Sarcin/Ricin (S/R), was studied to address three related questions: First, how do the stabilities of structured RNA 3D motifs, composed of non-Watson–Crick (non-WC) basepairs, compare to WC-paired helices of similar length and sequence? Second, what are the effects on the stabilities of such motifs of isosteric and non-isosteric base substitutions in the non-WC pairs? And third, is there selection for particular base combinations in non-WC basepairs, depending on the temperature regime to which an organism adapts? A survey of large and small subunit rRNAs from organisms adapted to different temperatures revealed the presence of systematic sequence variations at many non-WC paired sites of S/R motifs. UV melting analysis and enzymatic digestion assays of oligonucleotides containing the motif suggest that more stable motifs tend to be more rigid. We further found that the base substitutions at non-Watson–Crick pairing sites can significantly affect the thermodynamic stabilities of S/R motifs and these effects are highly context specific indicating the importance of base-stacking and base-phosphate interactions on motif stability. This study highlights the significance of non-canonical base pairs and their contributions to modulating the stability and flexibility of RNA molecules.

Affinity and Structural Analysis of the U1A RNA Recognition Motif with Engineered Methionines to Improve Experimental Phasing

RNA plays a central role in all organisms and can fold into complex structures to orchestrate function. Visualization of such structures often requires crystallization, which can be a bottleneck in the structure-determination process. To promote crystallization, an RNA-recognition motif (RRM) of the U1A spliceosomal protein has been co-opted as a crystallization module. Specifically, the U1-snRNA hairpin II (hpII) single-stranded loop recognized by U1A can be transplanted into an RNA target to promote crystal contacts and to attain phase information via molecular replacement or anomalous diffraction methods using selenomethionine. Herein, we produced the F37M/F77M mutant of U1A to augment the phasing capability of this powerful crystallization module. Selenomethionine-substituted U1A(F37M/F77M) retains high affinity for hpII (K _D of 59.7 ± 11.4 nM). The 2.20 Å resolution crystal structure reveals that the mutated sidechains make new S-π interactions in the hydrophobic core and are useful for single-wavelength anomalous diffraction. Crystals were also attained of U1A(F37M/F77M) in complex with a bacterial preQ₁-II riboswitch. The F34M/F37M/F77M mutant was introduced similarly into a lab-evolved U1A variant (TBP6.9) that recognizes the internal bulged loop of HIV-1 TAR RNA. We envision that this short RNA sequence can be placed into non-essential duplex regions to promote crystallization and phasing of target RNAs. We show that selenomethionine-substituted TBP6.9(F34M/F37M/F77M) binds a TAR variant wherein the apical loop was replaced with a GNRA tetraloop (K _D of 69.8 ± 2.9 nM), laying the groundwork for use of TBP6.9(F34M/F37M/F77M) as a crystallization module. These new tools are available to the research community.

Elucidation of Catalytic Strategies of Small Nucleolytic Ribozymes From Comparative Analysis of Active Sites

A number of small, self-cleaving ribozyme classes have been identified including the hammerhead, hairpin, hepatitis delta virus (HDV), Varkud satellite (VS), glmS, twister, hatchet, pistol, and twister sister ribozymes. Within the active sites of these ribozymes, myriad functional groups contribute to catalysis. There has been extensive structure-function analysis of individual ribozymes, but the extent to which catalytic devices are shared across different ribozyme classes is unclear. As such, emergent catalytic principles for ribozymes may await discovery. Identification of conserved catalytic devices can deepen our understanding of RNA catalysis specifically and of enzymic catalysis generally. To probe similarities and differences amongst ribozyme classes, active sites from more than 80 high-resolution crystal structures of self-cleaving ribozymes were compared computationally. We identify commonalities amongst ribozyme classes pertaining to four classic catalytic devices: deprotonation of the 2'OH nucleophile (γ), neutralization of the non-bridging oxygens of the scissile phosphate (β), neutralization of the O5' leaving group (δ), and in-line nucleophilic attack (α). In addition, we uncover conservation of two catalytic devices, each of which centers on the activation of the 2'OH nucleophile by a guanine: one to acidify the 2'OH by hydrogen bond donation to it (γ') and one to acidify the 2'OH by releasing it from non-productive interactions by competitive hydrogen bonding (γ''). Our findings reveal that the amidine functionalities of G, A, and C are especially important for these strategies, and help explain absence of U at ribozyme active sites. The identified γ' and γ'' catalytic strategies help unify the catalytic strategies shared amongst catalytic RNAs and may be important for large ribozymes, as well as protein enzymes that act on nucleic acids.

Intrinsic Base-Pair Rearrangement in the Hairpin Ribozyme Directs RNA Conformational Sampling and Tertiary Interface Formation

Dynamic fluctuations in RNA structure enable conformational changes that are required for catalysis and recognition. In the hairpin ribozyme, the catalytically active structure is formed as an intricate tertiary interface between two RNA internal loops. Substantial alterations in the structure of each loop are observed upon interface formation, or docking. The very slow on-rate for this relatively tight interaction has led us to hypothesize a double conformational capture mechanism for RNA-RNA recognition. We used extensive molecular dynamics simulations to assess conformational sampling in the undocked form of the loop domain containing the scissile phosphate (loop A). We observed several major accessible conformations with distinctive patterns of hydrogen bonding and base stacking interactions in the active-site internal loop. Several important conformational features characteristic of the docked state were observed in well-populated substates, consistent with the kinetic sampling of docking-competent states by isolated loop A. Our observations suggest a hybrid or multistage binding mechanism, in which initial conformational selection of a docking-competent state is followed by induced-fit adjustment to an in-line, chemically reactive state only after formation of the initial complex with loop B.

The Quick and the Dead: A Guide to Fast Phasing of Small Ribozyme and Riboswitch Crystal Structures

Ribozymes and riboswitches are examples of non-protein-coding (nc)RNA molecules that achieve biological activity by adopting complex three-dimensional folds. Visualization of such molecules at near-atomic resolution can enhance our understanding of how chemical groups are organized spatially, thereby providing novel insight into function. This approach has its challenges, which mainly entail sample crystallization followed by the application of empirical, structure-determination methods that often include experimental "phasing" of X-ray diffraction data. A paucity of high-quality crystals or a low symmetry space group are factors that demand rapid assessment of phasing potential during an ongoing experiment in order to assure a successful outcome. Here we describe the process of evaluating the anomalous signal-to-noise as a prelude to single wavelength or multiwavelength anomalous diffraction (SAD or MAD) phasing. Test cases include an autolytic 62-mer RNA enzyme known as the hairpin ribozyme, and a 33-mer riboswitch that binds the modified guanine metabolite preQ1. The crystals were derivatized with iridium (III) hexammine and osmium (III) pentaammine triflate, respectively. Each data set was then subjected to the XPREP and SHELX programs to assess the anomalous signal-to-noise and to locate the heavy-atom substructure. Subsequent noise filtering was conducted in SHELXE or RESOLVE. The methods described are applicable to the rapid phasing of RNA X-ray diffraction data, and contrast the efficacy of in-house X-rays with those attainable from synchrotron-radiation sources in terms of the potential to plan for and execute an experimental structure determination.

8-Azapurines as isosteric purine fluorescent probes for nucleic acid and enzymatic research

The 8-azapurines, and their 7-deaza and 9-deaza congeners, represent a unique class of isosteric (isomorphic) analogues of the natural purines, frequently capable of substituting for the latter in many biochemical processes. Particularly interesting is their propensity to exhibit pH-dependent room-temperature fluorescence in aqueous medium, and in non-polar media. We herein review the physico-chemical properties of this class of compounds, with particular emphasis on the fluorescence emission properties of their neutral and/or ionic species, which has led to their widespread use as fluorescent probes in enzymology, including enzymes involved in purine metabolism, agonists/antagonists of adenosine receptors, mechanisms of catalytic RNAs, RNA editing, etc. They are also exceptionally useful fluorescent probes for analytical and clinical applications in crude cell homogenates.

Comprehensive analysis of mutations in the hepatitis delta virus genome based on full-length sequencing in a nationwide cohort study and evolutionary pattern during disease progression

Delta hepatitis, caused by co-infection or super-infection of hepatitis D virus (HDV) in hepatitis B virus (HBV) -infected patients, is the most severe form of chronic hepatitis, often progressing to liver cirrhosis and liver failure. Although 15 million individuals are affected worldwide, molecular data on the HDV genome and its proteins, small and large delta antigen (S-/L-HDAg), are limited. We therefore conducted a nationwide study in HBV-HDV-infected patients from Iran and successfully amplified 38 HDV full genomes and 44 L-HDAg sequences from 34 individuals. Phylogenetic analyses of full-length HDV and L-HDAg isolates revealed that all strains clustered with genotype 1 and showed high genotypic distances to HDV genotypes 2 to 8, with a maximal distance to genotype 3. Longitudinal analyses in individual patients indicated a reverse evolutionary trend, especially in L-HDAg amino acid composition, over time. Besides multiple sequence variations in the hypervariable region of HDV, nucleotide substitutions preferentially occurred in the stabilizing P4 domain of the HDV ribozyme. A high rate of single amino acid changes was detected in structural parts of L-HDAg, whereas its post-translational modification sites were highly conserved. Interestingly, several non-synonymous mutations were positively selected that affected immunogenic epitopes of L-HDAg towards CD8 T-cell- and B-cell-driven immune responses. Hence, our comprehensive molecular analysis comprising a nationwide cohort revealed phylogenetic relationships and provided insight into viral evolution within individual hosts. Moreover, preferential areas of frequent mutations in the HDV ribozyme and antigen protein were determined in this study.

Comparison of ab Initio, DFT, and Semiempirical QM/MM Approaches for Description of Catalytic Mechanism of Hairpin Ribozyme

We have analyzed the capability of state-of-the-art multiscale computational approaches to provide atomic-resolution electronic structure insights into possible catalytic scenarios of the hairpin ribozyme by evaluating potential and free energy surfaces of the reactions by various hybrid QM/MM methods. The hairpin ribozyme is a unique catalytic RNA that achieves rate acceleration similar to other small self-cleaving ribozymes but without direct metal ion participation. Guanine 8 (G8) and adenine 38 (A38) have been identified as the catalytically essential nucleobases. However, their exact catalytic roles are still being investigated. In line with the available experimental data, we considered two reaction scenarios involving protonated A38H(+) as a general acid which is further assisted by either canonical G8 or deprotonated G8(-) forms. We used the spin-component scaled Møller-Plesset (SCS-MP2) method at the complete basis set limit as the reference method. The semiempirical AM1/d-PhoT and SCC-DFTBPR methods provided acceptable activation barriers with respect to the SCS-MP2 data but predicted significantly different reaction pathways. DFT functionals (BLYP and MPW1K) yielded the same reaction pathway as the SCS-MP2 method. The activation barriers were slightly underestimated by the GGA BLYP functional, although with accuracy comparable to the semiempirical methods. The SCS-MP2 method and hybrid MPW1K functional gave activation barriers that were closest to those derived from experimentally measured rate constants.

NMR localization of divalent cations at the active site of the Neurospora VS ribozyme provides insights into RNA-metal-ion interactions

Metal cations represent key elements of RNA structure and function. In the Neurospora VS ribozyme, metal cations play diverse roles; they are important for substrate recognition, formation of the active site, and shifting the pKa's of two key nucleobases that contribute to the general acid-base mechanism. Recently, we determined the NMR structure of the A730 loop of the VS ribozyme active site (SLVI) that contributes the general acid (A756) in the enzymatic mechanism of the cleavage reaction. Our studies showed that magnesium (Mg(2+)) ions are essential to stabilize the formation of the S-turn motif within the A730 loop that exposes the A756 nucleobase for catalysis. In this article, we extend these NMR investigations by precisely mapping the Mg(2+)-ion binding sites using manganese-induced paramagnetic relaxation enhancement and cadmium-induced chemical-shift perturbation of phosphorothioate RNAs. These experiments identify five Mg(2+)-ion binding sites within SLVI. Four Mg(2+) ions in SLVI are associated with known RNA structural motifs, including the G-U wobble pair and the GNRA tetraloop, and our studies reveal novel insights about Mg(2+) ion binding to these RNA motifs. Interestingly, one Mg(2+) ion is specifically associated with the S-turn motif, confirming its structural role in the folding of the A730 loop. This Mg(2+) ion is likely important for formation of the active site and may play an indirect role in catalysis.

Intermolecular domain docking in the hairpin ribozyme: metal dependence, binding kinetics and catalysis

The hairpin ribozyme is a prototype small, self-cleaving RNA motif. It exists naturally as a four-way RNA junction containing two internal loops on adjoining arms. These two loops interact in a cation-driven docking step prior to chemical catalysis to form a tightly integrated structure, with dramatic changes occurring in the conformation of each loop upon docking. We investigate the thermodynamics and kinetics of the docking process using constructs in which loop A and loop B reside on separate molecules. Using a novel CD difference assay to isolate the effects of metal ions linked to domain docking, we find the intermolecular docking process to be driven by sub-millimolar concentrations of the exchange-inert Co(NH₃)₆³⁺. RNA self-cleavage requires binding of lower-affinity ions with greater apparent cooperativity than the docking process itself, implying that, even in the absence of direct coordination to RNA, metal ions play a catalytic role in hairpin ribozyme function beyond simply driving loop-loop docking. Surface plasmon resonance assays reveal remarkably slow molecular association, given the relatively tight loop-loop interaction. This observation is consistent with a “double conformational capture” model in which only collisions between loop A and loop B molecules that are simultaneously in minor, docking-competent conformations are productive for binding.

Crystallographic analysis of small ribozymes and riboswitches

Ribozymes and riboswitches are RNA motifs that accelerate biological reactions and regulate gene expression in response to metabolite recognition, respectively. These RNA molecules gain functionality via complex folding that cannot be predicted a priori, and thus requires high-resolution three-dimensional structure determination to locate key functional attributes. Herein, we present an overview of the methods used to determine small RNA structures with an emphasis on RNA preparation, crystallization, and structure refinement. We draw upon examples from our own research in the analysis of the leadzyme ribozyme, the hairpin ribozyme, a class I preQ(1) riboswitch, and variants of a larger class II preQ(1) riboswitch. The methods presented provide a guide for comparable investigations of noncoding RNA molecules including a 48-solution, "first choice" RNA crystal screen compiled from our prior successes with commercially available screens.

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

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

Catalysis by the nucleolytic ribozymes

The nucleolytic ribozymes use general acid-base catalysis to contribute significantly to their rate enhancement. The VS (Varkud satellite) ribozyme uses a guanine and an adenine nucleobase as general base and acid respectively in the cleavage reaction. The hairpin ribozyme is probably closely similar, while the remaining nucleolytic ribozymes provide some interesting contrasts.

NMR structure of the A730 loop of the Neurospora VS ribozyme: insights into the formation of the active site

The Neurospora VS ribozyme is a small nucleolytic ribozyme with unique primary, secondary and global tertiary structures, which displays mechanistic similarities to the hairpin ribozyme. Here, we determined the high-resolution NMR structure of a stem–loop VI fragment containing the A730 internal loop, which forms part of the active site. In the presence of magnesium ions, the A730 loop adopts a structure that is consistent with existing biochemical data and most likely reflects its conformation in the VS ribozyme prior to docking with the cleavage site internal loop. Interestingly, the A730 loop adopts an S-turn motif that is also present in loop B within the hairpin ribozyme active site. The S-turn appears necessary to expose the Watson–Crick edge of a catalytically important residue (A756) so that it can fulfill its role in catalysis. The A730 loop and the cleavage site loop of the VS ribozyme display structural similarities to internal loops found in the active site of the hairpin ribozyme. These similarities provided a rationale to build a model of the VS ribozyme active site based on the crystal structure of the hairpin ribozyme.

Use of the spliceosomal protein U1A to facilitate crystallization and structure determination of complex RNAs

The structure determination of complex RNA molecules such as ribozymes, riboswitches and aptamers by X-ray crystallography hinges on the preparation of well-ordered crystals. Success usually results from molecular engineering to facilitate crystallization. An approach that has resulted in 10 new RNA structures in the past decade is the use of the U1A crystallization module. In this approach, the cognate site for the U1A spliceosomal protein is introduced into a functionally dispensable location in the RNA of interest, and the RNA is cocrystallized with the basic RNA-binding protein. In addition to facilitating crystallization, the presence of U1A can be useful for de novo phase determination. In this paper, some general considerations for the use of this approach to RNA crystallization are presented, and specifics of the application of the U1A module to the crystallization of the hairpin ribozyme and the tetracycline aptamer are reviewed.

Extensive molecular dynamics simulations showing that canonical G8 and protonated A38H+ forms are most consistent with crystal structures of hairpin ribozyme

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.

Molecular dynamics and quantum mechanics of RNA: conformational and chemical change we can believe in

Structure and dynamics are both critical to RNA’s vital functions in biology. Numerous techniques can elucidate the structural dynamics of RNA, but computational approaches based on experimental data arguably hold the promise of providing the most detail. In this Account, we highlight areas wherein molecular dynamics (MD) and quantum mechanical (QM) techniques are applied to RNA, particularly in relation to complementary experimental studies.

We have expanded on atomic-resolution crystal structures of RNAs in functionally relevant states by applying explicit solvent MD simulations to explore their dynamics and conformational changes on the submicrosecond time scale. MD relies on simplified atomistic, pairwise additive interaction potentials (force fields). Because of limited sampling, due to the finite accessible simulation time scale and the approximated force field, high-quality starting structures are required.

Despite their imperfection, we find that currently available force fields empower MD to provide meaningful and predictive information on RNA dynamics around a crystallographically defined energy minimum. The performance of force fields can be estimated by precise QM calculations on small model systems. Such calculations agree reasonably well with the Cornell et al. AMBER force field, particularly for stacking and hydrogen-bonding interactions. A final verification of any force field is accomplished by simulations of complex nucleic acid structures.

The performance of the Cornell et al. AMBER force field generally corresponds well with and augments experimental data, but one notable exception could be the capping loops of double-helical stems. In addition, the performance of pairwise additive force fields is obviously unsatisfactory for inclusion of divalent cations, because their interactions lead to major polarization and charge-transfer effects neglected by the force field. Neglect of polarization also limits, albeit to a lesser extent, the description accuracy of other contributions, such as interactions with monovalent ions, conformational flexibility of the anionic sugar−phosphate backbone, hydrogen bonding, and solute polarization by solvent. Still, despite limitations, MD simulations are a valid tool for analyzing the structural dynamics of existing experimental structures. Careful analysis of MD simulations can identify problematic aspects of an experimental RNA structure, unveil structural characteristics masked by experimental constraints, reveal functionally significant stochastic fluctuations, evaluate the structural role of base ionization, and predict structurally and potentially functionally important details of the solvent behavior, including the presence of tightly bound water molecules. Moreover, combining classical MD simulations with QM calculations in hybrid QM/MM approaches helps in the assessment of the plausibility of chemical mechanisms of catalytic RNAs (ribozymes).

In contrast, the reliable prediction of structure from sequence information is beyond the applicability of MD tools. The ultimate utility of computational studies in understanding RNA function thus requires that the results are neither blindly accepted nor flatly rejected, but rather considered in the context of all available experimental data, with great care given to assessing limitations through the available starting structures, force field approximations, and sampling limitations. The examples given in this Account showcase how the judicious use of basic MD simulations has already served as a powerful tool to help evaluate the role of structural dynamics in biological function of RNA.

Thermodynamic analysis of ligand binding and ligand binding-induced tertiary structure formation by the thiamine pyrophosphate riboswitch

The thi-box riboswitch regulates gene expression in response to the intracellular concentration of thiamine pyrophosphate (TPP) in archaea, bacteria, and eukarya. To complement previous biochemical, genetic, and structural studies of this phylogenetically widespread RNA domain, we have characterized its interaction with TPP by isothermal titration calorimetry. This shows that TPP binding is highly dependent on Mg(2+) concentration. The dissociation constant decreases from approximately 200 nM at 0.5 mM Mg(2+) concentration to approximately 9 nM at 2.5 mM Mg(2+) concentration. Binding is enthalpically driven, but the unfavorable entropy of binding decreases as Mg(2+) concentration rises, suggesting that divalent cations serve to pre-organize the RNA. Mutagenesis, biochemical analysis, and a new crystal structure of the riboswitch suggest that a critical element that participates in organizing the riboswitch structure is the tertiary interaction formed between the P3 and L5 regions. This tertiary contact is distant from the TPP binding site, but calorimetric analysis reveals that even subtle mutations in L5 can have readily detectable effects on TPP binding. The thermodynamic signatures of these mutations, namely decreased favorable enthalpy of binding and small effects on entropy of binding, are consistent with the P3-L5 association contributing allosterically to TPP-induced compaction of the RNA.

Exploring ribozyme conformational changes with X-ray crystallography

Relating three-dimensional fold to function is a central challenge in RNA structural biology. Toward this goal, X-ray crystallography has long been considered the "gold standard" for structure determinations at atomic resolution, although NMR spectroscopy has become a powerhouse in this arena as well. In the area of dynamics, NMR remains the dominant technique to probe the magnitude and timescales of molecular motion. Although the latter area remains largely unassailable by conventional crystallographic methods, inroads have been made on proteins using Laue radiation on timescales of ms to ns. Proposed 'fourth generation' radiation sources, such as free-electron X-ray lasers, promise ps- to fs-timescale resolution, and credible evidence is emerging that supports the feasibility of single molecule imaging. At present however, the preponderance of RNA structural information has been derived from timescale and motion insensitive crystallographic techniques. Importantly, developments in computing, automation and high-flux synchrotron sources have propelled the rapidity of 'conventional' RNA crystal structure determinations to timeframes of hours once a suitable set of phases is obtained. With a sufficient number of crystal structures, it is possible to create a structural ensemble that can provide insight into global and local molecular motion characteristics that are relevant to biological function. Here we describe techniques to explore conformational changes in the hairpin ribozyme, a representative non-protein-coding RNA catalyst. The approaches discussed include: (i) construct choice and design using prior knowledge to improve X-ray diffraction; (ii) recognition of long-range conformational changes and (iii) use of single-base or single-atom changes to create ensembles. The methods are broadly applicable to other RNA systems.

Single-atom imino substitutions at A9 and A10 reveal distinct effects on the fold and function of the hairpin ribozyme catalytic core

The hairpin ribozyme cleaves a phosphodiester bond within a cognate substrate. Structural and biochemical data indicate the conserved A9 and A10 bases reside close to the scissile bond but make distinct contributions to catalysis. To investigate these residues, we replaced the imino moiety of each base with N1-deazaadenosine. This single-atom change resulted in an 8-fold loss in k(obs) for A9 and displacement of the base from the active site; no effects were observed for A10. We propose that the imino moiety of A9 promotes a key water-mediated contact that favors transition-state formation, which suggests an enhanced chemical repertoire for RNA.

Comparative enzymology and structural biology of RNA self-cleavage

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

The structural basis for recognition of the PreQ0 metabolite by an unusually small riboswitch aptamer domain

Riboswitches are RNA elements that control gene expression through metabolite binding. The preQ(1) riboswitch exhibits the smallest known ligand-binding domain and is of interest for its economical organization and high affinity interactions with guanine-derived metabolites required to confer tRNA wobbling. Here we present the crystal structure of a preQ(1) aptamer domain in complex with its precursor metabolite preQ(0). The structure is highly compact with a core that features a stem capped by a well organized decaloop. The metabolite is recognized within a deep pocket via Watson-Crick pairing with C15. Additional hydrogen bonds are made to invariant bases U6 and A29. The ligand-bound state confers continuous helical stacking throughout the core fold, thus providing a platform to promote Watson-Crick base pairing between C9 of the decaloop and the first base of the ribosome-binding site, G33. The structure offers insight into the mode of ribosome-binding site sequestration by a minimal RNA fold stabilized by metabolite binding and has implications for understanding the molecular basis by which bacterial genes are regulated.

Molecular dynamics suggest multifunctionality of an adenine imino group in acid-base catalysis of the hairpin ribozyme

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

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

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

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

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

Hammerhead redux: does the new structure fit the old biochemical data?

The cleavage rates of 78 hammerhead ribozymes containing structurally conservative chemical modifications were collected from the literature and compared to the recently determined crystal structure of the Schistosoma mansoni hammerhead. With only a few exceptions, the biochemical data were consistent with the structure, indicating that the new structure closely resembles the transition state of the reaction. Since all the biochemical data were collected on minimal hammerheads that have a very different structure, the minimal hammerhead must be dynamic and occasionally adopt the quite different extended structure in order to cleave.

Focus on function: single molecule RNA enzymology

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.

A comparison of vanadate to a 2'-5' linkage at the active site of a small ribozyme suggests a role for water in transition-state stabilization

The potential for water to participate in RNA catalyzed reactions has been the topic of several recent studies. Here, we report crystals of a minimal, hinged hairpin ribozyme in complex with the transition-state analog vanadate at 2.05 A resolution. Waters are present in the active site and are discussed in light of existing views of catalytic strategies employed by the hairpin ribozyme. A second structure harboring a 2',5'-phosphodiester linkage at the site of cleavage was also solved at 2.35 A resolution and corroborates the assignment of active site waters in the structure containing vanadate. A comparison of the two structures reveals that the 2',5' structure adopts a conformation that resembles the reaction intermediate in terms of (1) the positioning of its nonbridging oxygens and (2) the covalent attachment of the 2'-O nucleophile with the scissile G+1 phosphorus. The 2',5'-linked structure was then overlaid with scissile bonds of other small ribozymes including the glmS metabolite-sensing riboswitch and the hammerhead ribozyme, and suggests the potential of the 2',5' linkage to elicit a reaction-intermediate conformation without the need to form metalloenzyme complexes. The hairpin ribozyme structures presented here also suggest how water molecules bound at each of the nonbridging oxygens of G+1 may electrostatically stabilize the transition state in a manner that supplements nucleobase functional groups. Such coordination has not been reported for small ribozymes, but is consistent with the structures of protein enzymes. Overall, this work establishes significant parallels between the RNA and protein enzyme worlds.

A posteriori design of crystal contacts to improve the X-ray diffraction properties of a small RNA enzyme

The hairpin ribozyme is a small catalytic RNA comprising two helix-loop-helix domains linked by a four-way helical junction (4WJ). In its most basic form, each domain can be formed independently and reconstituted without a 4WJ to yield an active enzyme. The production of such minimal junctionless hairpin ribozymes is achievable by chemical synthesis, which has allowed structures to be determined for numerous nucleotide variants. However, abasic and other destabilizing core modifications hinder crystallization. This investigation describes the use of a dangling 5'-U to form an intermolecular U.U mismatch, as well as the use of synthetic linkers to tether the loop A and B domains, including (i) a three-carbon propyl linker (C3L) and (ii) a nine-atom triethylene glycol linker (S9L). Both linker constructs demonstrated similar enzymatic activity, but S9L constructs yielded crystals that diffracted to 2.65 A resolution or better. In contrast, C3L variants diffracted to 3.35 A and exhibited a 15 A expansion of the c axis. Crystal packing of the C3L construct showed a paucity of 6(1) contacts, which comprise numerous backbone to 2'-OH hydrogen bonds in junctionless and S9L complexes. Significantly, the crystal packing in minimal structures mimics stabilizing features observed in the 4WJ hairpin ribozyme structure. The results demonstrate how knowledge-based design can be used to improve diffraction and overcome otherwise destabilizing defects.

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

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

Functional Analysis of Hairpin Ribozyme Active Site Architecture

The hairpin ribozyme is a small catalytic motif found in plant satellite RNAs where it catalyzes a reversible self-cleavage reaction during processing of replication intermediates. Crystallographic studies of hairpin ribozymes have provided high resolution views of the RNA functional groups that comprise the active site and stimulated biochemical studies that probed the contributions of nucleobase functional groups to catalytic chemistry. The dramatic loss of activity that results from perturbation of active site architecture points to the importance of positioning and orientation in catalytic rate acceleration. The current study focuses on the network of noncovalent interactions that align nucleophilic and leaving group oxygens in the orientation required for the S(N)2-type reaction mechanism and orient the active site nucleobases near the reactive phosphate to facilitate catalytic chemistry. Nucleotide modifications that alter or eliminate individual hydrogen bonding partners had different effects on the activation barrier to catalysis, the stability of ribozyme complexes in the ground state, and the internal equilibrium between cleavage and ligation of bound products. Furthermore, substitution of hydrogen bond donors and acceptors with seemingly equivalent pairs sometimes had very different functional consequences. These biochemical analyses augment high resolution structural information to provide insights into the functional significance of active site architecture.

Water in the active site of an all-RNA hairpin ribozyme and effects of Gua8 base variants on the geometry of phosphoryl transfer

The hairpin ribozyme requires functional group contributions from G8 to assist in phosphodiester bond cleavage. Previously, replacement of G8 by a series of nucleobase variants showed little effect on interdomain docking, but a 3-250-fold effect on catalysis. To identify G8 features that contribute to catalysis within the hairpin ribozyme active site, structures for five base variants were determined by X-ray crystallography in a resolution range between 2.3 and 2.7 A. For comparison, a native all-RNA "G8" hairpin ribozyme structure was refined to 2.05 A resolution. The native structure revealed a scissile bond angle (tau) of 158 degrees, which is close to the requisite 180 degrees "in-line" geometry. Mutations G8(inosine), G8(diaminopurine), G8(aminopurine), G8(adenosine), and G8(uridine) folded properly, but exhibited nonideal scissile bond geometries (tau ranging from 118 degrees to 93 degrees) that paralleled their diminished solution activities. A superposition ensemble of all structures, including a previously described hairpin ribozyme-vanadate complex, indicated the scissile bond can adopt a variety of conformations resulting from perturbation of the chemical environment and provided a rationale for how the exocyclic amine of nucleobase 8 promotes productive, in-line geometry. Changes at position 8 also caused variations in the A-1 sugar pucker. In this regard, variants A8 and U8 appeared to represent nonproductive ground states in which their 2'-OH groups mimicked the pro-R, nonbridging oxygen of the vanadate transition-state complex. Finally, the results indicated that ordered water molecules bind near the 2'-hydroxyl of A-1, lending support to the hypothesis that solvent may play an important role in the reaction.