Solution structure of a DNA double helix incorporating four consecutive non-Watson-Crick base-pairs

A series of DNA 21-mers containing a variety of the 4 x 4 internal loop sequence 5'-CAAG-3'/3'-ACGT-5' were studied using nuclear magnetic resonance (NMR) methodology and distance geometry (DG)/molecular dynamics (MD) approaches. Such oligomers exhibit excellent resolution in the NMR spectra and reveal many unusual NOEs (nuclear Overhauser effect) that allow for the detailed characterization of a DNA hairpin incorporating a track of four different non-Watson-Crick base-pairs in the stem. These include a wobble C.A base-pair, a sheared A.C base-pair, a sheared A.G base-pair, and a wobble G.T base-pair. Significantly different twisting angles were observed between the base-pairs in internal loop that results with excellent intra-strand and inter-strand base stacking within the four consecutive mismatches and the surrounding canonical base-pairs. This explains why it melts at 52 degrees C even though five out of ten base-pairs in the stem adopt non-Watson-Crick pairs. However, the 4 x 4 internal loop still fits into a B-DNA double helix very well without significant change in the backbone torsion angles; only zeta torsion angles between the tandem sheared base-pairs are changed to a great extent from the gauche(-) domain to the trans domain to accommodate the cross-strand base stacking in the internal loop. The observation that several consecutive non-canonical base-pairs can stably co-exist with Watson-Crick base-pairs greatly increases the limited repertoire of irregular DNA folds and reveals the possibility for unusual structural formation in the functionally important genomic regions that have potential to become single-stranded.     

Quadruple intercalated G-6 stack: a possible motif in the fold-back structure of the Drosophila centromeric dodeca-satellite?

The purine-rich strand d(GTACGGGACCGA)(n) of the Drosophila centromeric dodeca-satellite sequence is highly conserved and was found to form stable fold-back structures in which the homopurine 5'-GGGA-3' sequence was determined to play a crucial role. Here, we report the stable formation of the d(GGGA)(2) motif in the stem of a DNA hairpin closed by a single-residue d(ACC) loop. Similar to the zipper-like d(GGA)(2) motif observed in the human centromeric (TGGAA)(n) sequence, the central four guanosine bases in the d(GGGA)(2) motif do not pair, but interdigitate to form an elongated zipper-like quadruple-intercalated G-6 stack bracketed by sheared G.A base-pairs. Comparison between the current d(GGGA)(2) structure and the published crystal d(GAAA)(2) structure implies that the alignment of the unpaired purine bases plays an important role in determining the minor groove width of the purine-rich d(GPuPuA)(2) motif. Similarity between the zipper-like motifs possibly present in the Drosophila centromeric dodeca-satellite sequence and in the human centromeric (TGGAA)(n) sequence led us to propose that these special zipper-like motifs may constitute common cores in organizing eukaryotic centromeres.     

The thermal stability of DNA fragments with tandem mismatches at a d(CXYG).d(CY'X'G) site.

Temperature-Gradient Gel Electrophoresis (TGGE) was employed to determine the thermal stabilities of 28 DNA fragments, 373 bp long, with two adjacent mismatched base pairs, and eight DNAs with Watson-Crick base pairs at the same positions. Heteroduplex DNAs containing two adjacent mismatches were formed by melting and reannealing pairs of homologous 373 bp DNA fragments differing by two adjacent base pairs. Product DNAs were separated based on their thermal stability by parallel and perpendicular TGGE. The polyacrylamide gel contained 3.36 M urea and 19.2 % formamide to lower the DNA melting temperatures. The order of stability was determined in the sequence context d(CXYG).d(CY'X'G) where X.X' and Y.Y" represent the mismatched or Watson-Crick base pairs. The identity of the mismatched bases and their stacking interactions influence DNA stability. Mobility transition melting temperatures (T u) of the DNAs with adjacent mismatches were 1.0-3.6 degrees C (+/-0.2 degree C) lower than the homoduplex DNA with the d(CCAG).d(CTGG) sequence. Two adjacent G.A pairs, d(CGAG).d(CGAG), created a more stable DNA than DNAs with Watson-Crick A.T pairs at the same sites. The d(GA).d(GA) sequence is estimated to be 0.4 (+/-30%) kcal/mol more stable in free energy than d(AA).d(TT) base pairs. This result confirms the unusual stability of the d(GA).d(GA) sequence previously observed in DNA oligomers. All other DNAs with adjacent mismatched base pairs were less stable than Watson-Crick homoduplex DNAs. Their relative stabilities followed an order expected from previous results on single mismatches. Two homoduplex DNAs with identical nearest neighbor sequences but different next-nearest neighbor sequences had a small but reproducible difference in T u value. This result indicates that sequence dependent next neighbor stacking interactions influence DNA stability.     

