A-DNA accommodates adducts derived from diol epoxides of polycyclic aromatic hydrocarbons bound in a "side-stacking" mode

Molecular modeling of the major benzo[a]pyrene N2-dG adduct in cases where mutagenesis results are known in double stranded DNA

The potent mutagen/carcinogen benzo[a]pyrene (B[a]P) is metabolically activated to (+)-anti-B[a]PDE, which induces a full spectrum of mutations (e.g. GC-->TA, GC-->AT, etc.). One hypothesis for this complexity is that different mutations are induced by different conformations of its major adduct [+ta]-B[a]P-N2-dG when bypassed during DNA replication (probably by different DNA polymerases). Previous molecular modeling studies suggested that B[a]P-N2-dG adducts can in principle adopt at least 16 potential conformational classes in ds-DNA. Herein we report on molecular modeling studies with the eight conformations most likely to be relevant to base substitution mutagenesis in 10 cases where mutagenesis has been studied in ds-DNA plasmids in E. coli with B[a]P-N2-dG adducts of differing stereoisomers and DNA sequence contexts, as well as in five cases where the conformation is known by NMR. Of the approximately 11,000 structures generated in this study, the computed lowest energy structures are reported for 120 cases (i.e. eight conformations and 15 examples), and their conformations compared. Of the eight conformations, four are virtually always computed to be high in energy. The remaining four lower energy conformations include two with the BP moiety in the minor groove (designated: BPmi5 and BPmi3), and two base-displaced conformations, one with the dG moiety in the major groove (designated: Gma5) and one with the dG in the minor groove (designated: Gmi3). Interestingly, these four are the only conformations that have been observed for B[a]P-N2-dG adducts in NMR studies. Independent of sequence contexts and adduct stereochemistry, BPmi5 structures tend to look reasonably similar, as do BPmi3 structures, while the base-displaced structures Gma5 and BPmi3 tend to show greater variability in structure. A correlation was sought between modeling and mutagenesis results in the case of the low energy conformations BPmi5, BPmi3, Gma5 and Gma3. Plots of log[(G-->T)/(G-->A)] versus energy[(conformation X)-(conformation Y)] were constructed for all six pairwise combinations of these four conformations, and the only plot giving a straight line involved Gma5 and Gmi3. While this finding is striking, its significance is unclear (as discussed).

Toward an understanding of the role of DNA adduct conformation in defining mutagenic mechanism based on studies of the major adduct (formed at N(2)-dG) of the potent environmental carcinogen, benzo[a]pyrene

The process of carcinogenesis is initiated by mutagenesis, which often involves replication past damaged DNA. One question - what exactly is a DNA polymerase seeing when it incorrectly copies a damaged DNA base (e.g., inserting dATP opposite a dG adduct)? - has not been answered in any case. Herein, we reflect on this question, principally by considering the mutagenicity of one activated form of benzo[a]pyrene, (+)-anti-B[a]PDE, and its major adduct [+ta]-B[a]P-N(2)-dG. In previous work, [+ta]-B[a]P-N(2)-dG was shown to be capable of inducing>95% G-->T mutations in one sequence context (5'-TGC), and approximately 95% G-->A mutations in another (5'-AGA). This raises the question - how can a single chemical entity induce different mutations depending upon DNA sequence context? Our current working hypothesis is that adduct conformational complexity causes adduct mutational complexity, where DNA sequence context can affect the former, thereby influencing the latter. Evidence supporting this hypothesis was discussed recently (Seo et al., Mutation Res. [in press]). Assuming this hypothesis is correct (at least in some cases), one goal is to consider what these mutagenic conformations might be. Based on molecular modeling studies, 16 possible conformations for [+ta]-B[a]P-N(2)-dG are proposed. A correlation between molecular modeling and mutagenesis work suggests a hypothesis (Hypothesis 3): a base displaced conformation with the dG moiety of the adduct in the major vs. minor groove gives G-->T vs. G-->A mutations, respectively. (Hypothesis 4, which is a generalized version of Hypothesis 3, is also proposed, and can potentially rationalize aspects of both [+ta]-B[a]P-N(2)-dG and AP-site mutagenesis, as well as the so-called "A-rule".) Finally, there is a discussion of how conformational complexity might explain some unusual mutagenesis results that suggest [+ta]-B[a]P-N(2)-dG can become trapped in different conformations, and why we think it makes sense to interpret adduct mutagenesis results by modeling ds-DNA (at least in some cases), even though the mutagenic event must occur at a ss/ds-DNA junction in the presence of a DNA polymerase.

Solution structure of an oligodeoxynucleotide containing the human N-ras codon 12 sequence refined from 1H NMR using molecular dynamics restrained by nuclear Overhauser effects

