Solution structure of loop A from the hairpin ribozyme from tobacco ringspot virus satellite

Insights into the A-C Mismatch Conformational Ensemble in Duplex DNA and its Role in Genetic Processes through a Structure-based Review

Knowing the conformational ensembles formed by mismatches is crucial for understanding how they are generated and repaired and how they contribute to genomic instability. Here, we review structural and energetic studies of the A-C mismatch in duplex DNA and use the information to identify critical conformational states in its ensemble and their significance in genetic processes. In the 1970s, Topal and Fresco proposed the A-C wobble stabilized by two hydrogen bonds, one requiring protonation of adenine-N1. Subsequent NMR and X-ray crystallography studies showed that the protonated A-C wobble was in dynamic equilibrium with a neutral inverted wobble. The mismatch was shown to destabilize duplex DNA in a sequence- and pH-dependent manner by 2.4-3.8 kcal/mol and to have an apparent pKa ranging between 7.2 and 7.7. The A-C mismatch conformational repertoire expanded as structures were determined for damaged and protein-bound DNA. These structures included Watson-Crick-like conformations forming through tautomerization of the bases that drive replication errors, the reverse wobble forming through rotation of the entire nucleotide proposed to increase the fidelity of DNA replication, and the Hoogsteen base-pair forming through the flipping of the adenine base which explained the unusual specificity of DNA polymerases that bypass DNA damage. Thus, the A-C mismatch ensemble encompasses various conformational states that can be selectively stabilized in response to environmental changes such as pH shifts, intermolecular interactions, and chemical modifications, and these adaptations facilitate critical biological processes. This review also highlights the utility of existing 3D structures to build ensemble models for nucleic acid motifs.

Cross-correlated relaxation rates provide facile exchange signature in selectively labeled RNA

Gerhard Wagner has made numerous contributions to NMR spectroscopy, particularly his developments in the field of spin-relaxation stand out in directly mapping the spectral density functions of proteins. He and his group developed experimental techniques to reveal the importance of dynamics to protein biological function and drug discovery. On his 75th birthday, we take this opportunity to highlight how some of those seminal ideas developed for proteins are being extended to RNAs. The role of dynamics in the structure and function of RNA has been a major interest in drug design and therapeutics. Here we present the use of cross-correlated relaxation rates (η_xy) from anti-TROSY (R_2α) and TROSY (R_2β) to rapidly obtain qualitative information about the chemical exchange taking place within the bacterial and human A-site RNA system while reducing the sets of relaxation experiments required to map dynamics. We show that η_xy correlates with the order parameter which gives information on how flexible or rigid a residue is. We further show R_2β/η_xy can rapidly be used to probe chemical exchange as seen from its agreement with R_ex. In addition, we report the ability of R_2β/η_xy to determine chemical exchange taking place within the bacterial A-site RNA during structural transitions at pH 6.2 and 6.5. Finally, comparison of the R_2β/η_xy ratios indicates bacterial A-site has greater R_2β/η_xy values for G19 (1.34 s^-1), A20 (1.38 s^-1), U23 (1.63 s^-1) and C24 (1.51 s^-1) than human A-site [A19 (0.76 s^-1), A20 (1.01 s^-1), U23 (0.74 s^-1) and C24 (0.71 s^-1)]. Taken together, we have shown that the chemical exchange can quickly be analyzed for RNA systems from cross-correlated relaxation rates.

Alternate RNA Structures

RNA molecules fold into complex three-dimensional structures that sample alternate conformations ranging from minor differences in tertiary structure dynamics to major differences in secondary structure. This allows them to form entirely different substructures with each population potentially giving rise to a distinct biological outcome. The substructures can be partitioned along an existing energy landscape given a particular static cellular cue or can be shifted in response to dynamic cues such as ligand binding. We review a few key examples of RNA molecules that sample alternate conformations and how these are capitalized on for control of critical regulatory functions.

The effect of adenine protonation on RNA phosphodiester backbone bond cleavage elucidated by deaza-nucleobase modifications and mass spectrometry

The catalytic strategies of small self-cleaving ribozymes often involve interactions between nucleobases and the ribonucleic acid (RNA) backbone. Here we show that multiply protonated, gaseous RNA has an intrinsic preference for the formation of ionic hydrogen bonds between adenine protonated at N3 and the phosphodiester backbone moiety on its 5'-side that facilitates preferential phosphodiester backbone bond cleavage upon vibrational excitation by low-energy collisionally activated dissociation. Removal of the basic N3 site by deaza-modification of adenine was found to abrogate preferential phosphodiester backbone bond cleavage. No such effects were observed for N1 or N7 of adenine. Importantly, we found that the pH of the solution used for generation of the multiply protonated, gaseous RNA ions by electrospray ionization affects phosphodiester backbone bond cleavage next to adenine, which implies that the protonation patterns in solution are at least in part preserved during and after transfer into the gas phase. Our study suggests that interactions between protonated adenine and phosphodiester moieties of RNA may play a more important mechanistic role in biological processes than considered until now.

The Role of Proton Transfer on Mutations

Hydrogen bonds play a critical role in nucleobase studies as they encode genes, map protein structures, provide stability to the base pairs, and are involved in spontaneous and induced mutations. Proton transfer mechanism is a critical phenomenon that is related to the acid-base characteristics of the nucleobases in Watson-Crick base pairs. The energetic and dynamical behavior of the proton can be depicted from these characteristics and their adjustment to the water molecules or the surrounding ions. Further, new pathways open up in which protonated nucleobases are generated by proton transfer from the ionized water molecules and elimination of a hydroxyl radical in this review, the analysis will be focused on understanding the mechanism of untargeted mutations in canonical, wobble, Hoogsteen pairs, and mutagenic tautomers through the non-covalent interactions. Further, rare tautomer formation through the single proton transfer (SPT) and the double proton transfer (DPT), quantum tunneling in nucleobases, radiation-induced bystander effects, role of water in proton transfer (PT) reactions, PT in anticancer drugs-DNA interaction, displacement and oriental polarization, possible models for mutations in DNA, genome instability, and role of proton transfer using kinetic parameters for RNA will be discussed.

