Mapping the binding site of aflatoxin B1 in DNA: molecular modeling of the binding sites for the N(7)-guanine adduct of aflatoxin B1 in different DNA sequences.

Aflatoxin B1 (AFB1), a potent mutagen and carcinogen, forms an adduct exclusively at the N(7) position of guanine, but the structure of this adduct in double stranded DNA is not known. Molecular modeling (using the program, PSFRODO) in conjunction with molecular mechanical calculation (using the program, AMBER) are used to assess the binding modes available to this AFB1 adduct. Two modes appear reasonable; in one the AFB1 moiety is intercalated between the base pair containing the adducted guanine and the adjacent base pair on the 5'-side in reference to the adducted guanine, while in the second it is bound externally in the major groove of DNA. Rotational flexibility appears feasible in the latter providing four, potential binding sites. Molecular modeling reveals that the binding sites around the reactive guanine in different sequences are not uniformly compatible for interaction with AFB1. As the sequence is changed, one particular external binding site would be expected to give a pattern of reactivities that is reasonably consistent with the observed sequence specificity of binding that AFB1 shows in its reaction with DNA (Benasutti, M., Ejadi, S., Whitlow, M. D. and Loechler, E. L. (1988) Biochemistry 27, 472-481). The AFB1 moiety is face-stacked in the major groove with its long axis approximately perpendicular to the helix axis. Favorable interactions are formed between exocyclic amino groups that project into the major groove on cytosines and adenines surrounding the reactive guanine, and oxygens in AFB1; unfavorable interactions involve van der Waals contacts between the methyl group on thymine and the AFB1 moiety. "Some of the sequence specificity of binding data can be rationalized more readily if it is assumed that 5'-GG-3' sequences adopt an A-DNA structure." Based upon molecular modeling/potential energy minimization calculation, it is difficult to predict how reactivity would change in different DNA sequences in the case of the intercalative binding mode; however, several arguments suggest that intercalation might not be favored. From these considerations a model of the structure for the transition state in reaction of AFB1 with DNA is proposed involving one particular external binding site.

Mapping the binding site of aflatoxin B1 in DNA: systematic analysis of the reactivity of aflatoxin B1 with guanines in different DNA sequences.

The mutagenic and carcinogenic chemical aflatoxin B1 (AFB1) reacts almost exclusively at the N(7)-position of guanine following activation to its reactive form, the 8,9-epoxide (AFB1 oxide). In general N(7)-guanine adducts yield DNA strand breaks when heated in base, a property that serves as the basis for the Maxam-Gilbert DNA sequencing reaction specific for guanine. Using DNA sequencing methods, other workers have shown that AFB1 oxide gives strand breaks at positions of guanines; however, the guanine bands varied in intensity. This phenomenon has been used to infer that AFB1 oxide prefers to react with guanines in some sequence contexts more than in others and has been referred to as "sequence specificity of binding". Herein, data on the reaction of AFB1 oxide with several synthetic DNA polymers with different sequences are presented, and (following hydrolysis) adduct levels are determined by high-pressure liquid chromatography. These results reveal that for AFB1 oxide (1) the N(7)-guanine adduct is the major adduct found in all of the DNA polymers, (2) adduct levels vary in different sequences, and, thus, sequence specificity is also observed by this more direct method, and (3) the intensity of bands in DNA sequencing gels is likely to reflect adduct levels formed at the N(7)-position of guanine. Knowing this, a reinvestigation of the reactivity of guanines in different DNA sequences using DNA sequencing methods was undertaken. The reactivities of 190 guanines were determined quantitatively and considered in a pentanucleotide context, 5'-WXGYZ-3', where the central, underlined G represents the reactive guanine and W, X, Y, and Z can be any of the nucleotide bases. Methods are developed to determine that the X (5'-side) base and the Y (3'-side) base are most influential in determining guanine reactivity. The influence of the bases in the 5'-position (X) is 5'-G (1.0) greater than C (0.8) greater than A (0.3) greater than T (0.2), while the influence of the bases in the 3'-position (Y) is 3'-G (1.0) greater than T (0.8) greater than C (0.4) greater than A (0.3). These rules in conjunction with molecular modeling studies (to be published elsewhere) were used to assess the binding sites that might be utilized by AFB1 oxide in its reaction with DNA.

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.