Sequence and structural conservation in RNA ribose zippers

A comprehensive survey of long-range tertiary interactions and motifs in non-coding RNA structures

Understanding the 3D structure of RNA is key to understanding RNA function. RNA 3D structure is modular and can be seen as a composition of building blocks of various sizes called tertiary motifs. Currently, long-range motifs formed between distant loops and helical regions are largely less studied than the local motifs determined by the RNA secondary structure. We surveyed long-range tertiary interactions and motifs in a non-redundant set of non-coding RNA 3D structures. A new dataset of annotated LOng-RAnge RNA 3D modules (LORA) was built using an approach that does not rely on the automatic annotations of non-canonical interactions. An original algorithm, ARTEM, was developed for annotation-, sequence- and topology-independent superposition of two arbitrary RNA 3D modules. The proposed methods allowed us to identify and describe the most common long-range RNA tertiary motifs. Along with the prevalent canonical A-minor interactions, a large number of previously undescribed staple interactions were observed. The most frequent long-range motifs were found to belong to three main motif families: planar staples, tilted staples, and helical packing motifs.

Ancestry of RNA/RNA interaction regions within segmented ribosomes

The ribosome is the universally conserved ribozyme that translates DNA coded instructions into proteins with the assistance of other RNA molecules, including transfer and messenger RNAs. Of particular interest is the segmentation phenomena, which is found in trypanosomatids and other protists. In these organisms, the large subunit ribosomal RNA is assembled from multiple smaller RNAs. This phenomenon posits several challenges to the folding and stabilization of such ribosomes to retain functionality and efficiency. In earlier studies, RNA/protein interactions were suggested to fully compensate for the fragmentation. Recently, several conserved RNA/RNA interaction regions were described in the cryo-EM structures of segmented ribosomes from trypanosomatids. These regions also seemed to aid in the folding and stabilization of such ribosomes, even before the ribosomal proteins start their association. In the present study, the existence of conserved RNA/RNA interaction regions shared between trypanosomatid and Euglena gracilis segmented ribosomes was confirmed, despite differences in segmentation patterns. Analysis of the crystallographic structures of unsegmented ribosomes from other Eukaryotes, Bacteria, and Archaea allowed us to estimate the relative age of highly conserved RNA/RNA interaction regions. These results strongly suggest that common interaction regions likely date far back into the ribosomes of the last common ancestor. Results also revealed that single hydrogen bonds are overwhelmingly facilitated by the 2'OH, a distinctive RNA feature. This supports the notion that RNA predates DNA and places some constraints on alternative nucleic acids proposals.

Cryo neutron crystallography demonstrates influence of RNA 2'-OH orientation on conformation, sugar pucker and water structure

The ribose 2′-hydroxyl is the key chemical difference between RNA and DNA and primary source of their divergent structural and functional characteristics. Macromolecular X-ray diffraction experiments typically do not reveal the positions of hydrogen atoms. Thus, standard crystallography cannot determine 2′-OH orientation (H2′-C2′-O2′-HO2′ torsion angle) and its potential roles in sculpting the RNA backbone and the expansive fold space. Here, we report the first neutron crystal structure of an RNA, the Escherichia coli rRNA Sarcin-Ricin Loop (SRL). 2′-OD orientations were established for all 27 residues and revealed O-D bonds pointing toward backbone (O3′, 13 observations), nucleobase (11) or sugar (3). Most riboses in the SRL stem region show a 2′-OD backbone-orientation. GAGA-tetraloop riboses display a 2′-OD base-orientation. An atypical C2′-endo sugar pucker is strictly correlated with a 2′-OD sugar-orientation. Neutrons reveal the strong preference of the 2′-OH to donate in H-bonds and that 2′-OH orientation affects both backbone geometry and ribose pucker. We discuss 2′-OH and water molecule orientations in the SRL neutron structure and compare with results from a solution phase 10 μs MD simulation. We demonstrate that joint cryo-neutron/X-ray crystallography offers an all-in-one approach to determine the complete structural properties of RNA, i.e. geometry, conformation, protonation state and hydration structure.

Automatic Learning of Hydrogen-Bond Fixes in the AMBER RNA Force Field

The capability of current force fields to reproduce RNA structural dynamics is limited. Several methods have been developed to take advantage of experimental data in order to enforce agreement with experiments. Here, we extend an existing framework which allows arbitrarily chosen force-field correction terms to be fitted by quantification of the discrepancy between observables back-calculated from simulation and corresponding experiments. We apply a robust regularization protocol to avoid overfitting and additionally introduce and compare a number of different regularization strategies, namely, L1, L2, Kish size, relative Kish size, and relative entropy penalties. The training set includes a GACC tetramer as well as more challenging systems, namely, gcGAGAgc and gcUUCGgc RNA tetraloops. Specific intramolecular hydrogen bonds in the AMBER RNA force field are corrected with automatically determined parameters that we call gHBfix_opt. A validation involving a separate simulation of a system present in the training set (gcUUCGgc) and new systems not seen during training (CAAU and UUUU tetramers) displays improvements regarding the native population of the tetraloop as well as good agreement with NMR experiments for tetramers when using the new parameters. Then, we simulate folded RNAs (a kink–turn and L1 stalk rRNA) including hydrogen bond types not sufficiently present in the training set. This allows a final modification of the parameter set which is named gHBfix21 and is suggested to be applicable to a wider range of RNA systems.

Nonstandard RNA/RNA interactions likely enhance folding and stability of segmented ribosomes

The ribosome is the molecular factory that catalyzes all coded protein synthesis in extant organisms. Eukaryotic ribosomes are typically assembled out of four rRNAs; namely, 5S, 5.8S, 18S, and 28S. However, the 28S rRNA of some trypanosomatid organisms has been found to be segmented into six independent rRNAs of different sizes. The two largest segments have multiple sites where they jointly form stems comprised of standard base pairs that can hold them together. However, such regions of interaction are not observed among the four smaller RNAs. Early reports suggested that trypanosomatid segmented ribosome assembly was essentially achieved thanks to their association with rProteins. However, examination of cryo-EM ribosomal structures from Trypanosoma brucei, Leishmania donovani, and Trypanosoma cruzi reveals several long-range nonstandard RNA/RNA interactions. Most of these interactions are clusters of individual hydrogen bonds and so are not readily predictable. However, taken as a whole, they represent significant stabilizing energy that likely facilitates rRNA assembly and the overall stability of the segmented ribosomes. In the context of origin of life studies, the current results provide a better understanding of the true nature of RNA sequence space and what might be possible without an RNA replicase.

