Characterisation of hydrogen bonding networks in RNAs via magic angle spinning solid state NMR spectroscopy

It is demonstrated that the spatial proximity of ¹H nuclei in hydrogen bonded base-pairs in RNAs can be conveniently mapped via magic angle spinning solid state NMR experiments involving proton spin diffusion driven chemical shift correlation of low gamma nuclei such as the imino and amino nitrogens of nucleic acid bases. As different canonical and non-canonical base-pairing schemes encountered in nucleic acids are characterised by topologically different networks of proton dipolar couplings, different base-pairing schemes lead to characteristic cross-peak intensity patterns in such correlation spectra. The method was employed in a study of a 100 kDa RNA composed of 97 CUG repeats, or (CUG)₉₇ that has been implicated in the neuromuscular disease myotonic dystrophy. ¹⁵N–¹⁵N chemical shift correlation studies confirm the presence of Watson–Crick GC base pairs in (CUG)₉₇.     

Stacking geometry for non‐canonical G:U wobble base pair containing dinucleotide sequences in RNA: dispersion‐corrected DFT‐D study

Emergence of thousands of crystal structures of noncoding RNA molecules indicates its structural and functional diversity. RNA function is based upon a large variety of structural elements which are specifically assembled in the folded molecules. Along with the canonical Watson‐Crick base pairs, different orientations of the bases to form hydrogen‐bonded non‐canonical base pairs have also been observed in the available RNA structures. Frequencies of occurrences of different non‐canonical base pairs in RNA indicate their important role to maintain overall structure and functions of RNA. There are several reports on geometry and energetic stabilities of these non‐canonical base pairs. However, their stacking geometry and stacking stability with the neighboring base pairs are not well studied. Among the different non‐canonical base pairs, the G:U wobble base pair (G:U W:WC) is most frequently observed in the RNA double helices. Using quantum chemical method and available experimental data set we have studied the stacking geometry of G:U W:WC base pair containing dinucleotide sequences in roll‐slide parameters hyperspace for different values of twist. This study indicates that the G:U W:WC base pair can stack well with the canonical base pairs giving rise to large interaction energy. The overall preferred stacking geometry in terms of roll, twist and slide for the eleven possible dinucleotide sequences is seen to be quite dependent on their sequences. © 2015 Wiley Periodicals, Inc. Biopolymers 103: 328–338, 2015.     

RNA structure and dynamics: A base pairing perspective

RNA is now known to possess various structural, regulatory and enzymatic functions for survival of cellular organisms. Functional RNA structures are generally created by three-dimensional organization of small structural motifs, formed by base pairing between self-complementary sequences from different parts of the RNA chain. In addition to the canonical Watson–Crick or wobble base pairs, several non-canonical base pairs are found to be crucial to the structural organization of RNA molecules. They appear within different structural motifs and are found to stabilize the molecule through long-range intra-molecular interactions between basic structural motifs like double helices and loops. These base pairs also impart functional variation to the minor groove of A-form RNA helices, thus forming anchoring site for metabolites and ligands. Non-canonical base pairs are formed by edge-to-edge hydrogen bonding interactions between the bases. A large number of theoretical studies have been done to detect and analyze these non-canonical base pairs within crystal or NMR derived structures of different functional RNA. Theoretical studies of these isolated base pairs using ab initio quantum chemical methods as well as molecular dynamics simulations of larger fragments have also established that many of these non-canonical base pairs are as stable as the canonical Watson–Crick base pairs. This review focuses on the various structural aspects of non-canonical base pairs in the organization of RNA molecules and the possible applications of these base pairs in predicting RNA structures with more accuracy.     

Ion-induced folding of a kink turn that departs from the conventional sequence

