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.     

Solution structure of RNA duplexes containing alternating CG base pairs: NMR study of r(CGCGCG)2 and 2'-O-Me(CGCGCG)2 under low salt conditions.

Structures of r(CGCGCG)2 and 2'-O-Me(CGCGCG)2 have been determined by NMR spectroscopy under low salt conditions. All protons and phosphorus nuclei resonances have been assigned. Signals of H5'/5" have been assigned stereospecifically. All 3JH,H and 3JP,H coupling constants have been measured. The structures were determined and refined using an iterative relaxation matrix procedure (IRMA) and the restrained MD simulation. Both duplexes form half-turn, right-handed helices with several conformational features which deviate significantly from a canonical A-RNA structure. Duplexes are characterised as having C3'-endo sugar pucker, very low base-pair rise and high helical twist and inclination angles. Helices are overwound with <10 bp per turn. There is limited inter-strand guanine stacking for CG steps. Within CG steps of both duplexes, the planes of the inter-strand cytosines are not parallel while guanines are almost parallel. For the GC steps this pattern is reversed. The 2'-O-methyl groups are spatially close to the 5'-hydrogens of neighbouring residues from the 3'-side and are directed towards the minor groove of 2'-O-Me(CGCGCG)2 forming a hydrophobic layer. Solution structures of both duplexes are similar; the effect of 2'-O-methylation on the parent RNA structure is small. This suggests that intrinsic properties imposed by alternating CG base pairs govern the overall conformation of both duplexes.     

Crystal structure of an RNA duplex r(G GCGC CC)2 with non-adjacent G*U base pairs

