NMR structure of the active conformation of the Varkud satellite ribozyme cleavage site

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.

The Varkud Satellite Ribozyme: A Thirty-Year Journey through Biochemistry, Crystallography, and Computation

The discovery of catalytic RNAs or ribozymes introduced a new class of enzymes to biology. In addition to their increasingly important roles in modern life, ribozymes are key players in the RNA World hypothesis, which posits that life started or flourished with RNA supporting both genetic and enzymatic functions. Therefore, investigations into the mechanisms of ribozyme function provide an exciting opportunity to examine the foundational principles of biological catalysis. Ribozymes are also attractive model systems to investigate the relationship between structure and function in RNA. Endonucleolytic ribozymes represent the largest class of catalytic RNA, of which the Varkud satellite (VS) ribozyme is structurally the most complex. The last ribozyme to be discovered by accident, the VS ribozyme had eluded structural determination for over two decades. When we solved the first crystal structures of the VS ribozyme, an extensive body of biochemical and biophysical data had accumulated over the years with which we could evaluate the functional relevance of the structure. Conversely, the structures provided a new perspective from which to reexamine the functional data and test new hypotheses. The VS ribozyme is organized in a modular fashion where independently folding domains assemble into the active conformation of the ribozyme via three-way junctions. Structures of the VS ribozyme in complex with its substrate at different stages of activation enabled us to map the structural reorganization of the substrate that must precede catalysis. In addition to defining the global architecture of the RNA, the essential interactions between the substrate and catalytic domains, and the rearrangements in the substrate prior to catalysis, these structures provided detailed snapshots of the ribozyme active site, revealing potential catalytic interactions. High resolution structures of the active site bolstered the view that the catalytic mechanism involved nucleobase-mediated general acid-base catalysis and uncovered additional catalytic interactions between the cleavage site and catalytic residues. Informed by the crystal structures of the VS ribozyme, an integrated experimental and computational approach identified the key players and essential interactions that define the active site of the ribozyme. This confluence of biochemical, structural, and computational studies revealed the catalytic mechanism of the ribozyme at unprecedented detail. Additionally, comparative analyses of the active site structures of the VS ribozyme and other nucleic acid-based endoribonucleases revealed common architectural motifs and strikingly similar catalytic strategies. In this Account, we document the progress of VS ribozyme research starting from its discovery and extending to the elucidation of its detailed catalytic mechanism 30 years later.

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.

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.

Unraveling Mg2+-RNA binding with atomistic molecular dynamics

Interaction with divalent cations is of paramount importance for RNA structural stability and function. We report here a detailed molecular dynamics study of all the possible binding sites for Mg2+ on an RNA duplex, including both direct (inner sphere) and indirect (outer sphere) binding. In order to tackle sampling issues, we develop a modified version of bias-exchange metadynamics, which allows us to simultaneously compute affinities with previously unreported statistical accuracy. Results correctly reproduce trends observed in crystallographic databases. Based on this, we simulate a carefully chosen set of models that allows us to quantify the effects of competition with monovalent cations, RNA flexibility, and RNA hybridization. Our simulations reproduce the decrease and increase of Mg2+ affinity due to ion competition and hybridization, respectively, and predict that RNA flexibility has a site-dependent effect. This suggests a nontrivial interplay between RNA conformational entropy and divalent cation binding.

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.

Chemistry and Biology of Self-Cleaving Ribozymes

Self-cleaving ribozymes were discovered 30 years ago, but their biological distribution and catalytic mechanisms are only beginning to be defined. Each ribozyme family is defined by a distinct structure, with unique active sites accelerating the same transesterification reaction across the families. Biochemical studies show that general acid-base catalysis is the most common mechanism of self-cleavage, but metal ions and metabolites can be used as cofactors. Ribozymes have been discovered in highly diverse genomic contexts throughout nature, from viroids to vertebrates. Their biological roles include self-scission during rolling-circle replication of RNA genomes, co-transcriptional processing of retrotransposons, and metabolite-dependent gene expression regulation in bacteria. Other examples, including highly conserved mammalian ribozymes, suggest that many new biological roles are yet to be discovered.

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.

Linear free energy relationships in RNA transesterification: theoretical models to aid experimental interpretations