Kink turns (k-turns) are important structural motifs that create a sharp axial bend in RNA. Most conform to a consensus in which a three-nucleotide bulge is followed by consecutive G•A and A•G base pairs, and when these G•A pairs are modified in vitro this generally leads to a failure to adopt the k-turn conformation. Kt-23 in the 30S ribosomal subunit of Thermus thermophilus is a rare exception in which the bulge-distal A•G pair is replaced by a non-Watson–Crick A•U pair. In the context of the ribosome, Kt-23 adopts a completely conventional k-turn geometry. We show here that this sequence is induced to fold into a k-turn structure in an isolated RNA duplex by Mg²⁺ or Na⁺ ions. Therefore, the Kt-23 is intrinsically stable despite lacking the key A•G pair; its formation requires neither tertiary interactions nor protein binding. Moreover, the Kt-23 k-turn is stabilized by the same critical hydrogen-bonding interactions within the core of the structure that are found in more conventional sequences such as the near-consensus Kt-7. T. thermophilus Kt-23 has two further non-Watson–Crick base pairs within the non-canonical helix, three and four nucleotides from the bulge, and we find that the nature of these pairs influences the ability of the RNA to adopt k-turn conformation, although the base pair adjacent to the A•U pair is more important than the other.     

Identification and assignment of base pairs in three helical stems of Torulopsis utilis ribosomal 5S RNA and its RNase T1 and RNase T2 cleaved fragments via 500-MHz proton homonuclear overhauser enhancements

Imino proton resonances in the downfield region (10-14 ppm) of the 500-MHz 1H NMR spectrum of Torulopsis utilis 5S RNA are identified (A X U, G X C, or G X U) and assigned to base pairs in helices I, IV, and V via analysis of homonuclear Overhauser enhancements (NOE) from intact T. utilis 5S RNA, its RNase T1 and RNase T2 digested fragments, and a second yeast (Saccharomyces cerevisiae) 5S RNA whose nucleotide sequence differs at only six residues from that of T. utilis 5S RNA. The near-identical chemical shifts and NOE behavior of most of the common peaks from these four RNAs strongly suggest that helices I, IV, and V retain the same conformation after RNase digestion and that both T. utilis and S. cerevisiae 5S RNAs share a common secondary and tertiary structure. Of the four G X U base pairs identified in the intact 5S RNA, two are assigned to the terminal stem (helix I) and the other two to helices IV and V. Seven of the nine base pairs of the terminal stem have been assigned. Our experimental demonstration of a G X U base pair in helix V supports the 5S RNA secondary structural model of Luehrsen and Fox [Luehrsen, K. R., & Fox, G.E. (1981) Proc. Natl. Acad. Sci. U.S.A. 78, 2150-2154]. Finally, the base-pair proton peak assigned to the terminal G X U in helix V of the RNase T2 cleaved fragment is shifted downfield from that in the intact 5S RNA, suggesting that helices I and V may be coaxial in intact T. utilis 5S RNA.     

The G x U wobble base pair. A fundamental building block of RNA structure crucial to RNA function in diverse biological systems

The G x U wobble base pair is a fundamental unit of RNA secondary structure that is present in nearly every class of RNA from organisms of all three phylogenetic domains. It has comparable thermodynamic stability to Watson-Crick base pairs and is nearly isomorphic to them. Therefore, it often substitutes for G x C or A x U base pairs. The G x U wobble base pair also has unique chemical, structural, dynamic and ligand-binding properties, which can only be partially mimicked by Watson-Crick base pairs or other mispairs. These features mark sites containing G x U pairs for recognition by proteins and other RNAs and allow the wobble pair to play essential functional roles in a remarkably wide range of biological processes.     

Morpholino spin-labeling for base-pair sequencing of a 3'-terminal RNA stem by proton homonuclear Overhauser enhancements: yeast ribosomal 5S RNA