Sheared-type G(anti).C(syn) base-pair: a unique d(GXC) loop closure motif

Stable DNA loop structures closed by a novel G.C base-pair have been determined for the single-residue d(GXC) loops (X=A, T, G or C) in low-salt solution by high-resolution nuclear magnetic resonance (NMR) techniques. The closing G.C base-pair in these loops is not of the canonical Watson-Crick type, but adopts instead a unique sheared-type (trans Watson-Crick/sugar-edge) pairing, like those occurring in the sheared mismatched G.A or A.C base-pair, to draw the two opposite strands together. The cytidine residue in the closing base-pair is transformed into the rare syn domain to form two H-bonds with the guanine base and to prevent the steric clash between the G 2NH(2) and the C H-5 protons. Besides, the sugar pucker of the syn cytidine is still located in the regular C2'-endo domain, unlike the C3'-endo domain adopted for the pyrimidines of the out-of-alternation left-handed Z-DNA structure. The facile formation of the compact d(GXC) loops closed by a unique sheared-type G(anti).C(syn) base-pair demonstrates the great potential of the single-stranded d(GXC) triplet repeats to fold into stable hairpins.     

Site-resolved energetics in DNA triple helices containing G*TA and T*CG triads

Recognition of specific sites in double-helical DNA by triplex-forming oligonucleotides has been limited until recently to sites containing homopurine-homopyrimidine sequences. G*TA and T*CG triads, in which TA and CG base pairs are specifically recognized by guanine or by thymine, have now extended this recognition code to DNA target sites of mixed base sequences. In the present work, we have obtained a characterization of the stabilities of G*TA and T*CG triads, and of the effects of these triads upon canonical triads, in triple-helical DNA. The three DNA triplexes investigated are formed by the folding of the 31-mers d(GAAXAGGT(5)CCTYTTCT(5)CTTZTCC) with X = G, T, or C, Y = C, A, or G, and Z = C, G, or T. We have measured the exchange rates of imino protons in each triad of the three triplexes using nuclear magnetic resonance spectroscopy. The exchange rates are used to map the local free energy of structural stabilization in each triplex. The results indicate that the stability of Watson-Crick base pairs in the G*TA and T*CG triads is comparable to that of Watson-Crick base pairs in canonical triads. The presence of G*TA and T*CG triads, however, destabilizes neighboring canonical triads, two or three positions removed from the G*TA/T*CG site. Moreover, the long-range destabilizing effects induced by the T*CG triad are larger than those induced by the G*TA triad. These findings reveal the molecular basis for the lower overall stability of G*TA- and T*CG-containing triplexes.     

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.     

Unusual DNA duplex and hairpin motifs

Single-stranded DNA or double-stranded DNA has the potential to adopt a wide variety of unusual duplex and hairpin motifs in the presence (trans) or absence (cis) of ligands. Several principles for the formation of those unusual structures have been established through the observation of a number of recurring structural motifs associated with different sequences. These include: (i) internal loops of consecutive mismatches can occur in a B-DNA duplex when sheared base pairs are adjacent to each other to confer extensive cross- and intra-strand base stacking; (ii) interdigitated (zipper-like) duplex structures form instead when sheared G·A base pairs are separated by one or two pairs of purine·purine mismatches; (iii) stacking is not restricted to base, deoxyribose also exhibits the potential to do so; (iv) canonical G·C or A·T base pairs are flexible enough to exhibit considerable changes from the regular H-bonded conformation. The paired bases become stacked when bracketed by sheared G·A base pairs, or become extruded out and perpendicular to their neighboring bases in the presence of interacting drugs; (v) the purine-rich and pyrimidine-rich loop structures are notably different in nature. The purine-rich loops form compact triloop structures closed by a sheared G·A, A·A, A·C or sheared-like G_anti·C_syn base pair that is stacked by a single residue. On the other hand, the pyrimidine-rich loops with a thymidine in the first position exhibit no base pairing but are characterized by the folding of the thymidine residue into the minor groove to form a compact loop structure. Identification of such diverse duplex or hairpin motifs greatly enlarges the repertoire for unusual DNA structural formation.     

