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.     

Structural alignment of the (+)-trans-anti-benzo[a]pyrene-dG adduct positioned opposite dC at a DNA template-primer junction

This study reports on the solution conformation of the covalent (+)-trans-anti-[BP]dG adduct (derived from the binding of the highly mutagenic and tumorigenic (+)-anti-benzo[a]pyrene diol epoxide to the N2 of deoxyguanosine) positioned opposite dC at a junctional site in the d(A1-A2-C3-[BP]G4-C5- T6-A7-C8-C9-A10-T11-C12-C13).d(G14-G15-A16-T17-+ ++G18-G19-T20-A21-G22-C23) 13/10-mer DNA sequence. The 13-mer represents the template strand containing the junction [BP]dG4 lesion while the complementary 10-mer models a primer strand which extends upto and is complementary to the modified dG4 residue. The solution conformation has been determined by initially incorporating intramolecular and intermolecular proton-proton distances defined by lower and upper bounds deduced from NOESY spectra as restraints in molecular mechanics computations in torsion angle space and subsequently through restrained molecular dynamics calculations based on a NOE distance and intensity refinement protocol. The duplex segment retains a minimally perturbed B-DNA conformation with all base pairs, including the junctional [BP]dG4.dC23 pair, in Watson-Crick hydrogen-bonded alignments. The pyrenyl ring is not stacked over the adjacent dC5.dG22 base pair but is positioned on the minor groove-side of the [BP]dG moiety and directed toward the 5'-end of the template strand. The pyrenyl ring stacks over the base of the non-adjacent dA2 residue in one direction and the sugar ring of dC23 in the other direction. The solution structure of the (+)-trans-anti-[BP]dG adduct opposite dC in the 13/10-mer in which the modified deoxyguanosine adopts an anti glycosidic torsion angle (this study) is in striking contrast to the structure of the same (+)-trans-anti-[BP]dG moiety in a 13/9-mer of the same sequence but without the dC23 residue positioned opposite the adduct site [Cosman, M., et al. (1995) Biochemistry 34, 15334-15350]. For the latter case, the aromatic portion of the BP residue stacks over the adjacent dC5.dG22 base pair, the modified deoxyguanosine adopts a syn glycosidic torsion angle and is displaced toward the major groove direction. Insights into the factors that affect the sequence and context dependent conformations of stereoisomeric [BP]dG lesions have emerged following comparison of these two structures with the minor groove conformations of the same (+)-trans-anti-[BP]dG lesion in the fully complementary 11-mer duplex [Cosman, M., et al. (1992) Proc. Natl. Acad. Sci. U.S.A. 89, 1914-1918] and in the base displaced-intercalative conformation of the 11/10-mer deletion duplex containing a -1 deletion site opposite the lesion [Cosman, M., et al. (1994) Biochemistry 33, 11507-11517]. The contributing factors where applicable include Watson-Crick base pairing at the site of the lesion, positioning of the carcinogen within the floor of the minor groove, and the tendency of the bulky hydrophobic aromatic BP residue to assume stacked or intercalative conformations.     

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.     

Structure of a DNA duplex containing a single 2'-O-methyl-beta-D-araT: combined use of NMR, restrained molecular dynamics, and full relaxation matrix refinement

Two-dimensional 1H NMR spectroscopy was used to determine the solution structure of the double-stranded DNA oligonucleotide d(5'-CGCATATAGCC-3'): d(5'-GGCTAXATGCG-3'), where X is 1-(2-O-methyl-beta-D-arabinofuranosyl)thymine. The structure determination was based on a total relaxation matrix analysis of NOESY cross-peak intensities using the MARDIGRAS program. The improved RANDMARDI procedure was used during the calculations to include the experimental "noise" in the NOESY spectra. The NOE-derived distance restraints were applied in restrained molecular dynamics calculations. Twenty final structures each were generated for the modified DNA duplex from both A-form and B-form DNA starting structures. The root-mean-square deviation of the coordinates for the 40 structures was 0.82 A. The duplex adopts a normal B-DNA-type helix, and the spectra as well as the structure show that the modified nucleotide X adopts a C2'-endo (S) sugar conformation. There are no significant changes in the helix originating from the modified nucleotide. The CH3O group on X is directed toward the major groove, and there seems to be free space for further modifications at this position.     

