Computer-generated graphic models of the N2-substituted deoxyguanosine adducts of 2-acetylaminofluorene and benzo[a]pyrene and the O6-substituted deoxyguanosine adduct of 1-naphthylamine in the DNA double helix

Molecular Dynamics and In Vitro Quantification of Safrole DNA Adducts Reveal DNA Adduct Persistence Due to Limited DNA Distortion Resulting in Inefficient Repair

The formation and repair of N²-(trans-isosafrol-3'-yl)-2'-deoxyguanosine (S-3'-N²-dG) DNA adduct derived from the spice and herbal alkenylbenzene constituent safrole were investigated. DNA adduct formation and repair were studied in vitro and using molecular dynamics (MD) simulations. DNA adduct formation was quantified using liquid chromatography-mass spectrometry (LCMS) in wild type and NER (nucleotide excision repair) deficient CHO cells and also in HepaRG cells and primary rat hepatocytes after different periods of repair following exposure to safrole or 1'-hydroxysafrole (1'-OH safrole). The slower repair of the DNA adducts found in NER deficient cells compared to that in CHO wild type cells indicates a role for NER in repair of S-3'-N²-dG DNA adducts. However, DNA repair in liver cell models appeared to be limited, with over 90% of the adducts remaining even after 24 or 48 h recovery. In our further studies, MD simulations indicated that S-3'-N²-dG adduct formation causes only subtle changes in the DNA structure, potentially explaining inefficient activation of NER. Inefficiency of NER mediated repair of S-3'-N²-dG adducts points at persistence and potential bioaccumulation of safrole DNA adducts upon daily dietary exposure.

DNA–carcinogen interaction: covalent DNA-adducts of benzo(a)pyrene 7, 8-dihydrodiol 9, 10-epoxides studied by biochemical and biophysical techniques

Exposure to various chemicals, either due to occupation or lifestyle, is considered to be a major contributing factor to tumour formation in man (Higginson, 1969; Doll and Peto, 1981). An important and prevalent class of potent carcinogenic compounds present in he environment is polycyclic aromatic hydrocarbons (PAHs), which are found in various petroleum and combustion products derived from heat and power generation and motor vehicle exhausts (Baum, 1978). Furthermore, since PAHs are generally formed by pyrolysis of organic matters such as tobacco smoking and certain procedures of food preparation, the PAH exposure to humans is extensive.

DNA–carcinogen interaction: covalent DNA-adducts of benzo(a)pyrene 7, 8-dihydrodiol 9, 10-epoxides studied by biochemical and biophysical techniques

Exposure to various chemicals, either due to occupation or lifestyle, is considered to be a major contributing factor to tumour formation in man (Higginson, 1969; Doll & Peto, 1981). An important and prevalent class of potent carcinogenic compounds present in the environment is polycyclic aromatic hydrocarbons (PAHs), which are found in various petroleum and combustion products derived from heat and power generation and motor vehicle exhausts (Baum, 1978). Furthermore, since PAHs are generally formed by pyrolysis of organic matters such as tobacco smoking and certain procedures of food preparation, the PAH exposure to humans is extensive

DNA–carcinogen interaction: covalent DNA-adducts of benzo(a)pyrene 7, 8-dihydrodiol 9, 10-epoxides studied by biochemical and biophysical techniques

Exposure to various chemicals, either due to occupation or lifestyle, is considered to be a major contributing factor to tumour formation in man (Higginson, 1969; Doll & Peto, 1981). An important and prevalent class of potent carcinogenic compounds present in the environment is polycyclic aromatic hydrocarbons (PAHs), which are found in various petroleum and combustion products derived from heat and power generation and motor vehicle exhausts (Baum, 1978). Furthermore, since PAHs are generally formed by pyrolysis of organic matters such as tobacco smoking and certain procedures of food preparation, the PAH exposure to humans is extensive.

DNA–carcinogen interaction: covalent DNA-adducts of benzo(a)pyrene 7, 8-dihydrodiol 9, 10-epoxides studied by biochemical and biophysical techniques

Exposure to various chemicals, either due to occupation or lifestyle, is considered to be a major contributing factor to tumour formation in man (Higginson, 1969; Doll and Peto, 1981). An important and prevalent class of potent carcinogeniccompounds present in the environment is polycyclic aromatic hydrocarbons(PAHs), which are found in various petroleum and combustion products derived from heat and power generation and motor vehicle exhausts (Baum, 1978). Furthermore, since PAHs are generally formed by pyrolysis of organic matters such as tobacco smoking and certain procedures of food preparation, the PAH exposure to humans is extensive.

Free-energy perturbation methods to study structure and energetics of DNA adducts: results for the major N2-dG adduct of benzo[a]pyrene in two conformations and different sequence contexts

The potent mutagen/carcinogen benzo[a]pyrene (B[a]P) is activated to (+)-anti-B[a]PDE, which induces a variety of mutations (e.g., G --> T, G --> A, etc.) via its major adduct [+ta]-B[a]P-N2-dG. One hypothesis is that adducts (such as [+ta]-B[a]P-N2-dG) induce different mutations via different conformations, probably when replicated by different lesion-bypass DNA polymerases (DNAPs). We showed that Escherichia coli DNAP V was responsible for G --> T mutations with [+ta]-B[a]P-N2-dG in a 5'-TGT sequence (Yin et al., (2004) DNA Repair 3, 323), so we wish to study conformations of this adduct/sequence context by molecular modeling. The development of a CHARMM-based molecular dynamics (MD) simulations protocol with free-energy calculations in the presence of solvent and counterions is described. A representative base-pairing and base-displaced conformation of [+ta]-B[a]P-N2-dG in the 5'-TGT sequence are used: (1) BPmi5, which has the B[a]P moiety in the minor groove pointing toward the base on the 5'-side of the adduct, and (2) Gma5, which has the B[a]P moiety stacked with the surrounding base pairs and the dG moiety displaced into the major groove. The MD output structures are reasonable when compared to known NMR structures. Changes in DNA sequence context dramatically affect the biological consequences (e.g., mutagenesis) of [+ta]-B[a]P-N2-dG. Consequently, we also developed a MD-based free-energy perturbation (FEP) protocol to study DNA sequence changes. FEP involves the gradual "fading-out" of atoms in a starting structure (A) and "fading-in" of atoms in a final structure (B), which allows a realistic assessment of the energetic and structural changes when two structures A and B are closely related. Two DNA sequence changes are described: (1) 5'-TGT --> 5'-TGG, which involves two steps [T:A --> T:C --> G:C], and (2) 5'-TGT --> 5'-TGC, which involves three steps [T:A --> T:2AP --> C:2AP --> C:G], where 2AP (2-aminopurine) is included, because T:2AP and C:2AP retain more-or-less normal pairing orientations between complementary bases. FEP is also used to evaluate the impact that a 5'-TGT to 5'-UGT sequence change might have on mutagenesis with [+ta]-B[a]P-N2-dG. In summary, we developed (1) a CHARMM-based molecular dynamics (MD) simulations protocol with free-energy calculations in the presence of solvent and counterions to study B[a]P-N2-dG adducts in DNA duplexes, and (2) a MD-based free-energy perturbation (FEP) protocol to study DNA sequence context changes around B[a]P-N2-dG adducts.

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).

Molecular modeling of four stereoisomers of the major B[a]PDE adduct (at N(2)-dG) in five cases where the structure is known from NMR studies: molecular modeling is consistent with NMR results

The potent mutagen/carcinogen benzo[a]pyrene (B[a]P) is metabolically activated to (+)-anti-B[a]PDE, which is known to induce a variety 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-N(2)-dG when bypassed during DNA replication (perhaps by different DNA polymerases). Our previous molecular modeling studies have suggested that conformational complexity might be extensive in that B[a]P-N(2)-dG adducts appeared capable of adopting at least sixteen potential conformational classes in ds-DNA [e.g., Kozack and Loechler (1999) Carcinogenesis 21, 1953], although only eight seemed likely to be relevant to base substitution mutagenesis. Such molecular modeling studies are only likely to be valuable for the interpretation of mutagenesis results if global minimum energy conformations for adducts are found and if the differences in the energies of these different conformations can be computed reasonably accurately. One approach to assessing the reliability of our molecular modeling techniques is considered herein. Using a five-step molecular modeling protocol, which importantly included a molecular dynamics version of simulated annealing, eight conformations are studied in each of five cases. (The five cases are listed below, and were chosen because in each case the preferred solution conformation is known from a NMR study.) Of the eight conformations studied, the one computed to be lowest in energy is the same conformation as the one observed by NMR in four of the five cases: 5'-CGC sequence with [+ta]-, [-ta]-, and [+ca]-B[a]P-N(2)-dG, and 5'-TGC sequence with [+ta]-B[a]P-N(2)-dG. In the fifth case (5'-CGC sequence with [-ca]-B[a]P-N(2)-dG), the known NMR conformation is computed to be second lowest in energy, but it is within approximately 1.7 kcal of the computed lowest energy conformation. These results suggest that molecular modeling is surprisingly accurate in computing lowest energy conformations and that it should be useful in assessing the relative energies of different conformations. This is especially important given that currently molecular modeling is the only means available to study the energetics of minor conformations of DNA adducts.

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.

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.

Molecular mechanisms of mutagenesis by aromatic amines and amides

The mutagenic properties of 2-acetylaminofluorene-derived DNA adducts, including N-(deoxyguanosin-8-yl)-2-acetylaminofluorene, N-(deoxyguanosin-8-yl)-2-aminofluorene, N-(deoxyguanosin-N2-yl)-2-acetylaminofluorene, and several minor oxidation products have been explored, using site-specific techniques. Oligodeoxynucleotides containing a single AAF-derived DNA adduct were prepared by postsynthetic modification and used as templates in primer extension reactions catalyzed by bacterial and mammalian DNA polymerases. Base substitutions and deletions occurring during DNA synthesis were quantified. dG-C8-AAF promoted one- and two-base deletions and small amounts of incorporation of dCMP, dAMP, and/or dTMP opposite the lesion in reactions catalyzed by the 3'-->5' exonuclease-free Klenow fragment of DNA polymerase 1 (exo-) and polymerase alpha. dG-C8-AF did not miscode in reactions catalyzed by exo-; however, base misincorporation and deletions were observed in reactions with pol alpha. dG-N2-AAF promoted small amounts of dAMP incorporation in reactions catalyzed by exo-. The miscoding potential of minor oxidation products of dG-C8-AF was much higher than that of other adducts. Steady-state kinetics were used to measure frequencies of nucleotide insertion opposite the lesion and chain extension from the 3' terminus. Kinetic data were consistent with the results of primer extension studies. A mutation 'hot spot' was constructed and the influence of sequence context on the frequency of deletions generated by dG-C8-AAF was explored systematically in reactions catalyzed by exo-. Based on our results with aminofluorene DNA adducts, we propose a general mechanism for frameshift deletion mutagenesis. Site-specific methods also were used to establish the mutagenic potential of AAF-derived DNA adducts in mammalian cells. dG-C8-AAF and dG-C8-AF exhibited similar mutagenic specificities, predicting the occurrence of G-->T transversions and G-->A transitions in mammalian cells.

Molecular characterization of mutation and comparison of mutation profiles in the hprt gene of Chinese hamster ovary cells treated with benzo[a]pyrene trans-7,8-diol-anti-9,10-epoxide, 1-nitrobenzo[a]pyrene trans-7,8-diol-anti-9,10-epoxide, and 3-nitrobenzo[a]pyrene trans-7,8- diol-anti-9,10-epoxide

Both 1- and 3-nitrobenzo[a]pyrene (nitro-BaP) are environmental contaminants, potent mutagens in Salmonella, and moderate mutagens in Chinese hamster ovary (CHO) cells. The mutagenicity of their oxidized metabolites,trans-7,8-dihydroxy-anti-9,10-epoxy-7,8,9,10-epoxy -7,8,9,10-tetrahydro-1-nitrobenzo[a]pyrene (1-nitro-BaP-DE) and trans-7,8-dihydroxy-anti-9,10-epoxy-7,8,9,10-tetrahydro-3-nitrobenzo[a]- pyrene (3-nitro-BaPDE), together with trans-7,8-dihydroxy-anti-9, 10-ep- oxy-7,8,9,10-tetrahydrobenzo[a]pyrene (BaP-DE), was determined in CHO-K1 cells, and the resulting mutations at the hprt locus were characterized by polymerase chain reaction (PCR) amplification of reverse-transcribed hprt mRNA, followed by DNA sequence analysis. The mutant frequencies, in mutants/10(6) clonable cells, at 30 and 100 ng/ml, were BaP-DE, 248 and 456; 1-nitro-BaP-DE, 68 and 260; 3-nitro-BaP-DE, 81 and 232, respectively. In general, the three diolepoxides exhibited similar mutational spectra: 1) 64% (23/36 sequenced mutants) of BaP-DE, 53% (19/36) of 1-nitro-BaP-DE, and 64% (23/36) of 3-nitro-BaP-DE mutants resulted from simple base pair substitution, with the predominant mutation being G-->T transversion; 2) 90%, 100%, and 100% of mutations at G:C had the mutated dG on the nontranscribed DNA strand; and 3) about one quarter of the mutants produced by each mutagen had one or more PCR products with partial or complete exon deletions. The mutagens induced few frameshifts or complex mutations. Among the differences in mutational specificity for the three diolepoxides, the proportion of substituted dGs with 3' purines was significant (P < 0.05) for BaP-DE (16/19, 84%) and 3-nitro-BaP-DE (17/20, 85%), but not significant for 1-nitro-BaP-DE-induced mutants (11/17, 65%, P > 0.05). Also, high proportions of BaP-DE and 3-nitro-BaP-DE base pair substitutions at G:C occurred in DNA sequence contexts of 5'-GG-3', 5'-GGA-3', and 5'-TGGA-3', while the proportions of 1-nitro-BaP-DE mutants in these contexts were often lower. The results indicate that nitro substitution at C1 or C3 of BaP-DE reduces mutational potency in CHO cells and appears to have only subtle effects upon the mutational pattern in the hprt gene.

