The food mutagen 2-amino-3,8-dimethylimidazo[4,5-f]quinoxaline: a conformational analysis of its major DNA adduct and comparison with the 2-amino-3-methylimidazo[4,5-f]quinoline adduct

The heterocyclic amine 2-amino-3,8-dimethylimidazo[4,5-f]quinoxaline (MeIQx) is one of a group of heterocyclic amine carcinogens that exists in cooked meat and fish. It causes mutations in bacterial and mammalian assays and induces tumors in mammals. MeIQx is converted within cells to a reactive derivative which forms a major covalent adduct at carbon-8 of guanine in DNA. This adduct may alter the DNA conformation at critical stages of the replicative process, and cause mutations which initiate the carcinogenic process. Atomic resolution structures of the MeIQx-damaged DNA are not yet available experimentally. We have carried out an extensive molecular mechanics/energy minimization search to locate feasible structures for the major MeIQx adduct in DNA, using the sequence d(5'-C1-G2-C3-G4[IQ]-C5-G6-C7-3').d(5'-G8-C9-G10-C11-G12-C13-G14-3') with MeIQx modification at G4. We have created 1152 starting conformations which uniformly sampled each of the three flexible torsion angles that govern the MeIQx-DNA orientation at 15 degrees intervals, and minimized their energy. A mixture of conformations was generated, which were separated into families according to the position of the ring system of the carcinogenic amine: major groove, minor groove, and base-displaced-intercalated. While a generally similar mixture had been generated previously for the related carcinogen 2-amino-3-methylimidazo[4,5-f]quinoline (IQ) [Wu, X., et al. (1999) Chem. Res. Toxicol. 12, 895-905], differences were found which could be rationalized in terms of the additional methyl group in the MeIQx.     

No effects of chlorophyllin on IQ (2-amino-3-methylimidazo[4,5-f]-quinoline)-genotoxicity and -DNA adduct formation in Drosophila

Previously we demonstrated that chlorophyllin suppressed the genotoxicities of many carcinogens. However, the genotoxicity of IQ (2-amino-3-methylimidazo[4,5-f]quinoline), a carcinogenic heterocyclic amine, was not suppressed in Drosophila. On the contrary, it has been reported that chrolophyllin suppressed the genotoxicity of IQ in rodents, rainbow trout and Salmonella. We demonstrated that the chlorophyllin-induced suppression of MeIQx (2-amino-3,8-dimethylimidazo[4,5-f]quinoxaline)-genotoxicity was associated with a decrease in MeIQx-DNA adduct formation in Drosophila larval DNA. MeIQx represents another type of heterocyclic amine and is similar to IQ in structure. In this study we utilized (32)P-postlabeling to examine whether chlorophyllin reduced IQ-DNA adduct formation in Drosophila DNA in the same way as MeIQx. The results revealed that the formation of IQ-DNA adducts was unaffected by treatment with chlorophyllin. This was consistent with the absence of any inhibitory effect on genotoxicity as observed in the Drosophila repair test. These results suggest that IQ-behavior in Drosophila is not affected by chlorophyllin, indicating that the process of IQ-DNA adduct formation followed by expression of genotoxicity in Drosophila may be different from that in other organisms.     

DNA cross-linking by azinomycin B: Monte Carlo simulations in the evaluation of sequence selectivity

A new set of charges specifically developed for biologically relevant N7-alkylated purine adducts have been implemented in the AMBER force field of the MacroModel package and applied to the conformational search of azinomycin B-DNA interactions. To perform a sequence dependent reactivity relationship study, four DNA triplets known to interact differently with the drug, 5'-GCT-3', 5'-GCC-3', 5'-GTC-3', and 5'-GTT-3', have been modeled in B-form and intercalative conformations. Monte Carlo simulations of all possible monoadducts and intercalative complexes have been carried out and analyzed using a filtering criterion that estimates the probability of covalent bond formation and covalent cross-linking. We observed a good correlation between existing experimental data and our computational estimations that validate the approach. The comparison of the conformational properties of the drug-DNA monoadducts and complexes confirms the most probable mechanism of action involving an initial aziridine and subsequent epoxide alkylation. The different hydrogen bond network in the monoadducts and in the intercalative complexes between the drug and the three base-pair receptor is the primary reason for the different cross-linking reactivity. In addition, steric hindrance of the major groove exposed methyl group of central thymine-based triplets plays an important role in the lack of the reactivity of these sequences. Synthetic work on the azinomycins and the information coming from this computational study will be important for the design of more potent or DNA sequence-selective agents based on the azinomycin skeleton.     

Conformations of complexes between pyrrolo[1,4]benzodiazepines and DNA segments