The impact of an exocyclic cytosine adduct on DNA duplex properties: significant thermodynamic consequences despite modest lesion-induced structural alterations

The exocyclic base adduct 3,N4-deoxyethenocytosine (epsilonC) is a common DNA lesion that can arise from carcinogen exposure and/or as a biproduct of cellular processes. We have examined the thermal and thermodynamic impact of this lesion on DNA duplex properties, as well as the structural alterations imparted by the lesion. For these studies, we used calorimetric and spectroscopic techniques to investigate a family of 13-mer DNA duplexes of the form (5'CGCATGNGTACGC3')x(3'GCGTACNCATGCG5'), where the central NxN base pair represents the four standard Watson-Crick base pairs (corresponding to four control duplexes), and where either one of the N bases has been replaced by epsilonC, yielding eight test duplexes. Studies on these 12 duplexes permit us to assess the impact of the epsilonC lesion as a function of sequence context. Our spectroscopic and calorimetric data allow us to reach the following conclusions: (i) The epsilonC lesion imparts a large penalty on duplex stability, with sequence context only modestly modulating the extent of this lesion-induced destabilization. This result contrasts with our recent studies of duplexes with abasic sites, where sequence context was found to be the predominant determinant of thermodynamic damage. (ii) For the epsilonC-containing duplexes, sequence context effects are most often observed in the enthalpic contribution to lesion-induced duplex destabilization. However, due to compensating entropies, the free energy changes associated with this lesion-induced duplex destablization are nearly independent of sequence context. (iii) Despite significant lesion-induced changes in duplex energetics, our spectroscopic probes detect only modest lesion-induced changes in duplex structure. In fact, the overall duplex maintains a global B-form conformation, in agreement with NMR structural data. We discuss possible interpretations of the apparent disparity between the severe thermodynamic and relatively mild structural impacts of the epsilonC lesion on duplex properties. We also note and discuss the implications of empirical correlations between biophysical and biological properties of lesion-containing duplexes.     

Solution structure of the aminofluorene-intercalated conformer of the syn [AF]-C8-dG adduct opposite a--2 deletion site in the NarI hot spot sequence context

This paper addresses structural issues related to the capacity of aminofluorene [AF] for frameshift mutations of the -2 type on C8 covalent adduct formation at the G3 site in the d(C-G1-G2-C-G3-C-C) NarI hot spot sequence. This problem has been approached from a combined NMR and relaxation matrix analysis computational structural study of the [AF]dG adduct in the d(C-G-G-C-[AF]G-C-C).d(G-G-C-C-G) sequence context at the 12/10-mer adduct level (designated [AF]dG.del(-2) 12/10-mer). The proton spectra of this system are of exceptional quality and are consistent with the formation of an AF-intercalated conformer with the modified guanine in a syn alignment displaced along with the 5'-flanking cytosine residue into the major groove. The solution structure has been determined by initially incorporating intramolecular and intermolecular proton-proton distances defined by lower and upper bound deduced from NOESY spectra as restraints in molecular mechanics computations in torsion angle space and subsequently refined through restrainted molecular dynamics calculations based on a NOE distance and intensity refinement protocol. Strikingly, the [AF]dG.del(-2) 12/10-mer duplex adopts only one of two potential AF-intercalation alignments for the [AF]dG adduct opposite the -2 deletion site in the NarI sequence context with the extrusion of the dC-[AF]dG step favored completely over extrusion of the [AF]dG-dC step at the lesion site. This polarity establishes that the structural perturbation extends 5' rather than 3' to the [AF]dG lesion site in the adduct duplex. This structure of the [AF]dG adduct opposite a -2 deletion site shows distinct differences with conclusions reported on the alignment of the related acetylaminofluorene [AAF]dG adduct opposite a -2 deletion site in the identical NarI sequence context [Milhe, C., Fuchs, R. P. P., and Lefevre, J. F. (1996) Eur. J. Biochem. 235, 120-127]. In that study, qualitative NMR data without computational analysis were employed to conclude that the extrusion at the lesion site occurs at the [AAF]dG-dC step for the AAF-intercalated conformer of the adduct duplex. The structure of the [AF]dG adduct opposite a -2 deletion site determined in our group provides molecular insights into the architecture of extended slipped mutagenic intermediates involving aromatic amine intercalation and base-displaced syn modified guanines in AF and, by analogy, AAF-induced mutagenesis in the NarI hot spot sequence context.     

