NMR studies of G:A mismatches in oligodeoxyribonucleotide duplexes modelled after ribozymes.

The structures of two oligodeoxyribonucleotide duplexes, the base sequences of which were modelled after both a hammerhead ribozyme and a small metalloribozyme, were studied by NMR. Both duplexes contain adjacent G:A mismatches; one has a PyGAPu:PyGAPu sequence and the other a PyGAPy:PuGAPu sequence. It is concluded on the basis of many characteristic NOEs that in both duplexes G:A base pairs are formed in the unique 'sheared' form, where an amino proton instead of an imino proton of G is involved in the hydrogen bonding, and G and A bases are arranged 'side by side' instead of 'head to head'. A photo-CIDNP experiment, which gives unique and independent information on the solvent accessibility of nucleotide bases, also supports G:A base pairing rather than a bulged-out structure of G and A residues. This is the first demonstration that not only the PyGAPu:PyGAPu sequence but also the PyGAPy:PuGAPu sequence can form the unique sheared G:A base pairs. Taking the previous studies on G:A mismatches into account, the idea is suggested that a PyGA:GAPu sequence is a minimum and essential element for the formation of the sheared G:A base pairs. The sheared G:A base pairs in the PyGAPu:PyGAPu sequence are suggested to be more stable than those in the PyGAPy:PuGAPu sequence. This is explained rationally by the idea proposed above.     

Unusual DNA conformation at low pH revealed by NMR: parallel-stranded DNA duplex with homo base pairs.

We have investigated the conformational potentials of several DNA oligonucleotides containing sequences related to 5'-CGA in neutral pH and low pH (< 5.0) conditions. One-dimensional proton NMR spectra show that d(CGATCG), d(TCGATCGA), and d(CGATCGATCG) exhibit new sets of resonances at low pH (approximately 3.8-4.4), when compared to those from the neutral pH samples. The low pH form and the neutral pH form are in slow equilibrium. Analyses of the data suggest that these sequences under low pH conditions adopt structures distinct from B-DNA. Two-dimensional nuclear Overhauser effect spectroscopy (2D NOESY) data from the DNA hexamer d(CGATCG) of the neutral and low pH samples were used to analyze their respective structures in solution. An iterative NOE spectral-driven refinement procedure, SPEDREF [Robinson, H., & Wang, A. H.-J. (1992) Biochemistry 31, 3524-3533], was used to show that the neutral pH structure is close to canonical B-DNA. In contrast, analysis of the low pH form using the 2D NOESY data suggests that its structure is consistent with a right-handed parallel-stranded (PS) double helix with symmetrical non-Watson-Crick (C+:C, G:G, A:A, T:T) homo base pairs. Supporting evidence for the PS helix includes the asymmetric inversion-recovery relaxation times associated with the two ends of the helix. The structure is favored by the 5'-CGA sequence in which the cytosines provide the C+:C pairing for the nucleation step and the GpA step is significantly stabilized by the interstrand G-A stacking interactions.(ABSTRACT TRUNCATED AT 250 WORDS)     

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.     

Effects of flanking G . C base pairs on internal Watson-Crick, G . U, and nonbonded base pairs within a short ribonucleic acid duplex

A series of pentaribonucleotides, ApGpXpGpU (where X identical to A, G, C, or U), was synthesized to investigate the effects of flanking G . C pairs on internal Watson-Crick, G . U, and nonbonded base pairs. Sequences ApGpApCpU (Tm = 26 degrees C) and ApGpCpCpU (Tm = 25 degrees C) were each found to form a duplex with non-base-paired internal residues that stacked with the rest of the sequence but were not looped out. ApGpGpCpU also forms a duplex (Tm = 30 degrees C) but with dangling terminal nonbonded adenosines rather than internal nonbonded guanosines. ApGpUpCpU prefers a stacked single-strand conformation. In addition, contribution to duplex stability from an internal A . U or G . C base pair is enhanced by 6 degrees C when flanked by G . C base pairs as compared to A . U base pairs. G . C base pairs flanking an internal G . U base pair were found to be more tolerant to the altered conformation of a G . U pair and result in an increase to stability comparable with that found for an internal A . U base pair.     

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.     

