Solution structure of the d(T-C-G-A) duplex at acidic pH. A parallel-stranded helix containing C+ .C, G.G and A.A pairs

The solution structure of the d(T-C-G-A) sequence at acidic pH has been determined by a combination of NMR and molecular dynamics calculations including NOE intensity based refinements. This sequence forms a right-handed parallel-stranded duplex with C+ .C (three hydrogen bonds along Watson-Crick edge), G.G (two symmetry related N2-H.. N3 hydrogen bonds) and A.A (two symmetry related N6-H..N7 hydrogen bonds) homo base-pair formation at acidic pH. The duplex is stabilized by intra-strand base stacking at the C2-G3 step and cross-strand base stacking at the G3-A4 step. The thymine residues on partner strands are directed towards each other and are positioned over the C+ .C base-pair. All four residues adopt anti glycosidic torsion angles and C2'-endo type sugar conformations in the parallel-stranded d(T-C-G-A) duplex which exhibits large changes in twist angles between adjacent steps along the duplex. This study rules out previously proposed models for the structure of the d(T-C-G-A) duplex at acidic pH and supports earlier structural contributions, which established that d(C-G) and d(C-G-A) containing sequences at acidic pH pair through parallel-stranded alignment. We have also monitored hydration patterns in the symmetry related grooves of the parallel-stranded d(T-C-G-A) duplex.     

Structure of a 16-mer RNA duplex r(GCAGACUUAAAUCUGC)2 with wobble C.A+ mismatches

The crystal structure of a 16-mer, the longest known RNA duplex, has been determined at 2.5 A resolution. The hexadecamer r(GCAGACUUAAAUCUGC) contains isolated C.A/A.C mismatches with two hydrogen bonds. The two hydrogen bonds in the mismatches suggests that N1 of A is protonated even though the crystallization was done at neutral pH. Therefore, the C.A mismatch is a C.A+ wobble similar to the G.U wobble. The two C.A+ pairs are isolated by four Watson-Crick pairs and flanked by five Watson-Crick base-pairs on either sides. Kinks/bends of 20 degrees are observed at the wobble sites. The Watson-Crick base-pair A5.U26 on the 5'-side of the first C6.A27(+) wobble has a twist angle of 27 degrees compared to the 3'-side U7.A28 pair of 36 degrees. The twist angles are reversed (37 degrees and 26 degrees) in the second A11(+).C22 wobble because of the approximate dyad in the molecule, the flanking base-pair sequences are A.U pairs. The wobbles expand the major groove to 7.1 A/7.3 A. The duplexes form helical columns and are tightly packed around the 31-screw axis. The minor grooves of adjacent columns in juxtaposition interact through the O2' atoms and the anionic phosphate oxygen atoms.     

Structural characterization of an N-acetyl-2-aminofluorene (AAF) modified DNA oligomer by NMR, energy minimization, and molecular dynamics

An N-acetyl-2-aminofluorene (AAF) modified deoxyoligonucleotide duplex, d(C1-C2-A3-C4-[AAF-G5]-C6-A7-C8-C9).d(G10-G11-T12-G13-C14-++ +G15-T16-G17-G18), was studied by one- and two-dimensional NMR spectroscopy. Eight of the nine complementary nucleotides form Watson-Crick base pairs, as shown by NOEs between the guanine imino proton and cytosine amino protons for G.C base pairs or by an NOE between the thymine imino proton and adenine H2 proton for A.T base pairs. The AAF-G5 and C14 bases show no evidence of complementary hydrogen bond formation to each other. The AAF-G5 base adopts a syn conformation, as indicated by NOEs between the G5 imino proton and the A3-H3' and A3-H2'/H2" protons and by NOEs between the fluorene-H1 proton of AAF and the G5-H1' or C6-H1' proton. The NOEs from the C4-H6 proton to C4 sugar protons are weak, and thus the glycosidic torsion angle in this nucleotide is not well defined by these NMR data. The remaining bases are in the anti conformation, as depicted by the relative magnitude of the H8/H6 to H2' NOEs when compared to the H8/H6 to H1' NOEs. The three base pairs on each end of the duplex exhibit NOEs characteristic of right-handed B-form DNA. Distance restraints obtained from NOESY data recorded at 32 degrees C using a 100-ms mixing time were used in conformational searches by molecular mechanics energy minimization studies. The final, unrestrained, minimum-energy conformation was then used as input for an unrestrained molecular dynamics simulation. Chemical exchange cross peaks are observed, and thus the AAF-9-mer exists in more than a single conformation on the NMR time scale. The NMR data, however, indicate the presence of a predominant conformation (> or = 70%). The structure of the predominant conformation of the AAF-9-mer shows stacking of the fluorene moiety on an adjacent base pair, exhibiting features of the base-displacement [Grunberger, D., Nelson, J. H., et al. (1970) Proc. Natl. Acad. Sci. U.S.A. 66, 488-494] and insertion-denaturation models [Fuchs, R.P.P., & Daune, M. (1971) FEBS Lett. 14, 206-208], while the distal ring of the fluorene moiety protrudes into the minor groove.     