A molecular mechanics and dynamics study of the minor adduct between DNA and the carcinogen 2-(acetylamino)fluorene (dG-N2-AAF)

Experimental studies involving the carcinogenic aromatic amine 2-(acetylamino)fluorene (AAF) have afforded two acetylated DNA adducts, the major one bound to C8 of guanine and a minor adduct bound to N2 of guanine. The minor adduct may be important in carcinogenesis because it persists, while the major adduct is rapidly repaired. Primer extension studies of the minor adduct have indicated that it blocks DNA synthesis, with some bypass and misincorporation of adenine opposite the lesion [Shibutani, S., and Grollman, A.P. (1993) Chem. Res. Toxicol. 6, 819-824]. No experimental structural information is available for this adduct. Extensive minimized potential energy searches involving thousands of trials and molecular dynamics simulations were used to study the conformation of this adduct in three sequences: I, d(C1-G2-C3-[AAF]G4-C5-G6-C7).d(G8-C9-G10-C11-G12-C13-G14+ ++); II, the sequence of Shibutani and Grollman, d(C1-T2-A3-[AAF]G4-T5-C6-A7).d(T8-G9-A10-C11-T12-A13-G14); and III, which is the same as II but with a mismatched adenine in position 11, opposite the lesion. AAF was located in the minor groove in the low-energy structures of all sequences. In the lowest energy form of the C3-[AAF]G4-C5 sequence I, the fluorenyl rings point in the 3' direction along the modified strand and the acetyl in the 5' direction. These orientations are reversed in the second lowest energy structure of this sequence, and the energy of this structure is 1.4 kcal/mol higher. Watson Crick hydrogen bonding is intact in both structures. In the two lowest energy structures of the A3-[AAF]G4-T5 sequence II, the AAF is also located in the minor groove with Watson-Crick hydrogen bonding intact. However, in the lowest energy form, the fluorenyl rings point in the 5' direction and the acetyl in the 3' direction. The energy of the structure with opposite orientation is 5.1 kcal/mol higher. In sequence III with adenine mismatched to the modified guanine, the lowest energy form also had the fluorenyl rings oriented 5' in the minor groove with intact Watson-Crick base pairing. However, the mispaired adenine adopts a syn orientation with Hoogsteen pairing to the modified guanine. These results suggest that the orientation of the AAF in the minor groove may be DNA sequence dependent. Mobile aspects of favored structures derived from molecular dynamics simulations with explicit solvent and salt support the essentially undistorting nature of this lesion, which is in harmony with its persistence in mammalian systems.     

Unusual nucleotide conformations in GNRA and UNCG type tetraloop hairpins: evidence from Raman markers assignments

High resolution NMR data on UNCG and GNRA tetraloops (where N is any of the four nucleotides and R is a purine) have shown that they contain ribonucleosides with unusual 2'-endo/anti and 3'-endo/syn conformations, in addition to the 3'-endo/anti ones which are regularly encountered in RNA chains. In the current study, Raman spectroscopy has been used to probe these nucleoside conformations and follow the order (hairpin) to disorder (random chain) structural transitions in aqueous phase in the 5-80 degreesC temperature range. Spectral evolution of GCAA and GAAA tetraloops, as formed in very short hairpins with only three G.C base pairs in their stems (T m >60 degreesC), are reported and compared with those previously published on UUCG and UACG tetraloops, for which the syn orientation of the terminal guanine as well as the 2'-endo/anti conformation of the third rC residue have been confirmed by means of vibrational marker bands. Raman data obtained as a function of temperature show that the first uracil in the UUCG tetraloop is stacked and the two middle residues (rU and rC) are in the 2'-endo/anti conformation, in agreement with the previously published NMR results. As far as the new data concerning the GNRA type tetraloops are concerned, they lead us to conclude that: (i) in both cases (GCAA and GAAA tetraloops) the adenine bases are stacked; (ii) the second rC residue in the GCAA tetraloop has a 3'-endo/anti conformation; (iii) the sugar pucker associated with the third rA residue in both tetraloops possibly undergoes a 3'-endo/2'-endo interconversion as predicted by NMR results; (iv) the stem adopts a regular A-form structure; (v) all other nucleosides of these two GNRA tetraloops possess the usual 3'-endo/anti conformation.     

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.     