Cooperativity and Interdependency between RNA Structure and RNA-RNA Interactions

Complex RNA–RNA interactions are increasingly known to play key roles in numerous biological processes from gene expression control to ribonucleoprotein granule formation. By contrast, the nature of these interactions and characteristics of their interfaces, especially those that involve partially or wholly structured RNAs, remain elusive. Herein, we discuss different modalities of RNA–RNA interactions with an emphasis on those that depend on secondary, tertiary, or quaternary structure. We dissect recently structurally elucidated RNA–RNA complexes including RNA triplexes, riboswitches, ribozymes, and reverse transcription complexes. These analyses highlight a reciprocal relationship that intimately links RNA structure formation with RNA–RNA interactions. The interactions not only shape and sculpt RNA structures but also are enabled and modulated by the structures they create. Understanding this two-way relationship between RNA structure and interactions provides mechanistic insights into the expanding repertoire of noncoding RNA functions, and may inform the design of novel therapeutics that target RNA structures or interactions.

Broken symmetry between RNA enantiomers in a crystal lattice

Explaining the origin of the homochirality of biological molecules requires a mechanism of disrupting the natural equilibrium between enantiomers and amplifying the initial imbalance to significant levels. Authors of existing models have sought an explanation in the parity-breaking weak nuclear force, in some selectively acting external factor, or in random fluctuations that subsequently became amplified by an autocatalytic process. We have obtained crystals in which l- and d-enantiomers of short RNA duplexes assemble in an asymmetric manner. These enantiomers make different lattice contacts and have different exposures to water and metal ions present in the crystal. Apparently, asymmetry between enantiomers can arise upon their mutual interactions and then propagate via crystallization. Asymmetric racemic compounds are worth considering as possible factors in symmetry breaking and enantioenrichment that took place in the early biosphere.

Features and Functions of the A-Minor Motif, the Most Common Motif in RNA Structure

A-minor motifs are RNA tertiary structure motifs that generally involve a canonical base pair and an adenine base forming hydrogen bonds with the minor groove of the base pair. Such motifs are among the most common tertiary interactions in known RNA structures, comparable in number with the non-canonical base pairs. They are often found in functionally important regions of non-coding RNAs and, in particular, play a central role in protein synthesis. Here, we review local variations of the A-minor geometry and discuss difficulties associated with their annotation, as well as various structural contexts and common A-minor co-motifs, and diverse functions of A-minors in various processes in a living cell.

Engineering Crystal Packing in RNA Structures I: Past and Future Strategies for Engineering RNA Packing in Crystals

X-ray crystallography remains a powerful method to gain atomistic insights into the catalytic and regulatory functions of RNA molecules. However, the technique requires the preparation of diffraction-quality crystals. This is often a resource- and time-consuming venture because RNA crystallization is hindered by the conformational heterogeneity of RNA, as well as the limited opportunities for stereospecific intermolecular interactions between RNA molecules. The limited success at crystallization explains in part the smaller number of RNA-only structures in the Protein Data Bank. Several approaches have been developed to aid the formation of well-ordered RNA crystals. The majority of these are construct-engineering techniques that aim to introduce crystal contacts to favor the formation of well-diffracting crystals. A typical example is the insertion of tetraloop-tetraloop receptor pairs into non-essential RNA segments to promote intermolecular association. Other methods of promoting crystallization involve chaperones and crystallization-friendly molecules that increase RNA stability and improve crystal packing. In this review, we discuss the various techniques that have been successfully used to facilitate crystal packing of RNA molecules, recent advances in construct engineering, and directions for future research in this vital aspect of RNA crystallography.

Annotation of the local context of the RNA secondary structure improves the classification and prediction of A-minors

Non-coding RNAs play a crucial role in various cellular processes in living organisms, and RNA functions heavily depend on molecule structures composed of stems, loops, and various tertiary motifs. Among those, the most frequent are A-minor interactions, which are often involved in the formation of more complex motifs such as kink-turns and pseudoknots. We present a novel classification of A-minors in terms of RNA secondary structure where each nucleotide of an A-minor is attributed to the stem or loop, and each pair of nucleotides is attributed to their relative position within the secondary structure. By analyzing classes of A-minors in known RNA structures, we found that the largest classes are mostly homogeneous and preferably localize with known A-minor co-motifs, e.g. tetraloop-tetraloop receptor and coaxial stacking. Detailed analysis of local A-minors within internal loops revealed a novel recurrent RNA tertiary motif, the across-bulged motif. Interestingly, the motif resembles the previously known GAAA/11nt motif but with the local adenines performing the role of the GAAA-tetraloop. By using machine learning, we show that particular classes of local A-minors can be predicted from sequence and secondary structure. The proposed classification is the first step toward automatic annotation of not only A-minors and their co-motifs but various types of RNA tertiary motifs as well.

Beyond ribose and phosphate: Selected nucleic acid modifications for structure-function investigations and therapeutic applications

Over the past 25 years, the acceleration of achievements in the development of oligonucleotide-based therapeutics has resulted in numerous new drugs making it to the market for the treatment of various diseases. Oligonucleotides with alterations to their scaffold, prepared with modified nucleosides and solid-phase synthesis, have yielded molecules with interesting biophysical properties that bind to their targets and are tolerated by the cellular machinery to elicit a therapeutic outcome. Structural techniques, such as crystallography, have provided insights to rationalize numerous properties including binding affinity, nuclease stability, and trends observed in the gene silencing. In this review, we discuss the chemistry, biophysical, and structural properties of a number of chemically modified oligonucleotides that have been explored for gene silencing.

Incorporating a Thiophosphate Modification into a Common RNA Tetraloop Motif Causes an Unanticipated Stability Boost