An intercalated and thermally stable FAPY adduct of aflatoxin B1 in a DNA duplex: structural refinement from 1H NMR

The structure of a formamidopyrimidine (FAPY) adduct arising from imidazole ring opening of the initially formed trans-8, 9-dihydro-8-(N7-guanyl)-9-hydroxyaflatoxin B1 adduct under basic conditions and positioned in the 5'-d(CTATFAPYGATTCA)-3'*5'-d(TGAATCATAG)-3' oligodeoxynucleotide was determined. The FAPY adduct may be a major progenitor of aflatoxin B1-induced mutations in DNA. The freshly prepared sample showed biphasic melting, with transitions at 28 and 56 degreesC. NMR initially showed multiple subspectra. Over a period of several days at 4 degreesC, the sample converted to a single species with a Tm of 56 degreesC, 15 degrees C greater than the unmodified duplex. The deoxyribose was in the beta configuration about the anomeric carbon, evidenced by NOEs between FAPYG5 H3', H2', H2", and H1'. FAPY formation resulted in the loss of the guanine H8 proton, and the introduction of the formyl proton, which showed NOEs to FAPYG5 H1' and A6 N6Ha. A total of 31 NOEs from AFB1 to DNA protons were observed, mostly to the 5'-neighboring base, T4 in the modified strand. Sequential NOEs were interrupted between T4 and FAPYG5 in the modified strand, between C16 and A17 in the complementary strand, and between T4 N3H and FAPYG5 N1H. An NOE between FAPYG5 N1H and C16 N4H showed intact hydrogen bonding at FAPYG5*C16. Upfield chemical shifts were observed for T4 H6 and A17 H8. Molecular dynamics calculations converged with pairwise rmsd differences of <0.9 A. The sixth root residual was 8.7 x 10(-2). The AFB1 moiety intercalated from the major groove between FAPYG5 and T4*A17, and stacked with T4 and FAPYG5 and partially stacked with A17. The base step between T4*A17 and FAPYG5*C16 was increased from 3.4 to 7 A. The duplex unwound by about 15 degrees. The FAPY formyl group was positioned to form a hydrogen bond with A6 N6Ha. Strong stacking involving the AFB1 moiety, and this hydrogen bond explains the thermal stabilization of four base pairs by this adduct, and may be a significant factor in its processing.     

Structures of oligonucleotides containing 5-(methoxymethyl)-2'-deoxyuridine determined by NMR spectroscopy

Base pairing of 5-(methoxymethyl)-2'-deoxyuridine (MMdU) opposite either adenine or guanine in a seven base pair oligonucleotide duplex has been studied by NMR spectroscopy. When paired with A, we observe that the MMdU.A base pair adopts Watson-Crick geometry. The methoxymethyl substituent is not held in a fixed conformation and may rotate around the C5-CH2 and CH2-O bonds. Examination of the potential energy as a function of rotation around these bonds indicates the presence of four low energy conformations. No hydrogen bonding is indicated for the methoxymethyl substituent, and the four potential minima result from reduced steric clash. For the MMdU.G base pair, the two bases adopt a wobble geometry which does not change with increasing solvent pH. Similarly, we find four low energy conformations for the methoxymethyl substituent in the major groove of the DNA helix.     