RNA cleavage transesterification is of fundamental reaction in biology that is catalyzed by both protein and RNA enzymes. In this work, a series of RNA transesterification model reactions with a wide range of leaving groups are investigated with density-functional calculations in an aqueous solvation environment in order to study linear free energy relationships (LFERs) and their connection to transition state structure and bonding. Overall, results obtained from the polarizable continuum solvation model with UAKS radii produce the best linear correlations and closest overall agreement with experimental results. Reactions with a poor leaving group are predicted to proceed via a stepwise mechanism with a late transition state that is rate controlling. As leaving group becomes more acidic and labile, the barriers of both early and late transition states decrease. LFERs for each transition state are computed, with the late transition state barrier showing greater sensitivity to leaving group pKa. For sufficiently enhanced leaving groups, the reaction mechanism transits to a concerted mechanism characterized by a single early transition state. Further linear relationships were derived for bond lengths and bond orders as a function of leaving group pKa and rate constant values that can be used for prediction. This work provides important benchmark linear free energy data that allows a molecular-level characterization of the structure and bonding of the transition states for this important class of phosphoryl transfer reactions. The relations reported herein can be used to aid in the interpretation of data obtained from experimental studies of non-catalytic and catalytic mechanisms.

RNAMotifScanX: a graph alignment approach for RNA structural motif identification

RNA structural motifs are recurrent three-dimensional (3D) components found in the RNA architecture. These RNA structural motifs play important structural or functional roles and usually exhibit highly conserved 3D geometries and base-interaction patterns. Analysis of the RNA 3D structures and elucidation of their molecular functions heavily rely on efficient and accurate identification of these motifs. However, efficient RNA structural motif search tools are lacking due to the high complexity of these motifs. In this work, we present RNAMotifScanX, a motif search tool based on a base-interaction graph alignment algorithm. This novel algorithm enables automatic identification of both partially and fully matched motif instances. RNAMotifScanX considers noncanonical base-pairing interactions, base-stacking interactions, and sequence conservation of the motifs, which leads to significantly improved sensitivity and specificity as compared with other state-of-the-art search tools. RNAMotifScanX also adopts a carefully designed branch-and-bound technique, which enables ultra-fast search of large kink-turn motifs against a 23S rRNA. The software package RNAMotifScanX is implemented using GNU C++, and is freely available from http://genome.ucf.edu/RNAMotifScanX.

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.

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

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

Structural insights into substrate recognition by the Neurospora Varkud satellite ribozyme: importance of U-turns at the kissing-loop junction

Substrate recognition by the Neurospora Varkud satellite ribozyme depends on the formation of a magnesium-dependent kissing-loop interaction between the stem-loop I (SLI) substrate and stem-loop V (SLV) of the catalytic domain. From mutagenesis studies, it has been established that this I/V kissing-loop interaction involves three Watson–Crick base pairs and is associated with a structural rearrangement of the SLI substrate that facilitates catalysis. Here, we report the NMR structural characterization of this I/V kissing-loop using isolated stem-loops. NMR studies were performed on different SLI/SLV complexes containing a common SLV and shiftable, preshifted, or double-stranded SLI variants. These studies confirm the presence of three Watson–Crick base pairs at the kissing-loop junction and provide evidence for the structural rearrangement of shiftable SLI variants upon SLV binding. NMR structure determination of an SLI/SLV complex demonstrates that both the SLI and SLV loops adopt U-turn structures, which facilitates intermolecular Watson–Crick base pairing. Several other interactions at the I/V interface, including base triples and base stacking, help create a continuously stacked structure. These NMR studies provide a structural basis to understand the stability of the I/V kissing-loop interaction and lead us to propose a kinetic model for substrate activation in the VS ribozyme.

Affinity purification of T7 RNA transcripts with homogeneous ends using ARiBo and CRISPR tags

