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

Non-covalent interactions of the carcinogen (+)-anti-BPDE with exon 1 of the human K-ras proto-oncogene

Investigating the complementary, but different, effects of physical (non-covalent) and chemical (covalent) mutagen-DNA and carcinogen-DNA interactions is important for understanding possible mechanisms of development and prevention of mutagenesis and carcinogenesis. A highly mutagenic and carcinogenic metabolite of the polycyclic aromatic hydrocarbon benzo[α]pyrene, namely (+)-anti-BPDE, is known to undergo both physical and chemical complexation with DNA. Previous studies of BPDE-DNA complex formation have focused on processes that require substantial structural reorganization, such as intercalation, and consequently relatively long time scales. However, some initial processes which occur within shorter time scales, such as external non-covalent binding, and which do not require major DNA structural reorganization have not been thoroughly investigated. A detailed computational study of such initial BPDE-DNA interactions is needed to elucidate the temporal and structural origins of the major covalent adduct, a promutagenic, which is known to exist in an external (+)-trans-anti-BPDE-N(2)-dGuanosine configuration. Accordingly, the initial stages of external non-covalent BPDE-DNA binding are studied in this work as well as their relationship to subsequent formation of the major, also external, covalent adduct. To study mechanisms that occur prior to extensive DNA structural reorganization, we present a first and detailed codon by codon computational study of the non-covalent interactions of (+)-anti-BPDE with DNA. In particular, due to its relevance to carcinogenesis, the interaction of (+)-anti-BPDE with exon 1 of the human K-ras gene has been studied. External solvent-exposed non-covalent binding sites have been found which may be precursors of the major external trans adduct and, importantly, are located in codons 12 and 13 of the K-ras gene which are known to be key mutation hotspots. In addition, our study explains and correctly predicts preferential (+)-anti-BPDE binding at minor groove guanosines. A subtle combination of van der Waals and hydrogen bonding interactions has been found to be a primary factor in preferentially positioning (+)-anti-BPDE toward the 5' position of a guanosine's strand, consistent with proton NMR observations for the major trans adduct, and at 5'-TGG-3' sequences which are known to yield high binding probability.

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.

Aflatoxin B1 and DNA adducts. Proposed model for surface noncovalent and covalent complexes with N(7) of guanine. II

An alternate model for surface noncovalent and surface covalent binding of aflatoxin B1 to N(7) of guanine in DNA is proposed. This model considers the out-of-plane motions of C(8) of aflatoxin B1 in those interactions. The covalent intercalated fit of aflatoxin B1 into DNA arises from steric adjustments made by DNA at the covalent intercalation site as well as local strain in the bond angles about N(7) of guanine and C(8) of aflatoxin B1. The bond angle about N(7) deviates modestly from the sp2 value toward the sp3 value. This study suggests that the surface covalent aflatoxin B1-DNA complex serves only a minor role in aflatoxin's precarcinogenic interaction with DNA and is a likely correctable error.

The structure of the aflatoxin B1-DNA adduct at N7 of guanine. Theoretical intercalation and covalent adduct models

Two theoretical models are proposed for the conformational structure of both intercalated and covalent adduct complexes of aflatoxin B1, designated AFB1, with N7 of guanine of DNA. The covalent adduct model requires the DNA to kink a minimum of 39 degrees about the covalent site of the C8 and N7 atoms comprising the bond of the covalent complex. The preference of AFB1 for specific G bases within a sequence of GC content followed that of experimental studies with the added feature that for binding to the third G base of a tetramer sequence from the 3'-end, the AFB1 displayed enhanced binding at the 3' site of the targeted guanosine. Binding of AFB1 to the second G base of a tetramer sequence from the 3'-end leads to preference for a 5' site of the targeted guanosine. Inhibition of AFB1's interaction with the targeted DNA in the presence of intercalated ethidium bromide is explained by these proposed models.

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

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

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.

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.

Sequence selectivity, a test of the nature of the covalent adduct formed between benzo[a]pyrene and DNA

A theoretical study is presented of the energetic and structural properties of covalent adducts of benzo[a]pyrene and a DNA fragment. Energy optimisation is performed with the use of minimiser with constraints and an advanced semiempirical energy formula. Three types of adducts are studied: an external complex with the benzopyrene located in the DNA minor groove and two types of intercalative complexes with the carcinogen situated on the 3' side and 5' side of the covalently bound guanine. For each of the adducts the effects of DNA base sequence are examined. It is shown that the results for the intercalative complex with the carcinogen situated on the 5' side of the modified guanine correlate with the experimentally determined sequence preference.

Kinetic flow dichroism study of conformational changes in supercoiled DNA induced by ethidium bromide and noncovalent and covalent binding of benz[a]pyrene diol epoxide