Structure-specific binding of [Co(phen)2(HPIP)] to a DNA duplex containing sheared G:A mismatch base pairs

The binding of a Co(III) complex to the decanucleotide d(CCGAATGAGG)₂ containing two pairs of G:A mismatches was studied by 2D-NMR, UV absorption, and molecular modeling. NMR investigations indicate that racemic [Co(phen)₂(HPIP)]Cl₃ [HPIP = 2-(2-hydroxyphenyl) imidazo [4,5-f][1,10] phenanthroline] binds the decanucleotide by intercalation: the HPIP ligand selectively inserts between the stacked bases from the minor groove at the terminal regions and from the major groove at the sheared region. Further, molecular modeling revealed that the recognition shows strong enantioselectivity: the Λ-isomer preferentially intercalates into the T₆G₇:A₅A₄ region from the DNA major groove, while Δ-isomer favors the terminal C₁C₂:G₁₀G₉ region and intercalates from the minor groove. Detailed energy analysis suggests that the steric interaction, especially the electrostatic effect, is the primary determinants of the recognition event. Melting experiments indicate that binding stabilizes the DNA duplex and increases the melting temperature by 9.5 °C. The intrinsic binding constant of the complex to the mismatched duplex was determined to be 3.5 × 105 M⁻¹.     

Structural equilibria in RNA as revealed by 19F NMR

We have incorporated 5-fluorouridine into several sites within a 19-mer RNA modelled on the translational operator of the MS2 bacteriophage. The 19F NMR spectra demonstrate the different chemical shifts of helical and loop fluorouridines of the hairpin secondary structure. Addition of salt gives rise to a species in which the loop fluorouridine gains the chemical shift of its helical counterparts, due to the formation of the alternative bi-molecular duplex form. This is supported by UV thermal melting behaviour which becomes highly dependent on the RNA concentration. Distinct 19F NMR signals for duplex and hairpin forms allow the duplex-hairpin equilibrium constant to be determined under a range of conditions, enabling thermodynamic characterisation and its salt dependence to be determined. Mg2+ also promotes duplex formation, but more strongly than Na+, such that at 25 degrees C, 10 mM MgCl2 has a comparable duplex-promoting effect to 300 mM NaCl. A similar effect is observed with Sr2+, but not Ca2+ or Ba2+. Additional hairpin species are observed in the presence of Na+ as well as Mg2+, Ca2+, Sr2+ and Ba2+ ions. The overall, ensemble average, hairpin conformation is therefore salt-dependent. Electrostatic considerations are thus involved in the balance between different hairpin conformers as well as the duplex-hairpin equilibrium. The data presented here demonstrate that 19F NMR is a powerful tool for the study of conformational heterogeneity in RNA, which is particularly important for probing the effects of metal ions on RNA structure. The thermodynamic characterisation of duplex-hairpin equilibria will also be valuable in the development of theoretical models of nucleic acid structure.     

Direct identification of NH...N hydrogen bonds in non-canonical base pairs of RNA by NMR spectroscopy.

It is shown that the recently developed quantitative J(NN)HNN-COSY experiment can be used for the direct identification of hydrogen bonds in non-canonical base pairs in RNA. Scalar(2h)J(NN)couplings across NH.N hydrogen bonds are observed in imino hydrogen bonded GA base pairs of the hpGA RNA molecule, which contains a tandem GA mismatch, and in the reverse Hoogsteen AU base pairs of the E-loop of Escherichia coli 5S rRNA. These scalar couplings correlate the imino donor(15)N nucleus of guanine or uridine with the acceptor N1 or N7 nucleus of adenine. The values of the corresponding(2h)J(NN)coupling constants are similar in size to those observed in Watson-Crick base pairs. The reverse Hoogsteen base pairs could be directly detected for the E-loop of E.coli 5S rRNA both in the free form and in a complex with the ribosomal protein L25. This supports the notion that the E-loop is a pre-folded RNA recognition site that is not subject to significant induced conformational changes. Since Watson-Crick GC and AU base pairs are also readily detected the HNN-COSY experiment provides a useful and sensitive tool for the rapid identification of RNA secondary structure elements.     

