Configuration of wobble base pairs having pyrimidines as anticodon wobble bases: significance for codon degeneracy

Degeneracy of the genetic code was attributed by Crick to imprecise hydrogen-bonded base-pairing at the wobble position during codon–anticodon pairing. The Crick wobble rules define but do not explain the RNA base pair combinations allowed at this position. We select six pyrimidine bases functioning as anticodon wobble bases (AWBs) to study their H-bonded pairing properties with the four major RNA bases using density functional theory at the B3LYP/6-31G(d,p) level. This is done to assess the extent to which the configuration of a solitary RNA wobble base pair may in itself determine specificity and degeneracy of the genetic code by allowing or disallowing the given base pair during codon–anticodon pairing. Calculated values of select configuration markers for the base pairs screen well between allowed and disallowed base pairs for most cases examined here, where the base pair width emerges as an important factor. A few allowed wobble pairs invoke the involvement of RNA nucleoside conformation, as well as involvement of the exocyclic substituent in H-bonding. This study, however, cannot explain the disallowed status of the Ura?Gua wobble pair on the basis of configuration alone. Explanation of the allowed status of the V?Ura pair requires further study on the mediatory role of water molecules. Apart from these two cases, these computational results are sufficient, on the basis of base pair configuration alone, to account for the specificity and degeneracy of the genetic code for all known cases of codon–anticodon pairing which involve the pyrimidine AWBs studied here.     

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.     

Structure of SL4 RNA from the HIV-1 packaging signal.

The NMR-based structure is described for an RNA model of stem-loop 4 (SL4) from the HIV-1 major packaging domain. The GAGA tetraloop adopts a conformation similar to the classic GNRA form, although there are differences in the details. The type II tandem G.U pairs have a combination of wobble and bifurcated hydrogen bonds where the uracil 2-carbonyl oxygen is hydrogen-bonded to both G,H1 and G,H2. There is the likelihood of a Na(+) ion coordinated to the four carbonyl oxygens in the major groove for these G.U pairs and perhaps to the N7 lone pairs of the G bases as well. A continuous stack of five bases extends over nearly the whole length of the stem to the base of the loop in the RNA 16mer: C15/U14/G13/G5/C6. There is no evidence for a terminal G.A pair; instead, G1 appears quite unrestrained, and A16 stacks on both C15 and G2. Residues G2 through G5 exhibit broadened resonances, especially G3 and U4, suggesting enhanced mobility for the 5'-side of the stem. The structure shows G2/G3/U4 stacking along the same strand, nearly isolated from interaction with the other bases. This is probably an important factor in the signal broadening and apparent mobility of these residues and the low stability of the 16mer hairpin against thermal denaturation.     

Structural insights of non-canonical U?U pair and Hoogsteen interaction probed with Se atom

Unlike DNA, in addition to the 2′-OH group, uracil nucleobase and its modifications play essential roles in structure and function diversities of non-coding RNAs. Non-canonical U•U base pair is ubiquitous in non-coding RNAs, which are highly diversified. However, it is not completely clear how uracil plays the diversifing roles. To investigate and compare the uracil in U-A and U•U base pairs, we have decided to probe them with a selenium atom by synthesizing the novel 4-Se-uridine (^SeU) phosphoramidite and Se-nucleobase-modified RNAs (^SeU-RNAs), where the exo-4-oxygen of uracil is replaced by selenium. Our crystal structure studies of U-A and U•U pairs reveal that the native and Se-derivatized structures are virtually identical, and both U-A and U•U pairs can accommodate large Se atoms. Our thermostability and crystal structure studies indicate that the weakened H-bonding in U-A pair may be compensated by the base stacking, and that the stacking of the trans-Hoogsteen U•U pairs may stabilize RNA duplex and its junction. Our result confirms that the hydrogen bond (O4^…H-C5) of the Hoogsteen pair is weak. Using the Se atom probe, our Se-functionalization studies reveal more insights into the U•U interaction and U-participation in structure and function diversification of nucleic acids.     