The dynamic conformational changes due to the noncovalent intercalative binding of ethidium bromide and racemic trans-7,8-dihydroxy-anti-9,10-epoxy-7,8,9,10-tetrahydrobenzo[a]pyrene (BPDE), and the covalent binding of BPDE to supercoiled phi X174 DNA, have been studied by gel electrophoresis and a novel application of a kinetic flow linear dichroism technique. The magnitude of the linear dichroism (delta A) of the DNA oriented in the flow gradient is sensitive to the hydrodynamic shape of the DNA molecule which is affected by the binding of the drug or the carcinogen BPDE. While the linear dichroism of ethidium bromide supercoiled DNA is time independent, the delta A spectra of BPDE-DNA reaction mixtures vary on time scales of minutes, which correspond to the reaction rate constant of BPDE to form 7,8,9,10-tetrahydroxytetrahydrobenzo[a]pyrene hydrolysis products and covalent DNA adducts. The rapid noncovalent intercalation of BPDE causes an initial large increase in delta A (up to 250%, corresponding to the dichroism observed with relaxed circular DNA), followed by a slower decrease in the linear dichroism signal. This decrease in delta A is attributed to the removal of intercalated diol epoxide molecules and the resulting reversible increase in the number of superhelical turns. The kinetic flow dichroism spectra indicate that the noncovalent BPDE-DNA complexes are intercalative in nature, while the covalent adducts are characterized by a very different conformation in which the long axes of the pyrenyl residues are oriented at a large angle with respect to the average orientation of the planes of the DNA bases.(ABSTRACT TRUNCATED AT 250 WORDS)

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.

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.

DNA modification by chemical carcinogens

The chemistry and molecular biology of DNA adducts is only one part of the carcinogenic process. Many other factors will determine whether a particular chemical will exert a carcinogenic effect. For example, the size of particles upon which a carcinogenic may be adsorbed will influence whether or not, and if so where, deposition within the lung will occur. The simultaneous exposure to several different agents may enhance or inhibit the metabolism of a chemical to its ultimate carcinogenic form (Rice et al., 1984; Smolarek and Baird, 1984). The ultimate carcinogenic metabolites may be influenced in their ability to react with DNA by a number of factors such as internal levels of detoxifying enzymes, the presence of other metabolic intermediates such as glutathione with which they could react either enzymatically or non-enzymatically, and the state of DNA which is probably most heavily influenced by whether or not the cell is undergoing replication or particular sequences being expressed. Replicating forks have been shown to be more extensively modified than other areas of DNA. Another critical factor which can influence the final outcome of the DNA damage is whether or not the modifications can be repaired. If this occurs with high fidelity and the cell has not previously undergone replication then the effect of the damage by the carcinogen is likely to be minimal. The major area in which progress is needed is an understanding of what this damage really does to the cell such that after an additional period of time, which may be as long as twenty or more years, these prior events are expressed and cell proliferation occurs. Clearly additional stimulatory factors, for example tumor promoting agents such as the phorbol esters or phenobarbital, are often needed. After such prolonged periods it seems likely that the DNA adducts would no longer be present. However, the way in which their earlier presence is remembered is not clear. Simple mutations do not explain all the characteristics of tumor progression and, when it occurs, regression. Even if a specific site mutation does occur then its expression must be under other types of control. Any explanation of the action of DNA modification at the molecular level also requires that account be taken of the diverse nature of the DNA adducts from simple modifications such as methylation to bulkier adducts such as benzo[a]pyrene, aflatoxin or aromatic amines.(ABSTRACT TRUNCATED AT 400 WORDS)

Stereoselective covalent binding of anti-benzo(a)pyrene diol epoxide to DNA conformation of enantiomer adducts

The conformation of adducts derived from the reactions and covalent binding of the (+) and (-) enantiomers of 7 beta, 8 alpha-dihydroxy-9 alpha, 10 alpha-epoxy-7,8,9,10-tetrahydrobenzo(a)pyrene (anti-BaPDE) with double-stranded calf thymus DNA in vitro were investigated utilizing the electric linear dichroism technique. The linear dichroism and absorption spectra of the covalent DNA complexes are interpreted in terms of a superposition of two types of binding sites. One of these conformations (site I) is a complex in which the plane of the pyrene residue is close to parallel (within 30 degrees) to the planes of the DNA bases (quasi-intercalation), while the other (site II) is an external binding site; this latter type of adduct is attributed to the covalent binding of anti-BaPDE to the exocyclic amino group of deoxyguanine (N2-dG), while site I adducts are attributed to the O6-deoxyguanine and N6-deoxyadenine adducts identified in the product analysis of P. Brookes and M.R. Osborne (Carcinogenesis (1982) 3, 1223-1226). Site II adducts are dominant (approximately 90% in the covalent complexes derived from the (+) enantiomer), but account for only 50 +/- 5% of the adducts in the case of the (-)-enantiomer. The orientation of site II complexes is different by 20 +/- 10 degrees in the adducts derived from the binding of the (+) and the (-) enantiomers to DNA, the long axis of the pyrene chromophore being oriented more parallel to the axis of the DNA helix in the case of the (+) enantiomer. These findings support the proposals by Brookes and Osborne that the difference in spatial orientation of the N2-dG adducts of (-)-anti-BaPDE together with their lower abundance may account for the lower biological activity of the (-) enantiomer. The external site II adducts, rather than site I adducts, appear to be correlated with the biological activity of these compounds.