An alternating sheared AA pair and elements of stability for a single sheared purine-purine pair flanked by sheared GA pairs in RNA

A previous NMR structure of the duplex 5'GGU GGA GGCU/PCCG AAG CCG5' revealed an unusually stable RNA internal loop with three consecutive sheared GA pairs. Here, we report NMR studies of two duplexes, 5'GGU GGA GGCU/PCCA AAG CCG5' (replacing the UG pair with a UA closing pair) and 5'GGU GAA GGCU/PCCG AAG CCG5' (replacing the middle GA pair with an AA pair). An unusually stable loop with three consecutive sheared GA pairs forms in the duplex 5'GGU GGA GGCU/PCCA AAG CCG5'. The structure contrasts with that reported for this loop in the crystal structure of the large ribosomal subunit of Deinococcus radiodurans [Harms, J., Schluenzen, F., Zarivach, R., Bashan, A., Gat, S., Agmon, I., Bartels, H., Franceschi, F., and Yonath, A. (2001) Cell 107, 679-688]. The middle AA pair in the duplex 5'GGU GAA GGCU/PCCG AAG CCG5' rapidly exchanges orientations, resulting in alternative base stacking and pseudosymmetry with exclusively sheared pairs. The U GAA G/G AAG C internal loop is 2.1 kcal/mol less stable than the U GGA G/G AAG C internal loop at 37 degrees C. Structural, energetic, and dynamic consequences upon functional group substitutions within related 3 x 3 and 3 x 6 internal loops are also reported.     

Effects of N2,N2-dimethylguanosine on RNA structure and stability: crystal structure of an RNA duplex with tandem m2 2G:A pairs

Methylation of the exocyclic amino group of guanine is a relatively common modification in rRNA and tRNA. Single methylation (N²-methylguanosine, m²G) is the second most frequently encountered nucleoside analog in Escherichia coli rRNAs. The most prominent case of dual methylation (N²,N²-dimethylguanosine, m²₂G) is found in the majority of eukaryotic tRNAs at base pair m²₂G26:A44. The latter modification eliminates the ability of the N² function to donate in hydrogen bonds and alters its pairing behavior, notably vis-à-vis C. Perhaps a less obvious consequence of the N²,N²-dimethyl modification is its role in controlling the pairing modes between G and A. We have determined the crystal structure of a 13-mer RNA duplex with central tandem m²₂G:A pairs. In the structure both pairs adopt an imino-hydrogen bonded, pseudo-Watson–Crick conformation. Thus, the sheared conformation frequently seen in tandem G:A pairs is avoided due to a potential steric clash between an N²-methyl group and the major groove edge of A. Additionally, for a series of G:A containing self-complementary RNAs we investigated how methylation affects competitive hairpin versus duplex formation based on UV melting profile analysis.     

Analysis of the functional role of a G.A sheared base pair by in vitro genetics

A classical genetic strategy has been combined with an in vitro selection method to search for functional interactions between the two domains of the hairpin ribozyme. G(21) is located within internal loop B; it is proposed to form a sheared base pair with A(43) across loop B and to bind a Mg(2+) ion. Both nucleotides are important for ribozyme function, and G.A sheared base pairs are a very widespread motif in structured RNA. We took advantage of its presence in the hairpin ribozyme to study its functional role. Pseudorevertants, in which the loss of G(21) was compensated by mutations at other positions, were isolated by in vitro selection. The vast majority of G(21) revertants contained substitutions within domain A, pointing to functional communication between specific sites within the two domains of the hairpin ribozyme. The possibility of a direct or redundant contacts is supported by electrophoretic mobility shift studies showing that a complex formed between domain B of the ribozyme and the substrate was disrupted and restored by base substitutions that have analogous effects on catalytic activity. The functional significance of this complex, the role of the nucleotides involved, and the basis for magnesium ion requirement is discussed.     

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