Unique tertiary and neighbor interactions determine conservation patterns of Cis Watson-Crick A/G base-pairs

X-ray, phylogenetic and quantum chemical analysis of molecular interactions and conservation patterns of cis Watson-Crick (W.C.) A/G base-pairs in 16S rRNA, 23S rRNA and other molecules was carried out. In these base-pairs, the A and G nucleotides interact with their W.C. edges with glycosidic bonds oriented cis relative to each other. The base-pair is stabilised by two hydrogen bonds, the C1'-C1' distance is enlarged and the G(N2) amino group is left unpaired. Quantum chemical calculations show that, in the absence of other interactions, the unpaired amino group is substantially non-planar due to its partial sp(3) pyramidalization, while the whole base-pair is internally propeller twisted and very flexible. The unique molecular properties of the cis W.C. A/G base-pairs make them distinct from other base-pairs. They occur mostly at the ends of canonical helices, where they serve as interfaces between the helix and other motifs. The cis W.C. A/G base-pairs play crucial roles in natural RNA structures with salient sequence conservation patterns. The key contribution to conservation is provided by the unpaired G(N2) amino group that is involved in a wide range of tertiary and neighbor contacts in the crystal structures. Many of them are oriented out of the plane of the guanine base and utilize the partial sp(3) pyramidalization of the G(N2). There is a lack of A/G to G/A covariation, which, except for the G(N2) position, would be entirely isosteric. On the contrary, there is a rather frequent occurrence of G/A to G/U covariation, as the G/U wobble base-pair has an unpaired amino group in the same position as the cis W.C. G/A base-pair. The cis W.C. A/G base-pairs are not conserved when there is no tertiary or neighbor interaction. Obtaining the proper picture of the interactions and phylogenetic patterns of the cis W.C. A/G base-pairs requires a detailed analysis of the relation between the molecular structures and the energetics of interactions at a level of single H-bonds and contacts.     

Refined crystal structure of an octanucleotide duplex with G . T mismatched base-pairs

Single crystal X-ray diffraction techniques have been used to determine the structure of the DNA octamer d(G-G-G-G-C-T-C-C) at a resolution of 2.25 A. The asymmetric unit consists of two strands coiled about each other to produce an A-type DNA helix. The double helix contains six G . C Watson-Crick base-pairs and two G . T mismatched base-pairs. The mismatches adopt a "wobble" type structure in which both bases retain their major tautomer forms. The double helix is able to accommodate this G . T pairing with little distortion of the overall helical conformation. Crystals of this octamer melt at a substantially lower temperature than do those of a related octamer also containing two G . T base-pairs. We attribute this destabilization to disruption of the hydration network around the mismatch site combined with changes in intermolecular packing. Full details are given of conformational parameters, base stacking, intermolecular contacts and hydration involving 52 solvent molecules.     

Inosine.adenine base pairs in a B-DNA duplex.

The structure of the synthetic deoxydodecamer d(C-G-C-I-A-A-T-T-A-G-C-G) has been determined by single crystal X-ray diffraction techniques at 2.5A resolution. The refinement converged with a crystallographic residual, R = 0.19 and the location of 64 solvent molecules. The sequence crystallises as a B-DNA helix with 10 Watson-Crick base-pairs (4 A.T. and 6 G.C) and 2 inosine.adenine (I.A) pairs. The present work shows that in the purine.purine base-pairs the adenine adopts syn orientation with respect to the furanose moiety while the inosine is in the trans (anti) orientation. Two hydrogen bonds link the I.A. base-pair, one between N-1(I) and N-7(A), the other between O-6(I) and N-6(A). This bulky purine.purine base-pair is incorporated in the double helix at two positions with little distortion of either local or global conformation. The pairing observed in this study is presented as a model for I.A base-pairs in RNA codon-anticodon interactions and may help explain the thermodynamic stability of inosine containing base-pairs. Conformational parameters and base stacking interactions are presented and where appropriate compared with those of the native compound, d(C-G-C-G-A-A-T-T-C-G-C-G) and with other studies of oligonucleotides containing purine.purine base-pairs.     