Solution structure of a DNA duplex containing the exocyclic lesion 3,N4-etheno-2'-deoxycytidine opposite 2'-deoxyguanosine

Vinyl chloride reacts with cellular DNA producing 3,N4-etheno-2'-deoxycytidine (epsilonC) along with other exocyclic adducts. The solution structure of an oligodeoxynucleotide duplex containing an epsilonC.dG base pair was determined by high-resolution NMR spectroscopy and molecular dynamics simulations. NMR data indicated that the duplex adopts a right-handed helical structure having all residues in anti orientation around the glycosylic torsion angle. The epsilonC adduct has a sugar pucker in the C3'-endo/C4'-exo region while the rest of the residues are in the C2'-endo/C3'-exo range. NOE interactions established Watson-Crick alignments for canonical base pairs of the duplex. The imino proton of the lesion-containing base pair resonated as a sharp signal that was resistant to water exchange, suggesting hydrogen bonding. Restrained molecular dynamics simulations generated three-dimensional models in excellent agreement with the spectroscopic data. The refined structures are slightly bent at the lesion site without major perturbations of the sugar-phosphate backbone. The adduct is displaced and shifted toward the major groove of the helix while its partner on the complementary strand remains stacked. The epsilonC(anti).dG(anti) base pair alignment is sheared and stabilized by the formation of hydrogen bonds. The biological implications of structures of epsilonC-containing DNA duplexes are discussed.     

FTIR and UV spectroscopy of parallel-stranded DNAs with mixed A*T/G*C sequences and their A*T/I*C analogues

The infrared spectra of parallel-stranded (ps) hairpin duplexes with mixed A*T/G*C composition and either isolated or sequential G*C pairs were studied in comparison with antiparallel-stranded (aps) duplexes and a corresponding set of molecules with hypoxanthine as a G base analogue lacking the exocyclic amino group. The ps duplexes showed the characteristic bands for the C2=O2 and C4=O4 stretching vibrations of thymine residues in trans-Watson-Crick A*T pairing at 1683 and 1668 cm-1. The latter band was superimposed on the stretching vibration of the free C6=O6 group of guanine. Substitution of guanine by hypoxanthine inhibited the formation of ps hairpin duplexes whatever the sequence, demonstrating that in the H-bonding between G and C the 2-NH2 group is necessary for stabilizing all of the investigated ps duplexes. This result is in agreement with a model of trans-Watson-Crick G*C base pairs with two H-bonds [N2H2(G)-N3(C) and N1H(G)-O2(C)]. However, trans-Watson-Crick A*T and G*C base pairs with two H-bonds are not isomorphous, which may explain the decreased stability of the ps, but not the aps, duplexes upon increasing the number of A*T/G*C steps. Molecular modeling studies performed on two of the ps duplexes reveal the existence of propeller twist for avoiding a clash between the N2(G) and N4(C) amino groups, and favorable stacking of sequential G*C base pairs. The optimized hairpin ps duplexes invariably incorporated G*C base pairs with two H-bonds, regardless of the initial structures adopted for the force field calculations.     

Crystal structure of an RNA octamer duplex r(CCCIUGGG)2 incorporating tandem I.U wobbles

The crystal structure of the RNA octamer duplex r(CCCIUGGG)2has been elucidated at 2.5 A resolution. The crystals belong to the space group P21and have unit cell constants a = 33.44 A, b = 43.41 A, c = 49.39 A and beta = 104.7 degrees with three independent duplexes (duplexes 1-3) in the asymmetric unit. The structure was solved by the molecular replacement method and refined to an Rwork/Rfree of 0.185/0.243 using 3765 reflections between 8.0 and 2.5 A. This is the first report of an RNA crystal structure incorporating I.U wobbles and three molecules in the asymmetric unit. Duplex 1 displays a kink of 24 degrees between the mismatch sites, while duplexes 2 and 3 have two kinks each of 19 degrees and 27 degrees, and 24 degrees and 29 degrees, respectively, on either side of the tandem mismatches. At the I.U/U.I mismatch steps, duplex 1 has a twist angle of 33.9 degrees, close to the average for all base pair steps, but duplexes 2 and 3 are underwound, with twist angles of 24.4 degrees and 26.5 degrees, respectively. The tandem I.U wobbles show intrastrand purine-pyrimidine stacking but exhibit interstrand purine-purine stacking with the flanking C.G pairs. The three independent duplexes are stacked non-coaxially in a head-to-tail fashion to form infinite pseudo-continuous helical columns which form intercolumn hydrogen bonding interactions through the 2'-hydroxyl groups where the minor grooves come together.     