Genetic toxicity of 2-acetylaminofluorene, 2-aminofluorene and some of their metabolites and model metabolites

2-Acetylaminofluorene and 2-aminofluorene are among the most intensively studied of all chemical mutagens and carcinogens. Fundamental research findings concerning the metabolism of 2-acetylaminofluorene to electrophilic derivatives, the interaction of these derivatives with DNA, and the carcinogenic and mutagenic responses that are associated with the resulting DNA damage have formed the foundation upon which much of genetic toxicity testing is based. The parent compounds and their proximate and ultimate mutagenic and carcinogenic derivatives have been evaluated in a variety of prokaryotic and eukaryotic assays for mutagenesis and DNA damage. The reactive derivatives are active in virtually all systems, while 2-acetylaminofluorene and 2-aminofluorene are active in most systems that provide adequate metabolic activation. Knowledge of the structures of the DNA adducts formed by 2-acetylaminofluorene and 2-aminofluorene, the effects of the adducts on DNA conformation and synthesis, adduct distribution in tissues, cells and DNA, and adduct repair have been used to develop hypotheses to understand the genotoxic and carcinogenic effects of these compounds. Molecular analysis of mutations produced in cell-free, bacterial, in vitro mammalian, and intact animal systems have recently been used to extend these hypotheses.