New information content in RNA base pairing deduced from quantitative analysis of high-resolution structures

Non-canonical base pairs play important roles in organizing the complex three-dimensional folding of RNA. Here, we outline methodology developed both to analyze the spatial patterns of interacting base pairs in known RNA structures and to reconstruct models from the collective experimental information. We focus attention on the structural context and deformability of the seven pairing patterns found in greatest abundance in the helical segments in a set of well-resolved crystal structures, including (i–ii) the canonical A·U and G·C Watson–Crick base pairs, (iii) the G·U wobble pair, (iv) the sheared G·A pair, (v) the A·U Hoogsteen pair, (vi) the U·U wobble pair, and (vii) the G·A Watson–Crick-like pair. The non-canonical pairs stand out from the canonical associations in terms of apparent deformability, spanning a broader range of conformational states as measured by the six rigid-body parameters used to describe the spatial arrangements of the interacting bases, the root-mean-square deviations of the base-pair atoms, and the fluctuations in hydrogen-bonding geometry. The deformabilties, the modes of base-pair deformation, and the preferred sites of occurrence depend on sequence. We also characterize the positioning and overlap of the base pairs with respect to the base pairs that stack immediately above and below them in double-helical fragments. We incorporate the observed positions of the bases, base pairs, and intervening phosphorus atoms in models to predict the effects of the non-canonical interactions on overall helical structure.     

Thermodynamics of RNA-RNA duplexes with 2- or 4-thiouridines: implications for antisense design and targeting a group I intron

Antisense compounds are designed to optimize selective hybridization of an exogenous oligonucleotide to a cellular target. Typically, Watson-Crick base pairing between the antisense compound and target provides the key recognition element. Uridine (U), however, not only stably base pairs with adenosine (A) but also with guanosine (G), thus reducing specificity. Studies of duplex formation by oligonucleotides with either an internal or a terminal 2- or 4-thiouridine (s(2)U or s(4)U) show that s(2)U can increase the stability of base pairing with A more than with G, while s(4)U can increase the stability of base pairing with G more than with A. The latter may be useful when binding can be enhanced by tertiary interactions with a s(4)U-G pair. To test the effects of s(2)U and s(4)U substitutions on tertiary interactions, binding to a group I intron ribozyme from mouse-derived Pneumocystis carinii was measured for the hexamers, r(AUGACU), r(AUGACs(2)U), and r(AUGACs(4)U), which mimic the 3' end of the 5' exon. The results suggest that at least one of the carbonyl groups of the 3' terminal U of r(AUGACU) is involved in tertiary interactions with the catalytic core of the ribozyme and/or thio groups change the orientation of a terminal U-G base pair. Thus thio substitutions may affect tertiary interactions. Studies of trans-splicing of 5' exon mimics to a truncated rRNA precursor, however, indicate that thio substitutions have negligible effects on overall reactivity. Therefore, modified bases can enhance the specificity of base pairing while retaining other activities and, thus, increase the specificity of antisense compounds targeting cellular RNA.     

Wobble base-pairing in codon-anticodon interactions: a theoretical modelling study

The Crick wobble hypothesis attributes the phenomenon of codon degeneracy to a certain impreciseness of pairing between the third base of the codon and the first base of the anticodon. This theoretical study investigates the pairing properties of some wobble bases, including both, observed and unobserved pairs. Some wobble base-pairs are predicted to follow the Watson-Crick pairs in configuration and pairing facility, while others deviate from this norm. The observed U:V pair is unique in that a pairing configuration may be suggested for it wherein the hydrogen-bonding involves the exocyclic 5-carboxymethoxy group of V. By comparing the theoretical data on the configurations of these pairs with the evidence for their existence/non-existence in nature, some guidelines emerge for differentiating between observed and unobserved base pairs on the basis of the pairing configuration.     

