The global structure of the VS ribozyme

Impact of Ion-Mixing Entropy on Orientational Preferences of DNA Helices: FRET Measurements and Computer Simulations

Biological processes require DNA and RNA helices to pack together in specific interhelical orientations. While electrostatic repulsion between backbone charges is expected to be maximized when helices are in parallel alignment, such orientations are commonplace in nature. To better understand how the repulsion is overcome, we used experimental and computational approaches to investigate how the orientational preferences of DNA helices depend on the concentration and valence of mobile cations. We used Förster resonance energy transfer (FRET) to probe the relative orientations of two 24-bp helices held together via a freely rotating PEG linker. At low cation concentrations, the helices preferred more "cross"-like orientations over those closer to parallel, and this preference was reduced with increasing salt concentrations. The results were in good quantitative agreement with Poisson-Boltzmann (PB) calculations for monovalent salt (Na+). However, PB underestimated the ability of mixtures of monovalent and divalent ions (Mg2+) to reduce the conformational preference. As a complementary approach, we performed all-atom molecular dynamics (MD) simulations and found better agreement with the experimental results. While MD and PB predict similar electrostatic forces, MD predicts a greater accumulation of Mg2+ in the ion atmosphere surrounding the DNA. Mg2+ occupancy is predicted to be greater in conformations close to the parallel orientation than in conformations close to the crossed orientation, enabling a greater release of Na+ ions and providing an entropic gain (one bound ion for two released). MD predicts an entropy gain larger than that of PB because of the increased Mg2+ occupancy. The entropy changes have a negligible effect at low Mg2+ concentrations because the free energies are dominated by electrostatic repulsion. However, as the Mg2+ concentration increases, charge screening is more effective and the mixing entropy produces readily detectable changes in packing preferences. Our results underline the importance of mixing entropy of counterions in nucleic acid interactions and provide a new understanding on the impact of a mixed ion atmosphere on the packing of DNA helices.

Revolving hexameric ATPases as asymmetric motors to translocate double-stranded DNA genome along one strand

DsDNA translocation through nanoscale pores is generally accomplished by ATPase biomotors. The discovery of the revolving dsDNA translocation mechanism, as opposed to rotation, in bacteriophage phi29 elucidated how ATPase motors move dsDNA. Revolution-driven, hexameric dsDNA motors have been reported in herpesvirus, bacterial FtsK, Streptomyces TraB, and T7 phage. This review explores the common relationship between their structure and mechanisms. Commonalities include moving along the 5'→3' strand, inchworm sequential action leading to an asymmetrical structure, channel chirality, channel size, and 3-step channel gating for controlling motion direction. The revolving mechanism and contact with one of the dsDNA strands addresses the historic controversy of dsDNA packaging using nicked, gapped, hybrid, or chemically modified DNA. These controversies surrounding dsDNA packaging activity using modified materials can be answered by whether the modification was introduced into the 3'→5' or 5'→3' strand. Perspectives concerning solutions to the controversy of motor structure and stoichiometry are also discussed.

An integrative NMR-SAXS approach for structural determination of large RNAs defines the substrate-free state of a trans-cleaving Neurospora Varkud Satellite ribozyme

The divide-and-conquer strategy is commonly used for protein structure determination, but its applications to high-resolution structure determination of RNAs have been limited. Here, we introduce an integrative approach based on the divide-and-conquer strategy that was undertaken to determine the solution structure of an RNA model system, the Neurospora VS ribozyme. NMR and SAXS studies were conducted on a minimal trans VS ribozyme as well as several isolated subdomains. A multi-step procedure was used for structure determination that first involved pairing refined NMR structures with SAXS data to obtain structural subensembles of the various subdomains. These subdomain structures were then assembled to build a large set of structural models of the ribozyme, which was subsequently filtered using SAXS data. The resulting NMR-SAXS structural ensemble shares several similarities with the reported crystal structures of the VS ribozyme. However, a local structural difference is observed that affects the global fold by shifting the relative orientation of the two three-way junctions. Thus, this finding highlights a global conformational change associated with substrate binding in the VS ribozyme that is likely critical for its enzymatic activity. Structural studies of other large RNAs should benefit from similar integrative approaches that allow conformational sampling of assembled fragments.

Self-cleaving ribozymes: substrate specificity and synthetic biology applications

Various self-cleaving ribozymes appearing in nature catalyze the sequence-specific intramolecular cleavage of RNA and can be engineered to catalyze cleavage of appropriate substrates in an intermolecular fashion, thus acting as true catalysts. The mechanisms of the small, self-cleaving ribozymes have been extensively studied and reviewed previously. Self-cleaving ribozymes can possess high catalytic activity and high substrate specificity; however, substrate specificity is also engineerable within the constraints of the ribozyme structure. While these ribozymes share a common fundamental catalytic mechanism, each ribozyme family has a unique overall architecture and active site organization, indicating that several distinct structures yield this chemical activity. The multitude of catalytic structures, combined with some flexibility in substrate specificity within each family, suggests that such catalytic RNAs, taken together, could access a wide variety of substrates. Here, we give an overview of 10 classes of self-cleaving ribozymes and capture what is understood about their substrate specificity and synthetic applications. Evolution of these ribozymes in an RNA world might be characterized by the emergence of a new ribozyme family followed by rapid adaptation or diversification for specific substrates.

Self-cleaving ribozymes have become important tools of synthetic biology. Here we summarize the substrate specificity and applications of the main classes of these ribozymes.