1H NMR study of the solution structure of the self-complementary dodecanucleotide d(TGCA)3.

The deoxyoligonucleotide d(TGCA)3 is a candidate for exhibiting unusual conformations. Its 1H NMR spectrum under low salt conditions has been obtained at 400 MHz and assigned using two-dimensional NMR techniques. The sugar puckers and glycosidic torsions have been determined by inspecting the relative intensities of the intranucleotide NOEs and COSY crosspeaks. At low electrolyte concentration (100 mM NaCl) the molecule exists as a right-handed duplex with twelve Watson-Crick base-pairs and deoxyribose moieties assuming the O1'-endo to C1'-exo pucker.     

Consecutive GA pairs stabilize medium-size RNA internal loops

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

A search for base analogs to enhance third-strand binding to 'inverted' target base pairs of triplexes in the pyrimidine/parallel motif.

Eight base analogs were tested as third strand residues in otherwise homopyrimidine strands opposite each of the 'direct' (A.T and G.C) and 'inverted' (T.A and C.G) Watson-Crick base pairs, using UV melting profiles to assess triplex stability. The target duplexes contained 20 A.T base pairs and a central test base pair X.Y, while the third strand contained 20 T residues and a central Z test base. Z included 5-bromo-uracil, 5-propynyluracil, 5-propynylcytosine, 5-methyl-cytosine, 5-bromocytosine, hypoxanthine, 2-amino-purine and 2,6-diaminopurine. Some of the base analogs enhanced third strand binding to the target duplex with one or other 'inverted' central base pair relative to the binding afforded by any of the canonical bases. Other analogs did the same for the duplexes with the 'direct' target pairs. The increasing order of triplex stabilization by these base analogs is: opposite the 'inverted' base pairs, for T.A, A < C < 5-pC < 5-pU < T < 5-BrC < 5-meC < 5-BrU < 2-AP < 2,6-DAP < Hy < G, for C.G, 2-AP < A < Hy < G < 5-pC < 5-BrC < 5-meC < C < 2,6-DAP < T < 5-BrU < 5-pU; opposite the 'direct' base pairs, for A.T, 2-AP < A < 5-meC < C < G < Hy < 2,6-DAP < 5-pU < T = 5-BrU < 5-BrC < 5-pC, for G.C, G < 2,6-DAP < 2-AP < A < Hy < T < 5-BrU < 5-pU < 5-pC < 5-BrC < C < 5-meC.     

An intramolecular quadruplex of (GGA)(4) triplet repeat DNA with a G:G:G:G tetrad and a G(:A):G(:A):G(:A):G heptad, and its dimeric interaction.

The structure of d(GGAGGAGGAGGA) containing four tandem repeats of a GGA triplet sequence has been determined under physiological K(+) conditions. d(GGAGGAGGAGGA) folds into an intramolecular quadruplex composed of a G:G:G:G tetrad and a G(:A):G(:A):G(:A):G heptad. Four G-G segments of d(GGAGGAGGAGGA) are aligned parallel with each other due to six successive turns of the main chain at each of the GGA and GAGG segments. Two quadruplexes form a dimer stabilized through a stacking interaction between the heptads of the two quadruplexes. Comparison of the structure of d(GGAGGAGGAGGA) with the reported structure of d(GGAGGAN) (N=G or T) containing two tandem repeats of the GGA triplet revealed that although the two structures resemble each other to some extent, the extension of the repeats of the GGA triplet leads to distinct structural differences: intramolecular quadruplex for 12-mer versus intermolecular quadruplex for 7-mer; heptad versus hexad in the quadruplex; and three sheared G:A base-pairs versus two sheared G:A base-pairs plus one A:A base-pair per quadruplex. It was also suggested that d(GGAGGAGGAGGA) forms a similar quadruplex under low salt concentration conditions. This is in contrast to the case of d(GGAGGAN) (N=G or T), which forms a duplex under low salt concentration conditions. On the basis of these results, the structure of naturally occurring GGA triplet repeat DNA is discussed.     