Stacking of Crick Wobble pair and Watson-Crick pair: stability rules of G-U pairs at ends of helical stems in tRNAs and the relation to codon-anticodon Wobble interaction.

The occurrence of the noncomplementary G-U base pair at the end of a helix is found to be governed by stacking interactions. As a rule, a G-U pair with G on the 5'-side of a Watson-Crick base pair exhibits strikingly greater stacking overlap with the Watson-Crick base pair than a G-U pair on the 3'-side of a Watson-Crick base pair. The former arrangement is expected to be more stable and indeed is observed 29 times out of 32 in the known transfer RNA molecules. In accordance with this rule, the major wobble base pairs G-U or I-U in codon-anticodon interactions have G or I on the 5'-side of the anticodon. Similarly, in initiator tRNAs, this rule is obeyed where now the G is the first letter of the codon (5'-side). In the situation where U is in the wobble position of the anticodon, it is usually substituted at C(5) andmay also have a 2-thio group and it can read one to four codons depending on its modifications. A G at the wobble position of the anticodon can recognize the two codons ending with U or C and modification of G (unless it is I) does not change its reading properties.     

G.T wobble base-pairing in Z-DNA at 1.0 A atomic resolution: the crystal structure of d(CGCGTG).

The DNA oligomer d(CGCGTG) crystallizes as a Z-DNA double helix containing two guanine-thymine base pair mismatches of the wobble type. The crystal diffracts to 1 A resolution and the structure has been solved and refined. At this resolution, a large amount of information is revealed about the organization of the water molecules in the lattice generally and more specifically around the wobble base pairs. By comparing this structure with the analogous high resolution structure of d(CGCGCG) we can visualize the structural changes as well as the reorganization of the solvent molecules associated with wobble base pairing. There is only a small distortion of the Z-DNA backbone resulting from introduction of the GT mismatched base pairs. The water molecules cluster around the wobble base pair taking up all of the hydrogen bonding capabilities of the bases due to wobble pairing. These bridging water molecules serve to stabilize the base-base interaction and, thus, may be generally important for base mispairing either in DNA or in RNA molecules.     

Demonstration of G . U wobble base pairs by Raman and IR spectroscopy

When guanine and uracil form hydrogen bonds in the pairing scheme first proposed by Crick one would expect that poly(A,G) will form an unperturbed double helix with poly U at room temperature in a dilute electrolyte solution (0.1 M NaCl). We have demonstrated by Raman- and IR-spectroscopy that the secondary structure of poly(A.G) · poly U is very similar to the structure of poly A · poly U; only the thermal stability of the double helix seems slightly lower than the stability of poly A · poly U, whereas the average helix length is unaffected by the dispersed G · U base pairs. From our input ratio of guanine and adenine we estimate that about every fourth base pair is a wobble pair.     

Hydrogen and hydration of DNA and RNA oligonucleotides

In this paper, hydrogen bonding interaction and hydration in crystal structures of both DNA and RNA oligonucleotides are discussed. Their roles in the formation and stabilization of oligonucleotides have been covered. Details of the Watson-Crick base pairs G.C and A.U in DNA and RNA are illustrated. The geometry of the wobble (mismatched) G.U base pairs and the cis and almost trans conformations of the mismatched U.U base pairs in RNA is described. The difference in hydration of the Watson-Crick base pairs G.C, A.U and the wobble G.U in different sequences of codon-anticodon interaction in double helical molecules are indicative of the effect of hydration. The hydration patterns of the phosphate, the 2'-hydroxyl groups, the water bridges linking the phosphate group, N7 (purine) and N4 of Cs or O4 of Us in the major groove, the water bridges between the 2'-hydroxyl group and N3 (purine) and O2 (pyrimidine) in the minor groove are discussed.     

The 5S rRNA loop E: chemical probing and phylogenetic data versus crystal structure.