Preparation and characterization in solution of oligonucleotides alkylated by activated carcinogenic polycyclic aromatic hydrocarbons

The effects of aralkylation of selected oligonucleotides by a bulky chemical carcinogen, 7,12-dimethylbenz(a)anthracene (after activation) have been studied. The aralkylation involves the base adenine, designated A* at the modification site, in the center of synthetic heptameric, nonameric and pentadecameric oligonucleotides; complementary strands lacking any modification were also synthesized. The products were studied by UV melting curves and CD spectral techniques. Duplex formation was modified by such aralkylation of a central base in the oligomers. The extent of duplex formation was found to depend on chain length as follows: no evidence was found for duplex formation of the heptamer d(GTCA*GAC) + d(GTCTGAC); the nonamer, d(GTGCA*ATCC) + d(GGATTGCAC), appears to form a duplex at high salt concentrations and reduced temperature; the pentadecamer, d(CCGCT-GCGA*TCCGGC) + d(GCCGGATCGCAGCGG), forms a duplex at low salt concentration and room temperature, but its melting temperature is lower than that of the nonalkylated parent system. CD-spectra for the duplexes formed by the nonamer or pentadecamer are indicative of a right-handed helical conformations. On phosphordiesterase digestion it appears that the aralkylated adenine and the base on its 5'-side act as "stops" for enzymatic digestion from either direction. We suggest, from model building, that this inhibition of phosphodiesterase activity is the result of the steric bulk and disposition of the polycyclic aromatic hydrocarbon. We further suggest that unusual base pairing (mismatching), such as A...A, which would lead to an AT transversion, may be favored by the bulkiness of the aromatic group.

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