Molecular dynamics simulation reveals conformational switching of water-mediated uracil-cytosine base-pairs in an RNA duplex

A 4 ns molecular dynamics simulation of an RNA duplex (r-GGACUUCGGUCC)(2 )in solution with Na+ and Cl- as counterions was performed. The X-ray structure of this duplex includes two water-mediated uracil-cytosine pairs. In contrast to the other base-pairs in the duplex the water-mediated pairs switch between different conformations. One conformation corresponds to the geometry of the water-mediated UC pairs in the duplex X-ray structure with water acting both as hydrogen-bond donor and acceptor. Another conformation is close to that of a water-mediated UC base-pair found in the X-ray structure of the 23 S rRNA sarcin/ricin domain. In this case the oxygen of the water molecule is linked to two-base donor sites. For a very short time also a direct UC base-pair and a further conformation that is similar to the one found in the RNA duplex structure but exhibits an increased H3(U)...N3(C) distance is observed. Water molecules with unusually long residence times are involved in the water-mediated conformations. These results indicate that the dynamic behaviour of the water-mediated UC base-pairs differs from that of the duplex Watson-Crick and non-canonical guanine-uracil pairs with two or three direct hydrogen bonds. The conformational variability and increased flexibility has to be taken into account when considering these base-pairs as RNA building blocks and as recognition motifs.     

The structure of guanosine-thymidine mismatches in B-DNA at 2.5-A resolution

The structure of the deoxyoligomer d(C-G-C-G-A-A-T-T-T-G-C-G) was determined at 2.5-A resolution by single crystal x-ray diffraction techniques. The final R factor is 18% with the location of 71 water molecules. The oligomer crystallizes in a B-DNA-type conformation, with two strands interacting to form a dodecamer duplex. The double helix consists of four A X T and six G X C Watson-Crick base pairs and two G X T mismatches. The G X T pairs adopt a "wobble" structure with the thymine projecting into the major groove and the guanine into the minor groove. The mispairs are accommodated in the normal double helix by small adjustments in the conformation of the sugar phosphate backbone. A comparison with the isomorphous parent compound containing only Watson-Crick base pairs shows that any changes in the structure induced by the presence of G X T mispairs are highly localized. The global conformation of the duplex is conserved. The G X T mismatch has already been studied by x-ray techniques in A and Z helices where similar results were found. The geometry of the mispair is essentially identical in all structures so far examined, irrespective of the DNA conformation. The hydration is also similar with solvent molecules bridging the functional groups of the bases via hydrogen bonds. Hydration may be an important factor in stabilizing G X T mismatches. A characteristic of Watson-Crick paired A X T and G X C bases is the pseudo 2-fold symmetry axis in the plane of the base pairs. The G X T wobble base pair is pronouncedly asymmetric. This asymmetry, coupled with the disposition of functional groups in the major and minor grooves, provides a number of features which may contribute to the recognition of the mismatch by repair enzymes.     

Temperature dependence of the Raman spectrum of DNA. II. Raman signatures of premelting and melting transitions of poly(dA).poly(dT) and comparison with poly(dA-dT).poly(dA-dT).

The temperature dependence of the Raman spectrum of poly(dA).poly(dT) (dA: deoxyadenosine; dT: thymidine), a model for DNA containing consecutive adenine.thymine (A.T) pairs, has been analyzed using a spectrometer of high spectral precision and sensitivity. Three temperature intervals are distinguished: (a) premelting (10 < t < 70 degrees C), in which the native double helix is structurally altered but not dissociated into single strands; (b) melting (70 < t < 80 degrees C), in which the duplex is dissociated into single strands; and (c) postmelting (80 < t degrees C), in which no significant structural change can be detected. The distinctive Raman difference signatures observed between 10 and 70 degrees C and between 70 and 80 degrees C are interpreted in terms of the structural changes specific to premelting and melting transitions, respectively. Premelting alters the low-temperature conformation of the deoxyribose-phosphate backbone and eliminates base hydrogen bonding that is distinct from canonical Watson-Crick hydrogen bonding; these premelting perturbations occur without disruption of base stacking. Conversely, melting eliminates canonical Watson-Crick pairing and base stacking. The results are compared with those reported previously on poly(dA-dT).poly(dA-dT), the DNA structure consisting of alternating A.T and T.A pairs (L. Movileanu, J. M. Benevides, and G. J. Thomas, Jr. Journal of Raman Spectroscopy, 1999, Vol. 30, pp. 637-649). Poly(dA).poly(dT) and poly(dA-dT).poly(dA-dT) exhibit strikingly dissimilar temperature-dependent Raman profiles prior to the onset of melting. However, the two duplexes exhibit very similar melting transitions, including the same Raman indicators of ruptured Watson-Crick pairing, base unstacking and collapse of backbone order. A detailed analysis of the data provides a comprehensive Raman assignment scheme for adenosine and thymidine residues of B-DNA, delineates Raman markers diagnostic of consecutive A.T and alternating A.T/T.A tracts of DNA, and identifies the distinct Raman difference signatures for premelting and melting transitions in the two types of sequences.     