A significant fraction of the bases in a folded, structured RNA molecule participate in noncanonical base pairing interactions, often in the context of internal loops or multi-helix junction loops. The appearance of each new high-resolution RNA structure provides welcome data to guide efforts to understand and predict RNA 3D structure, especially when the RNA in question is a functionally conserved molecule. The recent publication of the crystal structure of the "Loop E" region of bacterial 5S ribosomal RNA is such an event [Correll CC, Freeborn B, Moore PB, Steitz TA, 1997, Cell 91:705-712]. In addition to providing more examples of already established noncanonical base pairs, such as purine-purine sheared pairings, trans-Hoogsteen UA, and GU wobble pairs, the structure provides the first high-resolution views of two new purine-purine pairings and a new GU pairing. The goal of the present analysis is to expand the capabilities of both chemical probing and phylogenetic analysis to predict with greater accuracy the structures of RNA molecules. First, in light of existing chemical probing data, we investigate what lessons could be learned regarding the interpretation of this widely used method of RNA structure probing. Then we analyze the 3D structure with reference to molecular phylogeny data (assuming conservation of function) to discover what alternative base pairings are geometrically compatible with the structure. The comparisons between previous modeling efforts and crystal structures show that the intricate involvements of ions and water molecules in the maintenance of non-Watson-Crick pairs render the process of correctly identifying the interacting sites in such pairs treacherous, except in cases of trans-Hoogsteen A/U or sheared A/G pairs for the adenine N1 site. The phylogenetic analysis identifies A/A, A/C, A/U and C/A, C/C, and C/U pairings isosteric with sheared A/G, as well as A/A and A/C pairings isosteric with both G/U and G/G bifurcated pairings. Thus, each non-Watson-Crick pair could be characterized by a phylogenetic signature of variations between isosteric-like pairings. In addition to the conservative changes, which form a dictionary of pairings isosterically compatible with those observed in the crystal structure, concerted changes involving several base pairs also occur. The latter covariations may indicate transitions between related but distinctive motifs within the loop E of 5S ribosomal RNA.     

An innate twist between Crick’s wobble and Watson-Crick base pairs

Non-Watson-Crick pairs like the G·U wobble are frequent in RNA duplexes. Their geometric dissimilarity (nonisostericity) with the Watson-Crick base pairs and among themselves imparts structural variations decisive for biological functions. Through a novel circular representation of base pairs, a simple and general metric scheme for quantification of base-pair nonisostericity, in terms of residual twist and radial difference that can also envisage its mechanistic effect, is proposed. The scheme is exemplified by G·U and U·G wobble pairs, and their predicable local effects on helical twist angle are validated by MD simulations. New insights into a possible rationale for contextual occurrence of G·U and other non-WC pairs, as well as the influence of a G·U pair on other non-Watson-Crick pair neighborhood and RNA-protein interactions are obtained from analysis of crystal structure data. A few instances of RNA-protein interactions along the major groove are documented in addition to the well-recognized interaction of the G·U pair along the minor groove. The nonisostericity-mediated influence of wobble pairs for facilitating helical packing through long-range interactions in ribosomal RNAs is also reviewed.     

Thermodynamics and kinetics for base-pair opening in the P1 duplex of the Tetrahymena group I ribozyme

The thermodynamics and kinetics for base-pair opening of the P1 duplex of the Tetrahymena group I ribozyme were studied by NMR hydrogen exchange experiments. The apparent equilibrium constants for base pair opening were measured for most of the imino protons in the P1 duplex using the base catalysts NH3, HPO4(2-) or TRIS. These equilibrium constants were also measured for several modified P1 duplexes, and the C-2.G23 base pair was the most stable base pair in all the duplexes. The conserved U-1*G22 base pair is required for activity of the ribozyme and the data here show that this wobble base pair destabilizes neighboring base pairs on only one side of the wobble. A 2'-OMe modification on the U-3 residue stabilized its own base pair but had little effect on the neighboring base pairs. Three base pairs, U-1*G22, C-2*G23 and A2*U21 showed unusual equilibrium constants for opening and possible implications of the opening thermodynamics of these base pairs on the undocking rates of the P1 helix with catalytic core are discussed.     