Coupling between conformational dynamics and catalytic function at the active site of the lead-dependent ribozyme

In common with other self-cleaving RNAs, the lead-dependent ribozyme (leadzyme) undergoes dynamic fluctuations to a chemically activated conformation. We explored the connection between conformational dynamics and self-cleavage function in the leadzyme using a combination of NMR spin-relaxation analysis of ribose groups and conformational restriction via chemical modification. The functional studies were performed with a North-methanocarbacytidine modification that prevents fluctuations to C2'-endo conformations while maintaining an intact 2'-hydroxyl nucleophile. Spin-relaxation data demonstrate that the active-site Cyt-6 undergoes conformational exchange attributed to sampling of a minor C2'-endo state with an exchange lifetime on the order of microseconds to tens of microseconds. A conformationally restricted species in which the fluctuations to the minor species are interrupted shows a drastic decrease in self-cleavage activity. Taken together, these data indicate that dynamic sampling of a minor species at the active site of this ribozyme, and likely of related naturally occurring motifs, is strongly coupled to catalytic function. The combination of NMR dynamics analysis with functional probing via conformational restriction is a general methodology for dissecting dynamics-function relationships in RNA.

Adenine protonation enables cyclic-di-GMP binding to cyclic-GAMP sensing riboswitches

In certain structural or functional contexts, RNA structures can contain protonated nucleotides. However, a direct role for stably protonated nucleotides in ligand binding and ligand recognition has not yet been demonstrated unambiguously. Previous X-ray structures of c-GAMP binding riboswitch aptamer domains in complex with their near-cognate ligand c-di-GMP suggest that an adenine of the riboswitch either forms two hydrogen bonds to a G nucleotide of the ligand in the unusual enol tautomeric form or that the adenine in its N1 protonated form binds the G nucleotide of the ligand in its canonical keto tautomeric state. By using NMR spectroscopy we demonstrate that the c-GAMP riboswitches bind c-di-GMP using a stably protonated adenine in the ligand binding pocket. Thereby, we provide novel insights into the putative biological functions of protonated nucleotides in RNA, which in this case influence the ligand selectivity in a riboswitch.

A novel form of RNA double helix based on G·U and C·A+ wobble base pairing

Wobble base pairs are critical in various physiological functions and have been linked to local structural perturbations in double-helical structures of nucleic acids. We report a 1.38-Å resolution crystal structure of an antiparallel octadecamer RNA double helix in overall A conformation, which includes a unique, central stretch of six consecutive wobble base pairs (W helix) with two G·U and four rare C·A+ wobble pairs. Four adenines within the W helix are N1-protonated and wobble-base-paired with the opposing cytosine through two regular hydrogen bonds. Combined with the two G·U pairs, the C·A+ base pairs facilitate formation of a half turn of W-helical RNA flanked by six regular Watson-Crick base pairs in standard A conformation on either side. RNA melting experiments monitored by differential scanning calorimetry, UV and circular dichroism spectroscopy demonstrate that the RNA octadecamer undergoes a pH-induced structural transition which is consistent with the presence of a duplex with C·A+ base pairs at acidic pH. Our crystal structure provides a first glimpse of an RNA double helix based entirely on wobble base pairs with possible applications in RNA or DNA nanotechnology and pH biosensors.

Thermodynamics and kinetics of RNA tertiary structure formation in the junctionless hairpin ribozyme

The hairpin ribozyme consists of two RNA internal loops that interact to form the catalytically active structure. This docking transition is a rare example of intermolecular formation of RNA tertiary structure without coupling to helix annealing. We have used temperature-dependent surface plasmon resonance (SPR) to characterize the thermodynamics and kinetics of RNA tertiary structure formation for the junctionless form of the ribozyme, in which loops A and B reside on separate molecules. We find docking to be strongly enthalpy-driven and to be accompanied by substantial activation barriers for association and dissociation, consistent with the structural reorganization of both internal loops upon complex formation. Comparisons with the parallel analysis of a ribozyme variant carrying a 2'-O-methyl modification at the self-cleavage site and with published data in other systems reveal a surprising diversity of thermodynamic signatures, emphasizing the delicate balance of contributions to the free energy of formation of RNA tertiary structure.

Structural Basis for Substrate Helix Remodeling and Cleavage Loop Activation in the Varkud Satellite Ribozyme

The Varkud satellite (VS) ribozyme catalyzes site-specific RNA cleavage and ligation reactions. Recognition of the substrate involves a kissing loop interaction between the substrate and the catalytic domain of the ribozyme, resulting in a rearrangement of the substrate helix register into a so-called "shifted" conformation that is critical for substrate binding and activation. We report a 3.3 Å crystal structure of the complete ribozyme that reveals the active, shifted conformation of the substrate, docked into the catalytic domain of the ribozyme. Comparison to previous NMR structures of isolated, inactive substrates provides a physical description of substrate remodeling, and implicates roles for tertiary interactions in catalytic activation of the cleavage loop. Similarities to the hairpin ribozyme cleavage loop activation suggest general strategies to enhance fidelity in RNA folding and ribozyme cleavage.

Insights into RNA structure and dynamics from recent NMR and X-ray studies of the Neurospora Varkud satellite ribozyme

