GU Pairs Research Articles

Anionic G•U pairs in bacterial ribosomal rRNAs.

Wobble GU pairs (or G•U) occur frequently within double-stranded RNA helices interspersed between standard G=C and A-U Watson-Crick pairs. Another type of G•U pair interacting via their Watson-Crick edges has been observed in the A site of ribosome structures between a modified U34 in the tRNA anticodon triplet and G + 3 in the mRNA. In such pairs, the electronic structure of the U is changed with a negative charge on N3(U), resulting in two H-bonds between N1(G)…O4(U) and N2(G)…N3(U). Here, we report that such pairs occur in other highly conserved positions in ribosomal RNAs of bacteria in the absence of U modification. An anionic cis Watson-Crick G•G pair is also observed and well conserved in the small subunit. These pairs are observed in tightly folded regions.

Nearest neighbor rules for RNA helix folding thermodynamics: improved end effects.

Nearest neighbor parameters for estimating the folding stability of RNA secondary structures are in widespread use. For helices, current parameters penalize terminal AU base pairs relative to terminal GC base pairs. We curated an expanded database of helix stabilities determined by optical melting experiments. Analysis of the updated database shows that terminal penalties depend on the sequence identity of the adjacent penultimate base pair. New nearest neighbor parameters that include this additional sequence dependence accurately predict the measured values of 271 helices in an updated database with a correlation coefficient of 0.982. This refined understanding of helix ends facilitates fitting terms for base pair stacks with GU pairs. Prior parameter sets treated 5′GGUC3′ paired to 3′CUGG5′ separately from other 5′GU3′/3′UG5′ stacks. The improved understanding of helix end stability, however, makes the separate treatment unnecessary. Introduction of the additional terms was tested with three optical melting experiments. The average absolute difference between measured and predicted free energy changes at 37°C for these three duplexes containing terminal adjacent AU and GU pairs improved from 1.38 to 0.27 kcal/mol. This confirms the need for the additional sequence dependence in the model.

Evolutionary Analyses of Base-Pairing Interactions in DNA and RNA Secondary Structures.

Pairs of nucleotides within functional nucleic acid secondary structures often display evidence of coevolution that is consistent with the maintenance of base-pairing. Here, we introduce a sequence evolution model, MESSI (Modeling the Evolution of Secondary Structure Interactions), that infers coevolution associated with base-paired sites in DNA or RNA sequence alignments. MESSI can estimate coevolution while accounting for an unknown secondary structure. MESSI can also use graphics processing unit parallelism to increase computational speed. We used MESSI to infer coevolution associated with GC, AU (AT in DNA), GU (GT in DNA) pairs in noncoding RNA alignments, and in single-stranded RNA and DNA virus alignments. Estimates of GU pair coevolution were found to be higher at base-paired sites in single-stranded RNA viruses and noncoding RNAs than estimates of GT pair coevolution in single-stranded DNA viruses. A potential biophysical explanation is that GT pairs do not stabilize DNA secondary structures to the same extent that GU pairs do in RNA. Additionally, MESSI estimates the degrees of coevolution at individual base-paired sites in an alignment. These estimates were computed for a SHAPE-MaP-determined HIV-1 NL4-3 RNA secondary structure. We found that estimates of coevolution were more strongly correlated with experimentally determined SHAPE-MaP pairing scores than three nonevolutionary measures of base-pairing covariation. To assist researchers in prioritizing substructures with potential functionality, MESSI automatically ranks substructures by degrees of coevolution at base-paired sites within them. Such a ranking was created for an HIV-1 subtype B alignment, revealing an excess of top-ranking substructures that have been previously identified as having structure-related functional importance, among several uncharacterized top-ranking substructures.

The multiple flavors of GoU pairs in RNA.

Wobble GU pairs (or GoU) occur frequently within double‐stranded RNA helices interspersed within the standard G═C and A─U Watson‐Crick pairs. However, other types of GoU pairs interacting on their Watson‐Crick edges have been observed. The structural and functional roles of such alternative GoU pairs are surprisingly diverse and reflect the various pairings G and U can form by exploiting all the subtleties of their electronic configurations. Here, the structural characteristics of the GoU pairs are updated following the recent crystallographic structures of functional ribosomal complexes and the development in our understanding of ribosomal translation.

Nuclear Magnetic Resonance Reveals That GU Base Pairs Flanking Internal Loops Can Adopt Diverse Structures