A structural intermediate pre-organizes the add adenine riboswitch for ligand recognition

Riboswitches are RNA sequences that regulate gene expression by undergoing structural changes upon the specific binding of cellular metabolites. Crystal structures of purine-sensing riboswitches have revealed an intricate network of interactions surrounding the ligand in the bound complex. The mechanistic details about how the aptamer folding pathway is involved in the formation of the metabolite binding site have been previously shown to be highly important for the riboswitch regulatory activity. Here, a combination of single-molecule FRET and SHAPE assays have been used to characterize the folding pathway of the adenine riboswitch from Vibrio vulnificus. Experimental evidences suggest a folding process characterized by the presence of a structural intermediate involved in ligand recognition. This intermediate state acts as an open conformation to ensure ligand accessibility to the aptamer and folds into a structure nearly identical to the ligand-bound complex through a series of structural changes. This study demonstrates that the add riboswitch relies on the folding of a structural intermediate that pre-organizes the aptamer global structure and the ligand binding site to allow efficient metabolite sensing and riboswitch genetic regulation.

Confluence of theory and experiment reveals the catalytic mechanism of the Varkud satellite ribozyme

The Varkud satellite ribozyme catalyses site-specific RNA cleavage and ligation, and serves as an important model system to understand RNA catalysis. Here, we combine stereospecific phosphorothioate substitution, precision nucleobase mutation and linear free-energy relationship measurements with molecular dynamics, molecular solvation theory and ab initio quantum mechanical/molecular mechanical free-energy simulations to gain insight into the catalysis. Through this confluence of theory and experiment, we unify the existing body of structural and functional data to unveil the catalytic mechanism in unprecedented detail, including the degree of proton transfer in the transition state. Further, we provide evidence for a critical Mg²⁺ in the active site that interacts with the scissile phosphate and anchors the general base guanine in position for nucleophile activation. This novel role for Mg²⁺ adds to the diversity of known catalytic RNA strategies and unifies functional features observed in the Varkud satellite, hairpin and hammerhead ribozyme classes.

Ligand recognition and helical stacking formation are intimately linked in the SAM-I riboswitch regulatory mechanism

Riboswitches are noncoding mRNA elements that control gene expression by altering their structure upon metabolite binding. Although riboswitch crystal structures provide detailed information about RNA-ligand interactions, little knowledge has been gathered to understand how riboswitches modulate gene expression. Here, we study the molecular recognition mechanism of the S-adenosylmethionine SAM-I riboswitch by characterizing the formation of a helical stacking interaction involving the ligand-binding process. We show that ligand binding is intimately linked to the formation of the helical stacking, which is dependent on the presence of three conserved purine residues that are flanked by stacked helices. We also find that these residues are important for the formation of a crucial long-range base pair formed upon SAM binding. Together, our results lend strong support to a critical role for helical stacking in the folding pathway and suggest a particularly important function in the formation of the long-range base pair.

Mg2+ Binding Promotes SLV as a Scaffold in Varkud Satellite Ribozyme SLI-SLV Kissing Loop Junction

Though the structure of the substrate stem loop I (SLI)-stem loop V (SLV) kissing loop junction of the Varkud Satellite ribozyme has been experimentally characterized, the dynamics of this Mg²⁺-dependent loop-loop interaction have been elusive. Specifically, each hairpin loop contains a U-turn motif, but only SLV shows a conformational shift triggered by Mg²⁺ ion association. Here, we use molecular dynamics simulations to analyze the binding and dynamics of this kissing loop junction. We show that SLV acts as a scaffold, providing stability to the junction. Mg²⁺ ions associate with SLV when it is part of the junction in a manner similar to when it is unbound, but there is no specificity in Mg²⁺ binding for the SLI loop. This suggests that the entropic penalty of ordering the larger SLI is too high, allowing SLV to act as a scaffold for multiple substrate loop sequences.

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.

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.

The NMR structure of the II-III-VI three-way junction from the Neurospora VS ribozyme reveals a critical tertiary interaction and provides new insights into the global ribozyme structure

As part of an effort to structurally characterize the complete Neurospora VS ribozyme, NMR solution structures of several subdomains have been previously determined, including the internal loops of domains I and VI, the I/V kissing-loop interaction and the III-IV-V junction. Here, we expand this work by determining the NMR structure of a 62-nucleotide RNA (J236) that encompasses the VS ribozyme II-III-VI three-way junction and its adjoining stems. In addition, we localize Mg(2+)-binding sites within this structure using Mn(2+)-induced paramagnetic relaxation enhancement. The NMR structure of the J236 RNA displays a family C topology with a compact core stabilized by continuous stacking of stems II and III, a cis WC/WC G•A base pair, two base triples and two Mg(2+) ions. Moreover, it reveals a remote tertiary interaction between the adenine bulges of stems II and VI. Additional NMR studies demonstrate that both this bulge-bulge interaction and Mg(2+) ions are critical for the stable folding of the II-III-VI junction. The NMR structure of the J236 RNA is consistent with biochemical studies on the complete VS ribozyme, but not with biophysical studies performed with a minimal II-III-VI junction that does not contain the II-VI bulge-bulge interaction. Together with previous NMR studies, our findings provide important new insights into the three-dimensional architecture of this unique ribozyme.

Fitness Landscapes of Functional RNAs