The molecular mechanics program AMBER, assisted by CHEMLAB II, was used to model the covalent and noncovalent binding of anthramycin, tomaymycin, and neothramycin A to the hexanucleotide conformation. Structures covalently bonded at N2 of guanine gave excellent fits when placed in either direction in the minor groove. However, energy analysis showed a preference for the direction wherein the side chain points toward the 5' end of the covalently bound strand. This preference agrees with published NMR studies. Noncovalent binding of anthramycin in the minor groove near guanine gave good fits with almost no distortion in the helix, and the reactive center of the ligand was close enough to N2 for subsequent covalent bond formation. Anthramycin also gave a good noncovalent complex near adenine in the minor groove, but binding in the major groove had decreased dispersion attractions. Binding of tomaymycin was similar to that of anthramycin, although the smaller size of tomaymycin resulted in less binding energy. Neothramycin noncovalent binding was characterized by strong electrostatic interactions, partly involving the 3-OH group, and by part of the molecule lying outside the minor groove. AMBER was used for the exploratory design of an anthramycin analogue that theoretically would bind as well as anthramycin but not cause cardiotoxicity. A related study involving anthramycin, tomaymycin, and the pentanucleotide duplex d(AAGAA/TTCTT) was undertaken to evaluate further the ability of AMBER to predict sequence specificity. It indicated a preferred direction of binding toward 5' in the minor groove of the duplex, but rather weak interaction with the noncovalently bound strand. This prediction agreed with experiments on tomaymycin that showed separation of the duplex and alignment of the drug toward the 5' end of the covalently bound strand.     

Mutagenicity of the 1-nitropyrene-DNA adduct N-(deoxyguanosin-8-yl)-1-aminopyrene in mammalian cells

The mutagenesis of the major DNA adduct N-(deoxyguanosin-8-yl)-1-aminopyrene (C8-AP-dG) formed by 1-nitropyrene was compared with the analogous C8-dG adducts of 2-aminofluorene (AF) and N-acetyl-2-aminofluorene (AAF) in simian kidney (COS-7) cells. The DNA sequence chosen for this comparison contained 5'-CCATC GCTACC-3' that has been used for solution NMR investigations. The structural and conformational differences among these lesions are well-established [Patel, D. J., Mao, B., Gu, Z., Hingerty, B. E., Gorin, A., Basu, A. K., and Broyde,S. (1998) NMR solution structures of covalent aromatic amine-DNA adducts and their mutagenic relevance. Chem. Res. Toxicol. 11, 391- 407.]. Accordingly, we found a notable difference in the viability of the progeny, which showed that the AAF adduct was most toxic and that the AF adduct was least toxic, with the AP adduct exhibiting intermediate toxicity. However, analysis of the progeny showed that translesion synthesis was predominantly error-free. Only low-level mutations (<3%) were detected with G-->T as the dominant type of mutation by all three DNA adducts. When C8-AP-dG was evaluated in a repetitive 5'-CGC GCG-3' sequence, higher mutational frequency ( approximately 8%) was observed. Again, G-->T was the major type of mutations in simian kidney cells, even though in bacteria CpG deletions predominate in this sequence [Hilario, P., Yan, S., Hingerty, B. E., Broyde, S., and Basu, A. K. (2002) Comparative mutagenesis of the C8-guanine adducts of 1-nitropyrene,and 1,6- and 1,8-dinitropyrene in a CpG repeat sequence: A slipped frameshift intermediate model for dinucleotide deletion. J. Biol. Chem. 277, 45068- 45074.]. Mutagenesis of C8-AP-dG in a 12-mer containing the local DNA sequence around codon 273 of the p53 tumor suppressor gene, where the adduct was located at the second base of this codon, was also investigated. In this 5'-GTGC GTGTTTGT-3' site, the mutations were slightly lower but not very different from the progeny derived from the 5'-CGC GCG-3' sequence. However, the mutational frequency increased by more than 50% when the 5'-C to the adduct was replaced with a 5-methylcytosine (5-MeC). With a 5-MeC, the most notable change in mutation was the enhancement of G-->A, which occurred 2.5 times relative to a 5'-C. The C8-AP-dG adduct in codon 273 dodecamer sequence with a 5'-C or 5-MeC was also evaluated in human embryonic kidney (293T) cells. Similar to COS cells, targeted mutations doubled with a 5-MeC 5' to the adduct. Except for an increase in G-->C transversions, the results in 293T were similar to that in COS cells. We conclude that C8-AP-dG mutagenesis depends on the type of cell in which it is replicated, the neighboring DNA sequence, and the methylation status of the 5'-C.     

Structure of the 1,N2-ethenodeoxyguanosine adduct opposite cytosine in duplex DNA: Hoogsteen base pairing at pH 5.2