The minor groove of undistorted A-DNA provides a good binding site for planar, hydrophobic moieties such as unmetabolized polycyclic aromatic hydrocarbons (PAHs), and the base pairs at the ends of short oligodeoxynucleotide helices. It also accommodates the chief adduct derived from the metabolically activated form of the carcinogen benzo[a]pyrene. B-DNA lacks such a site. Computerized models have been generated for the major (N2-guanine-linked) adducts formed at this site by both + and - enantiomers of anti-benzo[a]pyrene-7,8-dihydrodiol-9,10-epoxide (anti-BPDE) with poly(dG).poly(dC) in the A-DNA conformation. The BPDE adducts lie in the shallow, relatively hydrophobic minor groove of the A-DNA after empirical potential energy minimization using the program AMBER. We term this binding mode "side-stacking." The side-stacked + anti-BPDE may constitute the chief carcinogenic lesion derived from benzo[a]pyrene.

Binding of enantiomers of trans-7,8-dihydroxy-anti-9,10-epoxy-7,8,9,10-tetrahydro-benzo[a]pyrene to polynucleotides

DNA covalent binding studies with enantiomers of trans-7,8-dihydroxy- anti-9,10-epoxy-7,8,9,10-tetrahydro-benzo[a]pyrene (anti-BPDE) have been carried out by means of spectroscopic techniques (UV, CD, and fluorescence). Synthetic polynucleotides are employed to investigate binding differences between the G.C and A.T base pairs and to elucidate the bases for the stereoselective covalent binding of DNA toward anti-BPDE. The results indicate that of all the polynucleotides studied, only poly(dA-dT).poly(dA-dT) exhibits predominant intercalative covalent binding towards (+)-anti-BPDE and suffers the least covalent modification. Only minor intercalative covalent contributions are found in alternating polymer poly(dA-dC).poly(dG-dT). These observations parallel the DNA physical binding results of anti-BPDE and its hydrolysis products. They support the hypothesis that intercalative covalent adducts derive from intercalative physical binding while the external covalent adducts derive from external bimolecular associations. In contrast to the A.T polymers, the guanine containing polymers exhibit pronounced reduction in covalent modification by (-)-anti-BPDE. The intercalative covalent binding mode becomes relatively more important in the adducts formed by the (-) enantiomer as a consequence of decreased external guanine binding. These findings are consistent with the guanine specificity, stereoselective covalent binding at dG, the absence of stereoselectivity at dA for anti-BPDE, and the enhanced binding heterogeneity for the (-) enantiomer as found in the native DNA studies. The possible sequence and/or conformational dependence of such stereoselective covalent binding is indicated by the opposite pyrenyl CD sign exhibited by (+)-anti-BPDE bound to polynucleotides with pyrimidine on one strand and purine on another vs. that bound to polymers containing alternating purine-pyrimidine sequences.

Electrofluorescence study of polycyclic hydrocarbon diol-epoxide binding to DNA