RNA thermodynamics play an important role in determining the two- and three-dimensional structures of RNA. Internal loops of the sequence 5'-GMNU/3'-UNMG are relatively unstable thermodynamically. Here, five duplexes with GU-flanked 2 × 2 nucleotide internal loops were structurally investigated to reveal determinants of their instability. The following internal loops were investigated: 5'-GCAU/3'-UACG, 5'-UUCG/3'-GCUU, 5'-GCUU/3'-UUCG, 5'-GUCU/3'-UCUG, and 5'-GCCU/3'-UCCG. Two-dimensional nuclear magnetic resonance spectra indicate the absence of GU wobble base pairing in 5'-GCUU/3'-UUCG, 5'-GUCU/3'-UCUG, and 5'-GCCU/3'-UCCG. The 5'-GCUU/3'-UUCG loop has an unusual conformation of the GU base pairs, in which U's O2 carbonyl forms a bifurcated hydrogen bond with G's amino and imino protons. The internal loop of 5'-GUCU/3'-UCUG displays a shifted configuration in which GC pairs flank a U-U pair and several U's are in fast exchange between positions inside and outside the helix. In contrast, 5'-GCAU/3'-UACG and 5'-UUCG/3'-GCUU both have the expected GU wobble base pairs flanking the internal loop. Evidently, GU base pairs flanking internal loops are more likely to display atypical structures relative to Watson-Crick base pairs flanking internal loops. This appears to be more likely when the G of the GU pair is 5' to the loop. Such unusual structures could serve as recognition elements for biological function and as benchmarks for structure prediction methods.

Comparative genomics and evolution of trans-activating RNAs in Class 2 CRISPR-Cas systems

ABSTRACTTrans-activating CRISPR (tracr) RNA is a distinct RNA species that interacts with the CRISPR (cr) RNA to form the dual guide (g) RNA in type II and subtype V-B CRISPR-Cas systems. The tracrRNA-crRNA interaction is essential for pre-crRNA processing as well as target recognition and cleavage. The tracrRNA consists of an antirepeat, which forms an imperfect hybrid with the repeat in the crRNA, and a distal region containing a Rho-independent terminator. Exhaustive comparative analysis of the sequences and predicted structures of the Class 2 CRISPR guide RNAs shows that all these guide RNAs share distinct structural features, in particular, the nexus stem-loop that separates the repeat-antirepeat hybrid from the distal portion of the tracrRNA and the conserved GU pair at that end of the hybrid. These structural constraints might ensure full exposure of the spacer for target recognition. Reconstruction of tracrRNA evolution for 4 tight bacterial groups demonstrates random drift of repeat-antirepeat complementarity within a window of hybrid stability that is, apparently, maintained by selection. An evolutionary scenario is proposed whereby tracrRNAs evolved on multiple occasions, via rearrangement of a CRISPR array to form the antirepeat in different locations with respect to the array. A functional tracrRNA would form if, in the new location, the antirepeat is flanked by sequences that meet the minimal requirements for a promoter and a Rho-independent terminator. Alternatively, or additionally, the antirepeat sequence could be occasionally ‘reset’ by recombination with a repeat, restoring the functionality of tracrRNAs that drift beyond the required minimal hybrid stability.

Surprising Sequence Effects on GU Closure of Symmetric 2 × 2 Nucleotide RNA Internal Loops.

GU base pairs are important RNA structural motifs and often close loops. Accurate prediction of RNA structures relies upon understanding the interactions determining structure. The thermodynamics of some 2 × 2 nucleotide internal loops closed by GU pairs are not well understood. Here, several self-complementary oligonucleotide sequences expected to form duplexes with 2 × 2 nucleotide internal loops closed by GU pairs were investigated. Surprisingly, nuclear magnetic resonance revealed that many of the sequences exist in equilibrium between hairpin and duplex conformations. This equilibrium is not observed with loops closed by Watson-Crick pairs. To measure the thermodynamics of some 2 × 2 nucleotide internal loops closed by GU pairs, non-self-complementary sequences that preclude formation of hairpins were designed. The measured thermodynamics indicate that some internal loops closed by GU pairs are unusually unstable. This instability accounts for the observed equilibria between duplex and hairpin conformations. Moreover, it suggests that future three-dimensional structures of loops closed by GU pairs may reveal interactions that unexpectedly destabilize folding.

RNA Secondary Structure Determination by NMR.