The exocyclic 1, N(2)-ethenodeoxyguanosine (1,N(2)-is an element of dG) adduct, arising from the reaction of vinyl halides and other vinyl monomers, including chloroacetaldehyde, and lipid peroxidation products with dG, was examined at pH 5.2 in the oligodeoxynucleotide duplex 5'-d(CGCATXGAATCC)-3'.5'-d(GGATTCCATGCG)-3' (X = 1,N(2)-is an element of dG). Previously, X(anti).C(anti) pairing was established in this duplex, containing the 5'-TXG-3' sequence context, at pH 8.6 [Shanmugam, G., Goodenough, A. K., Kozekov, I. D., Harris, T. M., Guengerich, F. P., Rizzo, C. J., and Stone, M. P. (2007) Chem. Res. Toxicol. 20, 1601- 1611]. At pH 5.2, the 1,N(2)-is an element of dG adduct decreased the thermal stability of the duplex by approximately 13 degrees C. The 1,N(2)-is an element of dG adduct rotated about the glycosyl bond from the anti to the syn conformation. This resulted in the observation of a strong nuclear Overhauser effect (NOE) between the imidazole proton of 1,N(2)-is an element of dG and the anomeric proton of the attached deoxyribose, accompanied by an NOE to the minor groove A(20) H2 proton from the complementary strand. The syn conformation of the glycosyl bond at 1,N(2)-is an element of dG placed the exocyclic etheno moiety into the major groove. This resulted in the observation of NOEs between the etheno protons and the major groove protons of the 5'-neighboring thymine. The 1,N(2)-is an element of dG adduct formed a Hoogsteen pair with the complementary cytosine, characterized by downfield shifts of the amino protons of the cytosine complementary to the exocyclic adduct. The pattern of chemical shift perturbations indicated that the lesion introduced a localized structural perturbation involving the modified base pair and its 3'- and 5'-neighbor base pairs. A second conformational equilibrium was observed, in which both the modified base pair and its 3'-neighboring G.C base pair formed tandem Hoogsteen pairs. The results support the conclusion that at neutral pH, in the 5'-TXG-3' sequence, the 1,N(2)-is an element of dG adduct exists as a blend of conformations in duplex DNA. These involve the interconversion of the glycosyl torsion angle between the anti and the syn conformations, occurring at an intermediate rate on the NMR time scale.     

Molecular dynamics of a food carcinogen-DNA adduct in a replicative DNA polymerase suggest hindered nucleotide incorporation and extension

2-Amino-1-methyl-6-phenylimidazo[4,5-b]pyridine (PhIP) is the most abundant of the carcinogenic heterocyclic aromatic amines in the human diet, and the major mutagenic effect of dietary PhIP is G-->T transversions. The major PhIP-derived DNA adduct is to C8 of guanine. We have investigated this adduct in a PhIP-induced mutational hotspot 5'-GGGA-3' of the Apc tumor suppressor gene, frequently mutated in mammalian colon tumors. We have carried out a molecular dynamics study to elucidate on a structural level nucleotide incorporation and extension opposite this major adduct during replication. The PhIP adduct was modeled into the ternary complex closed conformation of DNA polymerase RB69, at incorporation and extension positions, with normal cytosine or mismatched partner adenine. RB69 polymerase is a member of the B family as are most replicative eukaryotic DNA polymerases such as DNA polymerase alpha. These systems were subjected to molecular dynamics simulations with AMBER. Our results show that the adduct can reside on the major groove side of the modified DNA template opposite an incoming dCTP or dATP. In the case of the normal partner, disturbance to the active site is observed at the incorporation step, but there is less perturbance in the extension simulation. In the case of the mismatched partner, a less disturbed active site is observed during the incorporation step, but extension appears to be more difficult. Disturbances include adverse impacts on Watson-Crick hydrogen bonding in the nascent base pair, on the distance between the alpha-phosphate of the incoming dNTP and the primer terminus 3'-OH, and on critical protein interactions with the dNTP. However, in all of these cases, a near reaction ready distance (within 3.5 angstroms) between the 3'-terminal oxygen of the primer and the Palpha of the incoming nucleotide triphosphate is sampled occasionally (0.4-23.5% of the time). Thus, error-free bypass or the induction of a G-->T transversion mutation could occur at times and contribute to an extent to the mutagenic effect of PhIP. Polymerase stalling would be the more common outcome and in vivo could lead to switch to an error-prone bypass polymerase.     

Dual roles of glycosyl torsion angle conformation and stereochemical configuration in butadiene oxide-derived N1 beta-hydroxyalkyl deoxyinosine adducts: a structural perspective