Affinity purification of RNA using the ARiBo tag technology currently provides an ideal approach to quickly prepare RNA with 3' homogeneity. Here, we explored strategies to also ensure 5' homogeneity of affinity-purified RNAs. First, we systematically investigated the effect of starting nucleotides on the 5' heterogeneity of a small SLI RNA substrate from the Neurospora VS ribozyme purified from an SLI-ARiBo precursor. A series of 32 SLI RNA sequences with variations in the +1 to +3 region was produced from two T7 promoters (class III consensus and class II 2.5) using either the wild-type T7 RNA polymerase or the P266L mutant. Although the P266L mutant helps decrease the levels of 5'-sequence heterogeneity in several cases, significant levels of 5' heterogeneity (≥1.5%) remain for transcripts starting with GGG, GAG, GCG, GGC, AGG, AGA, AAA, ACA, AUA, AAC, ACC, AUC, and AAU. To provide a more general approach to purifying RNA with 5' homogeneity, we tested the suitability of using a small CRISPR RNA stem-loop at the 5' end of the SLI-ARiBo RNA. Interestingly, we found that complete cleavage of the 5'-CRISPR tag with the Cse3 endoribonuclease can be achieved quickly from CRISPR-SLI-ARiBo transcripts. With this procedure, it is possible to generate SLI-ARiBo RNAs starting with any of the four standard nucleotides (G, C, A, or U) involved in either a single- or a double-stranded structure. Moreover, the 5'-CRISPR-based strategy can be combined with affinity purification using the 3'-ARiBo tag for quick purification of RNA with both 5' and 3' homogeneity.

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.

Additional roles of a peripheral loop-loop interaction in the Neurospora VS ribozyme

Many RNAs contain tertiary interactions that contribute to folding the RNA into its functional 3D structure. In the VS ribozyme, a tertiary loop-loop kissing interaction involving stem-loops I and V is also required to rearrange the secondary structure of stem-loop I such that nucleotides at the base of stem I, which contains the cleavage-ligation site, can adopt the conformation required for activity. In the current work, we have used mutants that constitutively adopt the catalytically permissive conformation to search for additional roles of the kissing interaction in vitro. Using mutations that disrupt or restore the kissing interaction, we find that the kissing interaction contributes ~1000-fold enhancement to the rates of cleavage and ligation. Large Mg(2+)-dependent effects on equilibrium were also observed: in the presence of the kissing interaction cleavage is favored >10-fold at micromolar concentrations of Mg(2+); whereas ligation is favored >10-fold at millimolar concentrations of Mg(2+). In the absence of the kissing interaction cleavage exceeds ligation at all concentrations of Mg(2+). These data provide evidence that the kissing interaction strongly affects the observed cleavage and ligation rate constants and the cleavage-ligation equilibrium of the ribozyme.

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.

Structural characterization of naturally occurring RNA single mismatches

RNA is known to be involved in several cellular processes; however, it is only active when it is folded into its correct 3D conformation. The folding, bending and twisting of an RNA molecule is dependent upon the multitude of canonical and non-canonical secondary structure motifs. These motifs contribute to the structural complexity of RNA but also serve important integral biological functions, such as serving as recognition and binding sites for other biomolecules or small ligands. One of the most prevalent types of RNA secondary structure motifs are single mismatches, which occur when two canonical pairs are separated by a single non-canonical pair. To determine sequence–structure relationships and to identify structural patterns, we have systematically located, annotated and compared all available occurrences of the 30 most frequently occurring single mismatch-nearest neighbor sequence combinations found in experimentally determined 3D structures of RNA-containing molecules deposited into the Protein Data Bank. Hydrogen bonding, stacking and interaction of nucleotide edges for the mismatched and nearest neighbor base pairs are described and compared, allowing for the identification of several structural patterns. Such a database and comparison will allow researchers to gain insight into the structural features of unstudied sequences and to quickly look-up studied sequences.

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.

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.

Catalytic strategies of self-cleaving ribozymes