The structure of the synthetic deoxyoctamer d(GGIGCTCC) has been determined by single crystal X-ray diffraction techniques to a resolution of 1.7A. The sequence crystallises in space group P6(1), with unit cell dimensions a = b = 45.07, c = 45.49A. The refinement converged with a crystallographic residual R = 0.14 and the location of 81 solvent molecules. The octamer forms an A-DNA duplex with 6 Watson-Crick (G.C) base pairs and 2 inosine-thymine (I.T) pairs. Refinement of the structure shows it to be essentially isomorphous with that reported for d(GGGGCTCC) with the mispairs adopting a "wobble" conformation. Conformational parameters and base stacking interactions are compared to those for the native duplex d(GGGGCCCC) and other similar sequences. A rationale for the apparent increased crystal packing efficiency and lattice stability of the I.T octamer is given.     

19F NMR of 5-fluorouracil-substituted transfer RNA transcribed in vitro: resonance assignment of fluorouracil-guanine base pairs.

5-Fluorouracil is readily incorporated into active tRNA(Val) transcribed in vitro from a recombinant phagemid containing a synthetic E. coli tRNA(Val) gene. This tRNA has the expected sequence and a secondary and tertiary structure resembling that of native 5-fluorouracil-substituted tRNA(Val), as judged by 19F NMR spectroscopy. To assign resonances in the 19F spectrum, mutant phagemids were constructed having base changes in the tRNA gene. Replacement of fluorouracil in the T-stem with cytosine, converting a FU-G to a C-G base pair, results in the loss of one downfield peak in the 19F NMR spectrum of the mutant tRNA(Val). The spectra of other mutant tRNAs having guanine for adenine substitutions that convert FU-A to FU-G base pairs all have one resonance shifted 4.5 to 5 ppm downfield. These results allow assignment of several 19F resonances and demonstrate that the chemical shift of the 19F signal from base-paired 5-fluorouracil differs considerably between Watson-Crick and wobble geometry.     

Structure and recognition of sheared tandem G x A base pairs associated with human centromere DNA sequence at atomic resolution

G x A mismatched base pairs are frequently found in nucleic acids. Human centromere DNA sequences contain unusual repeating motifs, e.g. , (GAATG)n x (CATTC)n found in the human chromosome. The purine-rich strand of this repeating pentamer sequence forms duplex and hairpin structures with unusual stability. The high stability of these structures is contributed by the "sheared" G x A base pairs which present a novel recognition surface for ligands and proteins. We have solved the crystal structure, by the multiple-wavelength anomalous diffraction (MAD) method of d(CCGAATGAGG) in which the centromere core sequence motif GAATG is embedded. Three crystal forms were refined to near-atomic resolution. The structures reveal the detailed conformation of tandem G x A base pairs whose unique hydrogen-bonding surface has interesting interactions with bases, hydrated magnesium ions, cobalt(III)hexaammine, spermine, and water molecules. The results are relevant in understanding the structure associated with human centromere sequence in particular and G x A base pairs in nucleic acids (including RNA, like ribozyme) in general.     

Copper insertion facilitates water-soluble porphyrin binding to rA.rU and rA.dT base pairs in duplex RNA and RNA.DNA hybrids