The solution structure of the N1-[1-hydroxy-3-buten-2(R)-yl]-2'-deoxyinosine adduct arising from the alkylation of adenine N1 by butadiene epoxide (BDO), followed by deamination to deoxyinosine, was determined in the oligodeoxynucleotide 5'-d(CGGACXAGAAG)-3'.5'-d(CTTCTTGTCCG)-3'. This oligodeoxynucleotide contained the BDO adduct at the second position of codon 61 of the human N-ras protooncogene (underlined) and was named the ras61 R-N1-BDO-(61,2) adduct. 1H NMR revealed a weak C5 H1' to X6 H8 nuclear Overhauser effects (NOE), followed by an intense X6 H8 to X6 H1' NOE. Simultaneously, the X6 H8 to X6 H3' NOE was weak. The resonances arising from the T16 and T17 imino protons were not observed. 1H NOEs between the butadiene moiety and the DNA positioned the adduct in the major groove. Structural refinement based upon a total of 394 NOE-derived distance restraints and 151 torsion angle restraints yielded a structure in which the modified deoxyinosine was in the syn conformation about the glycosyl bond, with a glycosyl bond angle of 83 degrees , and T17, the complementary nucleotide, was stacked into the helix but not hydrogen bonded with the adducted inosine. The refined structure provides a plausible hypothesis as to why these N1 deoxyinosine adducts strongly code for the incorporation of dCTP during trans lesion DNA replication, irrespective of stereochemistry, both in Escherichia coli [Rodriguez, D. A., Kowalczyk, A., Ward, J. B. J., Harris, C. M., Harris, T. M., and Lloyd, R. S. (2001) Environ. Mol. Mutagen. 38, 292-296] and in mammalian cells [Kanuri, M., Nechev, L. N., Tamura, P. J., Harris, C. M., Harris, T. M., and Lloyd, R. S. (2002) Chem. Res. Toxicol. 15, 1572-1580]. Rotation of the N1 deoxyinosine adduct into the syn conformation may facilitate incorporation of dCTP via Hoogsteen type templating with deoxyinosine, generating A to G mutations. However, conformational differences between the R- and the S-N1-BDO-(61,2) adducts, involving the positioning of the butenyl moiety in the major groove of DNA, suggest that adduct stereochemistry plays a secondary role in modulating the biological response to these adducts.     

Orientation of the crotonaldehyde-derived N2-[3-Oxo-1(S)-methyl-propyl]-dG DNA adduct hinders interstrand cross-link formation in the 5'-CpG-3' sequence