Base-pair sequences for 5S and 5.8S RNAs are not readily extracted from proton homonuclear nuclear Overhauser enhancement (NOE) connectivity experiments alone, due to extensive peak overlap in the downfield (11-15 ppm) proton NMR spectrum. In this paper, we introduce a new method for base-pair proton peak assignment for ribosomal RNAs, based upon the distance-dependent broadening of the resonances of base-pair protons spatially proximal to a paramagnetic group. Introduction of a nitroxide spin-label covalently attached to the 3'-terminal ribose provides an unequivocal starting point for base-pair hydrogen-bond proton NMR assignment. Subsequent NOE connectivities then establish the base-pair sequence for the terminal stem of a 5S RNA. Periodate oxidation of yeast 5S RNA, followed by reaction with 4-amino-2,2,6,6-tetramethylpiperidinyl-1-oxy (TEMPO-NH2) and sodium borohydride reduction, produces yeast 5S RNA specifically labeled with a paramagnetic nitroxide group at the 3'-terminal ribose. Comparison of the 500-MHz 1H NMR spectra of native and 3'-terminal spin-labeled yeast 5S RNA serves to identify the terminal base pair (G1 . C120) and its adjacent base pair (G2 . U119) on the basis of their proximity to the 3'-terminal spin-label. From that starting point, we have then identified (G . C, A . U, or G . U) and sequenced eight of the nine base pairs in the terminal helix via primary and secondary NOE's.     

Secondary structure and stability of the selenocysteine insertion sequences (SECIS) for human thioredoxin reductase and glutathione peroxidase

We have used high resolution NMR and thermodynamics to characterize the secondary structure and stability of the selenocysteine insertion sequences (SECIS) of human glutathione peroxidase (58 nt) and thioredoxin reductase (51 nt). These sequences are members of the two classes of SECIS recently identified with two distinct structures capable of directing selenocysteine incorporation into proteins in eukaryotes. UV melting experiments showed a single cooperative and reversible transition for each RNA, which indicates the presence of stable secondary structures. Despite their large size, the RNAs gave well resolved NMR spectra for the exchangeable protons. Using NOESY, the imino protons as well as the cytosine amino protons of all of the Watson-Crick base pairs were assigned. In addition, a number of non-canonical base pairs including the wobble G.U pairs were identified. The interbase-pair NOEs allowed definition of the hydrogen-bonded structure of the oligonucleotides, providing an experimental model of the secondary structure of these elements. The derived secondary structures are consistent with several features of the predicted models, but with some important differences, especially regarding the conserved sequence motifs.     

Stacking geometry between two sheared Watson-Crick basepairs: Computational chemistry and bioinformatics based prediction

BackgroundMolecular modeling of RNA double helices is possible using most probable values of basepair parameters obtained from crystal structure database. The A:A w:wC non-canonical basepair, involving Watson-Crick edges of two Adenines in cis orientation, appears quite frequently in database. Bimodal distribution of its Shear, due to two different H-bonding schemes, introduces the confusion in assigning most the probable value. Its effect is pronounced when the A:A w:wC basepair stacks on Sheared wobble G:U W:WC basepairs.MethodsWe employed molecular dynamics simulations of three possible double helices with GAG, UAG and GAU sequence motifs at their centers and quantum chemical calculation for non-canonical A:A w:wC basepair stacked on G:U W:WC basepair.ResultsWe noticed stable structures of GAG motif with specifically negative Shear of the A:A basepair but stabilities of the other motifs were not found with A:A w:wC basepairing. Hybrid DFT-D and MP2 stacking energy analyses on dinucleotide step sequences, A:A w:wC::G:U W:WC and A:A w:wC::U:G W:WC reveal that viable orientation of A:A::G:U prefers one of the H-bonding modes with negative Shear, supported by crystal structure database. The A:A::U:G dinucleotide, however, prefers structure with only positive Shear.ConclusionsThe quantum chemical calculations explain why MD simulations of GAG sequence motif only appear stable. In the cases of the GAU and UAG motifs “tug of war” situation between positive and negative Shears of A:A w:wC basepair induces conformational plasticity.General significanceWe have projected comprehensive reason behind the promiscuous nature of A:A w:wC basepair which brings occasional structural plasticity.     

Solution Structure of a GAAG Tetraloop in Helix 6 of SRP RNA from Pyrococcus furiosus

The NMR structure of a 12-mer RNA derived from the helix 6 of SRP RNA from Pyrococcus furiosus, whose loop-closing base pair is U:G, was determined, and the structural and thermodynamic properties of the RNA were compared with those of a mutant RNA with the C:G closing base pair. Although the structures of the two RNAs are similar to each other and adopt the GNRR motif, the conformational stabilities are significantly different to each other. It was suggested that weaker stacking interaction of the GAAG loop with the U:G closing base pair in 12-mer RNA causes the lower conformational stability.     