Neomycin, spermine and hexaamminecobalt (III) share common structural motifs in converting B- to A-DNA

The (dG)n.(dC)n-containing 34mer DNA duplex [d(A2G15C15T2)]2 can be effectively converted from the B-DNA to the A-DNA conformation by neomycin, spermine and Co(NH3)6(3+). Conversion is demonstrated by a characteristic red shift in the circular dichroism spectra and dramatic NMR spectral changes in chemical shifts. Additional support comes from the substantially stronger CH6/GH8-H3'NOE intensities of the ligand-DNA complexes than those from the native DNA duplex. Such changes are consistent with a deoxyribose pucker transition from the predominate C2'-endo (S-type) to the C3'-endo (N-type). The changes for all three ligand-DNA complexes are identical, suggesting that those three complex cations share common structural motifs for the B- to A-DNA conversion. The A-DNA structure of the 4:1 complex of Co(NH3)6(3+)/d(ACCCGCGGGT) has been analyzed by NOE-restrained refinement. The structural basis of the transition may be related to the closeness of the two negatively charged sugar-phosphate backbones along the major groove in A-DNA, which can be effectively neutralized by the multivalent positively charged amine functions of these ligands. In addition, ligands like spermine or Co(NH3)6(3+) can adhere to guanine bases in the deep major groove of the double helix, as is evident from the significant direct NOE cross-peaks from the protons of Co(NH3)6(3+) to GH8, GH1 (imino) and CH4 (amino) protons. Our results point to future directions in preparing more potent derivatives of Co(NH3)6(3+) for RNA binding or the induction of A-DNA.     

An unusual DNA adduct derived from the powerfully mutagenic environmental contaminant 3-nitrobenzanthrone

The covalent binding of an N-hydroxy metabolite of the powerfully mutagenic 3-nitrobenzanthrone (NBA) to 2'-deoxyguanosine (dG) and calf thymus DNA has been investigated in vitro. The major adduct obtained from the reaction of the N-acetoxy-N-acetyl derivative (N-Aco-N-Ac-ABA) of 3-aminobenzanthrone (ABA) and dG was identified as N-acetyl-3-amino-2-(2'-deoxyguanosin-8-yl)benzanthrone (dG-N-Ac-ABA) by 1H NMR and mass spectroscopies as well as by the reaction of N-Aco-N-Ac-ABA with the double-stranded calf thymus DNA. The coupling with the dG moiety occurred exclusively at C-2 of benzanthrone (BA), suggesting a significant contribution of a resonance-stabilized arenium ion intermediate derived from BA to the production of this new type of adduct. The preferred conformation of the adduct has been shown to be syn by 1H and 13C NMR.     

Solution conformation of the (-)-trans-anti-[BP]dG adduct opposite a deletion site in a DNA duplex: intercalation of the covalently attached benzo[a]pyrene into the helix with base displacement of the modified deoxyguanosine into the minor groove