By the use of a novel electrofluorescence method, estimates have been made of the geometry of binding to DNA of racemic mixtures of the anti-diol-epoxide derivatives of three polycyclic hydrocarbon carcinogens. These anti-configurations bind in a manner consistent with the planar diol-epoxide ring's being inclined at approximately 50 degrees to the DNA axis. This is true for the derivatives of benzo(a)pyrene, benz(a)anthracene and 3-methylcholanthrene. This binding is thus different from the regular intercalative interaction associated with the native hydrocarbons. As the (+ anti)-diol-epoxides are thought to be the initiatory compounds for carcinogenesis, the common binding characteristics for the three hydrocarbons may be significant in understanding the molecular interactions precursive to cancer.

Catalysis of carcinogen-detoxification by DNA: comparison of enantiomeric diol epoxides

The interactions of the (+)- and (-)-enantiomers of 7r,8t-dihydroxy-9t,10t-oxy-7,8,9,10-tetrahydrobenzo[a]pyrene (BPDE-I) with purified DNA have been studied in vitro. These compounds are formed by cellular metabolism of the potent environmental carcinogen benzo[a]pyrene, and the (+)-enantiomer is thought to be the ultimate carcinogenic metabolite. Non-covalent, intercalative binding was measured spectrophotometrically, hydrolysis was measured spectrofluorometrically and covalent binding was detected by liquid scintillation counting. No significant differences were found in the association constants for intercalative binding or in the ability of DNA to catalyse the hydrolysis of the two enantiomers. Covalent DNA binding was 4.5-fold higher for the (+)-enantiomer. When DNA was pretreated with a molar equivalent of the (-)-enantiomer, its subsequent ability to enhance the rate of BPDE-I hydrolysis and to bind covalently to (+)-BPDE-I was unimpaired. This suggests that the participation of the DNA in the hydrolysis reaction does not alter the DNA and therefore that the rate-enhancement is true catalysis.

Effects of adduct formation on the biological activity of single- and double-stranded øX174 DNA, modified by N-acetoxy-N-acetyl-2-aminofluorene

In order to establish a good quantitative relationship between the number of acetylaminofluorene adducts and the extent of inactivation of DNA, single-stranded (ss) øX174 DNA and øX174 RF DNA were modified to various extents with 3H-labelled N-acetoxy-N-acetyl-2-aminofluorene (N-AcO-AAF) and subsequently transfected to Escherichia coli spheroplasts having different repair capabilities. Exponential survival curves were obtained. In the case of ssDNA about one adduct per molecule appears to be lethal. On the other hand only 1 out of 10.2 adducts is found to inactivate RF DNA if tested on wild-type E. coli. However, when assayed on strains deficient in excision repair 1 out of 2.3 adducts leads to inactivation of RF DNA. RecA-dependent postreplication repair only has little influence on these figures. Product analysis of the modified DNAs shows that in RF DNA at least 76% of the interaction products is N-(deoxyguanosin-8-yl)-N-acetyl-2-aminofluorene (dGuo-C8-AAF) and at least 6% and at most 12% is 3-(deoxyguanosin-N2-yl)-N-acetyl-2-aminofluorene (dGuo-N2-AAF). In ssDNA only dGuo-C8-AAF is formed. No apurinic sites could be detected in the modified DNAs. From these results it can be concluded that in RF DNA most of the dGuo-C8-AAF is removed by excision repair. The remaining damage, consisting probably both of dGuo-N2-AAF and unexcised dGuo-C8-AAF, inactivates RF DNA. Inactivation can be explained by a model which shows that only damage in the minus strand of RF DNA inhibits replication and/or transcription.

Interactions of molecules with nucleic acids. X. Covalent intercalative binding of the carcinogenic BPDE I(+) to kinked DNA

A theoretical model is proposed for the covalent binding of (+) 7 beta,8 alpha-dihydroxy-9 alpha,10 alpha-epoxy-7,8,9,10- tetrahydrobenzo[a]pyrene denoted by BPDE I(+), to N2 on guanine. The DNA must kink a minimum of 39 degrees to allow proper hybrid configurations about the C10 and N2 atoms involved in bond formation and to allow stacking of the pyrene moiety with the non-bonded adjacent base pair. Conservative (same sugar puckers and glycosidic angles as in B-DNA) and non-conservative (alternating sugar puckers as in intercalation sites) conformations are found and they are proposed structures in pathways connecting B-DNA, an intercalation site, and a kink site in the formation of a covalently intercalative bound adduct of BPDE I(+) to N2 on guanine. Stereographic projections are presented for (3') and (5') binding in the DNA. Experimental data for bending of DNA by BPDE, orientation of BPDE in DNA and unwinding of superhelical DNA is explained. The structure of a covalent intercalative complex is predicted to result from the reaction. Also, an anti----syn transition of guanine results in a structure which allows the DNA to resume its overall B-form. The only change is that guanine has been rotated by 200 degrees about its glycosidic bond so that the BPDE I(+) is bound in the major groove. The latter step may allow the DNA to be stored with an adduct which may produce an error in the genetic code.