Isoalloxazine derivatives promote photocleavage of natural RNAs at G.U base pairs embedded within helices.

We have recently shown that isoalloxazine derivatives are able to photocleave RNA specifically at G.U base pairs embedded within a helical stack. The reaction involves the selective molecular recognition of G.U base pairs by the isoalloxazine ring and the removal of one nucleoside downstream of the uracil residue. Divalent metal ions are absolutely required for cleavage. Here we extend our studies to complex natural RNA molecules with known secondary and tertiary structures, such as tRNAs and a group I intron (td). G.U pairs were cleaved in accordance with the phylogenetically and experimentally derived secondary and tertiary structures. Tandem G.U pairs or certain G.U pairs located at a helix extremity were not affected. These new cleavage data, together with the RNA crystal structure, allowed us to perform molecular dynamics simulations to provide a structural basis for the observed specificity. We present a stable structural model for the ternary complex of the G. U-containing helical stack, the isoalloxazine molecule and a metal ion. This model provides significant new insight into several aspects of the cleavage phenomenon, mechanism and specificity for G. U pairs. Our study shows that in large natural RNAs a secondary structure motif made of an unusual base pair can be recognized and cleaved with high specificity by a low molecular weight molecule. This photocleavage reaction thus opens up the possibility of probing the accessibility of G.U base pairs, which are endowed with specific structural and functional roles in numerous structured and catalytic RNAs and interactions of RNA with proteins, in folded RNAs.     

Diversity of base-pair conformations and their occurrence in rRNA structure and RNA structural motifs

In addition to the canonical base-pairs comprising the standard Watson-Crick (C:G and U:A) and wobble U:G conformations, an analysis of the base-pair types and conformations in the rRNAs in the high-resolution crystal structures of the Thermus thermophilus 30S and Haloarcula marismortui 50S ribosomal subunits has identified a wide variety of non-canonical base-pair types and conformations. However, the existing nomenclatures do not describe all of the observed non-canonical conformations or describe them with some ambiguity. Thus, a standardized system is required to classify all of these non-canonical conformations appropriately. Here, we propose a new, simple and systematic nomenclature that unambiguously classifies base-pair conformations occurring in base-pairs, base-triples and base-quadruples that are associated with secondary and tertiary interactions. This system is based on the topological arrangement of the two bases and glycosidic bonds in a given base-pair. Base-pairs in the internal positions of regular secondary structure helices usually form with canonical base-pair groups (C:G, U:A, and U:G) and canonical conformations (C:G WC, U:A WC, and U:G Wb). In contrast, non-helical base-pairs outside of regular structure helices usually have non-canonical base-pair groups and conformations. In addition, many non-helical base-pairs are involved in RNA motifs that form a defined set of non-canonical conformations. Thus, each rare non-canonical conformation may be functionally and structurally important. Finally, the topology-based isostericity of base-pair conformations can rationalize base-pair exchanges in the evolution of RNA molecules.     

An intronic splicing enhancer binds U1 snRNPs to enhance splicing and select 5' splice sites

Intronic G triplets are frequently located adjacent to 5' splice sites in vertebrate pre-mRNAs and have been correlated with splicing efficiency and specificity via a mechanism that activates upstream 5' splice sites in exons containing duplicated sites (26). Using an intron dependent upon G triplets for maximal activity and 5' splice site specificity, we determined that these elements bind U1 snRNPs via base pairing with U1 RNA. This interaction is novel in that it uses nucleotides 8 to 10 of U1 RNA and is independent of nucleotides 1 to 7. In vivo functionality of base pairing was documented by restoring activity and specificity to mutated G triplets through compensating U1 RNA mutations. We suggest that the G-rich region near vertebrate 5' splice sites promotes accurate splice site recognition by recruiting the U1 snRNP.     

Thermodynamic insights into 2-thiouridine-enhanced RNA hybridization