A combined NMR-computational approach was employed to determine the solution structure of the (-)-trans-anti-[BP]dG adduct positioned opposite a -1 deletion site in the d(C1-C2-A3-T4-C5- [BP]G6-C7-T8-A9-C10-C11).d(G12-G13-T14-A15-G1 6-G17-A18-T19-G20-G21) sequence context. The (-)-trans-anti-[BP]dG moiety is derived from the binding of the (-)-anti-benzo[a]pyrene diol epoxide [(-)-anti-BPDE] to N2 of dG6 and has a 10R absolute configuration at the [BP]dG linkage site. The exchangeable and non-exchangeable protons of the benzo[a]pyrenyl moiety and the nucleic acid were assigned following analysis of two-dimensional NMR data sets in H2O and D2O solution. The solution conformation has been determined by incorporating intramolecular and intermolecular proton-proton distances defined by lower and upper bounds deduced from NOESY spectra as restraints in molecular mechanics computations in torsion angle space followed by restrained molecular dynamics calculations based on a NOE distance and intensity refinement protocol. Our structural studies establish that the aromatic BP ring system intercalates into the helix opposite the deletion site, while the modified deoxyguanosine residue is displaced into the minor groove with its face parallel to the helix axis. The intercalation site is wedge-shaped and the BP aromatic ring system stacks over intact flanking Watson-Crick dG.dC base pairs. The modified deoxyguanosine stacks over the minor groove face of the sugar ring of the 5'-flanking dC5 residue. The BP moiety is positioned with the benzylic ring oriented toward the minor groove and the distal pyrenyl aromatic ring directed toward the major groove. This conformation strikingly contrasts with the corresponding structure in the full duplex with the same 10R (-)-trans-anti-[BP]dG lesion positioned opposite a complementary dC residue [de los Santos et al. (1992) Biochemistry 31, 5245-5252); in this case the aromatic BP ring system is located in the minor groove, and there is no disruption of the [BP]dG.dC Watson-Crick base pairing alignment. The intercalation-base displacement features of the 10R (-)-trans-anti-[BP]dG adduct opposite a deletion site have features in common to those of the 10S (+)-trans-anti-[BP]dG adduct opposite a deletion site previously reported by Cosman et al. [(1994)(Biochemistry 33, 11507-11517], except that there is a nearly 180 degrees rotation of the BP residue about the axis of the helix at the base-displaced intercalation site and the modified deoxyguanosine is positioned in the opposite groove. In the 10S adduct, the benzylic ring is in the major groove and the aromatic ring systems point toward the minor groove. This work extends the theme of opposite orientations of adducts derived from chiral pairs of (+)- and (-)-anti-BPDE enantiomers; both 10S and 10R adducts can be positioned with opposite orientations either in the minor groove or at base displaced intercalation sites, depending on the presence or absence of the partner dC base in the complementary strand.     

Programmed cell death of identified peptidergic neurons involved in ecdysis behavior in the Moth, Manduca sexta

Molecular dynamics simulation in explicit solvent and continuum solvent models are applied to investigate the relative stability of A- and B-form helices for two DNA sequences, dA10-dT10 and dG10-dC10 in three structural forms. One structural form is based on an unrestrained molecular dynamics (MD) trajectory starting from a canonical B-DNA structure, the second is based on a MD trajectory starting in a canonical B-DNA structure with the sugars constrained to be C2'-endo and the third simulation started from a canonical A-DNA structure with the sugars constrained to C3'-endo puckers. For the energetic analysis, structures were taken as snapshots from nanosecond length molecular dynamics simulations computed in a consistent fashion in explicit solvent, applying the particle mesh Ewald method and the Cornell et al. force field. The electrostatic contributions to solvation free energies are computed using both a finite-difference Poisson-Boltzmann model and a pairwise Generalized Born model. The non-electrostatic contributions to the solvation free energies are estimated with a solvent accessible surface area dependent term. To estimate the gas phase component of the relative free energy between the various structures, the mean solute internal energies (determined with the Cornell et al. molecular mechanics potential including all pairwise interactions within the solute) and estimates of the solute entropy (using a harmonic approximation) were used. Consistent with experiment, the polyG-polyC (GC) structures are found to be much more A-phillic than the polyA-polyT (AT) structures, the latter being quite A-phobic. The dominant energy components responsible for this difference comes from the internal and van der Waal energies. A perhaps less appreciated difference between the GC and AT rich sequences is suggested by the calculated salt dependence which demonstrates a significantly enhanced ability to drive GC rich sequences towards an A-form structure compared to AT rich sequences. In addition to being A-phobic, the AT structure also has a noticably larger helical repeat than GC and other mixed sequence duplexes, consistent with experiment. Analysis of the average solvent density from the trajectories shows hydration patterns in qualitative agreement with experiment and previous theoretical treatments.     