Comparison of the solution structures of intramolecular DNA triple helices containing adjacent and non-adjacent CG.C+ triplets

The solution conformations of the intramolecular triple helices d(AGAAGA-X-TCTTCT-X-TC+TTC+T) and d(AAGGAA-X-TTCCTT-X-TTC+C+TT) (X = non-nucleotide linker) have been determined by NMR.1H NMR spectra in H2O showed that the third strand cytosine residues are fully paired with the guanine residues, each using two Hoogsteen hydrogen bonds. Determination of the13C chemical shifts of the cytosine C6 and C5 and their one-bond coupling constants (1 J CH) conclusively showed that the Hoogsteen cytosine residues are protonated at N3. The global conformations of the two molecules determined with >19 restraints per residue are very similar (RMSD = 0.96 A). However, some differences in local conformation and dynamics were observed for the central two base triplets of the two molecules. The C N3H were less labile in adjacent CG.C+triplets than in non-adjacent ones, indicating that the adjacent charge does not kinetically destabilize these triplets. The sugar conformations of the two adjacent cytosine residues were different and the 5'-residue was atypical of protonated cytosine. Hence, there are subtle effects of the interaction between two adjacent cytosine residues. The central two purines in each sequence showed non-standard backbone conformations, averaging between gamma approximately 60 degrees and gamma approximately 180 degrees. This may be related to the difference in the dependence of the thermodynamic stability on pH observed for these two sequences.     

RNA hydration: a detailed look

The crystal structure of the RNA duplex [r(CCCCGGGG)]2 has been refined to 1.46 A resolution with room temperature synchrotron diffraction data. This represents the highest resolution reported to date for an all-RNA oligonucleotide and is well beyond the best resolution ever achieved with an A-form DNA duplex. The analysis of the ordered hydration around the octamer duplex reveals conserved regular arrangements of water molecules in both grooves. In the major groove, all located first shell water molecules can be fitted into a pattern that is repeated through all eight base pairs, involves half the phosphate oxygens, and joins the two strands. In the minor groove, roughly across its narrowest dimension, tandem water molecules link the 2'-hydroxyl groups of adjacent nucleotides in base-pair steps in a similarly regular fashion. The structure provides evidence for an important role of the 2'-hydroxyl groups in the thermodynamic stabilization of RNA, beyond their known functions of locking the sugar pucker and mediating 3' --> 5' intrastrand O2'...O4' hydrogen bonds. The ribose 2'-hydroxyls lay the foundation for the enthalpic stability of the RNA relative to the DNA duplex, both as a scaffold for the water network in the minor groove and through their extensive individual hydration.     

1H NMR study of [d(GCGATCGC)]2 and its interaction with minor groove binding 4',6-diamidino-2-phenylindole

Two-dimensional NMR spectroscopy has been applied to study the solution binding of 4',6-diamidino-2-phenylindole (DAPI) to synthetic DNA duplex [d(GCGATCGC)]2. The structure of the complex at a molar ratio of 1:1 drug:duplex has been investigated. NMR results indicate that DAPI binds selectively in the minor groove of the DNA region containing only two A:T base pairs. The results disagree with conclusions drawn from footprinting experiments and show that the presence of the G3NH2 group in the minor groove does not prevent the binding. A molecular model is proposed that closely resembles the crystal structure previously published for the interaction of DAPI with the dodecamer [d(CGCGAATTCGCG)]2, containing four A:T base pairs in the binding site. In this model, DAPI lies in the minor groove, nearly isohelical, with its aromatic rings adjacent to H4' protons of T5 and C6 deoxyribose and the NH indole group oriented toward the DNA axis. The binding does not perturb the B-type conformation of the duplex, and the DNA oligomer conserves its 2-fold symmetry, indicating that fast exchange dynamics exist between the two stereochemically equivalent binding sites of the palindromic sequence. The binding constant and the exchange rate between free and bound species were also measured by NMR spectroscopy.     