Despite the large number of noncoding RNAs and their importance in several biological processes, our understanding of RNA structure and dynamics at atomic resolution is still very limited. Like many other RNAs, the Neurospora Varkud satellite (VS) ribozyme performs its functions through dynamic exchange of multiple conformational states. More specifically, the VS ribozyme recognizes and cleaves its stem-loop substrate via a mechanism that involves several structural transitions within its stem-loop substrate. The recent publications of high-resolution structures of the VS ribozyme, obtained by NMR spectroscopy and X-ray crystallography, offer an opportunity to integrate the data and closely examine the structural and dynamic properties of this model RNA system. Notably, these investigations provide a valuable example of the divide-and-conquer strategy for structural and dynamic characterization of a large RNA, based on NMR structures of several individual subdomains. The success of this divide-and-conquer approach reflects the modularity of RNA architecture and the great care taken in identifying the independently-folding modules. Together with previous biochemical and biophysical characterizations, the recent NMR and X-ray studies provide a coherent picture into how the VS ribozyme recognizes its stem-loop substrate. Such in-depth characterization of this RNA enzyme will serve as a model for future structural and engineering studies of dynamic RNAs and may be particularly useful in planning divide-and-conquer investigations. WIREs RNA 2017, 8:e1421. doi: 10.1002/wrna.1421 For further resources related to this article, please visit the WIREs website.

A Stably Protonated Adenine Nucleotide with a Highly Shifted pKa Value Stabilizes the Tertiary Structure of a GTP-Binding RNA Aptamer

RNA tertiary structure motifs are stabilized by a wide variety of hydrogen-bonding interactions. Protonated A and C nucleotides are normally not considered to be suitable building blocks for such motifs since their pK_a values are far from physiological pH. Here, we report the NMR solution structure of an in vitro selected GTP-binding RNA aptamer bound to GTP with an intricate tertiary structure. It contains a novel kind of base quartet stabilized by a protonated A residue. Owing to its unique structural environment in the base quartet, the pK_a value for the protonation of this A residue in the complex is shifted by more than 5 pH units compared to the pK_a for A nucleotides in single-stranded RNA. This is the largest pK_a shift for an A residue in structured nucleic acids reported so far, and similar in size to the largest pK_a shifts observed for amino acid side chains in proteins. Both RNA pre-folding and ligand binding contribute to the pK_a shift.

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.

Pistol ribozyme adopts a pseudoknot fold facilitating site-specific in-line cleavage