Structural features of the DNA hairpin d(ATCCTA-GTTA-TAGGAT): formation of a G-A base pair in the loop

The three-dimensional structure of the hairpin formed by d(ATCCTA-GTTA-TAGGAT) has been determined by means of two-dimensional NMR studies, distance geometry and molecular dynamics calculations. The first and the last residues of the tetraloop of this hairpin form a sheared G-A base pair on top of the six Watson-Crick base pairs in the stem. The glycosidic torsion angles of the guanine and adenine residues in the G-A base pair reside in the anti and high- anti domain ( approximately -60 degrees ) respectively. Several dihedral angles in the loop adopt non-standard values to accommodate this base pair. The first and second residue in the loop are stacked in a more or less normal helical fashion; the fourth loop residue also stacks upon the stem, while the third residue is directed away from the loop region. The loop structure can be classified as a so-called type-I loop, in which the bases at the 5'-end of the loop stack in a continuous fashion. In this situation, loop stability is unlikely to depend heavily on the nature of the unpaired bases in the loop. Moreover, the present study indicates that the influence of the polarity of a closing A.T pair is much less significant than that of a closing C.G base pair.     

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.     

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.     

Synthesis, miscoding specificity, and thermodynamic stability of oligodeoxynucleotide containing 8-methyl-2'-deoxyguanosine

8-Methyl-2'-deoxyguanosine (8-MedG) was synthesized by reacting dG under the methyl radical generating system and incorporated into oligodeoxynucleotides using phosphoramidite techniques. The site-specifically modified oligodeoxynucleotide containing a single 8-MedG was then used as a template for primer extension reactions catalyzed by the 3' --> 5' exonuclease-free (exo-) Klenow fragment of Escherichia Coli DNA polymerase I and mammalian DNA polymerase alpha. Primer extension catalyzed by the exo- Klenow fragment readily passed the 8-MedG lesion in the template while that catalyzed by pol alpha was retarded opposite the lesion. The fully extended products formed during DNA synthesis were analyzed to quantify the miscoding specificities of 8-MedG. Both DNA polymerases incorporated primarily dCMP, the correct base opposite the lesion, along with small amounts of incorporation of dGMP and dAMP. In addition, two-base deletion was observed only when the exo- Klenow fragment was used. The thermodynamic stability of 8-MedG in the duplex was also studied. The duplex containing 8-MedG:dG was more thermally and thermodynamically stable than that of dG:dG. The duplex containing 8-MedG:dA was more thermodynamically stable than that of dG:dA. We conclude that 8-MedG is a miscoding lesion and capable of generating G --> C and G --> T transversions and deletion in cells.     

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.     

NMR studies of the effects of the 5'-phosphate group on conformational properties of 5-methylaminomethyluridine found in the first position of the anticodon of Escherichia coli tRNA(Arg)4

5-Methylaminomethyluridine (mnm5U) exists in the first position of the anticodon (position 34) of Escherichia coli tRNA4Arg for codons AGA/AGG. In the present study, the temperature dependence of the ribose-puckering equilibrium of pmnm5U was analyzed by proton NMR spectroscopy. Thus, the enthalpy difference (delta H) between the C2'-endo and C3'-endo forms was obtained at 0.65 kcal.mol-1. By comparison of the delta H values of pU and pmnm5U, the 5-substitution was found to increase the relative stability of the C3'-endo form over the C2'-endo form significantly (by 0.56 kcal.mol-1). Furthermore, this conformational "rigidity" was concluded to depend on the 5'-phosphate group, because nucleoside U exhibits only a negligible change in the ribose-puckering equilibrium upon the 5-methylaminomethyl substitution. Further NMR analyses and molecular dynamics calculations revealed that interactions between the 5-methylaminomethyl and 5'-phosphate groups of pmnm5U restrict the conformation about the glycosidic bond to a low anti form, enhancing steric repulsion between the 2-carbonyl and 2'-hydroxyl groups in the C2'-endo form. This intrinsic conformational rigidity of the mnm5U residue in position 34 may contribute to the correct codon recognition.