The conformation of the crotonaldehyde-derived N(2)-[3-oxo-1(S)-methyl-propyl]-dG adduct in the oligodeoxynucleotide 5'-d(G(1)C(2)T(3)A(4)G(5)C(6)X(7)A(8)G(9)T(10)C(11)C(12))-3'.5'-d(G(13)G(14)A(15)C(16)T(17)C(18)G(19)C(20)T(21)A(22)G(23)C(2)(4))-3', where X = N(2)-[3-oxo-1(S)-methyl-propyl]-dG, is reported. This adduct arises from opening of the cyclic N(2)-(S-alpha-CH(3)-gamma-OH-1,N(2)-propano-2')-dG adduct when placed opposite dC in duplex DNA. This oligodeoxynucleotide contains the 5'-CpG-3' sequence in which the N(2)-(R-alpha-CH(3)-gamma-OH-1,N(2)-propano-2')-dG but not the N(2)-(S-alpha-CH(3)-gamma-OH-1,N(2)-propano-2')-dG adduct preferentially formed an interstrand carbinolamine cross-link [Kozekov, I. D., Nechev, L. V., Moseley, M. S., Harris, C. M., Rizzo, C. J., Stone, M. P., and Harris, T. M. (2003) J. Am. Chem. Soc. 125, 50-61; Cho, Y.-J., Wang, H., Kozekov, I. D., Kurtz, A. J., Jacob, J., Voehler, M., Smith, J., Harris, T. M., Lloyd, R. S., Rizzo, C. J., and Stone, M. P. (2006) Chem. Res. Toxicol. 19, 195-208]. Analysis of (1)H NOE data, chemical shift perturbations, and deoxyribose pseudorotations and backbone torsion angles suggested the presence of a stable and ordered DNA conformation at pH 9.3 and 30 degrees C, with minimal conformational perturbation. The spectral line widths of the adduct protons were comparable to those of the oligodeoxynucleotide, suggesting that the correlation times of these protons were similar to those of the overall duplex. The crotonaldehydic-derived methyl protons showed NOEs in the 5'-direction to C(18) H1', G(19) H1', and G(19) H4' in the complementary strand of the duplex. The aldehyde proton of the adduct exhibited NOEs in the 3'-direction to A(8) H1' and A(8) H4' in the modified strand. All of these NOEs involved DNA protons facing the minor groove. Molecular dynamics calculations, restrained by distances and torsion angles derived from the NMR data, revealed that within the minor groove, the aldehyde of the N(2)-[3-oxo-1(S)-methyl-propyl]-dG adduct oriented in the 3'-direction, while the 1(S) methyl group oriented in the 5'-direction. This positioned the aldehyde distal to the G(19) exocyclic amine and provided a rationale as to why the N(2)-(S-alpha-CH(3)-gamma-OH-1,N(2)-propano-2')-dG adduct generated interstrand cross-links less efficiently than did the N(2)-(R-alpha-CH(3)-gamma-OH-1,N(2)-propano-2')-dG adduct.     

Inhibition of 2-amino-3-methylimidazo[4,5-f]quinoline-DNA adduct formation in CDF1 mice by heat-altered derivatives of linoleic acid.

Grilled ground beef contains a number of carcinogens, including aminoimidazoazaarenes, such as 2-amino-3-methylimidazo[4,5-f]quinoline (IQ), as well as anticarcinogenic substances, such as heat-generated derivatives of linoleic acid (CLA). In the present study, CLA was administered by gavage every other day to young adult CDF1 mice for a period of 45 days (50 microliters/48 hr for days 1-24 and 100 microliters/48 hr for days 25-45), using trioctanoin as a control. On day 46 all animals received a single oral dose (50 mg/kg) of IQ and tissues were collected 24 hr later. Tissue DNA was purified and analysed for IQ-DNA adducts by 32P-postlabelling assays. Compared with controls, CLA treatment caused a 43.1 and 31.8% inhibition of adduct formation in the livers of male and female mice, respectively. In the lung and large intestine CLA had a 74.2 and 39.4% inhibitory effect, respectively, in the female only, whereas there was no effect in the stomach or small intestine of either sex. In the kidneys of females, CLA treatment inhibited IQ-DNA adduct formation almost completely (95.2%), whereas in the kidneys of males CLA had no effect. It is concluded that CLA inhibits IQ-DNA adduct formation in certain IQ target organs (liver and lung) and non-target organs (large intestine, kidney), but is inactive in other target organs (stomach) and non-target organs (small intestine) of the CDF1 mouse.     

Structure of an oligodeoxynucleotide containing a 1,N(2)-propanodeoxyguanosine adduct positioned in a palindrome derived from the Salmonella typhimurium hisD3052 gene: Hoogsteen pairing at pH 5.2

The structure of the 1,N(2)-Propanodeoxyguanosine (PdG) adduct was determined at pH 5.2 in the oligodeoxynucleotide duplex 5'-d(CGCGGTXTCCGCG)3'.5'-d(CGCGGACACCGCG)-3' (X = PdG). This sequence, referred to as the -TXT- sequence, is contained within the Salmonella typhimurium hisD3052 gene and contains a palindrome, representing a potential hotspot for frameshift mutagenesis. PdG provides a model for the primary adduct induced in DNA by malondialdehyde, the 3-(2'-deoxy-beta-D-erythro-pentofuranosyl)pyrimido[1,2-a]-purin-10(3H)-one (M(1)G) lesion. The solution structure was refined by molecular dynamics calculations restrained by a combination of NMR-derived distances and dihedral angles, using a simulated annealing protocol. PdG introduced a localized perturbation into the sequence at base pair X(7).C(20), which was pH-dependent. At neutral pH, conformational exchange resulted in spectral line broadening, and it was not possible to determine the structure. A stable structure was observed at pH 5.2 in which PdG rotated about the glycosyl bond into the syn conformation. This placed the exocyclic moiety into the major groove of the duplex. PdG formed a protonated Hoogsteen pair with nucleotide C(20) in the complementary strand. The pseudorotation of the deoxyribose at C(20) was altered to an approximately equal blend of C2'-endo and C3'-endo structures. However, these made little difference in the overall structure of the modified oligodeoxynucleotide. The structure was compared to that of PdG in the 5'-d(CGCXCGGCATG)-3'.5'-(CATGCCGCGCG)-3' sequence (the -CXC- sequence) at pH 5.8 [Singh, U. S., Moe, J. G., Reddy, G. R., Weisenseel, J. P., Marnett, L. J., and Stone, M. P. (1993) Chem. Res. Toxicol. 6, 825-836]. A sequence effect was observed. When PdG was placed into the -TXT- sequence at low pH, the structural perturbation was limited to the X(7).C(20) base pair. In contrast, when PdG was placed into the -CXC- sequence at low pH, both the modified base pair and its 3'-neighbor base pair were disrupted. The results are discussed in the context of differential outcomes for site-specific mutagenesis and replication bypass experiments when PdG was placed in the -TXT- and -CXC- sequences, respectively.     

Characterization of DNA adducts formed in vitro by reaction of N-hydroxy-2-amino-3-methylimidazo[4,5-f]quinoline and N-hydroxy-2-amino-3,8-dimethylimidazo[4,5-f]quinoxaline at the C-8 and N2 atoms of guanine.

The covalent binding of the carcinogenic N-hydroxy metabolites of 2-amino-3-methylimidazo-[4,5-f]quinoline (IQ) and 2-amino-3,8-dimethylimidazo[4,5-f]quinoxaline (MeIQx) to deoxynucleosides and DNA was investigated in vitro. Two major adducts were formed by the reaction of the N-acetoxy derivatives of IQ and MeIQx with deoxyguanosine (dG); however, no adducts were formed with deoxycytidine, deoxyadenosine, or thymidine. From proton NMR and mass spectroscopic characterization the adducts were identified as 5-(deoxyguanosin-N2-yl)-2-amino-3-methylimidazo[4,5-f]quinoline (dG-N2-IQ),N-(deoxyguanosin-8-yl)-2-amino-3-methylimidazo-[4,5-f]q uinoline (dG-C8-IQ), 5-(deoxyguanosin-N2-yl)-2-amino-3,8-dimethylimidazo[4,5-f]qu inoxaline (dG-N2-MeIQx), and N-(deoxyguanosin-8-yl)-2-amino-3,8-dimethylimidazo[4,5-f]qui noxaline (dG-C8-MeIQx). The level of dG-C8 adducts was approximately 8-10 times greater than the amount of dG-N2 adducts formed from the reaction of dG with the N-acetoxy derivatives of IQ and MeIQx. The C-8-substituted dG adduct was also the major adduct formed from reactions of DNA with N-acetoxy-IQ and N-acetoxy-MeIQx. Approximately 60-80% of the bound carcinogens were recovered from DNA as dG-C8 adducts upon enzymatic digestion. The dG-N2 adducts also were detected and accounted for approximately 4% of the bound IQ and 10% of the bound MeIQx. These results suggest that the relative contributions of the nitrenium and carbenium ion resonance forms as well as DNA macromolecular structure are major determinants for DNA adduct substitution sites. Investigations on adduct conformation of 1H NMR spectroscopy revealed that the anti form is preferred for the dG-N2 adducts of IQ and MeIQx, while the syn form is preferred for the dG-C8 adducts. The possible role of these adducts in the initiation of carcinogenesis is discussed.     

Computer simulation of the binding of naphthyridinomycin and cyanocycline A to DNA

Cyanocycline A was found to have a pKa of 6.6. Protonation of N14 was established by 1H NMR spectroscopy. In strongly acidic solution the oxazolidine ring opened irreversibly. A model was derived for the binding of naphthyridinomycin and cyanocycline A to the hexanucleotide duplex d(ATGCAT)2, by using the molecular mechanics and dynamics modules of AMBER 3.0. It involved protonation on the oxazolidine-ring nitrogen, reduction of the quinone ring to a hydroquinone, formation of an iminium ion with loss of the C7 substituent, noncovalent binding in the minor groove with the hydroquinone ring in the 3'-direction from guanine, and covalent binding to the 2-amino group of this guanine with C7 adopting the R configuration. This model is consistent with the experimental evidence on the DNA binding of these drugs. An alternative binding mode based on opening of the oxazolidine ring and alkylation at C3a also was feasible according to molecular mechanics calculations. The geometry of naphthyridinomycin does not permit interstrand cross-linking involving both C3a and C7, but formation of a cross-link to protein appears possible. When the covalent naphthyridinomycin-d(ATGCAT)2 models were refined in the presence of water and counterions, the models with the most favorable net binding enthalpies were the same as those produced by simulation in vacuum. Qualitative estimates of the relative entropy changes resulting from adduct formation were based on the number of ordered (hydrogen bonded) water molecules released from d(ATGCAT)2 and from the drug. In all cases but one, d(ATGCAT)2 loses five water molecules. It loses six in the C3a covalent model with 5',S geometry. Naphthyridinomycin hydroquinone loses up to two water molecules, depending on the particular adduct. The 3',R model was again favored for the C7 covalent adduct. Among the C3a covalent models, the one with 5',R geometry lost the second most water molecules, but it had the best binding enthalpy.     

Effect of CYP1A2 deficiency on heterocyclic amine DNA adduct levels in mice.

The contribution of CYP1A2 to the formation of DNA adducts of the cooked meat-derived heterocyclic amines (HCAs) 2-amino-3-methylimidazo[4,5-f]quinoline (IQ) and 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine (PhIP) was examined in CYP1A2-null (knock-out, KO) and wild-type (WT) mice. IQ (25 mg and 75 mg/kg) and PhIP (150 mg/kg) were administered by gavage to mice and DNA adduct levels in liver, kidney, mammary gland and colon were examined by the 32P-postlabeling assay. Three hours after either dose of IQ, adducts levels in liver and kidney of KO mice were 20-30% of the levels in WT mice, a difference that was statistically significant (Student's t-test, P < 0.05). In the colon, adduct levels in KO mice were significantly lower than in the WT mice only at the lowest dose of IQ (1.6+/-0.6 vs 4.6+/-0.7, respectively, relative adduct labeling (RAL) x 10(8), mean+/-S.E.M., n = 3-5 mice). In the mammary gland, however, there was no difference in IQ-DNA adduct levels in KO and WT mice at either dose of IQ. Three hours after dosing with PhIP, PhIP-DNA adduct levels were statistically significantly lower in KO mice than in WT mice in all tissues examined. PhIP-DNA adducts in liver and kidney of WT mice were 9.9+/-1.1 and 22.5+/-6.9, respectively, whereas no PhIP-DNA adducts were detected in either organ of KO mice (limit of detection, 1.4-2.8 x 10(9)). PhIP-DNA adduct levels in mammary gland and colon of WT mice were 47.1+/-9.5 and 58.0+/-21.7, respectively, but accordingly only 3.8+/-0.7 and 5.4+/-0.9 in KO mice. The findings indicate that CYP1A2, responsible for IQ and PhIP N-hydroxylation, the first step in the metabolic action, significantly effects DNA adduct formation in vivo. However, the data raise the possibility that other cytochromes P450 as well as other pathways of activation potentially contribute to DNA adduct formation in specific organs, depending on the HCA substrate.     

A molecular model for DNA cross-linking by the antitumor agent azinomycin B

A computational model for the covalent interstrand DNA cross-linking of the antitumor agent azinomycin B is reported and is based on Monte Carlo simulations of the four possible monoalkylation species and an examination of the low energy conformations of the cross-linked agent. The model was developed using a suitably modified version of the AMBER* force field with the experimentally determined triplet DNA target sequence 5'-d(GCT)-3' in both the native B-form and containing a preformed intercalation site.     

Fluorescence probing of aminofluorene-induced conformational heterogeneity in DNA duplexes

Fluorescence spectroscopy was used to study carcinogen-induced conformational heterogeneity in DNA duplexes. The fluorophore 2-aminopurine (AP) was incorporated adjacent (5') to the lesion (G*) in eight different DNA duplexes [d(5'-CTTCT PG* NCCTC-3'):d(5'-GAGGN XTAGAAG-3'), G* = FAF adduct, P = AP, N = G, A, C, T, and X = C, A] modified by FAF [ N-(2'-deoxyguanosin-8-yl)-7-fluoro-2-aminofluorene], a fluorine-tagged model DNA adduct derived from the potent carcinogen 2-aminofluorene. Steady-state measurements showed that fluorescence intensity and Stern-Volmer constants ( Ksv) derived from acrylamide quenching experiments decreased for all carcinogen-modified duplexes relative to the controls, which suggests greater AP stacking in the duplex upon adduct formation. Conformation-specific stacking of AP with the neighboring adduct was evidenced by a sequence-dependent variation in fluorescence intensity, position of emission maximum, degree of emission quenching by acrylamide, and temperature-dependent spectral changes. The magnitude of stacking was in the order of FAF residue in base-displaced stacked (S) > minor groove wedged (W) > major groove B type (B). This work represents a novel utility of AP in probing adduct-induced conformational heterogeneities in DNA duplexes.     

Stereospecific structural perturbations arising from adenine N(6) butadiene triol adducts in duplex DNA

Butadiene is oxidized in vivo to form stereoisomeric butadiene diol epoxides (BDE). These react with adenine N(6) in DNA yielding stereoisomeric N(6)-(2,3,4-trihydroxybutyl)-2'-deoxyadenosyl (BDT) adducts. When replicated in Escherichia coli, the (2R,3R)-N(6)-(2,3,4-trihydroxybutyl)-2'-deoxyadenosyl adduct yielded low levels of A-->G mutations whereas the (2S,3S)-N(6)-(2,3,4-trihydroxybutyl)-2'-deoxyadenosyl butadiene triol adduct yielded low levels of A-->C mutations [Carmical, J. R., Nechev, L. V., Harris, C. M., Harris, T. M., and Lloyd, R. S. (2000) Environ. Mol. Mutagen. 35, 48-56]. Accordingly, the structure of the (2R,3R)-N(6)-(2,3,4-trihydroxybutyl)-2'-deoxyadenosyl adduct at position X(6) in d(CGGACXAGAAG).d(CTTCTTGTCCG), the ras61 R,R-BDT-(61,2) adduct, was compared to the corresponding structure for the (2S,3S)-N(6)-(2,3,4-trihydroxybutyl)-2'-deoxyadenosyl adduct in the same sequence, the ras61 S,S-BDT-(61,2) adduct. Both the R,R-BDT-(61,2) and S,S-BDT-(61,2) adducts are oriented in the major groove of the DNA, accompanied by modest structural perturbations. However, structural refinement of the two adducts using a simulated annealing restrained molecular dynamics (rMD) approach suggests stereospecific differences in hydrogen bonding between the hydroxyl groups located at the beta- and gamma-carbons of the BDT moiety, and T(17) O(4) of the modified base pair X(6).T(17). The rMD calculations predict hydrogen bond formation between the gamma-OH and the T(17) O(4) in the R,R-BDT-(61,2) adduct whereas in the S,S-BDT-(61,2) adduct, hydrogen bond formation is predicted between the beta-OH and the T(17) O(4). This difference positions the two adducts differently in the major groove. This may account for the differential mutagenicity of the two adducts and suggests that the two adducts may interact differentially with other DNA processing enzymes. With respect to mutagenesis in E. coli, the minimal perturbation of DNA induced by both major groove adducts correlates with their facile bypass by three E. coli DNA polymerases in vitro and may account for their weak mutagenicity [Carmical, J. R., Nechev, L. V., Harris, C. M., Harris, T. M., and Lloyd, R. S. (2000) Environ. Mol. Mutagen. 35, 48-56].     

Structural refinement of the 8,9-dihydro-8-(N7-guanyl)-9-hydroxy-aflatoxin B(1) adduct in a 5'-Cp(AFB)G-3' sequence

The structure of the cationic 8,9-dihydro-8-(N7-guanyl)-9-hydroxy-aflatoxin B(1) adduct embedded in a 5'-CpG-3' sequence context and paired with deoxycytosine in the oligodeoxynucleotide d(ACATC(AFB)GATCT) x d(AGATCGATGT) was refined using molecular dynamics calculations restrained by NOE data and dihedral angle restraints obtained from NMR data. The aflatoxin moiety intercalated above the 5' face of the modified guanine. It stacked between C(5) x G(16) and (AFB)G(6) x C(15). The AFB(1) H5, OCH(3), and methylene protons faced into the minor groove, with the methylene protons oriented between the C(15) and G(16) nucleobases. The aflatoxin B(1) H6a, H8, H9, and H9a protons faced the major groove, with H6a and H9a pointing toward the 5' direction from the lesion site. The refined structure was compared to the structure of the aflatoxin B(1) adduct embedded in a 5'-ATGCAT-3' sequence in the oligodeoxynucleotide d(TAT(AFB)GCATA)(2) [Jones, W. R., Johnston, D. S., and Stone, M. P. (1998) Chem. Res. Toxicol.11, 873-881]. The structure of the intercalated aflatoxin B(1) lesion in the ATC(AFB)GAT sequence is similar to its structure in the d(AT(AFB)GCAT) sequence. This is consistent with a mechanism in which the precovalent intercalation of aflatoxin-8,9-exo-epoxide on the 5' face of guanine places the epoxide in close proximity and in the proper orientation to the N7 position of guanine, thus facilitating an S(N)2 reaction. The data provides additional insight into the nature of the disruption of the B-DNA duplex induced by aflatoxin B(1) intercalation. Overall, the results suggest that sequence contributes a minor role in modulating the structure of the cationic guanine N7 AFB(1) lesion in duplex DNA. On the other hand, structural differences are observed when the correctly paired structure is compared to the structure of the cationic AFB(1) adduct mispaired with dA [Giri, I., Johnston, D. S., and Stone, M. P. (2002) Biochemistry 41, 5462-5472].     

Conformational differences of the C8-deoxyguanosine adduct of 2-amino-3-methylimidazo[4,5-f]quinoline (IQ) within the NarI recognition sequence

2-amino-3-methylimidazo[4,5-f]quinoline (IQ) is a highly mutagenic heterocyclic amine found in cooked meats. The major DNA adduct of IQ is at the C8-position of dGuo. We have previously reported the incorporation of the C8-IQ adduct into oligonucleotides, namely, the G1-position of codon 12 of the N-ras oncogene sequence (G1G2T) and the G3-position of the NarI recognition sequence (G1G2CG3CC) (Elmquist et al. (2004) J. Am. Chem. Soc. 126, 11189-11201). Ultraviolet spectroscopy and circular dichroism studies indicated that the conformation of the adduct in the two oligonucleotides was different, and they were assigned as groove-bound and base-displaced intercalated, respectively. The conformation of the latter was subsequently confirmed through NMR and restrained molecular dynamics studies (Wang et al. (2006) J. Am. Chem. Soc. 128, 10085-10095). We report here the incorporation of the C8-IQ adduct into the G1- and G2-positions of the NarI sequence. A complete analysis of the UV, CD, and NMR chemical shift data for the IQ protons are consistent with the IQ adduct adopting a minor groove-bound conformation at the G1- and G2-positions of the NarI sequence. To further correlate the spectroscopic data with the adduct conformation, the C8-aminofluorene (AF) adduct of dGuo was also incorporated into the NarI sequence; previous NMR studies demonstrated that the AF-modified oligonucleotides were in a sequence-dependent conformational exchange between major groove-bound and base-displaced intercalated conformations. The spectroscopic data for the IQ- and AF-modified oligonucleotides are compared. The sequence-dependent conformational preferences are likely to play a key role in the repair and mutagenicity of C8-arylamine adducts.     

Molecular modeling of the intrastrand guanine-guanine DNA adducts produced by cisplatin and oxaliplatin.

Intrastrand DNA adducts formed by cisplatin and oxaliplatin were modeled with molecular mechanics minimization and restrained molecular dynamics simulations in a comparative study. A reasonable set of force field parameters for the Pt atom were refined by using the available cisplatinated DNA crystal structure as a guide. This crystal structure was also used as the starting structure for the simulations. Analysis of the resulting structures indicated that the covalent effects of oxaliplatin coordination on DNA structure were very similar to those of cisplatin. The most prominent difference between the two structures resulted from the presence of the 1,2-diaminocyclohexane ring in the oxaliplatin adduct. The modeling indicated that this ring protrudes directly outward into, and fills much of, the narrowed major groove of the bound DNA, forming a markedly altered and less polar major groove in the area of the adduct. The differences in the structure of the adducts produced by cisplatin and oxaliplatin are consistent with the observation that they are differentially recognized by the DNA mismatch repair system.