Dynamic programming methods for predicting RNA secondary structure often use thermodynamics and experimental restraints and/or constraints to limit folding space. Chemical mapping results typically restrain certain nucleotides not to be in AU or GC pairs. Two-dimensional nuclear magnetic resonance (NMR) spectra can reveal the order of AU, GC, and GU pairs in double helixes. This chapter describes a program, NMR-assisted prediction of secondary structure and chemical shifts (NAPSS-CS), that constrains possible secondary structures on the basis of the NMR determined order and 5'-3' direction of AU, GC, and GU pairs in helixes. NAPSS-CS minimally requires input of the order of base pairs as determined from nuclear Overhauser effect spectroscopy (NOESY) of imino protons. The program deduces the 5'-3' direction of the base pairs if certain chemical shifts are also input. Secondary structures predicted by the program provide assignments of input chemical shifts to particular nucleotides in the sequence, thus facilitating an important step for determination of the three dimensional structure by NMR. The method is particularly useful for revealing pseudoknots and an example is provided. The method may also allow determination of secondary structures when a sequence folds into two structures that exchange slowly.

Nuclear Magnetic Resonance-Assisted Prediction of Secondary Structure for RNA: Incorporation of Direction-Dependent Chemical Shift Constraints.

Knowledge of RNA structure is necessary to determine structure–function relationships and to facilitate design of potential therapeutics. RNA secondary structure prediction can be improved by applying constraints from nuclear magnetic resonance (NMR) experiments to a dynamic programming algorithm. Imino proton walks from NOESY spectra reveal double-stranded regions. Chemical shifts of protons in GH1, UH3, and UH5 of GU pairs, UH3, UH5, and AH2 of AU pairs, and GH1 of GC pairs were analyzed to identify constraints for the 5′ to 3′ directionality of base pairs in helices. The 5′ to 3′ directionality constraints were incorporated into an NMR-assisted prediction of secondary structure (NAPSS-CS) program. When it was tested on 18 structures, including nine pseudoknots, the sensitivity and positive predictive value were improved relative to those of three unrestrained programs. The prediction accuracy for the pseudoknots improved the most. The program also facilitates assignment of chemical shifts to individual nucleotides, a necessary step for determining three-dimensional structure.

Testing the nearest neighbor model for canonical RNA base pairs: revision of GU parameters.

Thermodynamic parameters for GU pairs are important for predicting the secondary structures of RNA and for finding genomic sequences that code for structured RNA. Optical melting curves were measured for 29 RNA duplexes with GU pairs to improve nearest neighbor parameters for predicting stabilities of helixes. The updated model eliminates a prior penalty assumed for terminal GU pairs. Six additional duplexes with the 5′GG/3′UU motif were added to the single representation in the previous database. This revises the ΔG°37 for the 5′GG/3′UU motif from an unfavorable 0.5 kcal/mol to a favorable −0.2 kcal/mol. Similarly, the ΔG°37 for the 5′UG/3′GU motif changes from 0.3 to −0.6 kcal/mol. The correlation coefficients between predicted and experimental ΔG°37, ΔH°, and ΔS° for the expanded database are 0.95, 0.89, and 0.87, respectively. The results should improve predictions of RNA secondary structure.

Consecutive Terminal GU Pairs Stabilize RNA Helices

Consecutive GU pairs at the ends of RNA helices provide significant thermodynamic stability between -1.0 and -3.8 kcal/mol at 37 °C, which is equivalent to approximately 2 orders of magnitude in the value of a binding constant. The thermodynamic stabilities of GU pairs depend on the sequence, stacking orientation, and position in the helix. In contrast to GU pairs in the middle of a helix that may be destabilizing, all consecutive terminal GU pairs contribute favorable thermodynamic stability. This work presents measured thermodynamic stabilities for 30 duplexes containing two, three, or four consecutive GU pairs at the ends of RNA helices and a model to predict the thermodynamic stabilities of terminal GU pairs. Imino proton NMR spectra show that the terminal GU nucleotides form hydrogen-bonded pairs. Different orientations of terminal GU pairs can have different conformations with equivalent thermodynamic stabilities. These new data and prediction model will help improve RNA secondary structure prediction, identification of miRNA target sequences with GU pairs, and efforts to understand the fundamental physical forces directing RNA structure and energetics.

The Influence of Anticodon-Codon Interactions and Modified Bases on Codon Usage Bias in Bacteria