The notion of fitness landscapes, a map between genotype and fitness, was proposed more than 80 years ago. For most of this time data was only available for a few alleles, and thus we had only a restricted view of the whole fitness landscape. Recently, advances in genetics and molecular biology allow a more detailed view of them. Here we review experimental and theoretical studies of fitness landscapes of functional RNAs, especially aptamers and ribozymes. We find that RNA structures can be divided into critical structures, connecting structures, neutral structures and forbidden structures. Such characterisation, coupled with theoretical sequence-to-structure predictions, allows us to construct the whole fitness landscape. Fitness landscapes then can be used to study evolution, and in our case the development of the RNA world.

Nuclear magnetic resonance structure of the III-IV-V three-way junction from the Varkud satellite ribozyme and identification of magnesium-binding sites using paramagnetic relaxation enhancement

The VS ribozyme is a catalytic RNA found within some natural isolates of Neurospora that is being used as a model system to improve our understanding of RNA structure, catalysis, and engineering. The catalytic domain contains five helical domains (SLII-SLVI) that are organized by two three-way junctions. The III-IV-V junction is required for high-affinity binding of the substrate domain (SLI) through formation of a kissing loop interaction with SLV. Here, we determine the high-resolution nuclear magnetic resonance (NMR) structure of a 47-nucleotide RNA containing the III-IV-V junction (J345). The J345 RNA adopts a Y-shaped fold typical of the family C three-way junctions, with coaxial stacking between stems III and IV and an acute angle between stems III and V. The NMR structure reveals that the core of the III-IV-V junction contains four stacked base triples, a U-turn motif, a cross-strand stacking interaction, an A-minor interaction, and a ribose zipper. In addition, the NMR structure shows that the cCUUGg tetraloop used to stabilize stem IV adopts a novel RNA tetraloop fold, different from the known gCUUGc tetraloop structure. Using Mn(2+)-induced paramagnetic relaxation enhancement, we identify six Mg(2+)-binding sites within J345, including one associated with the cCUUGg tetraloop and two with the junction core. The NMR structure of J345 likely represents the conformation of the III-IV-V junction in the context of the active VS ribozyme and suggests that this junction functions as a dynamic hinge that contributes to substrate recognition and catalysis. Moreover, this study highlights a new role for family C three-way junctions in long-range tertiary interactions.

A remarkably stable kissing-loop interaction defines substrate recognition by the Neurospora Varkud Satellite ribozyme

Kissing loops are tertiary structure elements that often play key roles in functional RNAs. In the Neurospora VS ribozyme, a kissing-loop interaction between the stem-loop I (SLI) substrate and stem-loop V (SLV) of the catalytic domain is known to play an important role in substrate recognition. In addition, this I/V kissing-loop interaction is associated with a helix shift in SLI that activates the substrate for catalysis. To better understand the role of this kissing-loop interaction in substrate recognition and activation by the VS ribozyme, we performed a thermodynamic characterization by isothermal titration calorimetry using isolated SLI and SLV stem-loops. We demonstrate that preshifted SLI variants have higher affinity for SLV than shiftable SLI variants, with an energetic cost of 1.8-3 kcal/mol for the helix shift in SLI. The affinity of the preshifted SLI for SLV is remarkably high, the interaction being more stable by 7-8 kcal/mol than predicted for a comparable duplex containing three Watson-Crick base pairs. The structural basis of this remarkable stability is discussed in light of previous NMR studies. Comparative thermodynamic studies reveal that kissing-loop complexes containing 6-7 Watson-Crick base pairs are as stable as predicted from comparable RNA duplexes; however, those with 2-3 Watson-Crick base pairs are more stable than predicted. Interestingly, the stability of SLI/ribozyme complexes is similar to that of SLI/SLV complexes. Thus, the I/V kissing loop interaction represents the predominant energetic contribution to substrate recognition by the trans-cleaving VS ribozyme.

RNA conformational changes analyzed by comparative gel electrophoresis

The study of biologically relevant native RNA structures is important to understand their cellular function(s). Native gel electrophoresis provides information about such native structures in solution as a function of experimental conditions. The application of native gel electrophoresis in a comparative manner allows to obtain precise information on relative angles subtended between given pair of stems in an RNA molecule. By adapting this approach, it is possible to obtain very specific structural information such as the amplitude of dihedral angles and helical rotation. As an example, we will describe how native gel electrophoresis can be used to study the folding of the S-adenosylmethionine (SAM) sensing riboswitch.

Mechanisms of RNA catalysis

Ribozymes are RNA molecules that act as chemical catalysts. In contemporary cells, most known ribozymes carry out phosphoryl transfer reactions. The nucleolytic ribozymes comprise a class of five structurally-distinct species that bring about site-specific cleavage by nucleophilic attack of the 2'-O on the adjacent 3'-P to form a cyclic 2',3'-phosphate. In general, they will also catalyse the reverse reaction. As a class, all these ribozymes appear to use general acid-base catalysis to accelerate these reactions by about a million-fold. In the Varkud satellite ribozyme, we have shown that the cleavage reaction is catalysed by guanine and adenine nucleobases acting as general base and acid, respectively. The hairpin ribozyme most probably uses a closely similar mechanism. Guanine nucleobases appear to be a common choice of general base, but the general acid is more variable. By contrast, the larger ribozymes such as the self-splicing introns and RNase P act as metalloenzymes.

Helix-length compensation studies reveal the adaptability of the VS ribozyme architecture