Computer modelling studies of the covalent interactions between DNA and the enantiomers of anti-7,8-diol,9,10-epoxy-benzo[a]pyrene

The molecular structures of adducts between the + and - enantiomers of 7,8-diol 9,10-epoxy benzo[a]pyrene and a double-stranded model for DNA, have been examined by empirical energy calculations. Low-energy structures were only obtained for A form, and not B form DNA. Both + and - adducts are of approximately equal energy. Some structural differences in the orientation of the BP chromophore in the two adducts were found.

Binding of pyrene to DNA, base sequence specificity and its implication

Solubilization as well as spectral studies of pyrene in natural DNA and synthetic deoxypolynucleotide solutions at neutral pH reveal at least two binding modes. Sites I are predominant in native DNA and in poly(dA-dT): poly(dA-dT) whereas sites II are found with denatured DNA and other polynucleotides such as poly(dA):poly(dT) and three different types of guanine containing copolymers which solubilize pyrene to a lesser extent. Spectral comparison with the covalent adducts of trans-7,8-dihydroxy-anti-9,10-epoxy-7,8,9,10- tetrahydro-benzo(a)pyrene (anti-BPDE) and the physical complexes of its tetraols lead to the suggestion of a base sequence specific binding model for this carcinogenic metabolite to account for the puzzling fact that although its physical binding is predominantly intercalative, the covalent adducts appear not to be intercalated. It is speculated that in neutral solutions, intercalation may have little, if any, to do with the chemical lesion of this metabolite to the guanine base of the DNA and may, on the contrary, provide an efficient pathway for detoxification.

AAF linked to the guanine amino group: a B-Z junction

Minimized conformational potential energy calculations have been performed for AAF linked to dCpdG at the guanine amino group. This is a model for the minor AAF adduct observed in DNA, whose conformational influence has been difficult to ascertain. A global minimum energy conformation was computed with torsion angles like those of the dCpdG residue of Z-DNA. This conformation was incorporated into a larger polymer model at a B-Z junction, with the carcinogen residing in the groove in the Z direction. Local minimum energy conformations of the B type were also computed. In addition, two minima were found with fluorenecytidine stacking. These results suggest that existing B-Z junctions may be vulnerable to modification by AAF at the guanine amino group, or that such junctions may be induced by the carcinogen if the base sequence is appropriate. Otherwise the carcinogen can be located in the minor groove of the B helix (5, 10, 11) or covalently intercalated (13-15).

The formation of acetylaminofluorene adducts in poly(dC-dG) and poly(dA-dT) on reaction with N-acetoxy-2-acetylaminofluorene and the effect of such modification upon the polymers as templates for DNA polymerases

N-Acetoxy-2-acetylaminofluorene (AcO-AAF) reacts with the alternating DNA-like polynucleotides poly(dC-dG) and poly(dA-dT) in vitro to give adducts of the guanine and adenine bases similar to those reported to be formed in DNA. A previously unobserved guanine adduct was detected in the poly(dC-dG). Using a double-labelled [U-14C-dG, 8-3H-G]-poly(dC-dG) we show that this adduct does not involve the 7- or 8-positions of the guanine. Similarly a thymine adduct of unknown structure was observed in poly(dA-dT). Modification of the polymers with AcO-AAF inhibits their capacity to act as templates for Escherichia coli DNA polymerase I and mammalian DNA polymerase alpha although the binding of the polymerases to the polynucleotides is unaffected. Such modification also leads to an increase in the levels of non-complementary nucleotides incorporated into newly synthesised DNA.

Mutagenicity and genotoxicity of nitroarenes. All nitro-containing chemicals were not created equal

The nitrated polycyclic aromatic hydrocarbons constitute a group of chemicals of environmental concern which display a broad spectrum of mutagenic, genotoxic and carcinogenic properties. Some members of the group are the most potent direct-acting bacterial mutagens while others exhibit low levels of potencies which require metabolic activation mixtures. Bacterial mutagenicity is dependent upon reduction of the nitro function. In mammalian cell systems the genetic and genotoxic effects of these nitrated chemicals include the induction of unscheduled DNA synthesis, sister-chromatid exchanges, chromosomal aberrations, gene mutations and cell transformation. The qualitative as well as quantitative expression of these effects is dependent upon the species and tissue of origin as well as culture history of the cell which in turn determine their enzymic capabilities and the conversion of these nitroarenes to ultimate mutagens and genotoxicants. In eukaryotic cells the following bioactivation pathways have been recognized: (a) reduction of the nitro moiety, (b) ring oxidation (the nature of which is influenced by the nitro function) followed by reduction of the nitro group, and (c) ring oxidation without concomitant reduction of the nitro moiety.

A comparison of the carcinogen-DNA adducts formed in rat liver in vivo after administration of single or multiple doses of N-methyl-4-aminoazobenzene