Differential DNA recognition by the enantiomers of 1-Rh(MGP)2 phi: a combination of shape selection and direct readout

The enantiomers of the symmetric metallointercalator complex 1-Rh(MGP)2phi5+ [MGP = 4-(guanidylmethyl)-1,10-phenanthroline; phi = phenanthrenequinone diimine] bound to DNA decamer duplexes containing their respective 6 bp recognition sequences have been investigated using 1H NMR. Shape selection due to the chirality of the metal center and hydrogen-bonding contacts of ancillary guanidinium groups to 3'-G N7 atoms define the recognition by complexes which bind by intercalation to duplex DNA. The titration of Lambda-Rh into the self-complementary decamer containing the recognition sequence (5'-GACATATGTC-3', L1) resulted in one symmetric bound conformation observed in the 1H NMR spectrum, indicating that the DNA duplex retains its symmetry in the presence of the metal complex. Upfield chemical shifts of duplex imino protons and the disruption of the NOE base-sugar contacts defined the central T5-A6 intercalation site. The downfield shift of the G8 imino proton supports the conclusion that the pendant guanidinium arms make simultaneous H-bonding contacts to the N7 atoms of 3'-G8 bases on either side of the site. A variable-temperature study of a partially titrated sample (2:3 Lambda-Rh/L1) showed the exchange rate (kobs) at 298 K to be 68 s-1 and the activation barrier to exchange (DeltaG of association) to be 2.7 kcal/mol, a value comparable to the stacking energy of one base step. The results presented coupled with biochemical data are therefore consistent with binding models in which Lambda-1-Rh(MGP)2phi5+ (Lambda-Rh) traps the recognition site 5'-CATATG-3' in an unwound state, permitting intercalation centrally and hydrogen bonding to guanines at the first and sixth base pair positions. The data suggest a different model of binding and recognition by Delta-Rh. The titration of Delta-Rh into a DNA decamer containing the 6 bp recognition site (D1, 5'-CGCATCTGAC-3'; D2, 5'-GTCAGATGCG-3') resulted in two, distinct conformers, in slow exchange on the NMR time scale. The rate of exchange between the two conformers (kobs) at 298 K is 37 s-1, most likely due to partial dissociation between binding modes. The slower rate relative to Lambda-Rh association reflects the relative rigidity of the D1 and/or D2 sequence in comparison to L1. NOE cross-peaks between the intercalating phi ligand and protons of T5-C6, as well as the upfield shifts observed for imino protons at this step, serve to define the central T5-C6 step as the single site of intercalation. The downfield shift of the 3'-G imino protons indicates the complex makes hydrogen bond contacts with these bases. The complex, which is too small to span a 6 bp B-form DNA sequence, nonetheless makes major groove contacts with 3'-G bases to either side of the site. Notably, both 3'-guanine bases are necessary to impart site specificity and slow dissociation kinetics with the 5'-CATCTG-3' site, as evidenced by the extremely exchange-broadened two-dimensional NOESY spectra of Delta-Rh bound to modified duplexes containing N7-deazaguanine at either G8 or G18; the loss of one major groove contact completely abolishes specificity for 5'-CATCTG-3'. DNA chemical shifts upon binding and intermolecular NOE contacts therefore support a model in which Delta-Rh intercalates in one of two canted binding conformations. Within this model, each intercalation mode allows one guanidinium-guanine hydrogen bond at a time, while bringing the other arm close to the phosphate backbone.     

Triplex formation at physiological pH by oligonucleotides incorporating 5-Me-dC-(N4-spermine)