The field of small self-cleaving nucleolytic ribozymes has been invigorated by the recent discovery of the twister, twister-sister, pistol and hatchet ribozymes. We report the crystal structure of a pistol ribozyme termed env25, which adopts a compact tertiary architecture stabilized by an embedded pseudoknot fold. The G-U cleavage site adopts a splayed-apart conformation with in-line alignment of the modeled 2'-O of G for attack on the adjacent to-be-cleaved P-O5' bond. Highly conserved residues G40 (N1 position) and A32 (N3 and 2'-OH positions) are aligned to act as a general base and a general acid, respectively, to accelerate cleavage chemistry, with their roles confirmed by cleavage assays on variants, and an increased pKa of 4.7 for A32. Our structure of the pistol ribozyme defined how the overall and local topologies dictate the in-line alignment at the G-U cleavage site, with cleavage assays on variants revealing key residues that participate in acid-base-catalyzed cleavage chemistry.

Structure, dynamics, and mechanism of the lead-dependent ribozyme

Leadzyme is a small catalytic RNA that was identified by in vitro selection for Pb2+-dependent cleavage from a tRNA library. Leadzyme employs a unique two-step Pb2+-specific mechanism to cleave within its active site. NMR and crystal structures of the active site revealed different folding patterns, but neither features the in-line alignment for attack by the 2'-OH nucleophilic group. These experimentally determined structures most likely represent ground states and are catalytically inactive. There are significant dynamics of the active site and the motif samples multiple conformations at the ground states. Various metal ion binding sites have been identified, including one that may be occupied by a catalytic Pb2+. Based on functional group analysis, a computational model of the transition state has been proposed. This model features a unique base triple that is consistent with sequence and functional group requirements for catalysis. This structure is likely only populated transiently, but imposing appropriate conformational constraints may significantly stabilize this state thereby promoting catalysis. Other ions may inhibit the cleavage by competing for the Pb2+ binding site, or by stabilizing the ground state thereby suppressing its transition to the catalytically active conformation. Some rare earth ions can enhance the reaction via an unknown mechanism. Because of its unique chemistry and dynamic behavior, leadzyme can continue to serve as an excellent model system for teaching us RNA biology and chemistry.

Crystal structure of the Varkud satellite ribozyme

Varkud Satellite (VS) ribozyme mediates rolling circle replication of a plasmid found in the Neurospora mitochondria. We report crystal structures of this ribozyme at 3.1Å resolution, revealing an intertwined dimer formed by an exchange of substrate helices. Within each protomer, an arrangement of three-way helical junctions organizes seven helices into a global fold that creates a docking site for the substrate helix of the other protomer, resulting in the formation of two active sites in trans. This mode of RNA-RNA association resembles the process of domain swapping in proteins and has implications for RNA regulation and evolution. Within each active site, adenine and guanine nucleobases abut the scissile phosphate, poised to serve direct roles in catalysis. Similarities to the active sites of the hairpin and hammerhead ribozymes highlight the functional significance of active site features, underscore the ability of RNA to access functional architectures from distant regions of sequence space, and suggest convergent evolution.

Structural features of a 3' splice site in influenza a

Influenza A is an RNA virus with a genome of eight negative sense segments. Segment 7 mRNA contains a 3′ splice site for alternative splicing to encode the essential M2 protein. On the basis of sequence alignment and chemical mapping experiments, the secondary structure surrounding the 3′ splice site has an internal loop, adenine bulge, and hairpin loop when it is in the hairpin conformation that exposes the 3′ splice site. We report structural features of a three-dimensional model of the hairpin derived from nuclear magnetic resonance spectra and simulated annealing with restrained molecular dynamics. Additional insight was provided by modeling based on ¹H chemical shifts. The internal loop containing the 3′ splice site has a dynamic guanosine and a stable imino (cis Watson–Crick/Watson–Crick) GA pair. The adenine bulge also appears to be dynamic with the A either stacked in the stem or forming a base triple with a Watson–Crick GC pair. The hairpin loop is a GAAA tetraloop closed by an AC pair.

The role of N7 protonation of guanine in determining the structure, stability and function of RNA base pairs

Ab initio computations and bioinformatics studies reveal that stabilization of some important RNA structural motifs might involve N7 protonation of guanine.

Building a stable RNA U-turn with a protonated cytidine

The U-turn is a classical three-dimensional RNA folding motif first identified in the anticodon and T-loops of tRNAs. It also occurs frequently as a building block in other functional RNA structures in many different sequence and structural contexts. U-turns induce sharp changes in the direction of the RNA backbone and often conform to the 3-nt consensus sequence 5'-UNR-3' (N = any nucleotide, R = purine). The canonical U-turn motif is stabilized by a hydrogen bond between the N3 imino group of the U residue and the 3' phosphate group of the R residue as well as a hydrogen bond between the 2'-hydroxyl group of the uridine and the N7 nitrogen of the R residue. Here, we demonstrate that a protonated cytidine can functionally and structurally replace the uridine at the first position of the canonical U-turn motif in the apical loop of the neomycin riboswitch. Using NMR spectroscopy, we directly show that the N3 imino group of the protonated cytidine forms a hydrogen bond with the backbone phosphate 3' from the third nucleotide of the U-turn analogously to the imino group of the uridine in the canonical motif. In addition, we compare the stability of the hydrogen bonds in the mutant U-turn motif to the wild type and describe the NMR signature of the C+-phosphate interaction. Our results have implications for the prediction of RNA structural motifs and suggest simple approaches for the experimental identification of hydrogen bonds between protonated C-imino groups and the phosphate backbone.

Hierarchy of RNA functional dynamics

RNA dynamics play a fundamental role in many cellular functions. However, there is no general framework to describe these complex processes, which typically consist of many structural maneuvers that occur over timescales ranging from picoseconds to seconds. Here, we classify RNA dynamics into distinct modes representing transitions between basins on a hierarchical free-energy landscape. These transitions include large-scale secondary-structural transitions at >0.1-s timescales, base-pair/tertiary dynamics at microsecond-to-millisecond timescales, stacking dynamics at timescales ranging from nanoseconds to microseconds, and other "jittering" motions at timescales ranging from picoseconds to nanoseconds. We review various modes within these three different tiers, the different mechanisms by which they are used to regulate function, and how they can be coupled together to achieve greater functional complexity.

Cobalt(III)hexaammine-dependent photocrosslinks in the hairpin ribozyme

We have utilized the hairpin ribozyme, an RNA enzyme whose structure has been solved by high-resolution methods, to develop a new tool for mapping nucleobase-stacking interactions and potential metal-binding sites in RNA molecules. This tool involves the photoactivation of a specifically bound cobalt(III)hexaammine molecule at wavelengths corresponding to excitation of the metal ion complex only; no base excitation is involved. The photoexcitation initiates a process which strongly promotes the formation of a novel covalent bond or crosslink between one base (termed the "first base"), which is close in space to the excited cobalt(III)hexaammine complex, and another base upon which the first base is closely stacked. These crosslinked species can be isolated and sequenced; their activities can be analyzed to ensure that the crosslinked structures represent an active conformation of the molecule. We have shown that, as in electron transfer in DNA, several criteria must be met to result in the successful formation of these crosslinks. These include the appropriate oxidation potential of the first donor base, the stacking and close interaction of the two donor bases involved in the crosslink, and the binding of a specific cobalt(III)hexaammine molecule to the first donor base. Additionally, we have determined that this crosslinking is pH-sensitive, although the cause of this sensitivity remains unknown. This tool has proven useful in the past for the analysis of the hairpin ribozyme folded structure, and has been applied to identify potential metal-binding sites on the hairpin and extended hammerhead ribozymes.

Luminescence of [Ru(bpy)2(dppz)]2+ bound to RNA mismatches

The luminescence of rac-[Ru(bpy)2(dppz)](2+) (bpy = 2,2'-bipyridine and dppz = dipyrido[3,2-a:2',3'-c]phenazine) was explored in the presence of RNA oligonucleotides containing a single RNA mismatch (CA and GG) in order to develop a probe for RNA mismatches. While there is minimal luminescence of [Ru(bpy)2(dppz)](2+) in the presence of matched RNA due to weak binding, the luminescence is significantly enhanced in the presence of a single CA mismatch. The luminescence differential between CA mismatched and matched RNA is substantially higher compared to the DNA analogue, and therefore, [Ru(bpy)2(dppz)](2+) appears to be also a sensitive light switch probe for a CA mismatch in duplex RNA. Although the luminescence intensity is lower in the presence of RNA than DNA, Förster resonance energy transfer (FRET) between the donor ruthenium complex and FRET acceptor SYTO 61 is successfully exploited to amplify the luminescence in the presence of the mismatch. Luminescence and quenching studies with sodium iodide suggest that [Ru(bpy)2(dppz)](2+) binds to these mismatches via metalloinsertion from the minor groove. This work provides further evidence that metalloinsertion is a general binding mode of octahedral metal complexes to thermodynamically destabilized mismatches not only in DNA but also in RNA.

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.

General acid-base catalysis mediated by nucleobases in the hairpin ribozyme

The catalytic mechanism by which the hairpin ribozyme accelerates cleavage or ligation of the phosphodiester backbone of RNA has been incompletely understood. There is experimental evidence for an important role for an adenine (A38) and a guanine (G8), and it has been proposed that these act in general acid-base catalysis. In this work we show that a large reduction in cleavage rate on substitution of A38 by purine (A38P) can be reversed by replacement of the 5'-oxygen atom at the scissile phosphate by sulfur (5'-PS), which is a much better leaving group. This is consistent with A38 acting as the general acid in the unmodified ribozyme. The rate of cleavage of the 5'-PS substrate by the A38P ribozyme increases with pH log-linearly, indicative of a requirement for a deprotonated base with a relatively high pK(a). On substitution of G8 by diaminopurine, the 5'-PS substrate cleavage rate at first increases with pH and then remains at a plateau, exhibiting an apparent pK(a) consistent with this nucleotide acting in general base catalysis. Alternative explanations for the pH dependence of hairpin ribozyme reactivity are discussed, from which we conclude that general acid-base catalysis by A38 and G8 is the simplest and most probable explanation consistent with all the experimental data.

Structural characterization of the viral and cRNA panhandle motifs from the infectious salmon anemia virus

Infectious salmon anemia virus (ISAV) has emerged as a virus of great concern to the aquaculture industry since it can lead to highly contagious and lethal infections in farm-raised salmon populations. While little is known about the transcription/replication cycle of ISAV, initial evidence suggests that it follows molecular mechanisms similar to those found in other orthomyxoviruses, which include the highly pathogenic influenza A (inf A) virus. During the life cycle of orthomyxoviruses, a panhandle structure is formed by the pairing of the conserved 5' and 3' ends of each genomic RNA. This structural motif serves both as a promoter of the viral RNA (vRNA)-dependent RNA polymerase and as a regulatory element in the transcription/replication cycle. As a first step toward characterizing the structure of the ISAV panhandle, here we have determined the secondary structures of the vRNA and the cRNA panhandles on the basis of solution nuclear magnetic resonance (NMR) and thermal melting data. The vRNA panhandle is distinguished by three noncanonical U · G pairs and one U · U pair in two stem helices that are linked by a highly stacked internal loop. For the cRNA panhandle, a contiguous stem helix with a protonated C · A pair near the terminus and tandem downstream U · U pairs was found. The observed noncanonical base pairs and base stacking features of the ISAV RNA panhandle motif provide the first insight into structural features that may govern recognition by the viral RNA polymerase.

New insights into the fundamental role of topological constraints as a determinant of two-way junction conformation

Recent studies have shown that topological constraints encoded at the RNA secondary structure level involving basic steric and stereochemical forces can significantly restrict the orientations sampled by helices across two-way RNA junctions. Here, we formulate these topological constraints in greater quantitative detail and use this topological framework to rationalize long-standing but poorly understood observations regarding the basic behavior of RNA two-way junctions. Notably, we show that the asymmetric nature of the A-form helix and the finite length of a bulge provide a physical basis for the experimentally observed directionality and bulge-length amplitude dependence of bulge induced inter-helical bends. We also find that the topologically allowed space can be modulated by variations in sequence, particularly with the addition of non-canonical GU base pairs at the junction, and, surprisingly, by the length of the 5′ and 3′ helices. A survey of two-way RNA junctions in the protein data bank confirms that junction residues have a strong preference to adopt looped-in, non-canonically base-paired conformations, providing a route for extending our bulge-directed framework to internal loop motifs and implying a simplified link between secondary and tertiary structure. Finally, our results uncover a new simple mechanism for coupling junction-induced topological constraints with tertiary interactions.

glmS Riboswitch binding to the glucosamine-6-phosphate α-anomer shifts the pKa toward neutrality

The glmS riboswitch regulates gene expression through a self-cleavage activity. The reaction is catalyzed with the assistance of the metabolite cofactor glucosamine-6-phosphate (GlcN6P), whose amino group is proposed to serve as the general acid during the reaction. This reaction is pH-dependent with a pK(a) that is lower than the observed pK(a) for the amine of GlcN6P in solution. GlcN6P, like other pyranose sugars, undergoes spontaneous and rapid interconversion between the α and β anomers at the C1 position. Here we demonstrate by NMR that the Bacillus anthracis glmS riboswitch selectively binds the α-anomer of GlcN6P with a maximum binding affinity of 0.36 mM and that binding is pH-dependent. We also report that the anomeric ratio between α and β is pH-dependent and the pK(a)s of the two amines differ by 0.5 pH units, α being the higher of the two (pK(a)=8.3). The pH dependence of binding reveals a pK(a) of 6.7, suggesting that the glmS RNA reduces the pK(a) of the GlcN6P amine by 1.6 units in the ground state. We reevaluated previously obtained kinetic data and found the reaction pK(a) is 6.9, within error of the binding data. The data support a model where the reaction pK(a) corresponds to that of the GlcN6P amine. This observation has broader relevance for considering how the microenvironment of an RNA, despite its anionic character, can reduce the pK(a)s of functional groups for use in catalysis.

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

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

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.

Protonation of base pairs in RNA: context analysis and quantum chemical investigations of their geometries and stabilities

Base pairs involving protonated nucleobases play important roles in mediating global macromolecular conformational changes and in facilitation of catalysis in a variety of functional RNA molecules. Here we present our attempts at understanding the role of such base pairs by detecting possible protonated base pairs in the available RNA crystal structures using BPFind software, in their specific structural contexts, and by the characterization of their geometries, interaction energies, and stabilities using advanced quantum chemical computations. We report occurrences of 18 distinct protonated base pair combinations from a representative data set of RNA crystal structures and propose a theoretical model for one putative base pair combination. Optimization of base pair geometries was carried out at the B3LYP/cc-pVTZ level, and the BSSE corrected interaction energies were calculated at the MP2/aug-cc-pVDZ level of theory. The geometries for each of the base pairs were characterized in terms of H-bonding patterns observed, rmsd values observed on optimization, and base pair geometrical parameters. In addition, the intermolecular interaction in these complexes was also analyzed using Morokuma energy decomposition. The gas phase interaction energies of the base pairs range from -24 to -49 kcal/mol and reveal the dominance of Hartree-Fock component of interaction energy constituting 73% to 98% of the total interaction energy values. On the basis of our combined bioinformatics and quantum chemical analysis of different protonated base pairs, we suggest resolution of structural ambiguities and correlate their geometric and energetic features with their structural and functional roles. In addition, we also examine the suitability of specific base pairs as key elements in molecular switches and as nucleators for higher order structures such as base triplets and quartets.

The human mitochondrial tRNAMet: structure/function relationship of a unique modification in the decoding of unconventional codons

Human mitochondrial mRNAs utilize the universal AUG and the unconventional isoleucine AUA codons for methionine. In contrast to translation in the cytoplasm, human mitochondria use one tRNA, hmtRNA(Met)(CAU), to read AUG and AUA codons at both the peptidyl- (P-), and aminoacyl- (A-) sites of the ribosome. The hmtRNA(Met)(CAU) has a unique post-transcriptional modification, 5-formylcytidine, at the wobble position 34 (f(5)C(34)), and a cytidine substituting for the invariant uridine at position 33 of the canonical U-turn in tRNAs. The structure of the tRNA anticodon stem and loop domain (hmtASL(Met)(CAU)), determined by NMR restrained molecular modeling, revealed how the f(5)C(34) modification facilitates the decoding of AUA at the P- and the A-sites. The f(5)C(34) defined a reduced conformational space for the nucleoside, in what appears to have restricted the conformational dynamics of the anticodon bases of the modified hmtASL(Met)(CAU). The hmtASL(Met)(CAU) exhibited a C-turn conformation that has some characteristics of the U-turn motif. Codon binding studies with both Escherichia coli and bovine mitochondrial ribosomes revealed that the f(5)C(34) facilitates AUA binding in the A-site and suggested that the modification favorably alters the ASL binding kinetics. Mitochondrial translation by many organisms, including humans, sometimes initiates with the universal isoleucine codons AUU and AUC. The f(5)C(34) enabled P-site codon binding to these normally isoleucine codons. Thus, the physicochemical properties of this one modification, f(5)C(34), expand codon recognition from the traditional AUG to the non-traditional, synonymous codons AUU and AUC as well as AUA, in the reassignment of universal codons in the mitochondria.

Compact intermediates in RNA folding

Large noncoding RNAs fold into their biologically functional structures via compact yet disordered intermediates, which couple the stable secondary structure of the RNA with the emerging tertiary fold. The specificity of the collapse transition, which coincides with the assembly of helical domains, depends on RNA sequence and counterions. It determines the specificity of the folding pathways and the magnitude of the free energy barriers to the ensuing search for the native conformation. By coupling helix assembly with nascent tertiary interactions, compact folding intermediates in RNA also play a crucial role in ligand binding and RNA-protein recognition.

Nucleobase-mediated general acid-base catalysis in the Varkud satellite ribozyme

Existing evidence suggests that the Varkud satellite (VS) ribozyme accelerates the cleavage of a specific phosphodiester bond using general acid-base catalysis. The key functionalities are the nucleobases of adenine 756 in helix VI of the ribozyme, and guanine 638 in the substrate stem loop. This results in a bell-shaped dependence of reaction rate on pH, corresponding to groups with pK(a) = 5.2 and 8.4. However, it is not possible from those data to determine which nucleobase is the acid, and which the base. We have therefore made substrates in which the 5' oxygen of the scissile phosphate is replaced by sulfur. This labilizes the leaving group, removing the requirement for general acid catalysis. This substitution restores full activity to the highly impaired A756G ribozyme, consistent with general acid catalysis by A756 in the unmodified ribozyme. The pH dependence of the cleavage of the phosphorothiolate-modified substrates is consistent with general base catalysis by nucleobase at position 638. We conclude that cleavage of the substrate by the VS ribozyme is catalyzed by deprotonation of the 2'-O nucleophile by G638 and protonation of the 5'-O leaving group by A756.

Catalytic importance of a protonated adenosine in the hairpin ribozyme active site

The hairpin ribozyme accelerates the rate of phosphodiester transfer reactions by at least 5 orders of magnitude. To achieve this rate enhancement, the active site forms via a substrate helix docking event that constrains the scissile phosphate linkage and positions G8 and A38 for catalysis, both of which have been implicated as sites of proton transfer in general acid-base catalysis. To investigate the functional groups required for hairpin activity, we previously reported a series of nucleotide analogue interference mapping experiments [Ryder, S. P., et al. (2001) RNA 7, 1454-1463]. The critical functional groups implicated in those studies were largely consistent with subsequent X-ray crystal structures, but the lack of A38 interference with 8-azaadenosine (n(8)A), a pK(a) perturbed nucleotide analogue, argued against functional base ionization at this site. This is inconsistent with a transition state crystal structure and other biochemical studies. To address this discrepancy, we investigated the hairpin ribozyme with an expanded set of pK(a) perturbed adenosine analogues containing fluorine. A38 was the only site that showed persistent and strong interference with low pK(a) analogues across a variety of construct/substrate pairs. This interference pattern suggests that A38 base ionization is required for catalytic activity. The lack of n(8)A interference at A38, in spite of its reduced pK(a), likely results from n(8)A stabilization of the docked state, which requires an unusual syn glycosidic base conformation at A38 for active site assembly. The fluorinated adenosine analogues are better suited to identify sites of functional ionization in systems where structural rearrangements are closely coupled to catalytic steps. All pK(a) reduced analogues, including those of the previous study, produce selective interference at A38 when substrates are stably bound and docked, consistent with the importance of base ionization at this site.

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.

Driving forces for nucleic acid pK(a) shifting in an A(+).C wobble: effects of helix position, temperature, and ionic strength

Secondary structure plays critical roles in nucleic acid function. Mismatches in DNA can lead to mutation and disease, and some mismatches involve a protonated base. Among protonated mismatches, A(+).C wobble pairs form near physiological pH and have relatively minor effects on helix geometry, making them especially important in biology. Herein, we investigate effects of helix position, temperature, and ionic strength on pK(a) shifting in A(+).C wobble pairs in DNA. We observe that pK(a) shifting is favored by internal A(+).C wobbles, which have low cooperativities of folding and make large contributions to stability, and disfavored by external A(+).C wobbles, which have high folding cooperativities but make small contributions to stability. An inverse relationship between pK(a) shifting and temperature is also found, which supports a model in which protonation is enthalpically favored overall and entropically correlated with cooperativity of folding. We also observe greater pK(a) shifts as the ionic strength decreases, consistent with anticooperativity between proton binding and counterion-condensed monovalent cation. Under the most favorable temperature and ionic strength conditions tested, a pK(a) of 8.0 is observed for the A(+).C wobble pair, which represents an especially large shift ( approximately 4.5 pK(a) units) from the unperturbed pK(a) value of adenosine. This study suggests that protonated A(+).C wobble pairs exist in DNA under biologically relevant conditions, where they can drive conformational changes and affect replication and transcription fidelity.

The i-motif in the bcl-2 P1 promoter forms an unexpectedly stable structure with a unique 8:5:7 loop folding pattern

Transcriptional regulation of the bcl-2 proto-oncogene is highly complex, with the majority of transcription driven by the P1 promoter site and the interaction of multiple regulatory proteins. A guanine- and cytosine-rich (GC-rich) region directly upstream of the P1 site has been shown to be integral to bcl-2 promoter activity, as deletion or mutation of this region significantly increases transcription. This GC-rich element consists of six contiguous runs of guanines and cytosines that have the potential to adopt DNA secondary structures, the G-quadruplex and i-motif, respectively. Our laboratory has previously demonstrated that the polypurine-rich strand of the bcl-2 promoter can form a mixture of three different G-quadruplex structures. In this current study, we demonstrate that the complementary polypyrimidine-rich strand is capable of forming one major intramolecular i-motif DNA secondary structure with a transition pH of 6.6. Characterization of the i-motif folding pattern using mutational studies coupled with circular dichroic spectra and thermal stability analyses revealed an 8:5:7 loop conformation as the predominant structure at pH 6.1. The folding pattern was further supported by chemical footprinting with bromine. In addition, a novel assay involving the sequential incorporation of a fluorescent thymine analog at each thymine position provided evidence of a capping structure within the top loop region of the i-motif. The potential of the GC-rich element within the bcl-2 promoter region to form DNA secondary structures suggests that the transition from the B-DNA to non-B-DNA conformation may play an important role in bcl-2 transcriptional regulation. Furthermore, the two adjacent large lateral loops in the i-motif structure provide an unexpected opportunity for protein and small molecule recognition.

A CA(+) pair adjacent to a sheared GA or AA pair stabilizes size-symmetric RNA internal loops

RNA internal loops are often important sites for folding and function. Residues in internal loops can have pK_a values shifted close to neutral pH because of the local structural environment. A series of RNA internal loops were studied at different pH by UV absorbance versus temperature melting experiments and imino proton nuclear magnetic resonance (NMR). A stabilizing CA pair forms at pH 7 in the and nearest neighbors when the CA pair is the first noncanonical pair (loop-terminal pair) in 3 × 3 nucleotide and larger size-symmetric internal loops. These and nearest neighbors, with CA adjacent to a closing Watson−Crick pair, are further stabilized when the pH is lowered from 7 to 5.5. The results are consistent with a significantly larger fraction (from ∼20% at pH 7 to ∼90% at pH 5.5) of adenines being protonated at the N1 position to form stabilizing wobble CA⁺ pairs adjacent to a sheared GA or AA pair. The noncanonical pair adjacent to the GA pair in can either stabilize or destabilize the loop, consistent with the sequence-dependent thermodynamics of GA pairs. No significant pH-dependent stabilization is found for most of the other nearest neighbor combinations involving CA pairs (e.g., and ), which is consistent with the formation of various nonwobble pairs observed in different local sequence contexts in crystal and NMR structures. A revised free-energy model, including stabilization by wobble CA⁺ pairs, is derived for predicting stabilities of medium-size RNA internal loops.

Direct measurement of the ionization state of an essential guanine in the hairpin ribozyme

Active site guanines are critical for self-cleavage reactions of several ribozymes, but their precise functions in catalysis are unclear. To learn whether protonated or deprotonated forms of guanine predominate in the active site, microscopic pKa values were determined for ionization of 8-azaguanosine substituted for G8 in the active site of a fully functional hairpin ribozyme in order to determine microscopic pKa values for 8-azaguanine deprotonation from the pH dependence of fluorescence. Microscopic pKa values above 9 for deprotonation of 8-azaguanine in the active site were about 3 units higher than apparent pKa values determined from the pH dependence of self-cleavage kinetics. Thus, the increase in activity with increasing pH does not correlate with deprotonation of G8, and most of G8 is protonated at neutral pH. These results do not exclude a role in proton transfer, but a simple interpretation is that G8 functions in the protonated form, perhaps by donating hydrogen bonds.

A catalytic metal ion interacts with the cleavage Site G.U wobble in the HDV ribozyme

The HDV ribozyme self-cleaves by a chemical mechanism involving general acid-base catalysis to generate 2',3'-cyclic phosphate and 5'-hydroxyl termini. Biochemical studies from several laboratories have implicated C75 as the general acid and hydrated magnesium as the general base. We have previously shown that C75 has a pK(a) shifted >2 pH units toward neutrality [Gong, B., Chen, J. H., Chase, E., Chadalavada, D. M., Yajima, R., Golden, B. L., Bevilacqua, P. C., and Carey, P. R. (2007) J. Am. Chem. Soc. 129, 13335-13342], while in crystal structures, it is well-positioned for proton transfer. However, no evidence for a hydrated magnesium poised to serve as a general base in the reaction has been observed in high-resolution crystal structures of various reaction states and mutants. Herein, we use solution kinetic experiments and parallel Raman crystallographic studies to examine the effects of pH on the rate and Mg(2+) binding properties of wild-type and 7-deazaguanosine mutants of the HDV ribozyme. These data suggest that a previously unobserved hydrated magnesium ion interacts with N7 of the cleavage site G.U wobble base pair. Integrating this metal ion binding site with the available crystal structures provides a new three-dimensional model for the active site of the ribozyme that accommodates all available biochemical data and appears competent for catalysis. The position of this metal is consistent with a role of a magnesium-bound hydroxide as a general base as dictated by biochemical data.

Unconventional decoding of the AUA codon as methionine by mitochondrial tRNAMet with the anticodon f5CAU as revealed with a mitochondrial in vitro translation system

Mitochondrial (mt) tRNA^Met has the unusual modified nucleotide 5-formylcytidine (f⁵C) in the first position of the anticodon. This tRNA must translate both AUG and AUA as methionine. By constructing an in vitro translation system from bovine liver mitochondria, we examined the decoding properties of the native mt tRNA^Met carrying f⁵C in the anticodon compared to a transcript that lacks the modification. The native mt Met-tRNA could recognize both AUA and AUG codons as Met, but the corresponding synthetic tRNA^Met lacking f⁵C (anticodon CAU), recognized only the AUG codon in both the codon-dependent ribosomal binding and in vitro translation assays. Furthermore, the Escherichia coli elongator tRNA^Met_m with the anticodon ac⁴CAU (ac⁴C = 4-acetylcytidine) and the bovine cytoplasmic initiator tRNA^Met (anticodon CAU) translated only the AUG codon for Met on mt ribosome. The codon recognition patterns of these tRNAs were the same on E. coli ribosomes. These results demonstrate that the f⁵C modification in mt tRNA^Met plays a crucial role in decoding the nonuniversal AUA codon as Met, and that the genetic code variation is compensated by a change in the tRNA anticodon, not by a change in the ribosome. Base pairing models of f⁵C-G and f⁵C-A based on the chemical properties of f⁵C are presented.

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

The hepatitis delta virus (HDV) ribozyme occurs in the genomic and antigenomic strands of the HDV RNA and within mammalian transcriptomes. Previous kinetic studies suggested that a wobble pair (G*U or A(+)*C) is preferred at the cleavage site; however, the reasons for this are unclear. We conducted sequence comparisons, which indicated that while G*U is the most prevalent combination at the cleavage site, G-C occurs to a significant extent in genomic HDV isolates, and G*U, G-C, and A-U pairs are present in mammalian ribozymes. We analyzed the folding of genomic HDV ribozymes by free energy minimization and found that variants with purine-pyrimidine combinations at the cleavage site are predicted to form native structures while pyrimidine-purine combinations misfold, consistent with earlier kinetic data and sequence comparisons. To test whether the cleavage site base pair contributes to catalysis, we characterized the pH and Mg(2+)-dependence of reaction kinetics of fast-folding genomic HDV ribozymes with cleavage site base pair purine-pyrimidine combinations: G*U, A-U, G-C, and A(+)*C. Rates for these native-folding ribozymes displayed highly similar pH and Mg(2+) concentration dependencies, with the exception of the A(+)*C ribozyme, which deviated at high pH. None of the four ribozymes underwent miscleavage. These observations support the A(+)*C ribozyme as being more active with a wobble pair at the cleavage site than with no base pair at all. Overall, the data support a model in which the cleavage site base pair provides a structural role in catalysis and does not need to be a wobble pair.

Structures, kinetics, thermodynamics, and biological functions of RNA hairpins

Most RNA comprises one strand and therefore can fold back on itself to form complex structures. At the heart of these structures is the hairpin, which is composed of a stem having Watson-Crick base pairing and a loop wherein the backbone changes directionality. First, we review the structure of hairpins including diversity in the stem, loop, and closing base pair. The function of RNA hairpins in biology is discussed next, including roles for isolated hairpins, as well as hairpins in the context of complex tertiary structures. We describe the kinetics and thermodynamics of hairpin folding including models for hairpin folding, folding transition states, and the cooperativity of folding. Lastly, we discuss some ways in which hairpins can influence the folding and function of tertiary structures, both directly and indirectly. RNA hairpins provide a simple means of controlling gene expression that can be understood in the language of physical chemistry.

Recent advances in RNA structure determination by NMR

Despite recent advances in the solution of NMR structures of RNA and RNA-ligand complexes, the rate limiting step remains the gathering of a large number of NOE and torsion restraints. Additional sources of information for structure determination of larger RNA molecules have recently become available, and it is possible to supplement NOE and J-coupling data with the measurement of dipolar couplings and cross-correlated relaxation rates in high-resolution NMR spectroscopy.

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.