GNRA (N = A, C, G, or U; R = A or G) tetraloops are common RNA secondary structural motifs and feature a phosphate stacked atop a nucleobase. The rRNA sarcin/ricin loop (SRL) is capped by GApGA, and the phosphate p stacks on G. We recently found that regiospecific incorporation of a single dithiophosphate (PS2) but not a monothiophosphate (PSO) instead of phosphate in the backbone of RNA aptamers dramatically increases the binding affinity for their targets. In the RNA:thrombin complex, the key contribution to the 1000-fold tighter binding stems from an edge-on contact between PS2 and a phenylalanine ring. Here we investigated the consequences of replacing the SRL phosphate engaged in a face-on interaction with guanine with either PS2 or PSO for stability. We found that PS2···G and Rp-PSO···G contacts stabilize modified SRLs compared to the parent loop to unexpected levels: up to 6.3 °C in melting temperature T_m and -4.7 kcal/mol in ΔΔG°. Crystal structures demonstrate that the vertical distance to guanine for the closest sulfur is just 0.05 Å longer on average compared to that of oxygen despite the larger van der Waals radius of the former (1.80 Å for S vs 1.52 Å for O). The higher stability is enthalpy-based, and the negative charge as assessed by a neutral methylphosphonate modification plays only a minor role. Quantum mechanical/molecular mechanical calculations are supportive of favorable dispersion attraction interactions by sulfur making the dominant contribution. A stacking interaction between phosphate and guanine (SRL) or uracil (U-turn) is also found in newly classified RNA tetraloop families besides GNRA.

Molecular structure of a U•A-U-rich RNA triple helix with 11 consecutive base triples

Three-dimensional structures have been solved for several naturally occurring RNA triple helices, although all are limited to six or fewer consecutive base triples, hindering accurate estimation of global and local structural parameters. We present an X-ray crystal structure of a right-handed, U•A-U-rich RNA triple helix with 11 continuous base triples. Due to helical unwinding, the RNA triple helix spans an average of 12 base triples per turn. The double helix portion of the RNA triple helix is more similar to both the helical and base step structural parameters of A′-RNA rather than A-RNA. Its most striking features are its wide and deep major groove, a smaller inclination angle and all three strands favoring a C3′-endo sugar pucker. Despite the presence of a third strand, the diameter of an RNA triple helix remains nearly identical to those of DNA and RNA double helices. Contrary to our previous modeling predictions, this structure demonstrates that an RNA triple helix is not limited in length to six consecutive base triples and that longer RNA triple helices may exist in nature. Our structure provides a starting point to establish structural parameters of the so-called ‘ideal’ RNA triple helix, analogous to A-RNA and B-DNA double helices.

Nuclear Magnetic Resonance of Single-Stranded RNAs and DNAs of CAAU and UCAAUC as Benchmarks for Molecular Dynamics Simulations

RNA and DNA are rapidly emerging as targets for therapeutics and as potential frameworks for nanotechnology. Accurate methods for predicting and designing structures and dynamics of nucleic acids would accelerate progress in these and other applications. Suitable approximations for modeling nucleic acids are being developed but require validation against disparate experimental observations. Here, nuclear magnetic resonance spectra for RNA and DNA single strands, CAAU and UCAAUC, are used as benchmarks to test molecular dynamics simulations with AMBER force fields OL3 and ROC-RNA for RNA and BSC1 for DNA. A detailed scheme for making comparisons is also presented. The results reflect recent progress in approximations and reveal remaining challenges.

SAM-VI riboswitch structure and signature for ligand discrimination

Riboswitches are metabolite-sensing, conserved domains located in non-coding regions of mRNA that are central to regulation of gene expression. Here we report the first three-dimensional structure of the recently discovered S-adenosyl-L-methionine responsive SAM-VI riboswitch. SAM-VI adopts a unique fold and ligand pocket that are distinct from all other known SAM riboswitch classes. The ligand binds to the junctional region with its adenine tightly intercalated and Hoogsteen base-paired. Furthermore, we reveal the ligand discrimination mode of SAM-VI by additional X-ray structures of this riboswitch bound to S-adenosyl-L-homocysteine and a synthetic ligand mimic, in combination with isothermal titration calorimetry and fluorescence spectroscopy to explore binding thermodynamics and kinetics. The structure is further evaluated by analysis of ligand binding to SAM-VI mutants. It thus provides a thorough basis for developing synthetic SAM cofactors for applications in chemical and synthetic RNA biology.

Structure of an RNA helix with pyrimidine mismatches and cross-strand stacking

The structure of a 22-base-pair RNA helix with mismatched pyrimidine base pairs is reported. The helix contains two symmetry-related CUG sequences: a triplet-repeat motif implicated in myotonic dystrophy type 1. The CUG repeat contains a U-U mismatch sandwiched between Watson-Crick pairs. Additionally, the center of the helix contains a dimerized UUCG motif with tandem pyrimidine (U-C/C-U) mismatches flanked by U-G wobble pairs. This region of the structure is significantly different from previously observed structures that share the same sequence and neighboring base pairs. The tandem pyrimidine mismatches are unusual and display sheared, cross-strand stacking geometries that locally constrict the helical width, a type of stacking previously associated with purines in internal loops. Thus, pyrimidine-rich regions of RNA have a high degree of structural diversity.

The role of sequence in altering the unfolding pathway of an RNA pseudoknot: a steered molecular dynamics study

This work highlights a sequence dependent unfolding pathway of an RNA pseudoknot under force-induced pulling conditions.

InterRNA: a database of base interactions in RNA structures