Synthesis of 2'-O-methyl-RNAs incorporating a 3-deazaguanine, and UV melting and computational studies on its hybridization properties

2'-O-methyl-RNAs incorporating 3-deazaguanine (c3G) were synthesized by use of N,N-diphenylcarbamoyl and N,N-dimethylaminomethylene as its base protecting groups to suppress sheared-type 5'-GA-3'/5'-GA-3' tandem mismatched base pairing which requires the N3 atom. These modified RNAs hybridized more weakly with the complementary and single mismatch-containing RNAs than the unmodified RNAs. The T(m) experiments were performed to clarify the effects of replacement of the fifth G with c(3)G on stabilization of 2'-O-methyl-(5'-CGGCGAGGAG-3')/5'-CUCCGAGCCG-3' and 2'-O-methyl-(5'-CGGGGACGAG-3')/5'-CUCGGACCCG-3'duplexes, which form sheared-type and face-to-face type 5'-GA-3'/5'-GA-3' tandem mismatched base pairs, respectively. Consequently, this replacement led to more pronounced destabilization of the former duplex that needs the N3 atom for the sheared-type base pair than the latter that does not need it for the face-to-face type base pair. A similar tendency was observed for 2'-O-methyl-RNA/DNA duplexes. These results suggest that the N3 atom of G plays an important role in stabilization of the canonical G/C base pair as well as the base discrimination and its loss suppressed formation of the undesired sheared-type mismatched base pair. Computational studies based on ab initio calculations suggest that the weaker hydrogen bonding ability and larger dipole moment of c3G can be the origin of the lower T(m).     

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.     

Structural basis for the guanosine requirement of the hairpin ribozyme

To form a catalytically active complex, the essential nucleotides of the hairpin ribozyme, embedded within the internal loops of the two domains, must interact with one another. Little is known about the nature of these essential interdomain interactions. In the work presented here, we have used recent topographical constraints and other biochemical data in conjunction with molecular modeling (constraint-satisfaction program MC-SYM) to generate testable models of interdomain interactions. Visual analysis of the generated models has revealed a potential interdomain base pair between the conserved guanosine immediately downstream of the reactive phosphodiester (G(+1)) and C(25) within the large domain. We have tested this former model through activity assays, using all 16 combinations of bases at positions +1 and 25. When the standard ribozyme was used, catalytic activity was severely suppressed with substrates containing U(+1), C(+1), or A(+1). Similarly, mutations of the putative pairing partner (C(25) to A(25) or G(25)) reduce activity by several orders of magnitude. The U(25) substitution retains a significant level of activity, consistent with the possible formation of a G.U wobble pair. Strikingly, when combinations of Watson-Crick (or wobble) base pairs were introduced in these positions, catalytic activity was restored, strongly suggesting the existence of the proposed interaction. These results provide a structural basis for the guanosine requirement of this ribozyme and indicate that the hairpin ribozyme can now be engineered to cleave a wider range of RNA sequences.     

Conformational features of the four successive non-Watson-Crick base pairs in RNA duplex.

A tridecaribonucleotide, r(UGAGCUUCGGCUC) doesn't form hairpin or interior loop and forms a double helix of 12 base pairs including the four successive nonstandard base pairs, U.G-U.C-C.U-G.U, in the crystal. Non-Watson-Crick base pairs, G.U and U.C are nicely incorporated in RNA duplex maintaining the regular A-form backbone. There exist the good overlapping between base pairings, U.G and U.C, so as to stabilize the nonstandard base pair track. Hydrogen bond networks involving water molecules in the major and minor grooves to stabilize this mismatch base pairing array, are observed and its conformational features are described.