Wobble dG X dT pairing in right-handed DNA: solution conformation of the d(C-G-T-G-A-A-T-T-C-G-C-G) duplex deduced from distance geometry analysis of nuclear Overhauser effect spectra

We report below on features of the three-dimensional structure of the d(C-G-T-G-A-A-T-T-C-G-C-G) self-complementary duplex (designated 12-mer GT) containing symmetrical G X T mismatches in the interior of the helix. The majority of the base and sugar protons in the 12-mer GT duplex were assigned by two-dimensional nuclear Overhauser effect (NOESY) spectra in H2O and D2O solution. A set of 92 short (less than 4.5-A) proton-proton distances defined by lower and upper bounds for one symmetrical half of the 12-mer GT duplex were estimated from NOESY data sets recorded as a function of mixing time. These experimental distances combined with nucleotide bond length parameters were embedded into Cartesian space; several trial structures were refined to minimize bond geometry and van der Waals and chirality error. Confidence in this approach is based on the similarity of the refined structures for the solution conformation of the 12-mer GT duplex. The G and T bases pair through two imino-carbonyl hydrogen bonds, and stacking is maintained between the G X T wobble pair and adjacent Watson-Crick G X C pairs. The experimental distance information is restricted to base and sugar protons, and hence structural features such as base pair overlap, glycosidic torsion angles, and sugar pucker are well-defined by this combination of NMR and distance geometry methods. By contrast, we are unable to define the torsion angles about the bonds C3'-O3'-P-O5'-C5'-C4' in the backbone of the nucleic acid.     

Nuclear magnetic resonance spectroscopy and molecular modeling reveal that different hydrogen bonding patterns are possible for G.U pairs: one hydrogen bond for each G.U pair in r(GGCGUGCC)(2) and two for each G.U pair in r(GAGUGCUC)(2)

G.U pairs occur frequently and have many important biological functions. The stability of symmetric tandem G.U motifs depends both on the adjacent Watson-Crick base pairs, e.g., 5'G > 5'C, and the sequence of the G.U pairs, i.e., 5'-UG-3' > 5'-GU-3', where an underline represents a nucleotide in a G.U pair [Wu, M., McDowell, J. A., and Turner, D. H. (1995) Biochemistry 34, 3204-3211]. In particular, at 37 degrees C, the motif 5'-CGUG-3' is less stable by approximately 3 kcal/mol compared with other symmetric tandem G.U motifs with G-C as adjacent pairs: 5'-GGUC-3', 5'-GUGC-3', and 5'-CUGG-3'. The solution structures of r(GAGUGCUC)(2) and r(GGCGUGCC)(2) duplexes have been determined by NMR and restrained simulated annealing. The global geometry of both duplexes is close to A-form, with some distortions localized in the tandem G.U pair region. The striking discovery is that in r(GGCGUGCC)(2) each G.U pair apparently has only one hydrogen bond instead of the two expected for a canonical wobble pair. In the one-hydrogen-bond model, the distance between GO6 and UH3 is too far to form a hydrogen bond. In addition, the temperature dependence of the imino proton resonances is also consistent with the different number of hydrogen bonds in the G.U pair. To test the NMR models, U or G in various G.U pairs were individually replaced by N3-methyluridine or isoguanosine, respectively, thus eliminating the possibility of hydrogen bonding between GO6 and UH3. The results of thermal melting studies on duplexes with these substitutions support the NMR models.