Oligonucleotide (ODN) directed triplex formation has therapeutic importance and depends on Hoogsteen hydrogen bonds between a duplex DNA and a third strand. While T*A:T triplets are formed at neutral pH, C+*G:C are favoured at acidic pH. Herein it is shown that 18-mer ODN containing spermine conjugated to 5-Me-dC at N4 (1-5), form triplexes with complementary 24-mer duplex 8:9 at neutral pH (7.3, 100 mM NaCl). Under such conditions, control ODN's carrying dC (6) or 5-Me-dC (7) did not show any triple helix formation. Remarkably, the triplexes from spermine-conjugates (1-5) have foremost stability at neutral pH (7.1), unlike the behavior of normal ODN's where optimal stability is at acidic pH (5.5). These results have importance in designing oligonucleotides for antigene applications.     

Solution conformation of the N-(deoxyguanosin-8-yl)-1-aminopyrene ([AP]dG) adduct opposite dC in a DNA duplex

Combined NMR-molecular mechanics computational studies were undertaken on the C8-deoxyguanosine adduct formed by the carcinogen 1-nitropyrene embedded in the d(C5-[AP]G6-C7).d(G16-C17-G18) sequence context in a 11-mer duplex, with dC opposite the modified deoxyguanosine. The exchangeable and nonexchangeable protons of the aminopyrene moiety and the nucleic acid were assigned following analysis of two-dimensional NMR data sets in H2O and D2O solution. There was a general broadening of several proton resonances for the three nucleotide d(G16-C17-G18) segment positioned opposite the [AP]dG6 lesion site resulting in weaker NOEs involving these protons in the adduct duplex. The solution conformation of the [AP]dG.dC 11-mer duplex has been determined by incorporating intramolecular and intermolecular proton-proton distances defined by upper and lower bounds deduced from NOESY spectra as restraints in molecular mechanics computations in torsion angle space. The aminopyrene ring of [AP]dG6 is intercalated into the DNA helix between intact Watson-Crick dC5.dG18 and dC7.dG16 base pairs. The modified deoxyguanosine ring of [AP]dG6 is displaced into the major groove and stacks with the major groove edge of dC5 in the adduct duplex. Both carbon and proton chemical shift data for the sugar resonances of the modified deoxyguanosine residue are consistent with a syn glycosidic torsion angle for the [AP]dG6 residue. The dC17 base on the partner strand is displaced from the center of the helix toward the major groove as a consequence of the aminopyrene ring intercalation into the helix. This base-displaced intercalative structure of the [AP]dG.dC 11-mer duplex exhibits several unusually shifted proton resonances which can be accounted for by the ring current contributions of the deoxyguanosinyl and pyrenyl rings of the [AP]dG6 adduct. In summary, intercalation of the aminopyrene moiety is accompanied by displacement of both [AP]dG6 and the partner dC17 into the major groove in the [AP]dG.dC 11-mer duplex.     

Stochastic model for a wavelike exothermal reaction in condensed heterogeneous systems

The labile protons of two 32-base-pair, four-arm models of immobile Holliday junctions have been studied by two-dimensional 1H nuclear magnetic resonance (NMR) spectroscopy. Overlap of resonances in the imino-imino region of two-dimensional nuclear Overhauser enhancement (NOE) spectra necessitates the use of a multi-pathway approach for obtaining sequence-specific assignments wherein all possible NOE connectivities to the labile protons are utilized, including those from the 2H of adenine, 5CH3 of thymine, and 5H of cytosine. Resonance assignments are obtained for all slowly exchanging imino and cytosine amino protons. Base-pairing up to and including the junction point is found in all four arms of both Holliday junctions. Several cross-arm NOE connectivities are identified and can be used to infer the geometry of the helical stacking domains. The two Holliday junctions studied, which differ only by the exchange of two base pairs at the branch point, appear to have opposite arm stacking geometries. These assignments form an important part of the critical background for detailed NMR analysis of Holliday junction structure and dynamics.     

Dual recognition of double-stranded DNA by 2'-aminoethoxy-modified oligonucleotides: the solution structure of an intramolecular triplex obtained by NMR spectroscopy