[Structure: see text]. Five naturally occurring nucleolytic ribozymes have been identified: the hammerhead, hairpin, glmS, hepatitis delta virus (HDV), and Varkud satellite (VS) ribozymes. All of these RNA enzymes catalyze self-scission of the RNA backbone using a chemical mechanism equivalent to that of RNase A. RNase A uses four basic strategies to promote this reaction: geometric constraints, activation of the nucleophile, transition-state stabilization, and leaving group protonation. In this Account, we discuss the current thinking on how nucleolytic ribozymes harness RNase A's four sources of catalytic power. The geometry of the phosphodiester cleavage reaction constrains the nucleotides flanking the scissile phosphate so that they are unstacked from a canonical A-form helix and thus require alternative stabilization. Crystal structures and mutational analysis reveal that cross-strand base pairing, along with unconventional stacking and tertiary hydrogen-bonding interactions, work to stabilize the splayed conformation in nucleolytic ribozymes. Deprotonation of the 2'-OH nucleophile greatly increases its nucleophilicity in the strand scission reaction. Crystal structures of the hammerhead, hairpin, and glmS ribozymes reveal the N1 of a G residue within hydrogen-bonding distance of the 2'-OH. In each case, this residue has also been shown to be important for catalysis. In the HDV ribozyme, a hydrated magnesium has been implicated as the general base. Catalysis by the VS ribozyme requires both an A and a G, but the precise role of either has not been elucidated. Enzymes can lower the energy of a chemical reaction by binding more tightly to the transition state than to the ground states. Comparison of the hairpin ground- and transition-state mimic structures reveal greater hydrogen bonding to the transition-state mimic structure, suggesting transition-state stabilization as a possible catalytic strategy. However, the hydrogen-bonding pattern in the glmS ribozyme transition-state mimic structure and the ground-state structures are equivalent. Protonation of the 5'-O leaving group by a variety of functional groups can promote the cleavage reaction. In the HDV ribozyme, the general acid is a conserved C residue. In the hairpin ribozyme, a G residue has been implicated in protonation of the leaving group. An A in the hammerhead ribozyme probably plays a similar role. In the glmS ribozyme, an exogenous cofactor may provide the general acid. This diversity is in contrast to the relatively small number of functional groups that serve as a general base, where at least three of the nucleolytic ribozymes may use the N1 of a G.

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.

An important role of G638 in the cis-cleavage reaction of the Neurospora VS ribozyme revealed by a novel nucleotide analog incorporation method

We describe a chemical coupling procedure that allows joining of two RNAs, one of which contains a site-specific base analog substitution, in the absence of divalent ions. This method allows incorporation of nucleotide analogs at specific positions even into large, cis-cleaving ribozymes. Using this method we have studied the effects of substitution of G638 in the cleavage site loop of the VS ribozyme with a variety of purine analogs having different functional groups and pK(a) values. Cleavage rate versus pH profiles combined with kinetic solvent isotope experiments indicate an important role for G638 in proton transfer during the rate-limiting step of the cis-cleavage reaction.

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.

DNA mimicry by a high-affinity anti-NF-kappaB RNA aptamer

The binding of RNA molecules to proteins or other ligands can require extensive RNA folding to create an induced fit. Understanding the generality of this principle involves comparing structures of RNA before and after complex formation. Here we report the NMR solution structure of a 29-nt RNA aptamer whose crystal structure had previously been determined in complex with its transcription factor target, the p50₂ form of NF-κB. The RNA aptamer internal loop structure has pre-organized features that are also found in the complex, including non-canonical base pairing and cross-strand base stacking. Remarkably, the free RNA aptamer structure possesses a major groove that more closely resembles B-form DNA than RNA. Upon protein binding, changes in RNA structure include the kinking of the internal loop and distortion of the terminal tetraloop. Thus, complex formation involves both pre-formed and induced fit binding interactions. The high affinity of the NF-κB transcription factor for this RNA aptamer may largely be due to the structural pre-organization of the RNA that results in its ability to mimic DNA.

A comparative analysis of the triloops in all high-resolution RNA structures reveals sequence structure relationships

Despite an increasing number of experimentally determined RNA structures, the gap between the number of structures and that of RNA families is still growing. To overcome this limitation, efficient and reliable RNA modeling methodologies must be developed. In order to reach this goal, here, we show how triloop sequence-structure relationships have been inferred through a systematic analysis of all triloops found in available high-resolution structures. The structural annotation of all triloops allowed us to define discrete states of the triloop's conformational space, and therefore an explicit sequence-to-structure relation. The sequence-structure relationships inferred from this explicit relation are presented in a convenient modeling table that provides a limited set of possible three-dimensional structures given any triloop sequence. The table is indexed by the two nucleotides that form the triloop's flanking base pair, since they are shown to provide the most information about the triloop three-dimensional structures. We also report the observations in the X-ray crystallographic structures of important conformational variations, which we believe might be the result of RNA dynamic.

FR3D: finding local and composite recurrent structural motifs in RNA 3D structures