Compensatory mutations in RNA are generally regarded as those that maintain base pairing, and their identification forms the basis of phylogenetic predictions of RNA secondary structure. However, other types of compensatory mutations can provide higher-order structural and evolutionary information. Here, we present a helix-length compensation study for investigating structure–function relationships in RNA. The approach is demonstrated for stem-loop I and stem-loop V of the Neurospora VS ribozyme, which form a kissing–loop interaction important for substrate recognition. To rapidly characterize the substrate specificity (k_cat/K_M) of several substrate/ribozyme pairs, a procedure was established for simultaneous kinetic characterization of multiple substrates. Several active substrate/ribozyme pairs were identified, indicating the presence of limited substrate promiscuity for stem Ib variants and helix-length compensation between stems Ib and V. 3D models of the I/V interaction were generated that are compatible with the kinetic data. These models further illustrate the adaptability of the VS ribozyme architecture for substrate cleavage and provide global structural information on the I/V kissing–loop interaction. By exploring higher-order compensatory mutations in RNA our approach brings a deeper understanding of the adaptability of RNA structure, while opening new avenues for RNA research.

Metal-ion-dependent folding of a uranyl-specific DNAzyme: insight into function from fluorescence resonance energy transfer studies

Fluorescence resonance energy transfer (FRET) has been used to study the global folding of an uranyl (UO(2)(2+))-specific 39E DNAzyme in the presence of Mg(2+), Zn(2+), Pb(2+), or UO(2)(2+). At pH 5.5 and physiological ionic strength (100 mM Na(+)), two of the three stems in this DNAzyme folded into a compact structure in the presence of Mg(2+) or Zn(2+). However, no folding occurred in the presence of Pb(2+) or UO(2)(2+); this is analogous to the "lock-and-key" catalysis mode first observed in the Pb(2+)-specific 8-17 DNAzyme. However, Mg(2+) and Zn(2+) exert different effects on the 8-17 and 39E DNAzymes. Whereas Mg(2+) or Zn(2+)-dependent folding promoted 8-17 DNAzyme activity, the 39E DNAzyme folding induced by Mg(2+) or Zn(2+) inhibited UO(2)(2+)-specific activity. Group IIA series of metal ions (Mg(2+), Ca(2+), Sr(2+)) also caused global folding of the 39E DNAzyme, for which the apparent binding affinity between these metal ions and the DNAzyme decreases as the ionic radius of the metal ions increases. Because the ionic radius of Sr(2+) (1.12 Å) is comparable to that of Pb(2+) (1.20 Å), but contrary to Pb(2+), Sr(2+) induces the DNAzyme to fold under identical conditions, ionic size alone cannot account for the unique folding behaviors induced by Pb(2+) and UO(2)(2+). Under low ionic strength (30 mM Na(+)), all four metal ions (Mg(2+), Zn(2+), Pb(2+), and UO(2)(2+)), caused 39E DNAzyme folding, suggesting that metal ions can neutralize the negative charge of DNA-backbone phosphates in addition to playing specific catalytic roles. Mg(2+) at low (<2 mM) concentration promoted UO(2)(2+)-specific activity, whereas Mg(2+) at high (>2 mM) concentration inhibited the UO(2)(2+)-specific activity. Therefore, the lock-and-key mode of DNAzymes depends on ionic strength, and the 39E DNAzyme is in the lock-and-key mode only at ionic strengths of 100 mM or greater.

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.

Constitutive regulatory activity of an evolutionarily excluded riboswitch variant

The exquisite specificity of the adenine-responsive riboswitch toward its cognate metabolite has been shown to arise from the formation of a Watson-Crick interaction between the adenine ligand and residue U65. A recent crystal structure of a U65C adenine aptamer variant has provided a rationale for the phylogenetic conservation observed at position 39 for purine aptamers. The G39-C65 variant adopts a compact ligand-free structure in which G39 is accommodated by the ligand binding site and is base-paired to the cytosine at position 65. Here, we demonstrate using a combination of biochemical and biophysical techniques that the G39-C65 base pair not only severely impairs ligand binding but also disrupts the functioning of the riboswitch in vivo by constitutively activating gene expression. Folding studies using single-molecule FRET revealed that the G39-C65 variant displays a low level of dynamic heterogeneity, a feature reminiscent of ligand-bound wild-type complexes. A restricted conformational freedom together with an ability to significantly fold in monovalent ions are exclusive to the G39-C65 variant. This work provides a mechanistic framework to rationalize the evolutionary exclusion of certain nucleotide combinations in favor of sequences that preserve ligand binding and gene regulation functionalities.

Do the hairpin and VS ribozymes share a common catalytic mechanism based on general acid-base catalysis? A critical assessment of available experimental data

The active centers of the hairpin and VS ribozymes are both generated by the interaction of two internal loops, and both ribozymes use guanine and adenine nucleobases to accelerate cleavage and ligation reactions. The centers are topologically equivalent and the relative positioning of key elements the same. There is good evidence that the cleavage reaction of the VS ribozyme is catalyzed by the guanine (G638) acting as general base and the adenine (A756) as general acid. We now critically evaluate the experimental mechanistic evidence for the hairpin ribozyme. We conclude that all the available data are fully consistent with a major contribution to catalysis by general acid-base catalysis involving the adenine (A38) and guanine (G8). It appears that the two ribozymes are mechanistically equivalent.

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.

Metal ions: supporting actors in the playbook of small ribozymes

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