The solution structure of an intramolecular triple helical oligonucleotide has been solved by NMR. The third strand of the pyrimidine x purine x pyrimidine triplex is composed of 2'-aminoethoxy-modified riboses, whereas the remaining part of the nucleic acid is DNA. The structure around the aminoethoxy modification was obtained with the help of selective isotope labeling in conjunction with isotope-editing experiments. Dinucleotide steps and interstrand connectivities, as well as the complete backbone conformation of the triplex, were derived from J-couplings, NOEs, and 31P chemical shifts. The structure of this triplex, solved by distance geometry, explains the extraordinary stability and increase in rate of triplex formation induced by 2'-aminoethoxy-modified oligonucleotides: apart from the formation of seven base triples, a well-defined hydrogen-bonding network is formed across the Crick-Hoogsteen groove involving the amino protons of the aminoethoxy moieties and the phosphates of the purine strand of the DNA. The modified strand adopts a conformation which is close to an A-type helix, whereas the DNA duplex conformation is best described as an unwound B-type helix. The groove dimensions and helical parameters of the 2'-aminoethoxy-modified rY x dRdY triplex are surprisingly well conserved in comparison with DNA triplexes.     

Refined solution structure of 8,9-dihydro-8-(N7-guanyl)-9-hydroxyaflatoxin B1 opposite CpA in the complementary strand of an oligodeoxynucleotide duplex as determined by 1H NMR

The solution structure of d(CCATCAFBGATCC).d(GGATCAGATGG), containing the 8,9-dihydro-8-(N7-guanyl)-9-hydroxyaflatoxin B1 adduct, was refined using molecular dynamics restrained by NOE data obtained from 1H NMR. The modified guanosine was positioned opposite cytosine, while the aflatoxin moiety was positioned opposite adenosine in the complementary strand. Sequential 1H NOEs were interrupted between C5 and AFBG6, but intrastrand NOEs were traced through the aflatoxin moiety, via H6a of aflatoxin and H8 of the modified guanine. Opposite the lesion, the NOE between A16 H1' and G17 H8 was weak. A total of 43 NOEs were observed between DNA protons and aflatoxin protons. Molecular dynamics calculations restrained with 259 experimental and empirical distances, and using sp2 hybridization at AFBG6 N7, refined structures with pairwise rms differences < 0.85 A, excluding terminal base pairs. Relaxation matrix calculations yielded a sixth root rms difference between refined structures and NOE intensity data of 7.3 x 10(-2). The aflatoxin moiety intercalated on the 5'-face of the modified guanine. The extra adenine A16 was inserted between base pair AFBG6.C15 and the aflatoxin moiety. A 36 degree bending between the plane of base pair AFBG6.C15 and the plane of the aflatoxin moiety was predicted. The aflatoxin moiety stacked below the top domain of the oligodeoxynucleotide, which consisted of base pairs C1.G21, C2.G20, A3.T19, T4.A18, and C5.G17. The bottom domain consisted of base pairs AFBG6.C15, A7.T14, T8.A13, C9.G12, and C10.G11. The average winding angle between base pair C5.G17, the intercalated aflatoxin moiety, A16, and base pair AFBG6.C15 was reduced to 10 degrees. The preponderance of base pair substitutions in the aflatoxin B1 mutational spectrum, particularly G-->T transversions, suggests that the stability of this modified oligodeoxynucleotide, which models a templated +1 addition mutation, does not reliably predict the frequency of frame shifts.     

Thermodynamic, kinetic, and conformational properties of a parallel intermolecular DNA triplex containing 5' and 3' junctions