The crystal structure of a self-complementary RNA duplex r(GGGCGCUCC)2with non-adjacent G*U and U*G wobble pairs separated by four Watson-Crick base pairs has been determined to 2.5 A resolution. Crystals belong to the space group R3; a = 33.09 A,alpha = 87.30 degrees with a pseudodyad related duplex in the asymmetric unit. The structure was refined to a final Rworkof 17.5% and Rfreeof 24.0%. The duplexes stack head-to-tail forming infinite columns with virtually no twist at the junction steps. The 3'-terminal cytosine nucleosides are disordered and there are no electron densities, but the 3' penultimate phosphates are observed. As expected, the wobble pairs are displaced with guanine towards the minor groove and uracil towards the major groove. The largest twist angles (37.70 and 40.57 degrees ) are at steps G1*C17/G2*U16 and U7*G11/C8*G10, while the smallest twist angles (28.24 and 27.27 degrees ) are at G2*U16/G3*C15 and C6*G12/U7*G11 and conform to the pseudo-dyad symmetry of the duplex. The molecule has two unequal kinks (17 and 11 degrees ) at the wobble sites and a third kink at the central G5 site which may be attributed to trans alpha (O5'-P), trans gamma (C4'-C5') backbone conformations. The 2'-hydroxyl groups in the minor groove form inter-column hydrogen bonding, either directly or through water molecules.     

Sheared Aanti.Aanti base pairs in a destabilizing 2 x 2 internal loop: the NMR structure of 5'(rGGCAAGCCU)2

The 5'(rGGCAAGCCU)(2) duplex contains tandem A.A pairs. The three-dimensional structure of the 5'(rGGCAAGCCU)(2) duplex was modeled by molecular dynamics and energy minimization with NMR-derived distance and dihedral angle restraints. Although the 5'(rCAAG)(2) loop is thermodynamically destabilizing by 1.1 kcal/mol, the tandem A.A pairs adopt a predominant conformation: a sheared anti-anti (A.A trans Hoogsteen/Sugar-edge) alignment similar to that observed in the crystal structure of the P4-P6 domain of the Tetrahymena thermophila intron [Cate, J. H., Gooding, A. R., Podell, E., Zhou, K., Golden, B. L., Kundrot, C. E., Cech, T. R., and Doudna, J. A. (1996) Science 273, 1678-1685]. The NMR-derived structure of the 5'(rGGCAAGCCU)(2) duplex exhibits cross-strand hydrogen bonds from N3 of A4 to an amino hydrogen of A5 and from the 2' oxygen of the A4 sugar to the other amino hydrogen of A5. An intrastrand hydrogen bond is formed from the 2' OH hydrogen of A4 to O5' of A5. The cross-strand A5 bases are stacked. The Watson-Crick G-C regions are essentially A-form. The sheared anti-anti (A.A trans Hoogsteen/Sugar-edge) alignment provides potential contact sites for tertiary interactions and, therefore, is a possible target site for therapeutics. Thus, thermodynamically destabilizing internal loops can be preorganized for tertiary interactions or ligand binding.     

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.     

An NMR structural study of deaminated base pairs in DNA.

The structurally aberrant base pairs TG, UG and TI may occur in DNA as a consequence of deamination of 5-methylcytosine, cytosine and adenine respectively. Results of NMR spectroscopic studies are reported here for these deaminated base pairs in a model seven base pair long oligonucleotide duplex. We find that in all three cases, the DNA helix is a normal B form and both mispaired bases are intrahelical and hydrogen bonded with one another in a wobble geometry. Similarly, in all three cases, all sugars are found to be normal C2' endo in conformation. Symmetric structural perturbations are observed in the helix twist on the 3' side of the mispaired pyrimidine and on the 5' side of the mispaired purine. In all three cases, the amino group of the G residue on the 3' side of the mispaired pyrimidine shows hindered rotation. Although less thermodynamically stable than helices containing only Watson-Crick base pairs, these helices melt normally from the ends and not from the mispair outwards.     

The intrinsic structure and stability of out-of-alternation base pairs in Z-DNA.

Alternating pyrimidine-purine sequences typically form Z-DNA, with the pyrimidines in the anti and purines in the syn conformations. The observation that dC and dT nucleotides can also adopt the syn conformation (i.e. the nucleotides are out-of-alternation) extends the range of sequences that can convert to this left-handed form of DNA. Here, we study the effects of placing two adjacent d(G*C) base pairs as opposed to a single d(G*C) base pair or two d(A*T) base pairs out-of-alternation by comparing the structure of d(m5CGGCm5CG)2with the previously published structures of d(m5CGGGm5CG)*d(m5CGCCm5CG) and d(m5CGATm5CG)2. A high buckle and loss of stacking interactions are observed as intrinsic properties of the out-of-alternation base pairs regardless of sequence and the context of the dinucleotide. From solution titrations, we find that the destabilizing effect of out-of-alternation d(G*C) base pairs are identical whether these base pairs are adjacent or isolated. We can therefore conclude that it is these intrinsic distortions in the structure of the base pairs and not neighboring effects that account for the inability of out-of-alternation base pairs to adopt the left-handed Z conformation.     

Unnatural base pairs between 2-amino-6-(2-thienyl)purine and the complementary bases.

The unnatural base, 2-amino-6-(2-thienyl)purine (designated as s), instead of 2-amino-6-(N,N-dimethylamino)purine (designated as x), was designed in order to improve the specificity and efficiency of the base pairing with pyridin-2-one (designated as y). DNA fragments containing s were chemically synthesized, and the thermal stability and the enzymatic reactions involving the s-y pairing were examined. Thermal denaturation experiments showed that the DNA duplex (12-mer) containing the s-y pair was more stable than that containing the x-y pair. The incorporation of dyTP was also more advantageous to the s-y pairing than the x-y pairing in single-nucleotide insertion experiments using the Klenow fragment of Escherichia coli DNA polymerase I.     

Dissociation kinetics of 19 base paired oligonucleotide-DNA duplexes containing different single mismatched base pairs.

The dissociation kinetics of 19 base paired oligonucleotide-DNA duplex containing a various single mismatched base pair are studied on dried agarose gels. The kinetics of the dissociation are first order under our experimental conditions. The incorporation of a single mismatched base pair destabilizes the DNA duplexes to some extent, the amount depending on the nature of the mismatched base pair. G-T and G-A mismatches slightly destabilize a duplex, while A-A, T-T, C-T and C-A mismatches significantly destabilize it. The activation energy for the overall dissociation processes for these oligonucleotide-DNA duplexes containing 19 base pairs is 52 +/- 2 Kcal mol-1 as determined from the slope of Arrhenius plot.     

Noncanonical G(syn)-G(anti) base pairs stabilized by sulphate anions in two X-ray structures of the (GUGGUCUGAUGAGGCC) RNA duplex

The structures of two crystal forms of the RNA 16-mer with the sequence GUGGUCUGAUGAGGCC, grown in the presence of a high concentration of sulphate ions, have been determined using synchrotron radiation at 1.4- and 2.0-Å resolution. RNA with this sequence is known as one of the two strands of the noncleavable form of the hammerhead ribozyme. In both crystal structures, two G(syn)–G(anti) noncanonical base pairs are observed in the middle of a 14 base-pair (bp) duplex having 5′-dangling GU residues. Both structures contain sulphate anions interacting with the G–G bp stabilizing G in its syn conformation and bridging the two RNA strands. In both cases the interactions take place in the major groove, although the anions are accommodated within different helix geometries, most pronounced in the changing width of the major groove. In one structure, where a single sulphate spans both G–G pairs, the major groove is closed around the anion, while in the other structure, where each of the two G–G pairs is associated with a separate sulphate, the groove is open. This work provides the first examples of a G–G pair in syn-anti conformation, which minimizes the purine–purine clash in the center of the duplex, while utilizing its residual hydrogen bonding potential in specific interactions with sulphate anions.     

Crystal structure of d(GCGAAAGCT) containing a parallel-stranded duplex with homo base pairs and an anti-parallel duplex with Watson-Crick base pairs

A DNA fragment d(GCGAAAGCT), known to adopt a stable mini-hairpin structure in solution, has been crystallized in the space group I4₁22 with the unit-cell dimensions a = b = 53.4 Å and c = 54.0 Å, and the crystal structure has been determined at 2.5 Å resolution. The four nucleotide residues CGAA of the first half of the oligomer form a parallel duplex with another half through the homo base pairs, C₂:C₂⁺ (singly-protonated between the Watson– Crick sites), G₃:G₃ (between the minor groove sites), A₄:A₄ (between the major groove sites) and A₅:A₅ (between the Watson–Crick sites). The two strands remaining in the half of the parallel duplex are split away in different directions, and they pair in an anti-parallel B-form duplex with the second half extending from a neighboring parallel duplex, so that an infinite column is formed in a head-to-tail fashion along the c-axis. It seems that a hexa-ammine cobalt cation supports such a branched and bent conformation of the oligomer. One end of the parallel duplex is stacked on the corresponding end of the adjacent parallel duplex; between them, the guanine base of the first residue is stacked on the fourth ribose of another duplex.     

Synthesis and crystal structure of an octamer RNA r(guguuuac)/r(guaggcac) with G.G/U.U tandem wobble base pairs: comparison with other tandem G.U pairs

We have determined the crystal structure of the RNA octamer duplex r(guguuuac)/r(guaggcac) with a tandem wobble pair, G·G/U·U (motif III), to compare it with U·G/G·U (motif I) and G·U/U·G (motif II) and to better understand their relative stabilities. The crystal belongs to the rhombohedral space group R3. The hexagonal unit cell dimensions are a = b = 41.92 Å, c = 56.41 Å, and γ = 120°, with one duplex in the asymmetric unit. The structure was solved by the molecular replacement method at 1.9 Å resolution and refined to a final R factor of 19.9% and R_free of 23.3% for 2862 reflections in the resolution range 10.0–1.9 Å with F ≥ 2σ(F). The final model contains 335 atoms for the RNA duplex and 30 water molecules. The A-RNA stacks in the familiar head-to-tail fashion forming a pseudo-continuous helix. The uridine bases of the tandem U·G pairs have slipped towards the minor groove relative to the guanine bases and the uridine O2 atoms form bifurcated hydrogen bonds with the N1 and N2 of guanines. The N2 of guanine and O2 of uridine do not bridge the ‘locked’ water molecule in the minor groove, as in motifs I and II, but are bridged by water molecules in the major groove. A comparison of base stacking stabilities of motif III with motifs I and II confirms the result of thermodynamic studies, motif I > motif III > motif II.     

Hydration of the RNA duplex r(CGCAAAUUUGCG)2 determined by NMR.

The so-called spine of hydration in the minor groove of AnTn tracts in DNA is thought to stabilise the structure, and kinetically bound water detected in the minor groove of such DNA species by NMR has been attributed to a narrow minor groove [Liepinsh, E., Leupin, W. and Otting, G. (1994) Nucleic Acids Res. 22, 2249-2254]. We report here an NMR study of hydration of an RNA dodecamer which has a wide, shallow minor groove. Complete assignments of exchangeable protons, and a large number of non-exchangeable protons in r(CGCAAAUUUGCG)2 have been obtained. In addition, ribose C2'-OH resonances have been detected, which are probably involved in hydrogen bonds. Hydration at different sites in the dodecamer has been measured using ROESY and NOESY experiments at 11.75 and 14.1 T. Base protons in both the major and minor grooves are in contact with water, with effective correlation times for the interaction of approximately 0.5 ns, indicating weak hydration, in contrast to the hydration of adenine C2H in the homologous DNA sequence. NOEs to H1' in the minor groove are consistent with hydration water present that is not observed in the analogous DNA sequence. Hydration kinetics in nucleic acids may be determined by chemical factors such as hydrogen-bonding more than by simple conformational factors such as groove width.     

The effects of unnatural base pairs and mispairs on DNA duplex stability and solvation

In an effort to develop unnatural DNA base pairs we examined six pyridine-based nucleotides, d3MPy, d4MPy, d5MPy, d34DMPy, d35DMPy and d45DMPy. Each bears a pyridyl nucleobase scaffold but they are differentiated by methyl substitution, and were designed to vary both inter- and intra-strand packing within duplex DNA. The effects of the unnatural base pairs on duplex stability demonstrate that the pyridine scaffold may be optimized for stable and selective pairing, and identify one self pair, the pair formed between two d34DMPy nucleotides, which is virtually as stable as a dA:dT base pair in the same sequence context. In addition, we found that the incorporation of either the d34DMPy self pair or a single d34DMPy paired opposite a natural dA significantly increases oligonucleotide hybridization fidelity at other positions within the duplex. Hypersensitization of the duplex to mispairing appears to result from global and interdependent solvation effects mediated by the unnatural nucleotide(s) and the mispair. The results have important implications for our efforts to develop unnatural base pairs and suggest that the unnatural nucleotides might be developed as novel biotechnological tools, diagnostics, or therapeutics for applications where hybridization stringency is important.     

The effect of base sequence on the stability of RNA and DNA single base bulges

Forty-eight RNA duplexes were constructed that contained all common single base bulges at six different locations. The stabilities of the RNAs were determined by temperature gradient gel electrophoresis (TGGE). The relative stability of a single base bulge was dependent on both base identity and the nearest neighbor context. The single base bulges were placed into two categories. A bulged base with no identical neighboring base was defined as a Group I base bulge. Group II-bulged bases had at least one neighboring base identical to it. Group II bulges were generally more stable than Group I bulges in the same nearest neighbor environments. This indicates that position degeneracy of an unpaired base enhances stability. Differences in the mobility transition temperatures between the RNA fragments with bulges and the completely base-paired reference RNAs were related to free energy differences. Simple models for estimating the free energy contribution of single base bulges were evaluated from the free energy difference data. The contribution of a Group I bulge 5'-(XNZ)-3'.5'-(Z'-X')-3' where N is the unpaired base and X.X' and Z.Z' the neighboring base pairs, could be well-represented (+/-0.34 kcal/mol) by the equation, DeltaG((X)(N)()(Z))(.)((Z)(')(-)(X)(')()) = 3.11 + 0. 40DeltaG(s)()((XZ))(.)((Z)(')(X)(')()). DeltaG(s)()((XZ))(. )((Z)(')(X)(')()) is the stacking energy of the closing base pair doublet. By adding a constant term, delta = -0.3 kcal/mol, to the right side of the above equation, free energies of Group II bulges could also be predicted with the same accuracy. The term delta represents the stabilizing effect due to position degeneracy. A similar equation/model was applied to previous data from 32 DNA fragments with single base bulges. It predicted the free energy differences with a similar standard deviation.     

Solid-phase synthesis and high-resolution NMR studies of two synthetic double-helical RNA dodecamers: r(CGCGAAUUCGCG) and r(CGCGUAUACGCG)

Ten-micromole solid-phase RNA synthesis has been successfully performed on an automated nucleic acid synthesizer with coupling efficiencies up to 99%, using the tert-butyldimethylsilyl group to protect the 2'-hydroxyl. The tert-butyldimethylsilyl group was easily removed by tetrabutylammonium fluoride under conditions in which virtually no 2'- to 3'-isomerization was found to occur. By use of this approach, the self-complementary RNA dodecamers r(CGCGAAUUCGCG) and r(CGCGUAUACGCG) were synthesized on an automated nucleic acid synthesizer, purified by TLC, and studied by high-resolution NMR. Imino protons were assigned from one-dimensional nuclear Overhauser effects. The nonexchangeable base, H1', and H2' protons were assigned by the sequential NOESY connectivity method. The NOE data from these two oligomers were analyzed qualitatively and compared to the ideal A- and B-type helix models of Arnott et al. (1972a,b). The internucleotide H6/H8 NOEs to the preceding H1' in r(CGCGUAUACGCG) were found to be sequence-dependent and probably reflect the roll angles between adjacent bases. The internucleotide H6/H8 to H2' NOEs of these oligomers correspond very well to an A-type conformation, but the interstrand adenine H2 NOEs to the following H1' were much stronger than those predicted from the fiber model. These srong interstrand NOEs can be rationalized by base pair slide to favor more interstrand base overlap, as predicted by Callidine and Drew (1984).     

Nitrogen-15 chemical shifts in AT (adenine-thymine) and CG (cytosine-guanine) nucleic acid base pairs

This paper presents ab initio (DFT) calculations of the 15N chemical shifts in AT (Adenine-Thymine) and CG (Cytosine-Guanine) nucleic acid base pairs. Calculations were done on 14 AT and 18 CG base pairs using experimental (X-ray) geometries obtained from several DNA decamers. The calculated chemical shifts are compared with the experimental values in the pure bases and subjected to statistical analysis to explore their sensitivity to the local geometry and pair helix parameters. The results indicate that the 15N chemical shifts, isotropic and principal components are quite sensitive to small changes in the geometry of the pairs, but they do not correlate well with the helix pair parameters. From the statistical analysis, several linear correlations between structural parameters and chemical shifts emerge. These relationships may serve as a foundation to extract information on molecular structure from 15N chemical shift measurements.     

Buoyant density of DNA-Hoechst 33258 (bisbenzimide) complexes in CsCl gradients: Hoechst 33258 binds to single AT base pairs.

Buoyant density of DNA in CsCl gradients with Hoechst 33258 (bisbenzimide) was investigated as a function of guanine plus cytosine content of the DNA (%GC; in mole percent). A formula for calculating %GC from the refractive index (nD) of the isopycnic CsCl/Hoechst 33258 solution over the range of 0-75 %GC was established: %GC = 351762.28 X nD - 123778.66 X nD2 - 249789.47 (the coefficients must not be rounded off). The shape of this curve indicates that under these conditions, in contrast to dilute buffers, Hoechst 33258 binds to single AT base pairs on DNA. Resolution of DNA bands in CsCl/Hoechst 33258 gradients is 1.6 to 2.1 times better than comparative CsCl gradients without the dye. Potential application to %GC determination is discussed.