N-Methyl-4-aminoazobenzene (MAB) is believed to be metabolized in the liver to an electrophilic N-sulfonyloxy ester which binds covalently to cellular macromolecules, resulting in the induction of hepatic neoplasia. Previous in vivo studies in the rat detected only two hepatic MAB-DNA adducts, 3-(deoxyguanosin-N2-yl)-MAB(N2-dG) and N-(deoxyguanosin-8-yl)-MAB(C8-dG), which respectively accounted for 25% and 70% of the total MAB bound to DNA at 8 h after a single dose of the carcinogen. Subsequently, the C8-dG adduct was shown to be rapidly lost from the DNA while the N2-dG adduct was a persistent lesion. Since a single dose of MAB is not sufficient for complete carcinogenic activity, we sought to identify the MAB-DNA adducts present in rat liver after multiple oral doses of [3H]MAB. The MAB was administered by intubation at a level of 0.2 mmol/kg for 1, 3 or 4 doses and animals were sacrificed at 8 h after the last dose. Hepatic DNA was isolated by extraction and hydroxylapatite chromatography and was enzymatically hydrolyzed to MAB-mononucleoside adducts, which were quantitated by high pressure liquid chromatography (HPLC). After 3 doses, N2-dG, C8-dG, and an unknown adduct were detected. By 4 doses, these accounted for 51%, 25% and 23% of the total adducts. This data is consistent with rapid removal of the C8-dG derivative and the relative persistence of the N2-dG and the unknown adduct. The latter was shown to exhibit chromatographic and pH-dependent solvent partitioning properties that were identical to a product also present in DNA treated with the synthetic ultimate carcinogen, N-benzoyloxy-MAB. Analysis of this adduct by field desorption mass spectrometry (M+ = 460) and, after perdeuteromethylation, by electron impact mass spectrometry (M+ = 528; M-N(CH3)(CD3) = 481) indicated the structure to be a deoxyadenosin-N6-yl derivative substituted through an aromatic ring of MAB. Further analysis by 270 MHz 1H-NMR spectroscopy allowed complete assignment of the MAB and adenyl resonances and was uniquely consistent with a 3-(deoxyadenosin-N6-yl)-MAB structure. Since this persistent adduct is potentially mutagenic due to possible tautomeric equilibria between the N6-amino and N6-imino structures, it may represent an initiating lesion in MAB hepatocarcinogenesis.

Defective excision repair in a mutant of Micrococcus radiodurans hypermutable by some monofunctional alkylating agents

The lethal and mutagenic effects of methyl methanesulphonate (MMS), ethyl methanesulphonate (EMS), and N-methyl-N'-nitro-N-nitrosoguanidine (MNNG) can be dissociated in a mitomycin C (MTC)-sensitive mutant, strain 302, of Micrococcus radiodurans. As regards lethality 302 is extremely sensitive, compared with the wild type, to MTC and decarbamoyl MTC (DCMTC), slightly sensitive to EMS, MNNG, nitrous acid, 7-bromomethylbenz[alpha]anthracene (BrMBA), and N-acetoxy-N-2-acetylaminofluorene (AAAF), and resistant to MMS, hydroxylamine, and ICR 191G. As regards mutability it is, compared to the wild type, very sensitive to MMS, EMS, and MNNG, and slightly sensitive to hydroxylamine and nitrous acid but not to any other agent examined. Alkaline sucrose gradient studies indicate the 302 does not incise DNA containing BrMBA adducts, although it does incise DNA damaged by AAAF but probably not to the same extent as wild type. We put forward the hypothesis that the hypermutability of 302 is due to the non-removal of bases or nucleotides, modified in exocyclic positions, which have altered base-pairing capabilities, while lethality results from the non-removal of bases or nucleotides, also modified in exocyclic positions, which no longer form hydrogen-bonded base pairs.

Characterization of DNA adducts of the carcinogen N-methyl-4-aminoazobenzene in vitro and in vivo

Since the susceptibility of specific tissues to tumor formation has been correlated with the persistence of DNA-carcinogen adducts, the identity and persistence of DNA adducts formed from the hepatocarcinogen N-methyl-4-aminoazobenzene (MAB) has been determined. The synthetic ultimate carcinogen N-benzoyloxy-N-methyl-4-aminoazobenzene (N-BxO-MAB) was reacted in vitro with either calf thymus or rat liver DNA to yield approx. 1 bound residue per 1000 nucleotides. After enzymatic hydrolysis of the DNA and high pressure liquid chromatographic analysis, at least six MAB adducts were detected. Two of the products cochromatographed with MAB-DNA adducts formed in rat liver in vivo following oral administration of the precarcinogen MAB. These two adducts were identified by mass, UV and nuclear magnetic resonance (NMR) spectroscopy as N-(deoxyguanosin-8-yl)- and 3-(deoxyguanosin-N2-yl)-MAB. The former adduct was initially the predominant product in vivo, but it could not be detected 7 days following treatment. The latter adduct remained at a constant level for 14 days and therefore appears to be a persistent lesion.

Molecular models of induced DNA premutational damage and mutational pathways for the carcinogen 4-nitroquinoline 1-oxide and its metabolites