Most transfer RNAs (tRNAs) can translate more than one synonymous codon, and most codons can be translated by more than one isoacceptor tRNA. The rates of translation of synonymous codons are dependent on the concentrations of the tRNAs and on the rates of pairing of each anticodon-codon combination. Translational selection causes a significant bias in codon frequencies in highly expressed genes in most bacteria. By comparing codon frequencies in high and low-expression genes, we determine which codons are preferred for each amino acid in a large sample of bacterial genomes. We relate this to the number of copies of each tRNA gene in each genome. In two-codon families, preferred codons have Watson-Crick pairs (GC and AU) between the third codon base and the wobble base of the anticodon rather than GU pairs. This suggests that these combinations are more rapidly recognized by the ribosome. In contrast, in four-codon families, preferred codons do not correspond to Watson-Crick rules. In some cases, a wobble-U tRNA can pair with all four codons. In these cases, A and U codons are preferred over G and C. This indicates that the nonstandard UU combination appears to be translated surprisingly well. Differences in modified bases at the wobble position of the anticodon appear to be responsible for the differences in behavior of tRNAs in two- and four-codon families. We discuss the way changes in the bases in the anticodon influence both the speed and the accuracy of translation. The number of tRNA gene copies and the strength of translational selection correlate with the growth rate of the organism, as we would expect if the primary cause of translational selection in bacteria is the requirement to optimize the speed of protein production.

Cover Picture: Solution Structure of Domain 6 from a Self‐Splicing Group II Intron Ribozyme: A Mg2+ Binding Site is Located Close to the Stacked Branch Adenosine (ChemBioChem 3/2007)

The cover picture shows the NMR solution structure of domain 6 from a group II intron ribozyme derived from the intron ai5γ located in the mitochondrial cytochrome c oxidase I gene of Saccharomyces cerevisiae. As indicated at the top, these large molecular catalytic RNAs cut themselves out of the primary RNA transcript by a two‐step mechanism similar to the spliceosome. The 2′‐OH of a highly conserved adenosine within domain 6 serves as the nucleophile in the first step of splicing. Our solution structure shows for the first time that this branch adenosine resides within the helix, sandwiched between two adjacent GU wobble pairs. These GU pairs thereby form an excellent binding site for a Mg2+ in the major groove. Being an essential cofactor for group II intron catalysis, this novel Mg2+ is in close enough proximity to the catalytic 2′‐OH to exert its influence on catalysis. Further details on the structure and four additional Mg2+ binding sites within domain 6 can be found in the article by R. K. O. Sigel et al. on p. 306 ff.

NMR Structures of (rGCUGAGGCU)2 and (rGCGGAUGCU)2: Probing the Structural Features That Shape the Thermodynamic Stability of GA Pairs,

The NMR structures of [see text] and [see text] are reported. The internal loop, [see text], is about 2 kcal/mol more stable than [see text] at 37 degrees C. The duplexes assemble into similar global folds characterized by the formation of tandem sheared GA pairs. The different stabilities of the loops are accompanied by differences in the local structure of the closing GU pairs. In the [see text] internal loop, the GU pairs form canonical wobble configurations with two hydrogen bonds, whereas in [see text], the GU pairs form a single hydrogen bond involving the amino group, GH22, and the carbonyl group, UO4. This pairing is similar to the GU closing pair of the 690 hairpin loop found in E. coli 16S rRNA. The [see text] and [see text] structures reveal how the subtle interplay between stacking and hydrogen bonding determines sequence dependent conformation and thermodynamic stability. Thus, this work provides structural and thermodynamic benchmarks for theoreticians in the ongoing effort to understand the sequence dependence of RNA physicochemical properties.

Facilitating RNA Structure Prediction with Microarrays

Determining RNA secondary structure is important for understanding structure-function relationships and identifying potential drug targets. This paper reports the use of microarrays with heptamer 2'-O-methyl oligoribonucleotides to probe the secondary structure of an RNA and thereby improve the prediction of that secondary structure. When experimental constraints from hybridization results are added to a free-energy minimization algorithm, the prediction of the secondary structure of Escherichia coli 5S rRNA improves from 27 to 92% of the known canonical base pairs. Optimization of buffer conditions for hybridization and application of 2'-O-methyl-2-thiouridine to enhance binding and improve discrimination between AU and GU pairs are also described. The results suggest that probing RNA with oligonucleotide microarrays can facilitate determination of secondary structure.

Cleavage of dsRNAs hyper-edited by ADARs occurs at preferred editing sites