Folding of the lysine riboswitch: importance of peripheral elements for transcriptional regulation

The Bacillus subtilis lysC lysine riboswitch modulates its own gene expression upon lysine binding through a transcription attenuation mechanism. The riboswitch aptamer is organized around a single five-way junction that provides the scaffold for two long-range tertiary interactions (loop L2–loop L3 and helix P2–loop L4)—all of this for the creation of a specific lysine binding site. We have determined that the interaction P2–L4 is particularly important for the organization of the ligand-binding site and for the riboswitch transcription attenuation control. Moreover, we have observed that a folding synergy between L2–L3 and P2–L4 allows both interactions to fold at lower magnesium ion concentrations. The P2–L4 interaction is also critical for the close juxtaposition involving stems P1 and P5. This is facilitated by the presence of lysine, suggesting an active role of the ligand in the folding transition. We also show that a previously uncharacterized stem–loop located in the expression platform is highly important for the riboswitch activity. Thus, folding elements located in the aptamer and the expression platform both influence the lysine riboswitch gene regulation.

Structure of the three-way helical junction of the hepatitis C virus IRES element

The hepatitis C virus internal ribosome entry site (IRES) element contains a three-way junction that is important in the overall RNA conformation, and for its role in the internal initiation of translation. The junction also illustrates some important conformational principles in the folding of three-way helical junctions. It is formally a 3HS(4) junction, with the possibility of two alternative stacking conformers. However, in principle, the junction can also undergo two steps of branch migration that would form 2HS(1)HS(3) and 2HS(2)HS(2) junctions. Comparative gel electrophoresis and ensemble fluorescence resonance energy transfer (FRET) studies show that the junction is induced to fold by the presence of Mg(2+) ions in low micromolar concentrations, and suggest that the structure adopted is based on coaxial stacking of the two helices that do not terminate in a hairpin loop (i.e., helix IIId). Single-molecule FRET studies confirm this conclusion, and indicate that there is no minor conformer present based on an alternative choice of helical stacking partners. Moreover, analysis of single-molecule FRET data at an 8-msec resolution failed to reveal evidence for structural transitions. It seems probable that this junction adopts a single conformation as a unique and stable fold.

The ligand-free state of the TPP riboswitch: a partially folded RNA structure

Riboswitches are elements of mRNA that regulate gene expression by undergoing structural changes upon binding of small ligands. Although the structures of several riboswitches have been solved with their ligands bound, the ligand-free states of only a few riboswitches have been characterized. The ligand-free state is as important for the functionality of the riboswitch as the ligand-bound form, but the ligand-free state is often a partially folded structure of the RNA, with conformational heterogeneity that makes it particularly challenging to study. Here, we present models of the ligand-free state of a thiamine pyrophosphate riboswitch that are derived from a combination of complementary experimental and computational modeling approaches. We obtain a global picture of the molecule using small-angle X-ray scattering data and use an RNA structure modeling software, MC-Sym, to fit local structural details to these data on an atomic scale. We have used two different approaches to obtaining these models. Our first approach develops a model of the RNA from the structures of its constituent junction fragments in isolation. The second approach treats the RNA as a single entity, without bias from the structure of its individual constituents. We find that both approaches give similar models for the ligand-free form, but the ligand-bound models differ for the two approaches, and only the models from the second approach agree with the ligand-bound structure known previously from X-ray crystallography. Our models provide a picture of the conformational changes that may occur in the riboswitch upon binding of its ligand. Our results also demonstrate the power of combining experimental small-angle X-ray scattering data with theoretical structure prediction tools in the determination of RNA structures beyond riboswitches.

Comparative gel electrophoresis analysis of helical junctions in RNA

Comparative gel electrophoresis provides information on the relative angles subtended between helical arms at a branchpoint in RNA. It is based upon the comparison of electrophoretic mobility in polyacrylamide gels of species containing two long arms, with the remaining one(s) being significantly shorter. Although the method currently lacks a really well-established basis of physical theory, it is very powerful, yet simple to apply. It has had a number of significant successes in RNA, DNA and DNA-protein complexes, and in all cases to date the results have stood the test of time and eventual comparison with crystallographic analysis.

Analysis of RNA folding by native polyacrylamide gel electrophoresis

Polyacrylamide gel electrophoresis under native conditions (native PAGE) is a well-established and versatile method for probing nucleic acid conformation and nucleic acid-protein interactions. Native PAGE has been used to measure RNA folding equilibria and kinetics under a wide variety of conditions. Advantages of this method are its adaptability, absolute determination of reaction endpoints, and direct analysis of conformational hetereogeneity within a sample. Native PAGE is also useful for resolving ligand-induced structural changes.

Formation of an active site in trans by interaction of two complete Varkud Satellite ribozymes

The complete VS ribozyme comprises seven helical segments, connected by three three-way RNA junctions. In the presence of Mg(2+) ions, cleavage occurs within the internal loop of helix I. This requires the participation of a guanine (G638) within the helix I loop, and a remote adenine (A756) within an internal loop of helix VI. Previous structural studies have suggested that helix I docks into the fold of the remaining part of the ribozyme, bringing A756 and G638 close to the scissile phosphate to allow the cleavage reaction to proceed. We show here that while either A756C or G638A individually exhibit very low cleavage activity, a mixture of the two variants leads to cleavage of the A756C RNA, but not the G638A RNA. The rate of cleavage depends on the concentration of the VS G638A RNA, as expected for a bimolecular interaction. This regaining of cleavage activity by complementation indicates that helix I of one VS RNA can interact with another VS RNA molecule to generate a functional active site in trans.

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.