The structure of d(GGCAGGTGGTG).d(CACCACCTGCC), consisting of codons 11, 12 (underlined), and 13 of the human n-ras protooncogene, was refined from 1H NMR data. Patterns of internucleotide NOEs consistent with a B-form helix were observed for each strand. NOE intensities between purine H8 and H1' protons were small compared to intensities between cytosine H5 and H6 protons, indicative of glycosyl torsion angles in the anti range. Cross-peaks were observed between purine H8 and pyrimidine H5 and CH3 protons on adjacent bases in the direction of purine (5'-->3')pyrimidine, but not in the direction pyrimidine(5'-->3')purine. Watson-Crick hydrogen bonding between bases was intact. A total of 232 experimental distance restraints were obtained. Of these, 143 were intra-residue restraints and 89 were inter-residue restraints. A restrained molecular dynamics/simulated annealing approach yielded 6 MD structures calculated from a B-form starting structure and 6 MD structures from an A-form starting structure. These refined to an average pairwise rms difference of 0.92 angstrom, with maximum pairwise rmsd of 1.35 angstroms. Accuracy of emergent structures was assessed by relaxation matrix back-calculation. The sixth-root residual index of 7.0 x 10(-2) measured between the refined structures and the NOE intensity data suggested that the former were in reasonable agreement with the NOE data. The refined solution structures were in the B-family. Similar to the human n-ras codon 61 sequence [Feng, B., & Stone, M.P. (1995) Chem. Res. Toxicol. 8, 821-832], the ras12 sequence contained local variations in B-like conformation which did not confer large structural alterations upon the duplex, but perhaps modulated the reactivity of the first as compared to the second guanine in codon 12.

Solution conformation of the major adduct between the carcinogen (+)-anti-benzo[a]pyrene diol epoxide and DNA

We have synthesized, separated, and purified approximately 10 mg of a deoxyundecanucleotide duplex containing a single centrally positioned covalent adduct between (+)-anti-benzo[a]pyrene (BP) diol epoxide and the exocyclic amino group of guanosine. Excellent proton NMR spectra are observed for the (+)-trans-anti-BP diol epoxide-N2-dG adduct positioned opposite dC and flanked by G.C pairs in the d[C1-C2-A3-T4-C5-(BP)G6-C7-T8-A9-C10-C11].d[12- G13-T14-A15-G16-C17-G18-A19-T20-G 21-G22] duplex +ADdesignated (BP)G.C 11-mer+BD. We have determined the solution structure centered about the BP covalent adduct site in the (BP)G.C 11-mer duplex by incorporating intramolecular and intermolecular proton-proton distance bounds deduced from the NMR data sets as constraints in energy minimization computations. The BP ring is positioned in the minor groove and directed toward the 5' end of the modified strand. One face of the BP ring of (BP)G6 is stacked over the G18 and A19 sugar-phosphate backbone on the partner strand and the other face is exposed to solvent. A minimally perturbed B-DNA helix is observed for the d[T4-C5-(BP)G6-C7-T8].d[A15-G16-C17-G18-A19] segment centered about the adduct site with Watson-Crick alignment for both the (BP)G6.C17 pair and flanking G.C pairs. A widening of the minor groove at the adduct site is detected that accommodates the BP ring whose long axis makes an angle of approximately 45 degrees with the average direction of the DNA helix axis. Our study holds future promise for the characterization of other steroisomerically pure adducts of BP diol epoxides with DNA to elucidate the molecular basis of structure-activity relationships associated with the stereoisomer-dependent spectrum of mutational and carcinogenic activities.

Adduct-induced base-shifts: a mechanism by which the adducts of bulky carcinogens might induce mutations

Most carcinogens have been shown to be mutagens, and DNA adducts are formed when mutagenic/carcinogenic substances react with DNA. It is generally believed these adducts (or their derivatives) induce misreplication events that result in mutations. Many of the more potently mutagenic substances are bulky and three-dimensionally complex, such as the polycyclic aromatic hydrocarbons, aromatic amines, and aflatoxins; little is known about the mechanisms by which they induce mutations. Several theories exist and herein an additional mechanism is proposed by which bulky adducts might induce mutations at GC base pairs. Molecular modeling in conjunction with molecular mechanical calculation is used to assess if the mutagen/carcinogen moiety of the adduct might be able to shift the position of the base moiety of the adduct in such a way that misreplication events might be facilitated. This mechanism is referred to as adduct-induced base-shift, and two classes appeared possible; adduct-induced base-wobble and adduct-induced base-rotation. The latter has been proposed previously. By adduct-induced, base-wobble, the mutagen/carcinogen moiety of the adduct induces a shift in the position of the base moiety of the adduct with respect to the helix axis, which might facilitate mispairing events that are reminisent of non-Watson/Crick pairing that occurs at the wobble base of tRNA during translation. For example, in some guanine adducts, the guanine appears more thymine-like, which might facilitate G.A mispairing and thereby ultimately GC to TA transversion mutations. Adduct-induced base-rotation involves the rotation of the adducted base from the anti to the syn conformation and a variety of mispairing events might result.