The carcinogen 4-nitroquinoline 1-oxide (4NQO) and its metabolites undergo intercalative or covalent binding with DNA. Recent evidence indicates that the latter binding pattern is probably facilitated by an initial weaker intercalative interaction that can align potentially reactive sites on a 4NQO-metabolite and adjacent stacked bases. In the present study, we have proposed numerous possible covalent reaction products between 4NQO and its metabolites with DNA mini-helices based on chemical properties and key 'short-contacts' after energy-minimization in 21 different intercalative-like complexes. It is known from numerous experimental studies that 90% of the quinoline-bound DNAs in vivo involve guanine with the remaining 10% apparently involving adenine residues. The results of the present study suggest that this trend is not due to the greater affinity of the quinolines for guanine, but instead results from secondary processes involving the preferential formation of apurinic sites at aralkyl-adenine residues over that of aralkyl-guanine residues. In addition, observed mutational patterns can be rationalized in terms of the proposed reaction-products. The role of DNA repair mechanisms in the removal and correction of the different proposed reaction products are discussed. The binding pattern of several other aromatic carcinogens are similar to those depicted in the present work for the 4NQO-metabolites; hence the present study may be of some general significance.

DNA alkylation and unwinding induced by benzo[a]pyrene diol epoxide: modulation by ionic strength and superhelicity

Superhelical and partially relaxed DNAs of simian virus 40 were allowed to react in vitro with (+/-)-7 beta,8 alpha-dihydroxy-9 alpha,10 alpha-epoxy-7,8,9,10-tetrahydrobenzo[a]pyrene (BaP diol epoxide). The modified DNA contained N2 guanine and N6 adenine hydrocarbon adducts in the ratio 86:14. Superhelical simian virus 40 DNA was approximately 6% more susceptible to modification than was partially relaxed viral DNA. Counterions inhibited DNA alkylation by up to 90%, Mg2+ being 50-fold more effective than Na+. The sensitivity of covalent binding to helix stability is consistent with a reaction complex in which BaP diol epoxide is intercalated. The superhelical density of the modified DNA substrates was determined electrophoretically relative to partially relaxed standards, and an unwinding angle for the hydrocarbon adducts was calculated. The angle was dependent upon the superhelicity of the DNA molecule and ranged from 330 degrees to 30 degrees. These data indicate that the modified base pairs are disrupted and, in the presence of torsional strain, act as centers for the further denaturation of up to eight adjacent base pairs. In the absence of such strain the alkylation sites have an ordered structure, with the attached hydrocarbon probably oriented in the minor or major groove of the helix.

Modification of deoxyribonucleic acid by a diol epoxide of benzo[a]pyrene. Relation to deoxyribonucleic acid structure and conformation and effects on transfectional activity

The effects of secondary structure on DNA modification by (+/-)-7 beta, 9 alpha-dihydroxy-9 alpha, 10 alpha-epoxy-7,8,9,10-tetrahydrobenzol[a]pyrene [(+/-)BPDE I] were investigated. No differences in the total extent of (+/-) BPDE I binding to double- and single-stranded calf thymus DNA were found. High-performance liquid chromatography (LC) of the nucleoside adducts obtained from hydrolysates of native and denatured calf thymus, as well as from superhelical and linear plasmid DNA, indicated that in all cases the major adduct (60--80% of total adducts) was formed by reaction of the (+) enantiomer of BPDE I with the N-2 position of dG residues in the DNA. A minor adduct formed from the reaction of the (-) enantiomer with dG residues was also detected and was present in greater amounts in denautred DNA than in native DNA. Small amounts of BPDE I--dA and BPDE I--dC adducts were also detected in both the single- and double-stranded DNAs. Restriction enzyme analysis of BPDE I modified SV40 and phage lambda DNA provided evidence that the modification of DNA by this carcinogen is fairly random with respect to nucleotide sequence. Partial hydrolysis of modified plasmid DNA by the single-strand-specific S1 nuclease and LC analysis of the nucleoside adducts in the digested and undigested fractions of the DNA revealed no preferential excision by the S1 nuclease of the different BPDE I--deoxynucleoside adducts. Functional changes in BPDE I modified DNA were demonstrated. With increasing extents of modification, there was a decrease in the ability of plasmid DNA to transfect a receptive Escherichia coli strain to antibiotic resistance.

Rapid isolation of carcinogen-bound DNA and RNA by hydroxyapatite chromatography

Carcinogen-bound DNA and RNA are conveniently isolated by solvent extraction and hydroxyapatite (HAP) chromatography. Tissue is suspended in 8 M urea-0.24 M sodium phosphate-1% sodium dodecyl sulfate-10 mM EDTA, pH 6.8 (MUP-SDS-EDTA) and extracted with chloroform-isoamyl alcohol-phenol (24:1:25; CIP) to remove protein. RNA and DNA are separated by passing the aqueous solution through an HAP column; RNA is eluted with MUP, DNA with 0.48 M sodium phosphate, pH 6.8. Examples presented are: (1) calf thymus DNA that has been reacted with N-acetoxy-2-acetylaminofluorene (N-OAc-AAF), (2) isolated intact rat hepatocytes incubated with N-hydroxy-AAF and (3) livers from Sprague-Dawley rats treated with N-hydroxy-AAF.