Nucleobase modifications dramatically alter nucleic acid structure and thermodynamics. 2-thiouridine (s²U) is a modified nucleobase found in tRNAs and known to stabilize U:A base pairs and destabilize U:G wobble pairs. The recently reported crystal structures of s²U-containing RNA duplexes do not entirely explain the mechanisms responsible for the stabilizing effect of s²U or whether this effect is entropic or enthalpic in origin. We present here thermodynamic evaluations of duplex formation using ITC and UV thermal denaturation with RNA duplexes containing internal s²U:A and s²U:U pairs and their native counterparts. These results indicate that s²U stabilizes both duplexes. The stabilizing effect is entropic in origin and likely results from the s²U-induced preorganization of the single-stranded RNA prior to hybridization. The same preorganizing effect is likely responsible for structurally resolving the s²U:U pair-containing duplex into a single conformation with a well-defined H-bond geometry. We also evaluate the effect of s²U on single strand conformation using UV- and CD-monitored thermal denaturation and on nucleoside conformation using ¹H NMR spectroscopy, MD and umbrella sampling. These results provide insights into the effects that nucleobase modification has on RNA structure and thermodynamics and inform efforts toward improving both ribozyme-catalyzed and nonenzymatic RNA copying.     

Database of non-canonical base pairs found in known RNA structures

Atomic resolution RNA structures are being published at an increasing rate. It is common to find a modest number of non-canonical base pairs in these structures in addition to the usual Watson-Crick pairs. This database summarizes the occurrence of these rare base pairs in accordance with standard nomenclature. The database, http://prion.bchs.uh.edu/, contains information such as sequence context, sugar pucker conformation, anti / syn base conformations, chemical shift, p K (a)values, melting temperature and free energy. Of the 29 anticipated pairs with two or more hydrogen bonds, 20 have been encountered to date. In addition, four unexpected pairs with two hydrogen bonds have been reported bringing the total to 24. Single hydrogen bond versions of five of the expected geometries have been encountered among the single hydrogen bond interactions. In addition, 18 different types of base triplets have been encountered, each of which involves three to six hydrogen bonds. The vast majority of the rare base pairs are antiparallel with the bases in the anti configuration relative to the ribose. The most common are the GU wobble, the Sheared GA pair, the Reverse Hoogsteen pair and the GA imino pair.     

Comparative analysis of hairpin ribozyme structures and interference data

Great strides in understanding the molecular underpinnings of RNA catalysis have been achieved with advances in RNA structure determination by NMR spectroscopy and X-ray crystallography. Despite these successes the functional relevance of a given structure can only be assessed upon comparison with biochemical studies performed on functioning RNA molecules. The hairpin ribozyme presents an excellent case study for such a comparison. The active site is comprised of two stems each with an internal loop that forms a series of non-canonical base pairs. These loops dock into each other to create an active site for catalysis. Recently, three independent structures have been determined for this catalytic RNA, including two NMR structures of the isolated loop A and loop B stems and a high-resolution crystal structure of both loops in a docked conformation. These structures differ significantly both in their tertiary fold and the nature of the non-canonical base pairs formed within each loop. Several of the chemical groups required to achieve a functioning hairpin ribozyme have been determined by nucleotide analog interference mapping (NAIM). Here we compare the three hairpin structures with previously published NAIM data to assess the convergence between the structural and functional data. While there is significant disparity between the interference data and the individual NMR loop structures, there is almost complete congruity with the X-ray structure. The only significant differences cluster around an occluded pocket adjacent to the scissile phosphate. These local differences may suggest a role for these atoms in the transition state, either directly in chemistry or via a local structural rearrangement.     

Globular proteins,GU wobbling,and the evolution of the genetic code

Summary It has previously been shown that the formation of GU base pairs in RNA copying processes leads to an accumulation of G and U in both strands of the replicating RNA, which results in a non-random distribution of base triplets. In the present paper, this distribution is calculated, and, using the ²-test, a correlation between the distribution of triplets and the amino acid composition of the evolutionarily conservative interior regions of selected globular proteins is established.It is suggested that GU wobbling in early replication of RNA could have led to the observed amino acid composition of present-day protein interiors. If this hypothesis is correct, the GU wobbling must have been very extensive in the imprecisely replicating RNA, even reaching values close to the critical for stability of its double-helical structure. Implications of the hypothesis both for the evolution of the genetic code and of proteins are discussed.     