New methods are described for finding recurrent three-dimensional (3D) motifs in RNA atomic-resolution structures. Recurrent RNA 3D motifs are sets of RNA nucleotides with similar spatial arrangements. They can be local or composite. Local motifs comprise nucleotides that occur in the same hairpin or internal loop. Composite motifs comprise nucleotides belonging to three or more different RNA strand segments or molecules. We use a base-centered approach to construct efficient, yet exhaustive search procedures using geometric, symbolic, or mixed representations of RNA structure that we implement in a suite of MATLAB programs, "Find RNA 3D" (FR3D). The first modules of FR3D preprocess structure files to classify base-pair and -stacking interactions. Each base is represented geometrically by the position of its glycosidic nitrogen in 3D space and by the rotation matrix that describes its orientation with respect to a common frame. Base-pairing and base-stacking interactions are calculated from the base geometries and are represented symbolically according to the Leontis/Westhof basepairing classification, extended to include base-stacking. These data are stored and used to organize motif searches. For geometric searches, the user supplies the 3D structure of a query motif which FR3D uses to find and score geometrically similar candidate motifs, without regard to the sequential position of their nucleotides in the RNA chain or the identity of their bases. To score and rank candidate motifs, FR3D calculates a geometric discrepancy by rigidly rotating candidates to align optimally with the query motif and then comparing the relative orientations of the corresponding bases in the query and candidate motifs. Given the growing size of the RNA structure database, it is impossible to explicitly compute the discrepancy for all conceivable candidate motifs, even for motifs with less than ten nucleotides. The screening algorithm that we describe finds all candidate motifs whose geometric discrepancy with respect to the query motif falls below a user-specified cutoff discrepancy. This technique can be applied to RMSD searches. Candidate motifs identified geometrically may be further screened symbolically to identify those that contain particular basepair types or base-stacking arrangements or that conform to sequence continuity or nucleotide identity constraints. Purely symbolic searches for motifs containing user-defined sequence, continuity and interaction constraints have also been implemented. We demonstrate that FR3D finds all occurrences, both local and composite and with nucleotide substitutions, of sarcin/ricin and kink-turn motifs in the 23S and 5S ribosomal RNA 3D structures of the H. marismortui 50S ribosomal subunit and assigns the lowest discrepancy scores to bona fide examples of these motifs. The search algorithms have been optimized for speed to allow users to search the non-redundant RNA 3D structure database on a personal computer in a matter of minutes.

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.

The NMR structure of an internal loop from 23S ribosomal RNA differs from its structure in crystals of 50s ribosomal subunits

Internal loops play an important role in structure and folding of RNA and in recognition of RNA by other molecules such as proteins and ligands. An understanding of internal loops with propensities to form a particular structure will help predict RNA structure, recognition, and function. The structures of internal loops 5' 1009CUAAG1013 3'/3' 1168GAAGC1164 5' and 5' 998CUAAG1002 3'/3' 1157GAAGC1153 5' from helix 40 of the large subunit rRNA in Deinococcus radiodurans and Escherichia coli, respectively, are phylogenetically conserved, suggesting functional relevance. The energetics and NMR solution structure of the loop were determined in the duplex 5' 1GGCUAAGAC9 3'/3' 18CCGAAGCUG10 5'. The internal loop forms a different structure in solution and in the crystal structures of the ribosomal subunits. In particular, the crystal structures have a bulged out adenine at the equivalent of position A15 and a reverse Hoogsteen UA pair (trans Watson-Crick/Hoogsteen UA) at the equivalent of U4 and A14, whereas the solution structure has a single hydrogen bond UA pair (cis Watson-Crick/sugar edge A15U4) between U4 and A15 and a sheared AA pair (trans Hoogsteen/sugar edge A14A5) between A5 and A14. There is cross-strand stacking between A6 and A14 (A6/A14/A15 stacking pattern) in the NMR structure. All three structures have a sheared GA pair (trans Hoogsteen/sugar edge A6G13) at the equivalent of A6 and G13. The internal loop has contacts with ribosomal protein L20 and other parts of the RNA in the crystal structures. These contacts presumably provide the free energy to rearrange the base pairing in the loop. Evidently, molecular recognition of this internal loop involves induced fit binding, which could confer several advantages. The predicted thermodynamic stability of the loop agrees with the experimental value, even though the thermodynamic model assumes a Watson-Crick UA pair.

Consecutive GA pairs stabilize medium-size RNA internal loops