The binding of the copper(II) complex of water-soluble meso-tetrakis(N-methylpyridinium-4-yl)porphyrin (TMPyP) to double-helical polynucleotides has been studied by optical absorption, circular dichroism (CD), and resonance Raman spectroscopic methods. The target polymers were RNA and RNA.DNA hybrids consisting of rA.rU, rI.rC, rA.dT, and rI.dC base pairs. Relative to the metal-free H(2)TMPyP [Uno, T., Hamasaki, K., Tanigawa, M., and Shimabayashi, S. (1997) Inorg. Chem. 36, 1676-1683], CuTMPyP binds to poly(rA).poly(dT) and poly(rA).poly(rU) with a greatly increased binding constant. The external self-stacking of the porphyrin on the surface of the polymers was evident from the strong conservative-type induced CD signals. The signal intensity correlated almost linearly with the number of stacking sites on the polymer except for poly(rA).poly(dT), which showed extraordinarily strong CD signals. Thus, the bound porphyrin may impose an ordered architecture on the polymer surface, the stacking being facilitated by the more planar nature of the CuTMPyP than the nonmetal counterpart. Resonance Raman spectra of the stacked CuTMPyP were indistinguishable from those of the intercalated one with positive delta(Cbeta-H) and negative delta(Cm-Py) bending shifts, and hence the stacked porphyrins are suggested to adopt a similar structure to that of intercalated ones. Porphyrin flattening by copper insertion opens a new avenue for medical applications of porphyrins, blocking biological events related to RNA and hybrids in malignant cells.     

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.     

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.     

Thermodynamics of RNA duplexes with tandem mismatches containing a uracil-uracil pair flanked by C.G/G.C or G.C/A.U closing base pairs

The thermodynamics governing the denaturation of RNA duplexes containing 8 bp and a central tandem mismatch or 10 bp were evaluated using UV absorbance melting curves. Each of the eight tandem mismatches that were examined had one U-U pair adjacent to another noncanonical base pair. They were examined in two different RNA duplex environments, one with the tandem mismatch closed by G.C base pairs and the other with G.C and A.U closing base pairs. The free energy increments (Delta Gdegrees(loop)) of the 2 x 2 loops were positive, and showed relatively small differences between the two closing base pair environments. Assuming temperature-independent enthalpy changes for the transitions, (Delta Gdegrees(loop)) for the 2 x 2 loops varied from 0.9 to 1.9 kcal/mol in 1 M Na(+) at 37 degrees C. Most values were within 0.8 kcal/mol of previously estimated values; however, a few sequences differed by 1.2-2.0 kcal/mol. Single strands employed to form the RNA duplexes exhibited small noncooperative absorbance increases with temperature or transitions indicative of partial self-complementary duplexes. One strand formed a partial self-complementary duplex that was more stable than the tandem mismatch duplexes it formed. Transitions of the RNA duplexes were analyzed using equations that included the coupled equilibrium of self-complementary duplex and non-self-complementary duplex denaturation. The average heat capacity change (DeltaC(p)) associated with the transitions of two RNA duplexes was estimated by plotting DeltaH degrees and DeltaS degrees evaluated at different strand concentrations as a function of T(m) and ln T(m), respectively. The average DeltaC(p) was 70 +/- 5 cal K(-)(1) (mol of base pairs)(-)(1). Consideration of this heat capacity change reduced the free energy of formation at 37 degrees C of the 10 bp control RNA duplexes by 0.3-0.6 kcal/mol, which may increase Delta Gdegrees(loop) values by similar amounts.     

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.     

Solution conformation of a bulged adenosine base in an RNA duplex by relaxation matrix refinement

Bulges are common structural motifs in RNA secondary structure and are thought to play important roles in RNA-protein and RNA-drug interactions. Adenosine bases are the most commonly occurring unpaired base in double helical RNA secondary structures. The solution conformation and dynamics of a 25-nucleotide RNA duplex containing an unpaired adenosine, r(GGCAGAGUGCCGC): r(GCGGCACCUGCC) have been studied by NMR spectroscopy and MORASS iterative relaxation matrix structural refinement. The results show that the bulged adenosine residue stacks into the RNA duplex with little perturbation around the bulged region. Most of the bases in the RNA duplex adopt C(3)'-endo conformation, exhibiting the N-type sugar pucker as found in the A form helices. The sugars of the bulged residue and the 5' flanking residue to it are found to exhibit C(2)'-endo conformation. None of the residues are in syn conformation.