The interaction of the 11-mer oligodeoxypyrimidine d(TCTTCTUTCCT) with the 17 bp duplex d(CGCTAGAAGAAAGGACG).d(CGTCCUTTCTTCTAGCG) in forming an intermolecular DNA triplex has been examined in solution by surface plasmon resonance (SPR), UV thermal denaturation, circular dichroism (CD), and NMR methods. Thermodynamic data were also acquired for the shorter 15 bp target duplex d(CGCTAGAAGAAAGGA). d(TCCUTTCTTCTAGCG), which forms a 3' flush-ended parallel triplex. CD titrations at pH 5 gave a triplex --> (duplex + strand) dissociation constant Kd of 0.5 microM at 15 degreesC and approximately 2 microM at 25 degreesC for both the 11-15.15 and 11-17.17 systems, in agreement with analysis of the UV melting data and a direct calorimetric measurement. In contrast, the "apparent" Kd value determined by SPR was 10-20-fold smaller. The rate constant for dissociation (kd) of the third strand from the triplex was found to be approximately 0.0002 s-1 at 25 degreesC by SPR. The rate constant for exchange between the triplex and duplex states determined by NMR was approximately 2 s-1 at 40 degreesC. The dissociation kinetics measured by SPR are considerably underestimated, which largely accounts for the poor estimation of Kd using this technique. Extensive 1H NMR assignments were obtained for both the 17 bp DNA duplex and the triplex. Large changes in chemical shifts were observed in the purine strand of the host duplex, but only small shift changes were induced in the complementary pyrimidine strand. Dramatic differences in shifts were observed for the G and A residues, especially in the minor groove, consistent with only small, localized conformational changes in the underlying duplex. The magnitude of the shift changes decreased to baseline within one base of the 3' triplex-duplex junction and over two to three bases at the 5' junction. Chemical shift changes at the 5' junction suggest small conformational anomalies at this site. COSY and NOESY spectra indicate that the nucleotides are in the "S" domain in both the triplex and duplex states. These data rule out major conformation changes at the triplex-duplex boundaries. NOEs between pyrimidines in the third strand and those in the duplex showed proximity for these bases in the major groove, which could be ascribed to buckling of the Hoogsteen bases out of the plane of the Watson-Crick base pairs.     

Conformations of nicked and gapped DNA structures by NMR and molecular dynamic simulations in water

We have analyzed and compared the molecular structures and dynamics of DNA duplexes containing a nick or a gap of one nucleotide where the base in front of the gap is a guanine. The continuous strand has the sequence 5'(CAGAGTCXCTGGCTC) where the residue X is absent for the nick, 14-mer, and where it is a G residue for the gap. Duplexes were formed with the two corresponding 7-mers. Neither of these is phosphorylated adjacent at the nick site, but it is a good model for a single strand break. For the nick structure, the quantitative NMR data show that the global conformation is very close to canonical B-form DNA, but it displays enhanced local flexibility. For the gap structure, we observe only one species in which the extra G is well stacked into the helix. The two half-helices around this residue also show a B-form conformation. As with the nick duplex, the adjacent G imino protons show enhanced exchange with solvent. The gap does not close completely. Using distance constraints, MD calculations show that the nick conformation is very close to a duplex with no lesion but is indeed more flexible in the central part. The gapped structure shows two families of conformations. One is close to B-DNA, the other is significantly kinked at the gap which reduces the size of the cavity. We observe a spine of hydration within the cavities, similar, but of different geometry in the two cases.     

Extension of the range of DNA sequences available for triple helix formation: stabilization of mismatched triplexes by acridine-containing oligonucleotides

Triple helix formation usually requires an oligopyrimidine*oligopurine sequence in the target DNA. A triple helix is destabilized when the oligopyrimidine*oligopurine target contains one (or two) purine*pyrimidine base pair inversion(s). Such an imperfect target sequence can be recognized by a third strand oligonucleotide containing an internally incorporated acridine intercalator facing the inverted purine*pyrimidine base pair(s). The loss of triplex stability due to the mismatch is partially overcome. The stability of triplexes formed at perfect and imperfect target sequences was investigated by UV thermal denaturation experiments. The stabilization provided by an internally incorporated acridine third strand oligonucleotide depends on the sequences flanking the inverted base pair. For triplexes containing a single mismatch the highest stabilization is observed for an acridine or a propanediol tethered to an acridine on its 3'-side facing an inverted A*T base pair and for a cytosine with an acridine incorporated to its 3'-side or a guanine with an acridine at its 5'-side facing an inverted G*C base pair. Fluorescence studies provided evidence that the acridine was intercalated into the triplex. The target sequences containing a double base pair inversion which form very unstable triplexes can still be recognized by oligonucleotides provided they contain an appropriately incorporated acridine facing the double mismatch sites. Selectivity for an A*T base pair inversion was observed with an oligonucleotide containing an acridine incorporated at the mismatched site when this site is flanked by two T*A*T base triplets. These results show that the range of DNA base sequences available for triplex formation can be extended by using oligonucleotide intercalator conjugates.