The structure and folding of branched RNA analyzed by fluorescence resonance energy transfer

Fluorescence resonance energy transfer (FRET) is a spectroscopic means of obtaining distance information over a range up to ~80Å in solution. It is based on the dipolar coupling between the electronic transition moments of a donor and acceptor fluorophore attached at known positions on the RNA species of interest. It can be applied in ensembles of molecules, either by steady-state fluorescence or by lifetime measurements, but it is also very appropriate for single-molecule studies. In addition to the provision of distance information, recent studies have emphasized the orientation dependence of energy transfer.

Application of fluorescent measurements for characterization of riboswitch-ligand interactions

Riboswitches are recently discovered messenger RNA motifs involved in gene regulation. They modulate gene expression at various levels, such as transcription, translation, splicing, and mRNA degradation. Because riboswitches exhibit relatively complex structures, they are able to form highly complex ligand-binding sites, which enable the specific recognition of target metabolites in a complex cellular environment. Practically in all studied cases, riboswitches use ligand-induced conformational changes to control gene expression. To monitor the structural reorganization of riboswitches, we use the local fluorescent reporter 2-aminopurine (2AP), which is a structural analog of adenine. The 2AP fluorescence is strongly quenched when the fluorophore is involved in stacking interactions with surrounding bases, and can, therefore, be used to monitor local structural rearrangements. Here, we show specific examples in which 2AP fluorescence can be used to monitor structural changes in the aptamer domain of the S-adenosyl methionine (SAM) riboswitch and where it can be used as a ligand for the guanine riboswitch.

Analysis of branched nucleic acid structure using comparative gel electrophoresis

Electrophoresis in polyacrylamide gels provides a simple yet powerful means of analyzing the relative disposition of helical arms in branched nucleic acids. The electrophoretic mobility of DNA or RNA with a central discontinuity is determined by the angle subtended between the arms radiating from the branchpoint. In a multi-helical branchpoint, comparative gel electrophoresis can provide a relative measure of all the inter-helical angles and thus the shape and symmetry of the molecule. Using the long-short arm approach, the electrophoretic mobility of all the species with two helical arms that are longer than all others is compared. This can be done as a function of conditions, allowing the analysis of ion-dependent folding of branched DNA and RNA species. Notable successes for the technique include the four-way (Holliday) junction in DNA and helical junctions in functionally significant RNA species such as ribozymes. Many of these structures have subsequently been proved correct by crystallography or other methods, up to 10 years later in the case of the Holliday junction. Just as important, the technique has not failed to date. Comparative gel electrophoresis can provide a window on both fast and slow conformational equilibria such as conformer exchange in four-way DNA junctions. But perhaps the biggest test of the approach has been to deduce the structures of complexes of four-way DNA junctions with proteins. Two recent crystallographic structures show that the global structures were correctly deduced by electrophoresis, proving the worth of the method even in these rather complex systems. Comparative gel electrophoresis is a robust method for the analysis of branched nucleic acids and their complexes.

The complete VS ribozyme in solution studied by small-angle X-ray scattering

We have used small-angle X-ray solution scattering to obtain ab initio shape reconstructions of the complete VS ribozyme. The ribozyme occupies an electron density envelope with an irregular shape, into which helical sections have been fitted. The ribozyme is built around a core comprising a near-coaxial stack of three helices, organized by two three-way helical junctions. An additional three-way junction formed by an auxiliary helix directs the substrate stem-loop, juxtaposing the cleavage site with an internal loop to create the active complex. This is consistent with the current view of the probable mechanism of trans-esterification in which adenine and guanine nucleobases contributed by the interacting loops combine in general acid-base catalysis.

Single VS ribozyme molecules reveal dynamic and hierarchical folding toward catalysis

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

Probing RNA structural dynamics and function by fluorescence resonance energy transfer (FRET)

Biological function of RNA is often mediated by cyclic switching between several (meta-)stable arrangements of tertiary structure. Fluorophore labeling of RNA offers a unique view into these folding and conformational switching events, since a fluorescence signal is sensitive to its molecular environment and can be continuously monitored in real time to produce kinetic rate information. This unit focuses on the practical implications of using fluorescence resonance energy transfer (FRET) to probe RNA structural dynamics and function. FRET is a particularly powerful fluorescence technique since, in addition to kinetic data, it provides insights into the structural basis of a conformational rearrangement. Protocols describe how to postsynthetically label RNA for FRET and how to acquire and analyze FRET data. Support protocols describe methods for deprotecting synthetic RNA and for purifying RNA by gel electrophoresis and HPLC. Considerations for selecting appropriate RNA, fluorophores, and labeling strategies are discussed in detail in the commentary.

Role of SLV in SLI substrate recognition by the Neurospora VS ribozyme