Structural basis for a lethal mutation in U6 RNA

U6 RNA is essential for nuclear pre-mRNA splicing and has been implicated directly in catalysis of intron removal. The U80G mutation at the essential magnesium binding site of the U6 3' intramolecular stem-loop region (ISL) is lethal in yeast. To further understand the structure and function of the U6 ISL, we have investigated the structural basis for the lethal U80G mutation by NMR and optical spectroscopy. The NMR structure reveals that the U80G mutation causes a structural rearrangement within the ISL resulting in the formation of a new Watson-Crick base pair (C67 x G80), and disrupts a protonated C67 x A79 wobble pair that forms in the wild-type structure. Despite the structural change, the accessibility of the metal binding site is unperturbed, and cadmium titration produces similar phosphorus chemical shift changes for both the U80G mutant and wild-type RNAs. The thermodynamic stability of the U80G mutant is significantly increased (Delta Delta G(fold) = -3.6 +/- 1.9 kcal/mol), consistent with formation of the Watson-Crick pair. Our structural and thermodynamic data, in combination with previous genetic data, suggest that the lethal basis for the U80G mutation is stem-loop hyperstabilization. This hyperstabilization may prevent the U6 ISL melting and rearrangement necessary for association with U4.     

NCIR: a database of non-canonical interactions in known RNA structures

The secondary and tertiary structure of an RNA molecule typically includes a number of non-canonical base–base interactions. The known occurrences of these interactions are tabulated in the NCIR database, which can be accessed from http://prion.bchs.uh.edu/bp_type/. The number of examples is now over 1400, which is an increase of >700% since the database was first published. This dramatic increase reflects the addition of data from the recently published crystal structures of the 50S (2.4 Å) and 30S (3.0 Å) ribosomal subunits. In addition, non-canonical interactions observed in published crystal and NMR structures of tRNAs, group I introns, ribozymes, RNA aptamers and synthetic oligonucleotides are included. Properties associated with these interactions, such as sequence context, sugar pucker conformation, glycosidic angle conformation, melting temperature, chemical shift and free energy, are also reported when available. Out of the 29 anticipated pairs with at least two hydrogen bonds, 28 have been observed to date. In addition, several novel examples, not generally predicted, have also been encountered, bringing the total of such pairs to 36. Added to this list are a variety of single, bifurcated, triple and quadruple interactions. The most common non-canonical pairs are the sheared GA, GA imino, AU reverse Hoogsteen, and the GU and AC wobble pairs. The most frequent triple interaction connects N3 of an A with the amino of a G that is also involved in a standard Watson–Crick pair.     

Automated and assisted RNA resonance assignment using NMR chemical shift statistics

The three-dimensional structure determination of RNAs by NMR spectroscopy relies on chemical shift assignment, which still constitutes a bottleneck. In order to develop more efficient assignment strategies, we analysed relationships between sequence and ¹H and ¹³C chemical shifts. Statistics of resonances from regularly Watson–Crick base-paired RNA revealed highly characteristic chemical shift clusters. We developed two approaches using these statistics for chemical shift assignment of double-stranded RNA (dsRNA): a manual approach that yields starting points for resonance assignment and simplifies decision trees and an automated approach based on the recently introduced automated resonance assignment algorithm FLYA. Both strategies require only unlabeled RNAs and three 2D spectra for assigning the H2/C2, H5/C5, H6/C6, H8/C8 and H1′/C1′ chemical shifts. The manual approach proved to be efficient and robust when applied to the experimental data of RNAs with a size between 20 nt and 42 nt. The more advanced automated assignment approach was successfully applied to four stem-loop RNAs and a 42 nt siRNA, assigning 92–100% of the resonances from dsRNA regions correctly. This is the first automated approach for chemical shift assignment of non-exchangeable protons of RNA and their corresponding ¹³C resonances, which provides an important step toward automated structure determination of RNAs.