Internal loops in RNA are important for folding and function. Consecutive noncanonical pairs can form in internal loops having at least two nucleotides on each side. Thermodynamic and structural insights into such internal loops should improve approximations for their stabilities and predictions of secondary and three-dimensional structures. Most natural internal loops are purine rich. A series of oligoribonucleotides that form purine-rich internal loops of 5-10 nucleotides, including kink-turn loops, were studied by UV melting, exchangeable proton and phosphorus NMR. Three consecutive GA pairs with the motif 5' Y GGA/3' R AAG or GGA R 3'/AAG Y 5' (i.e., 5' GGA 3'/3' AAG 5' closed on at least one side with a CG, UA, or UG pair with Y representing C or U and R representing A or G) stabilize internal loops having 6-10 nucleotides. Certain motifs with two consecutive GA pairs are also stabilizing. In internal loops with three or more nucleotides on each side, the motif 5' U G/3' G A has stability similar to 5' C G/3' G A. A revised model for predicting stabilities of internal loops with 6-10 nucleotides is derived by multiple linear regression. Loops with 2 x 3 nucleotides are predicted well by a previous thermodynamic model.

Solution structure of an informationally complex high-affinity RNA aptamer to GTP

Higher-affinity RNA aptamers to GTP are more informationally complex than lower-affinity aptamers. Analog binding studies have shown that the additional information needed to improve affinity does not specify more interactions with the ligand. In light of those observations, we would like to understand the structural characteristics that enable complex aptamers to bind their ligands with higher affinity. Here we present the solution structure of the 41-nt Class I GTP aptamer (K(d) = 75 nM) as determined by NMR. The backbone of the aptamer forms a reverse-S that shapes the binding pocket. The ligand nucleobase stacks between purine platforms and makes hydrogen bonds with the edge of another base. Interestingly, the local modes of interaction for the Class I aptamer and an RNA aptamer that binds ATP with a K(d) of 6 microM are very much alike. The aptamers exhibit nearly identical levels of binding specificity and fraction of ligand sequestered from the solvent (81%-85%). However, the GTP aptamer is more informationally complex (approximately 45 vs. 35 bits) and has a larger recognition bulge (15 vs. 12 nucleotides) with many more stabilizing base-base interactions. Because the aptamers have similar modes of ligand binding, we conclude that the stabilizing structural elements in the Class I aptamer are responsible for much of the difference in K(d). These results are consistent with the hypothesis that increasing the number of intra-RNA interactions, rather than adding specific contacts to the ligand, is the simplest way to improve binding affinity.

The RNA Ontology Consortium: an open invitation to the RNA community

The aim of the RNA Ontology Consortium (ROC) is to create an integrated conceptual framework-an RNA Ontology (RO)-with a common, dynamic, controlled, and structured vocabulary to describe and characterize RNA sequences, secondary structures, three-dimensional structures, and dynamics pertaining to RNA function. The RO should produce tools for clear communication about RNA structure and function for multiple uses, including the integration of RNA electronic resources into the Semantic Web. These tools should allow the accurate description in computer-interpretable form of the coupling between RNA architecture, function, and evolution. The purposes for creating the RO are, therefore, (1) to integrate sequence and structural databases; (2) to allow different computational tools to interoperate; (3) to create powerful software tools that bring advanced computational methods to the bench scientist; and (4) to facilitate precise searches for all relevant information pertaining to RNA. For example, one initial objective of the ROC is to define, identify, and classify RNA structural motifs described in the literature or appearing in databases and to agree on a computer-interpretable definition for each of these motifs. To achieve these aims, the ROC will foster communication and promote collaboration among RNA scientists by coordinating frequent face-to-face workshops to discuss, debate, and resolve difficult conceptual issues. These meeting opportunities will create new directions at various levels of RNA research. The ROC will work closely with the PDB/NDB structural databases and the Gene, Sequence, and Open Biomedical Ontology Consortia to integrate the RO with existing biological ontologies to extend existing content while maintaining interoperability.

Statistical analysis of RNA backbone

Local conformation is an important determinant of RNA catalysis and binding. The analysis of RNA conformation is particularly difficult due to the large number of degrees of freedom (torsion angles) per residue. Proteins, by comparison, have many fewer degrees of freedom per residue. In this work, we use and extend classical tools from statistics and signal processing to search for clusters in RNA conformational space. Results are reported both for scalar analysis, where each torsion angle is separately studied, and for vectorial analysis, where several angles are simultaneously clustered. Adapting techniques from vector quantization and clustering to the RNA structure, we find torsion angle clusters and RNA conformational motifs. We validate the technique using well-known conformational motifs, showing that the simultaneous study of the total torsion angle space leads to results consistent with known motifs reported in the literature and also to the finding of new ones.

Application of NMR and EPR methods to the study of RNA