Substrate recognition by the VS ribozyme involves a magnesium-dependent loop/loop interaction between the SLI substrate and the SLV hairpin from the catalytic domain. Recent NMR studies of SLV demonstrated that magnesium ions stabilize a U-turn loop structure and trigger a conformational change for the extruded loop residue U700, suggesting a role for U700 in SLI recognition. Here, we kinetically characterized VS ribozyme mutants to evaluate the contribution of U700 and other SLV loop residues to SLI recognition. To help interpret the kinetic data, we structurally characterized the SLV mutants by NMR spectroscopy and generated a three-dimensional model of the SLI/SLV complex by homology modeling with MC-Sym. We demonstrated that the mutation of U700 by A, C, or G does not significantly affect ribozyme activity, whereas deletion of U700 dramatically impairs this activity. The U700 backbone is likely important for SLI recognition, but does not appear to be required for either the structural integrity of the SLV loop or for direct interactions with SLI. Thus, deletion of U700 may affect other aspects of SLI recognition, such as magnesium ion binding and SLV loop dynamics. As part of our NMR studies, we developed a convenient assay based on detection of unusual (31)P and (15)N N7 chemical shifts to probe the formation of U-turn structures in RNAs. Our model of the SLI/SLV complex, which is compatible with biochemical data, leads us to propose novel interactions at the loop I/loop V interface.

Structural inference of native and partially folded RNA by high-throughput contact mapping

The biological behaviors of ribozymes, riboswitches, and numerous other functional RNA molecules are critically dependent on their tertiary folding and their ability to sample multiple functional states. The conformational heterogeneity and partially folded nature of most of these states has rendered their characterization by high-resolution structural approaches difficult or even intractable. Here we introduce a method to rapidly infer the tertiary helical arrangements of large RNA molecules in their native and non-native solution states. Multiplexed hydroxyl radical (.OH) cleavage analysis (MOHCA) enables the high-throughput detection of numerous pairs of contacting residues via random incorporation of radical cleavage agents followed by two-dimensional gel electrophoresis. We validated this technology by recapitulating the unfolded and native states of a well studied model RNA, the P4-P6 domain of the Tetrahymena ribozyme, at subhelical resolution. We then applied MOHCA to a recently discovered third state of the P4-P6 RNA that is stabilized by high concentrations of monovalent salt and whose partial order precludes conventional techniques for structure determination. The three-dimensional portrait of a compact, non-native RNA state reveals a well ordered subset of native tertiary contacts, in contrast to the dynamic but otherwise similar molten globule states of proteins. With its applicability to nearly any solution state, we expect MOHCA to be a powerful tool for illuminating the many functional structures of large RNA molecules and RNA/protein complexes.

The tolerance to exchanges of the Watson Crick base pair in the hammerhead ribozyme core is determined by surrounding elements

Tertiary interacting elements are important features of functional RNA molecules, for example, in all small nucleolytic ribozymes. The recent crystal structure of a tertiary stabilized type I hammerhead ribozyme revealed a conventional Watson-Crick base pair in the catalytic core, formed between nucleotides C3 and G8. We show that any Watson-Crick base pair between these positions retains cleavage competence in two type III ribozymes. In the Arabidopsis thaliana sequence, only moderate differences in cleavage rates are observed for the different base pairs, while the peach latent mosaic viroid (PLMVd) ribozyme exhibits a preference for a pyrimidine at position 3 and a purine at position 8. To understand these differences, we created a series of chimeric ribozymes in which we swapped sequence elements that surround the catalytic core. The kinetic characterization of the resulting ribozymes revealed that the tertiary interacting loop sequences of the PLMVd ribozyme are sufficient to induce the preference for Y3-R8 base pairs in the A. thaliana hammerhead ribozyme. In contrast to this, only when the entire stem-loops I and II of the A. thaliana sequences are grafted on the PLMVd ribozyme is any Watson-Crick base pair similarly tolerated. The data provide evidence for a complex interplay of secondary and tertiary structure elements that lead, mediated by long-range effects, to an individual modulation of the local structure in the catalytic core of different hammerhead ribozymes.

A chemo-genetic approach for the study of nucleobase participation in nucleolytic ribozymes

A novel chemo-genetic approach for the analysis of general acid-base catalysis by nucleobases in ribozymes is reviewed. This involves substitution of a C-nucleoside with imidazole in place of a natural nucleobase. The Varkud satellite ribozyme in which the nucleobase at the critical 756 position has been replaced by imidazole is active in both cleavage and ligation reactions. Similarly, a modified hairpin ribozyme with the nucleobase at position 8 substituted by imidazole is active in cleavage and ligation reactions. Although the rates are lower than those of the natural ribozymes, they are significantly greater than other variants at these positions. The dependence of the hairpin ribozyme reaction rates on pH has been studied. Both cleavage and ligation reactions display a bell-shaped pH dependence, consistent with general acid-base catalysis involving the nucleotide at position 8.

A loop loop interaction and a K-turn motif located in the lysine aptamer domain are important for the riboswitch gene regulation control

The lysine riboswitch is associated to the lysC gene in Bacillus subtilis, and the binding of lysine modulates the RNA structure to allow the formation of an intrinsic terminator presumably involved in transcription attenuation. The complex secondary structure of the lysine riboswitch aptamer is organized around a five-way junction that undergoes structural changes upon ligand binding. Using single-round transcription assays, we show that a loop-loop interaction is important for lysine-induced termination of transcription. Moreover, upon close inspection of the secondary structure, we find that an unconventional kink-turn motif is present in one of the stems participating in the loop-loop interaction. We show that the K-turn adopts a pronounced kink and that it binds the K-turn-binding protein L7Ae of Archaeoglobus fulgidus in the low nanomolar range. The functional importance of this K-turn motif is revealed from single-round transcription assays, which show its importance for efficient transcription termination. This motif is essential for the loop-loop interaction, and consequently, for lysine binding. Taken together, our results depict for the first time the importance of a K-turn-dependent loop-loop interaction for the transcription regulation of a lysine riboswitch.