Long double-stranded RNAs (dsRNAs) may undergo covalent modification (hyper-editing) by adenosine deaminases that act on RNA (ADARs), whereby up to 50–60% of adenosine residues are converted to inosine. Previously, we have described a ribonuclease activity in various cell extracts that specifically targets dsRNAs hyper-edited by ADARs. Such a ribonuclease may play an important role in viral defense, or may alternatively be involved in down-regulation of other RNA duplexes. Cleavage of hyper-edited dsRNA occurs within sequences containing multiple IU pairs but not in duplexes that contain either isosteric GU pairs or Watson–Crick base pairs. Here, we describe experiments aimed at further characterizing cleavage of hyper-edited dsRNA. Using various inosine-containing dsRNAs we show that cleavage occurs preferentially at a site containing both IU and UI pairs, and that inclusion of even a single GU pair inhibits cleavage. We also show that cleavage occurs on both strands within a single dsRNA molecule and requires a 2′-OH group. Strikingly, we show that ADAR1, ADAR2 or dADAR all preferentially generate the preferred cleavage site when hyper-editing a long dsRNA.

Conformational Determinants of Tandem GU Mismatches in RNA: Insights from Molecular Dynamics Simulations and Quantum Mechanical Calculations

Structure and energetic properties of base pair mismatches in duplex RNA have been the focus of numerous investigations due to their role in many important biological functions. Such efforts have contributed to the development of models for secondary structure prediction of RNA, including the nearest-neighbor model. In RNA duplexes containing GU mismatches, 5'-GU-3' tandem mismatches have a different thermodynamic stability than 5'-UG-3' mismatches. In addition, 5'-GU-3' mismatches in some sequence contexts do not follow the nearest-neighbor model for stability. To characterize the underlying atomic forces that determine the structural and thermodynamic properties of GU tandem mismatches, molecular dynamics (MD) simulations were performed on a series of 5'-GU-3' and 5'-UG-3' duplexes in different sequence contexts. Overall, the MD-derived structural models agree well with experimental data, including local deviations in base step helicoidal parameters in the region of the GU mismatches and the model where duplex stability is associated with the pattern of GU hydrogen bonding. Further analysis of the simulations, validated by data from quantum mechanical calculations, suggests that the experimentally observed differences in thermodynamic stability are dominated by GG interstrand followed by GU intrastrand base stacking interactions that dictate the one versus two hydrogen bonding scenarios for the GU pairs. In addition, the inability of 5'-GU-3' mismatches in different sequence contexts to all fit into the nearest-neighbor model is indicated to be associated with interactions of the central four base pairs with the surrounding base pairs. The results emphasize the role of GG and GU stacking interactions on the structure and thermodynamics of GU mismatches in RNA.

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 A) and 30S (3.0 A) 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.

15N NMR of RNA Fragments Containing Specifically Labeled Tandem GA Pairs

15N and 13C NMR have been used to study three specifically labeled RNA fragments that include two tandem GU wobble pairs of different thermal stability and a tetraloop GU pair. The data are also compared to earlier 15N NMR work on an isolated, intrahelical GU pair. The results provide qualitative information on the relative contributions by stacking and hydrogen bonding to chemical shift changes at specific sites in a variety of GU pairs. Chemical shifts for the unpaired guanine amino groups in the isolated and both tandem wobble pairs are all significantly upfield of corresponding paired GC aminos. The guanine amino of the tetraloop is further downfield, as expected for a base paired hydrogen bond donor. Thus, 15N NMR can be diagnostic for the presence of base−base pairing in RNA. Chemical shifts of the guanine N1 in all four GU pairs show significant shielding effects, particularly those of the more stable tandem GU pair. Similar shielding is also seen in the 13C data for the C2 atom of the GU wobble pa...

15N NMR of RNA Fragments Containing Specifically Labeled GU and GC Pairs

15N and 13C NMR have been used to study three specifically labeled RNA fragments that include two tandem GU wobble pairs of different thermal stability and a tetraloop GU pair. The data are also compared to earlier 15N NMR work on an isolated, intrahelical GU pair. The results provide qualitative information on the relative contributions by stacking and hydrogen bonding to chemical shift changes at specific sites in a variety of GU pairs. Chemical shifts for the unpaired guanine amino groups in the isolated and both tandem wobble pairs are all significantly upfield of corresponding paired GC aminos. The guanine amino of the tetraloop is further downfield, as expected for a base paired hydrogen bond donor. Thus, 15N NMR can be diagnostic for the presence of base−base pairing in RNA. Chemical shifts of the guanine N1 in all four GU pairs show significant shielding effects, particularly those of the more stable tandem GU pair. Similar shielding is also seen in the 13C data for the C2 atom of the GU wobble pairs. Our results demonstrate that these adjacent ring atoms (N1 and C2) are in the shielding regions of neighboring bases, and that specific labeling can provide information on base stacking. Thus, where global structure is known, specific labeling can be invaluable as a complement to probe local interactions. It is likely, however, that the most important use of specific labeling will be with complex systems in which global structural information is incomplete.