The application of techniques based on magnetic resonance, specifically electron paramagnetic resonance (EPR) and nuclear magnetic resonance (NMR), has provided a wealth of new information on RNA structures, as well as insights into the dynamics and function of these important biomolecules. NMR spectroscopy is very successful for determining the solution structures of small RNA domains, aptamers and ribozymes, and exploring their intramolecular dynamics and interactions with ligands. EPR-based methods have been used to map local dynamic and structural features of RNA, to explore different modes of RNA-ligand interaction, to obtain long-range structural restraints and to probe metal-ion-binding sites.

The role of phosphate groups in the VS ribozyme-substrate interaction

The VS ribozyme trans-cleavage substrate interacts with the catalytic RNA via tertiary interactions. To study the role of phosphate groups in the ribozyme-substrate interaction, 18 modified substrates were synthesized, where an epimeric phosphorothioate replaces one of the phosphate diester linkages. Sites in the stem-loop substrate where phosphorothioate substitution impaired reaction cluster in two regions. The first site is the scissile phosphate diester linkage and nucleotides downstream of this and the second site is within the loop region. The addition of manganese ions caused recovery of the rate of reaction for phosphorothioate substitutions between A621 and A622 and U631 and C632, suggesting that these two phosphate groups may serve as ligands for two metal ions. In contrast, significant manganese rescue was not observed for the scissile phosphate diester linkage implying that electrophilic catalysis by metal ions is unlikely to contribute to VS ribozyme catalysis. In addition, an increase in the reaction rate of the unmodified VS ribozyme was observed when a mixture of magnesium and manganese ions acted as the cofactor. One possible explanation for this effect is that the cleavage reaction of the VS ribozyme is rate limited by a metal dependent docking of the substrate on the ribozyme.

The solution structure of the VS ribozyme active site loop reveals a dynamic "hot-spot"

The VS ribozyme is the largest ribozyme in its class and is also the least structurally characterized thus far. The current working model of the VS ribozyme locates the active site in stem-loop VI. The solution structure of this active site loop was determined using high resolution NMR spectroscopy. The structure reveals that the ground-state conformation of the active site differs significantly from that determined previously from chemical structure probing and mutational analysis of the ribozyme in its active conformation, which contains several looped out bases. In contrast, the base-pairing scheme found for the isolated loop contains three mismatched base-pairs: an A+-C, a G-U wobble, and a sheared G-A base-pair and no looped out bases. Dynamics observed within the active site loop provide insight into the mechanism by which the RNA can rearrange its secondary structure into an "activated" conformation prior to cleavage. These findings lend support to the idea that RNA secondary structure is more fluid than once believed and that a better understanding of structure and dynamic features of ribozymes is required to unravel the intricacies of their catalytic abilities.

Using an RNA secondary structure partition function to determine confidence in base pairs predicted by free energy minimization

A partition function calculation for RNA secondary structure is presented that uses a current set of nearest neighbor parameters for conformational free energy at 37 degrees C, including coaxial stacking. For a diverse database of RNA sequences, base pairs in the predicted minimum free energy structure that are predicted by the partition function to have high base pairing probability have a significantly higher positive predictive value for known base pairs. For example, the average positive predictive value, 65.8%, is increased to 91.0% when only base pairs with probability of 0.99 or above are considered. The quality of base pair predictions can also be increased by the addition of experimentally determined constraints, including enzymatic cleavage, flavin mono-nucleotide cleavage, and chemical modification. Predicted secondary structures can be color annotated to demonstrate pairs with high probability that are therefore well determined as compared to base pairs with lower probability of pairing.

The Varkud satellite ribozyme

The VS ribozyme is the largest nucleolytic ribozyme, for which there is no crystal structure to date. The ribozyme consists of five helical sections, organized by two three-way junctions. The global structure has been determined by solution methods, particularly FRET. The substrate stem-loop binds into a cleft formed between two helices, while making a loop-loop contact with another section of the ribozyme. The scissile phosphate makes a close contact with an internal loop (the A730 loop), the probable active site of the ribozyme. This loop contains a particularly critical nucleotide A756. Most changes to this nucleotide lead to three-orders of magnitude slower cleavage, and the Watson-Crick edge is especially important. NAIM experiments indicate that a protonated base is required at this position for the ligation reaction. A756 is thus a strong candidate for nucleobase participation in the catalytic chemistry.