A major component of RNA structure stabilization are the hydrogen bonded interactions between the base residues. The importance and biological relevance for large clusters of base interactions can be much more easily investigated when their occurrences have been systematically detected, catalogued and compared. In this paper, we describe the database InterRNA (INTERactions in RNA structures database—http://mfrlab.org/interrna/) that contains records of known RNA 3D motifs as well as records for clusters of bases that are interconnected by hydrogen bonds. The contents of the database were compiled from RNA structural annotations carried out by the NASSAM (http://mfrlab.org/grafss/nassam) and COGNAC (http://mfrlab.org/grafss/cognac) computer programs. An analysis of the database content and comparisons with the existing corpus of knowledge regarding RNA 3D motifs clearly show that InterRNA is able to provide an extension of the annotations for known motifs as well as able to provide novel interactions for further investigations.

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 structure-based mechanism for tRNA and retroviral RNA remodelling during primer annealing

To prime reverse transcription, retroviruses require annealing of a transfer RNA molecule to the U5 primer binding site (U5-PBS) region of the viral genome. The residues essential for primer annealing are initially locked in intramolecular interactions; hence, annealing requires the chaperone activity of the retroviral nucleocapsid (NC) protein to facilitate structural rearrangements. Here we show that, unlike classical chaperones, the Moloney murine leukaemia virus NC uses a unique mechanism for remodelling: it specifically targets multiple structured regions in both the U5-PBS and tRNA(Pro) primer that otherwise sequester residues necessary for annealing. This high-specificity and high-affinity binding by NC consequently liberates these sequestered residues--which are exactly complementary--for intermolecular interactions. Furthermore, NC utilizes a step-wise, entropy-driven mechanism to trigger both residue-specific destabilization and residue-specific release. Our structures of NC bound to U5-PBS and tRNA(Pro) reveal the structure-based mechanism for retroviral primer annealing and provide insights as to how ATP-independent chaperones can target specific RNAs amidst the cellular milieu of non-target RNAs.

Purine biosynthetic intermediate-containing ribose-phosphate polymers as evolutionary precursors to RNA

The RNA world hypothesis proposes that RNA once functioned as the principal genetic material and biological catalyst. However, RNA is a complex molecule made up of phosphate, ribose, and nucleobase moieties, and its evolution is unclear. Yakhnin has proposed a period of prebiotic chemical evolution prior to the advent of replication and Darwinian evolution, in which macromolecules containing polyols joined by phosphodiester linkages underwent spontaneous transesterification reactions with selection for stability. Although he proposes that the nucleobases were obtained during this stage from less stable macromolecules, the ultimate source of the nucleobases is not addressed. We propose that the purine nucleobases arose in situ from simpler precursors attached to a ribose-phosphate backbone, and that the weaker and less specific intra- and interstrand interactions between these precursors were the forerunners to the base pairing and base stacking interactions of the modern RNA nucleobases. Further, in line with Granick's hypothesis of biosynthetic pathways recapitulating evolution, we propose that these simpler precursors were the same or similar to intermediates of the modern de novo purine biosynthetic pathway. We propose that successive nucleobase precursors formed progressively stronger interactions that stabilized the ribose-phosphate polymer, and that the increased stability of the parent polymer drove the selection and further chemical evolution of the purine nucleobases. Such interactions may have included hydrogen bonding between ribose hydroxyls, hydrogen bonding between carbonyl oxygens and protonated amine side groups, the intra- and interstrand coordination of metal cations, and the stacking of imidazole rings. Five of the eleven steps of the modern de novo purine biosynthetic pathway have previously been shown to have alternative nonenzymatic syntheses, while a sixth step has also been proposed to occur nonenzymatically, supporting a prebiotic origin for the pathway.

Ribonucleotides as nucleotide excision repair substrates

The incorporation of ribonucleotides in DNA has attracted considerable notice in recent years, since the pool of ribonucleotides can exceed that of the deoxyribonucleotides by at least 10-20-fold, and single ribonucleotide incorporation by DNA polymerases appears to be a common event. Moreover ribonucleotides are potentially mutagenic and lead to genome instability. As a consequence, errantly incorporated ribonucleotides are rapidly repaired in a process dependent upon RNase H enzymes. On the other hand, global genomic nucleotide excision repair (NER) in prokaryotes and eukaryotes removes damage caused by covalent modifications that typically distort and destabilize DNA through the production of lesions derived from bulky chemical carcinogens, such as polycyclic aromatic hydrocarbon metabolites, or via crosslinking. However, a recent study challenges this lesion-recognition paradigm. The work of Vaisman et al. (2013) [34] reveals that even a single ribonucleotide embedded in a deoxyribonucleotide duplex is recognized by the bacterial NER machinery in vitro. In their report, the authors show that spontaneous mutagenesis promoted by a steric-gate pol V mutant increases in uvrA, uvrB, or uvrC strains lacking rnhB (encoding RNase HII) and to a greater extent in an NER-deficient strain lacking both RNase HI and RNase HII. Using purified UvrA, UvrB, and UvrC proteins in in vitro assays they show that despite causing little distortion, a single ribonucleotide embedded in a DNA duplex is recognized and doubly-incised by the NER complex. We present the hypothesis to explain the recognition and/or verification of this small lesion, that the critical 2'-OH of the ribonucleotide - with its unique electrostatic and hydrogen bonding properties - may act as a signal through interactions with amino acid residues of the prokaryotic NER complex that are not possible with DNA. Such a mechanism might also be relevant if it were demonstrated that the eukaryotic NER machinery likewise incises an embedded ribonucleotide in DNA.

RNA modularity for synthetic biology

RNA molecules are highly modular components that can be used in a variety of contexts for building new metabolic, regulatory and genetic circuits in cells. The majority of synthetic RNA systems to date predominately rely on two-dimensional modularity. However, a better understanding and integration of three-dimensional RNA modularity at structural and functional levels is critical to the development of more complex, functional bio-systems and molecular machines for synthetic biology applications.

Comparative sequence and structure analysis reveals the conservation and diversity of nucleotide positions and their associated tertiary interactions in the riboswitches

The tertiary motifs in complex RNA molecules play vital roles to either stabilize the formation of RNA 3D structure or to provide important biological functionality to the molecule. In order to better understand the roles of these tertiary motifs in riboswitches, we examined 11 representative riboswitch PDB structures for potential agreement of both motif occurrences and conservations. A total of 61 unique tertiary interactions were found in the reference structures. In addition to the expected common A-minor motifs and base-triples mainly involved in linking distant regions the riboswitch structures three highly conserved variants of A-minor interactions called G-minors were found in the SAM-I and FMN riboswitches where they appear to be involved in the recognition of the respective ligand’s functional groups. From our structural survey as well as corresponding structure and sequence alignments, the agreement between motif occurrences and conservations are very prominent across the representative riboswitches. Our analysis provide evidence that some of these tertiary interactions are essential components to form the structure where their sequence positions are conserved despite a high degree of diversity in other parts of the respective riboswitches sequences. This is indicative of a vital role for these tertiary interactions in determining the specific biological function of riboswitch.

An RNA folding motif: GNRA tetraloop-receptor interactions

Nearly two decades after Westhof and Michel first proposed that RNA tetraloops may interact with distal helices, tetraloop–receptor interactions have been recognized as ubiquitous elements of RNA tertiary structure. The unique architecture of GNRA tetraloops (N=any nucleotide, R=purine) enables interaction with a variety of receptors, e.g., helical minor grooves and asymmetric internal loops. The most common example of the latter is the GAAA tetraloop–11 nt tetraloop receptor motif. Biophysical characterization of this motif provided evidence for the modularity of RNA structure, with applications spanning improved crystallization methods to RNA tectonics. In this review, we identify and compare types of GNRA tetraloop–receptor interactions. Then we explore the abundance of structural, kinetic, and thermodynamic information on the frequently occurring and most widely studied GAAA tetraloop–11 nt receptor motif. Studies of this interaction have revealed powerful paradigms for structural assembly of RNA, as well as providing new insights into the roles of cations, transition states and protein chaperones in RNA folding pathways. However, further research will clearly be necessary to characterize other tetraloop–receptor and long-range tertiary binding interactions in detail – an important milestone in the quantitative prediction of free energy landscapes for RNA folding.

Propensities for loop structures of RNA & DNA backbones

RNA oligonucleotides exhibit a large tendency to bend and form a loop conformation which is a major motif contributing to their complex three-dimensional structure. This is in contrast to DNA molecules that predominantly form the double-helix structure. In this paper we investigate by molecular dynamics simulation, as well as, by its combination with the replica-exchange method, the propensity of RNA chains containing the GCUAA pentaloop to form spontaneously a hairpin conformation. The results were then compared with those of analogous hybrid oligonucleotides in which the ribose groups in the loop-region were substituted by deoxyriboses. We find that the RNA oligomers exhibit a marginal excess stability to form loop structures. The equilibrium constant for opening the loop to an extended conformation is twice as large in the hybrid than it is in the RNA chain. Analyses of the hydrogen bonds indicate that the excess stability for forming a hairpin is a result of hydrogen bonds the 2'-hydroxyls in the loop region form with other groups in the loop. Of these hydrogen bonds, the most important is the hydrogen bond donated from the 2'-OH at the first position of the loop to N7 of adenine at the forth position. RNA and DNA backbones are characterized by different backbone dihedral angles and sugar puckering that can potentially facilitate or hamper the hydrogen bonds involving the 2'-OH. Nevertheless, the sugar puckerings of all the pentaloop nucleotides were not significantly different between the two chains displaying the C3'-endo conformation characteristic to the A-form double helix. All of the other backbone dihedrals also did not show any considerable difference in the loop-region except of the δ-dihedral. In this case, the RNA loop exhibited bimodal distributions corresponding to, both, the RNA and DNA backbones, whereas the loop of the hybrid chain behaved mostly as that of a DNA backbone. Thus, it is possible that the behavior of the δ-dihedrals in the loop-region of the RNA adopts conformations that facilitate the intra-nucleotide hydrogen bondings of the 2'-hydroxyls, and consequently renders loop structures in RNA more stable.

Structure of an A-form RNA duplex obtained by degradation of 6S RNA in a crystallization droplet

In the course of a crystallographic study of a 132 nt variant of Aquifex aeolicus 6S RNA, a crystal structure of an A-form RNA duplex containing 12 base pairs was solved at a resolution of 2.6 Å. In fact, the RNA duplex is part of the 6S RNA and was obtained by accidental but precise degradation of the 6S RNA in a crystallization droplet. 6S RNA degradation was confirmed by microscopic observation of crystals and gel electrophoresis of crystallization droplets. The RNA oligomers obtained form regular A-form duplexes containing three GoU wobble-type base pairs, one of which engages in intermolecular contacts through a ribose-zipper motif at the crystal-packing interface.

The GA-minor submotif as a case study of RNA modularity, prediction, and design

Complex natural RNAs such as the ribosome, group I and group II introns, and RNase P exemplify the fact that three-dimensional (3D) RNA structures are highly modular and hierarchical in nature. Tertiary RNA folding typically takes advantage of a rather limited set of recurrent structural motifs that are responsible for controlling bends or stacks between adjacent helices. Herein, the GA minor and related structural motifs are presented as a case study to highlight several structural and folding principles, to gain further insight into the structural evolution of naturally occurring RNAs, as well as to assist the rational design of artificial RNAs.

Ribose 2'-Hydroxyl Groups Stabilize RNA Hairpin Structures Containing GCUAA Pentaloop

The chemical structure of RNA and DNA is very similar; however, the three-dimensional conformation of these two nucleic acids is very different. Whereas the DNA adopts a repetitive structure of a double-stranded helix, RNA is primarily single stranded with a complex three-dimensional structure in which the hairpin is the most common secondary structure. Apart from the difference between uracil and thymine, the difference in the chemical structure between RNA and DNA is the presence of a hydroxyl group at position 2' of the sugar (ribose) instead of a hydrogen (deoxyribose). In this paper, we present molecular dynamics simulations addressing the contribution of 2'-hydroxyls to the stability of a GCUAA pentaloop motif. The results indicate that the 2'-hydroxyls stabilize the hairpin conformation of the GCUAA pentaloop relative to an analogous oligonucleotide in which the ribose sugars in the loop region were substituted with deoxyriboses. The magnitude of the stabilization was found to be 23.8 ± 4.1 kJ/mol using an alchemical mutations free energy method and 4.2 ± 6.5 kJ/mol using potential of mean force calculations. The latter indicates that in addition to its larger thermodynamic stability the RNA hairpin is also kinetically more stable. We find that the excess stability is a result of intrahairpin hydrogen bonds in the loop region between the 2'-hydroxyls and sugars, bases, and phosphates. The hydrogen bonds with the sugars and phosphates involve predominantly interactions with adjacent nucleotides. However, the hydrogen bonds with the bases involve also interactions between groups on opposite sides of the loop or with the middle base of the loop and are therefore likely to contribute significantly to the stability of the loop. Of these hydrogen bonds, the most frequent is observed between the 2'-hydroxyl at the first position of the pentaloop with N6/N7 of adenine at the forth position, as well as between the 2'-hydroxyl at position -1 with N6 of adenine at the fifth position. Our results contribute to the notion that one of the important roles of the ribose sugars in RNA is to facilitate hairpin formation.

Structural variation and uniformity among tetraloop-receptor interactions and other loop-helix interactions in RNA crystal structures

Tetraloop-receptor interactions are prevalent structural units in RNAs, and include the GAAA/11-nt and GNRA-minor groove interactions. In this study, we have compiled a set of 78 nonredundant loop-helix interactions from X-ray crystal structures, and examined them for the extent of their sequence and structural variation. Of the 78 interactions in the set, only four were classical GAAA/11-nt motifs, while over half (48) were GNRA-minor groove interactions. The GNRA-minor groove interactions were not a homogeneous set, but were divided into five subclasses. The most predominant subclass is characterized by two triple base pair interactions in the minor groove, flanked by two ribose zipper contacts. This geometry may be considered the "standard" GNRA-minor groove interaction, while the other four subclasses are alternative ways to form interfaces between a minor groove and tetraloop. The remaining 26 structures in the set of 78 have loops interacting with mostly idiosyncratic receptors. Among the entire set, a number of sequence-structure correlations can be identified, which may be used as initial hypotheses in predicting three-dimensional structures from primary sequences. Conversely, other sequence patterns are not predictive; for example, GAAA loop sequences and GG/CC receptors bind to each other with three distinct geometries. Finally, we observe an example of structural evolution in group II introns, in which loop-receptor motifs are substituted for each other while maintaining the larger three-dimensional geometry. Overall, the study gives a more complete view of RNA loop-helix interactions that exist in nature.

The right angle (RA) motif: a prevalent ribosomal RNA structural pattern found in group I introns

The right angle (RA) motif, previously identified in the ribosome and used as a structural module for nano-construction, is a recurrent structural motif of 13 nucleotides that establishes a 90° bend between two adjacent helices. Comparative sequence analysis was used to explore the sequence space of the RA motif within ribosomal RNAs in order to define its canonical sequence space signature. We investigated the sequence constraints associated with the RA signature using several artificial self-assembly systems. Thermodynamic and topological investigations of sequence variants associated with the RA motif in both minimal and expanded structural contexts reveal that the presence of a helix at the 3' end of the RA motif increases the thermodynamic stability and rigidity of the resulting three-helix junction domain. A search for the RA in naturally occurring RNAs as well as its experimental characterization led to the identification of the RA in groups IC1 and ID intron ribozymes, where it is suggested to play an integral role in stabilizing peripheral structural domains. The present study exemplifies the need of empirical analysis of RNA structural motifs for facilitating the rational design and structure prediction of RNAs.

NASSAM: a server to search for and annotate tertiary interactions and motifs in three-dimensional structures of complex RNA molecules

Similarities in the 3D patterns of RNA base interactions or arrangements can provide insights into their functions and roles in stabilization of the RNA 3D structure. Nucleic Acids Search for Substructures and Motifs (NASSAM) is a graph theoretical program that can search for 3D patterns of base arrangements by representing the bases as pseudo-atoms. The geometric relationship of the pseudo-atoms to each other as a pattern can be represented as a labeled graph where the pseudo-atoms are the graph’s nodes while the edges are the inter-pseudo-atomic distances. The input files for NASSAM are PDB formatted 3D coordinates. This web server can be used to identify matches of base arrangement patterns in a query structure to annotated patterns that have been reported in the literature or that have possible functional and structural stabilization implications. The NASSAM program is freely accessible without any login requirement at http://mfrlab.org/grafss/nassam/.

The molecular interactions that stabilize RNA tertiary structure: RNA motifs, patterns, and networks

RNA molecules adopt specific three-dimensional structures critical to their function. Many essential metabolic processes, including protein synthesis and RNA splicing, are carried out by RNA molecules with elaborate tertiary structures (e.g. 3QIQ, right). Indeed, the ribosome and self-splicing introns are complex RNA machines. But even the coding regions in messenger RNAs and viral RNAs are flanked by highly structured untranslated regions, which provide regulatory information necessary for gene expression. RNA tertiary structure is defined as the three-dimensional arrangement of RNA building blocks, which include helical duplexes, triple-stranded structures, and other components that are held together through connections collectively termed RNA tertiary interactions. The structural diversity of these interactions is now a subject of intense investigation, involving the techniques of NMR, X-ray crystallography, chemical genetics, and phylogenetic analysis. At the same time, many investigators are using biophysical techniques to elucidate the driving forces for tertiary structure formation and the mechanisms for its stabilization. RNA tertiary folding is promoted by maximization of base stacking, much like the hydrophobic effect that drives protein folding. RNA folding also requires electrostatic stabilization, both through charge screening and site binding of metals, and it is enhanced by desolvation of the phosphate backbone. In this Account, we provide an overview of the features that specify and stabilize RNA tertiary structure. A major determinant for overall tertiary RNA architecture is local conformation in secondary-structure junctions, which are regions from which two or more duplexes project. At junctions and other structures, such as pseudoknots and kissing loops, adjacent helices stack on one another, and these coaxial stacks play a major role in dictating the overall architectural form of an RNA molecule. In addition to RNA junction topology, a second determinant for RNA tertiary structure is the formation of sequence-specific interactions. Networks of triple helices, tetraloop-receptor interactions, and other sequence-specific contacts establish the framework for the overall tertiary fold. The third determinant of tertiary structure is the formation of stabilizing stacking and backbone interactions, and many are not sequence specific. For example, ribose zippers allow 2'-hydroxyl groups on different RNA strands to form networks of interdigitated hydrogen bonds, serving to seal strands together and thereby stabilize adjacent substructures. These motifs often require monovalent and divalent cations, which can interact diffusely or through chelation to specific RNA functional groups. As we learn more about the components of RNA tertiary structure, we will be able to predict the structures of RNA molecules from their sequences, thereby obtaining key information about biological function. Understanding and predicting RNA structure is particularly important given the recent discovery that although most of our genome is transcribed into RNA molecules, few of them have a known function. The prevalence of RNA viruses and pathogens with RNA genomes makes RNA drug discovery an active area of research. Finally, knowledge of RNA structure will facilitate the engineering of supramolecular RNA structures, which can be used as nanomechanical components for new materials. But all of this promise depends on a better understanding of the RNA parts list, and how the pieces fit together.

A new way to see RNA

Unlike proteins, the RNA backbone has numerous degrees of freedom (eight, if one counts the sugar pucker), making RNA modeling, structure building and prediction a multidimensional problem of exceptionally high complexity. And yet RNA tertiary structures are not infinite in their structural morphology; rather, they are built from a limited set of discrete units. In order to reduce the dimensionality of the RNA backbone in a physically reasonable way, a shorthand notation was created that reduced the RNA backbone torsion angles to two (η and θ, analogous to φ and ψ in proteins). When these torsion angles are calculated for nucleotides in a crystallographic database and plotted against one another, one obtains a plot analogous to a Ramachandran plot (the η/θ plot), with highly populated and unpopulated regions. Nucleotides that occupy proximal positions on the plot have identical structures and are found in the same units of tertiary structure. In this review, we describe the statistical validation of the η/θ formalism and the exploration of features within the η/θ plot. We also describe the application of the η/θ formalism in RNA motif discovery, structural comparison, RNA structure building and tertiary structure prediction. More than a tool, however, the η/θ formalism has provided new insights into RNA structure itself, revealing its fundamental components and the factors underlying RNA architectural form.

Quantum chemical studies of nucleic acids: can we construct a bridge to the RNA structural biology and bioinformatics communities?

In this feature article, we provide a side-by-side introduction for two research fields: quantum chemical calculations of molecular interaction in nucleic acids and RNA structural bioinformatics. Our main aim is to demonstrate that these research areas, while largely separated in contemporary literature, have substantial potential to complement each other that could significantly contribute to our understanding of the exciting world of nucleic acids. We identify research questions amenable to the combined application of modern ab initio methods and bioinformatics analysis of experimental structures while also assessing the limitations of these approaches. The ultimate aim is to attain valuable physicochemical insights regarding the nature of the fundamental molecular interactions and how they shape RNA structures, dynamics, function, and evolution.

The tertiary structure of group II introns: implications for biological function and evolution

Group II introns are some of the largest ribozymes in nature, and they are a major source of information about RNA assembly and tertiary structural organization. These introns are of biological significance because they are self-splicing mobile elements that have migrated into diverse genomes and played a major role in the genomic organization and metabolism of most life forms. The tertiary structure of group II introns has been the subject of many phylogenetic, genetic, biochemical and biophysical investigations, all of which are consistent with the recent crystal structure of an intact group IIC intron from the alkaliphilic eubacterium Oceanobacillus iheyensis. The crystal structure reveals that catalytic intron domain V is enfolded within the other intronic domains through an elaborate network of diverse tertiary interactions. Within the folded core, DV adopts an activated conformation that readily binds catalytic metal ions and positions them in a manner appropriate for reaction with nucleic acid targets. The tertiary structure of the group II intron reveals new information on motifs for RNA architectural organization, mechanisms of group II intron catalysis, and the evolutionary relationships among RNA processing systems. Guided by the structure and the wealth of previous genetic and biochemical work, it is now possible to deduce the probable location of DVI and the site of additional domains that contribute to the function of the highly derived group IIB and IIA introns.

Protonation states of the key active site residues and structural dynamics of the glmS riboswitch as revealed by molecular dynamics

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

Quantifying the relationship between sequence and three-dimensional structure conservation in RNA

Background

In recent years, the number of available RNA structures has rapidly grown reflecting the increased interest on RNA biology. Similarly to the studies carried out two decades ago for proteins, which gave the fundamental grounds for developing comparative protein structure prediction methods, we are now able to quantify the relationship between sequence and structure conservation in RNA.

Results

Here we introduce an all-against-all sequence- and three-dimensional (3D) structure-based comparison of a representative set of RNA structures, which have allowed us to quantitatively confirm that: (i) there is a measurable relationship between sequence and structure conservation that weakens for alignments resulting in below 60% sequence identity, (ii) evolution tends to conserve more RNA structure than sequence, and (iii) there is a twilight zone for RNA homology detection.

Discussion

The computational analysis here presented quantitatively describes the relationship between sequence and structure for RNA molecules and defines a twilight zone region for detecting RNA homology. Our work could represent the theoretical basis and limitations for future developments in comparative RNA 3D structure prediction.

Molecular dynamics simulations suggest that RNA three-way junctions can act as flexible RNA structural elements in the ribosome

We present extensive explicit solvent molecular dynamics analysis of three RNA three-way junctions (3WJs) from the large ribosomal subunit: the 3WJ formed by Helices 90–92 (H90–H92) of 23S rRNA; the 3WJ formed by H42–H44 organizing the GTPase associated center (GAC) of 23S rRNA; and the 3WJ of 5S rRNA. H92 near the peptidyl transferase center binds the 3′-CCA end of amino-acylated tRNA. The GAC binds protein factors and stimulates GTP hydrolysis driving protein synthesis. The 5S rRNA binds the central protuberance and A-site finger (ASF) involved in bridges with the 30S subunit. The simulations reveal that all three 3WJs possess significant anisotropic hinge-like flexibility between their stacked stems and dynamics within the compact regions of their adjacent stems. The A-site 3WJ dynamics may facilitate accommodation of tRNA, while the 5S 3WJ flexibility appears to be essential for coordinated movements of ASF and 5S rRNA. The GAC 3WJ may support large-scale dynamics of the L7/L12-stalk region. The simulations reveal that H42–H44 rRNA segments are not fully relaxed and in the X-ray structures they are bent towards the large subunit. The bending may be related to L10 binding and is distributed between the 3WJ and the H42–H97 contact.

The RNA backbone plays a crucial role in mediating the intrinsic stability of the GpU dinucleotide platform and the GpUpA/GpA miniduplex

The side-by-side interactions of nucleobases contribute to the organization of RNA, forming the planar building blocks of helices and mediating chain folding. Dinucleotide platforms, formed by side-by-side pairing of adjacent bases, frequently anchor helices against loops. Surprisingly, GpU steps account for over half of the dinucleotide platforms observed in RNA-containing structures. Why GpU should stand out from other dinucleotides in this respect is not clear from the single well-characterized H-bond found between the guanine N2 and the uracil O4 groups. Here, we describe how an RNA-specific H-bond between O2′(G) and O2P(U) adds to the stability of the GpU platform. Moreover, we show how this pair of oxygen atoms forms an out-of-plane backbone ‘edge’ that is specifically recognized by a non-adjacent guanine in over 90% of the cases, leading to the formation of an asymmetric miniduplex consisting of ‘complementary’ GpUpA and GpA subunits. Together, these five nucleotides constitute the conserved core of the well-known loop-E motif. The backbone-mediated intrinsic stabilities of the GpU dinucleotide platform and the GpUpA/GpA miniduplex plausibly underlie observed evolutionary constraints on base identity. We propose that they may also provide a reason for the extreme conservation of GpU observed at most 5′-splice sites.

RNA structural motifs that entail hydrogen bonds involving sugar-phosphate backbone atoms of RNA

The growing number of high-resolution crystal structures of large RNA molecules provides much information for understanding the principles of structural organization of these complex molecules. Several in-depth analyses of nucleobase-centered RNA structural motifs and backbone conformations have been published based on this information, including a systematic classification of base pairs by Leontis and Westhof. However, hydrogen bonds involving sugar-phosphate backbone atoms of RNA have not been analyzed systematically until recently, although such hydrogen bonds appear to be common both in local and tertiary interactions. Here we review some backbone structural motifs discussed in the literature and analyze a set of eight high-resolution multi-domain RNA structures. The analyzed RNAs are highly structured: among 5372 nucleotides in this set, 89% are involved in at least one "long-range" RNA-RNA hydrogen bond, i.e., hydrogen bonds between atoms in the same residue or sequential residues are ignored. These long-range hydrogen bonds frequently use backbone atoms as hydrogen bond acceptors, i.e., OP1, OP2, O2', O3', O4', or O5', or as a donor (2'OH). A surprisingly large number of such hydrogen bonds are found, considering that neither single-stranded nor double-stranded regions will contain such hydrogen bonds unless additional interactions with other residues exist. Among 8327 long-range hydrogen bonds found in this set of structures, 2811, or about one-third, are hydrogen bonds entailing RNA backbone atoms; they involve 39% of all nucleotides in the structures. The majority of them (2111) are hydrogen bonds entailing ribose hydroxyl groups, which can be used either as a donor or an acceptor; they constitute 25% of all hydrogen bonds and involve 31% of all nucleotides. The phosphate oxygens OP1 or OP2 are used as hydrogen bond acceptors in 12% of all nucleotides, and the ribose ring oxygen O4' and phosphodiester oxygens O3' and O5' are used in 4%, 4%, and 1% of all nucleotides, respectively. Distributions of geometric parameters and some examples of such hydrogen bonds are presented in this report. A novel motif involving backbone hydrogen bonds, the ribose-phosphate zipper, is also identified.

Tertiary motifs revealed in analyses of higher-order RNA junctions

RNA junctions are secondary-structure elements formed when three or more helices come together. They are present in diverse RNA molecules with various fundamental functions in the cell. To better understand the intricate architecture of three-dimensional (3D) RNAs, we analyze currently solved 3D RNA junctions in terms of base-pair interactions and 3D configurations. First, we study base-pair interaction diagrams for solved RNA junctions with 5 to 10 helices and discuss common features. Second, we compare these higher-order junctions to those containing 3 or 4 helices and identify global motif patterns such as coaxial stacking and parallel and perpendicular helical configurations. These analyses show that higher-order junctions organize their helical components in parallel and helical configurations similar to lower-order junctions. Their sub-junctions also resemble local helical configurations found in three- and four-way junctions and are stabilized by similar long-range interaction preferences such as A-minor interactions. Furthermore, loop regions within junctions are high in adenine but low in cytosine, and in agreement with previous studies, we suggest that coaxial stacking between helices likely forms when the common single-stranded loop is small in size; however, other factors such as stacking interactions involving noncanonical base pairs and proteins can greatly determine or disrupt coaxial stacking. Finally, we introduce the ribo-base interactions: when combined with the along-groove packing motif, these ribo-base interactions form novel motifs involved in perpendicular helix-helix interactions. Overall, these analyses suggest recurrent tertiary motifs that stabilize junction architecture, pack helices, and help form helical configurations that occur as sub-elements of larger junction networks. The frequent occurrence of similar helical motifs suggest nature's finite and perhaps limited repertoire of RNA helical conformation preferences. More generally, studies of RNA junctions and tertiary building blocks can ultimately help in the difficult task of RNA 3D structure prediction.

Mapping metal-binding sites in the catalytic domain of bacterial RNase P RNA

Ribonuclease P (RNase P) is a ribonucleoprotein enzyme that contains a universally conserved, catalytically active RNA component. RNase P RNA requires divalent metal ions for folding, substrate binding, and catalysis. Despite recent advances in understanding the structure of RNase P RNA, no comprehensive analysis of metal-binding sites has been reported, in part due to the poor crystallization properties of this large RNA. We have developed an abbreviated yet still catalytic construct, Bst P7Delta RNA, which contains the catalytic domain of the bacterial RNase P RNA and has improved crystallization properties. We use this mutant RNA as well as the native RNA to map metal-binding sites in the catalytic core of the bacterial RNase P RNA, by anomalous scattering in diffraction analysis. The results provide insight into the interplay between RNA structure and focalization of metal ions, and a structural basis for some previous biochemical observations with RNase P. We use electrostatic calculations to extract the potential functional significance of these metal-binding sites with respect to binding Mg(2+). The results suggest that with at least one important exception of specific binding, these sites mainly map areas of diffuse association of magnesium ions.

Analysis of natural variants of the hepatitis C virus internal ribosome entry site reveals that primary sequence plays a key role in cap-independent translation

The HCV internal ribosome entry site (IRES) spans a region of ∼340 nt that encompasses most of the 5′ untranslated region (5′UTR) of the viral mRNA and the first 24–40 nt of the core-coding region. To investigate the implication of altering the primary sequence of the 5′UTR on IRES activity, naturally occurring variants of the 5′UTR were isolated from clinical samples and analyzed. The impact of the identified mutations on translation was evaluated in the context of RLuc/FLuc bicistronic RNAs. Results show that depending on their location within the RNA structure, these naturally occurring mutations cause a range of effects on IRES activity. However, mutations within subdomain IIId hinder HCV IRES-mediated translation. In an attempt to explain these data, the dynamic behavior of the subdomain IIId was analyzed by means of molecular dynamics (MD) simulations. Despite the loss of function, MD simulations predicted that mutant G266A/G268U possesses a structure similar to the wt-RNA. This prediction was validated by analyzing the secondary structure of the isolated IIId RNAs by circular dichroism spectroscopy in the presence or absence of Mg²⁺ ions. These data strongly suggest that the primary sequence of subdomain IIId plays a key role in HCV IRES-mediated translation.