A guanine nucleobase important for catalysis by the VS ribozyme

A guanine (G638) within the substrate loop of the VS ribozyme plays a critical role in the cleavage reaction. Replacement by any other nucleotide results in severe impairment of cleavage, yet folding of the substrate is not perturbed, and the variant substrates bind the ribozyme with similar affinity, acting as competitive inhibitors. Functional group substitution shows that the imino proton on the N1 is critical, suggesting a possible role in general acid-base catalysis, and this in accord with the pH dependence of the reaction rate for the natural and modified substrates. We propose a chemical mechanism for the ribozyme that involves general acid-base catalysis by the combination of the nucleobases of guanine 638 and adenine 756. This is closely similar to the probable mechanism of the hairpin ribozyme, and the active site arrangements for the two ribozymes appear topologically equivalent. This has probably arisen by convergent evolution.

Structural transitions and thermodynamics of a glycine-dependent riboswitch from Vibrio cholerae

Riboswitches are complex folded RNA domains found in noncoding regions of mRNA that regulate gene expression upon small molecule binding. Recently, Breaker and coworkers reported a tandem aptamer riboswitch (VCI-II) that binds glycine cooperatively. Here, we use hydroxyl radical footprinting and small-angle X-ray scattering (SAXS) to study the conformations of this tandem aptamer as a function of Mg(2+) and glycine concentration. We fit a simple three-state thermodynamic model that describes the energetic coupling between magnesium-induced folding and glycine binding. Furthermore, we characterize the structural conformations of each of the three states: In low salt with no magnesium present, the VCI-II construct has an extended overall conformation, presumably representing unfolded structures. Addition of millimolar concentrations of Mg(2+) in the absence of glycine leads to a significant compaction and partial folding as judged by hydroxyl radical protections. In the presence of millimolar Mg(2+) concentrations, the tandem aptamer binds glycine cooperatively. The glycine binding transition involves a further compaction, additional tertiary packing interactions and further uptake of magnesium ions relative to the state in high Mg(2+) but no glycine. Employing density reconstruction algorithms, we obtain low resolution 3-D structures for all three states from the SAXS measurements. These data provide a first glimpse into the structural conformations of the VCI-II aptamer, establish rigorous constraints for further modeling, and provide a framework for future mechanistic studies.

Structural methods for studying IRES function

Internal ribosome entry sites (IRESs) substitute RNA sequences for some or all of the canonical translation initiation protein factors. Therefore, an important component of understanding IRES function is a description of the three-dimensional structure of the IRES RNA underlying this mechanism. This includes determining the degree to which the RNA folds, the global RNA architecture, and higher resolution information when warranted. Knowledge of the RNA structural features guides ongoing mechanistic and functional studies. In this chapter, we present a roadmap to structurally characterize a folded RNA, beginning from initial studies to define the overall architecture and leading to high-resolution structural studies. The experimental strategy presented here is not unique to IRES RNAs but is adaptable to virtually any RNA of interest, although characterization of RNA-protein interactions requires additional methods. Because IRES RNAs have a specific function, we present specific ways in which the data are interpreted to gain insight into that function. We provide protocols for key experiments that are particularly useful for studying IRES RNA structure and that provide a framework onto which additional approaches are integrated. The protocols we present are solution hydroxyl radical probing, RNase T1 probing, native gel electrophoresis, sedimentation velocity analytical ultracentrifugation, and strategies to engineer RNA for crystallization and to obtain initial crystals.

Structural analysis of the catalytic core of human telomerase RNA by FRET and molecular modeling

Telomerase is the ribonucleoprotein reverse transcriptase involved in the maintenance of the telomeres, the termini of eukaryotic chromosomes. The RNA component of human telomerase (hTR) consists of 451 nucleotides with the 5' half folding into a highly conserved catalytic core comprising the template region and an adjacent pseudoknot domain (nucleotides 1-208). While the secondary structure of hTR is established, there is little understanding of its three-dimensional (3D) architecture. Here, we have used fluorescence resonance energy transfer (FRET) between fluorescently labelled peptide nucleic acids, hybridized to defined single stranded regions of full length hTR, to evaluate long-range distances. Using molecular modeling, the distance constraints derived by FRET were subsequently used, together with the known secondary structure, to generate a 3D model of the catalytic core of hTR. An overlay of a large set of models generated has provided a low-resolution structure (6.5-8.0 A) that can readily be refined as new structural information becomes available. A notable feature of the modeled structure is the positioning of the template adjacent to the pseudoknot, which brings a number of conserved nucleotides close in space.

Folding of the adenine riboswitch

The pbuE adenine riboswitch undergoes metal ion-dependent folding that involves a loop-loop interaction. Binding of 2-aminopurine to the aptamer domain strongly correlates with the ability of the loops to interact, and single-molecule FRET studies reveal that folding proceeds via a discrete intermediate. Folding occurs in the absence of adenine ligand, but ligand binding stabilizes the folded structure by increasing the folding rate and decreasing the unfolding rate, and it lowers the magnesium ion concentration required to promote the loop-loop interaction. Individual aptamer molecules exhibit great heterogeneity in folding and unfolding rates, but this is reduced in the presence of adenine. In the full riboswitch, the adenine binding domain fails to fold because of conformational competition by the terminator stem. Thus, riboswitch function should depend on the relative rates of ligand binding and the transcriptional process.