Minimized conformational potential energy calculations have been performed for the major adduct of N-2-acetylaminofluorene (AAF) and aminofluorene (AF) with dCpdG. The DNA backbone conformations that produce base displacement, about which no information was previously available, were ascertained. For the AAF adduct, five different combinations of the DNA backbone torsions ω', ω and ψ produce low-energy conformers in which the fluorene is coplanar with cytidine, and the guanine is syn and twisted nearly perpendicular to the AAF. For the deacetylated AF adduct, three additional backbone types yielded fluorene-cytidine stacked conformers. In these conformers the guanine was anti. The various modes of achieving base-displacement may occur under differing circumstances, such as extent of modification, base sequence, salt concentration, and the nature of the adduct (acetylated or deacetylated).

It is well known that CpG dinucleotide steps in DNA, which are highly methylated at the 5-position of cytosine (meC) in human tissues, exhibit a disproportionate number of mutations within certain codons of the p53 gene. There is ample published evidence indicating that the reactivity of guanine with anti-B[a]PDE (a metabolite of the environmental carcinogen benzo[a]pyrene) at CpG mutation hot spots is enhanced by the methylation of the cytosine residue flanking the target guanine residue on the 5′-side. In this work we demonstrate that such a methylation can also dramatically affect the conformational characteristics of an adduct derived from the reaction of one of the two enantiomers of anti-B[a]PDE with the exocyclic amino group of guanine ([BP]G adduct). A detailed NMR study indicates that the 10R (−)-trans-anti-[BP]G adduct undergoes a transition from a minor groove-binding alignment of the aromatic BP ring system in the unmethylated C-[BP]G sequence context, to an intercalative BP alignment with a concomitant displacement of the modified guanine residue into the minor groove in the methylated meC-[BP]G sequence context. By contrast, a minor groove-binding alignment was observed for the stereoisomeric 10S (+)-trans-anti-[BP]G adduct in both the C-[BP]G and meC-[BP]G sequence contexts. This remarkable conformational switch resulting from the presence of a single methyl group at the 5-position of the cytosine residue flanking the lesion on the 5′-side, is attributed to the hydrophobic effect of the methyl group that can stabilize intercalated adduct conformations in an adduct stereochemistry-dependent manner. Such conformational differences in methylated and unmethylated CpG sequences may be significant because of potential alterations in the cellular processing of the [BP]G adducts by DNA transcription, replication, and repair enzymes.

An extensive conformational analysis has been carried out for two diastereoisomeric pairs of model estrogen quinone-derived DNA adducts, N6-(2-hydroxyestron-6(alpha,beta)-yl)-2'-deoxyadenosine (2-OHE1-6(alpha,beta)-N6-dA) and N2-(2-hydroxyestron-6(alpha,beta)-yl)-2'-deoxyguanosine (2-OHE1-6(alpha,beta)-N2-dG), in a B-DNA duplex and at a primer-template junction in a pol alpha family DNA polymerase. In vitro primer extension studies in pol alpha [Terashima, I., et al. (1998) Biochemistry 37, 13807-13815] have shown that the adenine adducts can incorporate dT, together with a small proportion of the incorrect base dC opposite the lesion, and they block less strongly than the guanine adducts. We have carried out conformational searches with energy minimization for four DNA duplexes containing 2-OHE1-6alpha-N6-dA, 2-OHE1-6beta-N6-dA, 2-OHE1-6alpha-N2-dG, or 2-OHE1-6beta-N2-dG. Our searches revealed that the four-ring nonplanar 2-hydroxyestrone (2-OHE1) moiety strongly prefers to reside in the major groove of the adenine adducts or the minor groove of the guanine adducts in a B-DNA duplex, with stereochemistry-dependent orientational differences in each case. No low energy conformations involving intercalation of the 2-OHE1 moiety were located in the searches. This stems from the largely nonplanar, nonaromatic nature of the 2-OHE1 ring system and implies that the proclivity for such bulky, nonplanar adducts to reside at the DNA helix exterior is a plausible conformational feature of other structurally similar estrogen quinone-derived DNA adducts, independent of base sequence context. In addition, the adenine adduct isomers, located in the major groove, manifest serious disturbance to the Watson-Crick base pairs at and near the lesion site, suggesting repair susceptibility. Possible structures of these adducts in a pol alpha family polymerase were also investigated through molecular modeling. The results rationalized the experimental in vitro primer extension studies. In addition, poor accommodation of the beta-stereoisomers within the polymerase was noted, suggesting that these stereoisomers would be more prone to cause blockage. Stereochemistry-dependent differences in adduct orientation could be expected to produce different biochemical effects, as has been observed in adducts derived from polycyclic aromatic hydrocarbons.

The pyridyloxobutylating agents derived from metabolically activated tobacco-specific nitrosamines can covalently modify guanine bases in DNA at the O(6) position. The adduct formed, O(6)-[4-oxo-4-(3-pyridyl)butyl]guanine ([POB]dG), results in mutations that can lead to tumor formation, posing a significant cancer risk to humans exposed to tobacco smoke. A combined NMR-molecular mechanics computational approach was used to determine the solution structure of the [POB]dG adduct within an 11mer duplex sequence d(CCATAT-[POB]G-GCCC).d(GGGCCATATGG). In agreement with the NMR results, the POB ligand is located in the major groove, centered between the flanking 5'-side dT.dA and the 3'-side dG.dC base pairs and thus in the plane of the modified [POB]dG.dC base pair, which is displaced slightly into the minor groove. The modified base pair in the structure adopts wobble base pairing (hydrogen bonds between [POB]dG(N1) and dC(NH4) amino proton and between [POB]dG(NH2) amino proton and dC(N3)). A hydrogen bond appears to occur between the POB carbonyl oxygen and the partner dC's second amino proton. The modified guanine purine base, partner cytosine pyrimidine base, and POB pyridyl ring form a triplex via this unusual hydrogen-bonding pattern. The phosphodiester backbone twists at the lesion site, accounting for the unusual phosphorus chemical shift differences relative to those for the control DNA duplex. The helical distortions and wobble base pairing induced by the covalent binding of POB to the O(6)-position of dG help explain the significant decrease of 17.6 degrees C in melting temperature of the modified duplex relative to the unmodified control.

Many carcinogens exert their cancer-causing effects by reacting with DNA either directly or following metabolic activation, resulting in covalently linked combination molecules known as carcinogen-DNA adducts. The presence of such lesions in the genome increases the error frequency of the replication machinery, causing mutations that contribute to the initiation and progression of cancer. Cellular DNA repair pathways remove carcinogen adducts from DNA, thus averting the mutagenic potential of many DNA lesions by reducing their presence in the genome. Bulky DNA adducts, like those derived from a number of activated environmental carcinogens such as polycyclic aromatic hydrocarbons (PAHs), are primarily repaired by the nucleotide excision repair (NER) pathway. Transcription-coupled NER (TC-NER) preferentially removes lesions from the transcribed strand of actively expressed genes, and RNA polymerase II stalled at the lesion quite possibly initiates the pathway. Among the bulky DNA adducts that are subject to TC-NER are those resulting from the reaction of the metabolically activated PAH benzo[a]pyrene (BP) with DNA. The P450 mixed-function oxygenases convert BP into a number of reactive intermediates, including tumorigenic (+)- and non-tumorigenic (−)-anti-benzo[a]pyrene diol epoxide (BPDE) that react with DNA via trans epoxide opening to form (+)-trans-anti-[BP]-N2-dG ((+)-ta[BP]G) and (−)-trans-anti-[BP]-N2-dG ((−)-ta[BP]G), respectively. To test the effect of these lesions on RNA synthesis, in vitro transcription assays using human nuclear extracts were performed with DNA templates containing an RNAPII promoter and a stereochemically pure (+)- or (−)-ta[BP]G adduct on the transcribed or non-transcribed strand. Transcription past (+)- or (−)-ta[BP]G adducts was investigated in the same sequence context to examine stereochemical effects. The (+)-ta[BP]G adduct was investigated in two different local sequence contexts to determine if the surrounding bases influence the adduct's ability to block transcription. These experiments revealed that (+)- and (−)-ta[BP]G adducts on the transcribed strand of the DNA template block RNAPII in a sequence and stereochemistry-dependent manner; however, adducts on the non-transcribed strand do not block elongation significantly but may increase pausing at innate pause sites. In order to elucidate biologically influential differences between the (+)- and (−)-ta[BP]G structures, the DUPLEX program was used to carry out potential energy minimization searches at model transcription junctions. The lowest-energy minimum for the (+)-ta[BP]G adduct gives a structure in which the benzo[a]pyrenyl ring system resides in the minor groove of the heteroduplex region. In contrast, the lowest-energy minimum for a (−)-ta[BP]G adduct shows an orientation in which the benzo[a]pyrenyl group adopts a carcinogen/base-stacked conformation. These conformational preferences may contribute to the differential treatment of (+)- and (−)-ta[BP]G adducts by human RNAPII. In addition, while previous experiments showed that BPDE adducts cause T7RNAP to produce a ladder of truncated transcripts, RNAPII is blocked entirely at only one or two positions by the (+)- and (−)-ta[BP]G adducts, depending on sequence context. It is likely that these differences between the behaviors of T7RNAP and human RNAPII are a result of the structural characteristics of the enzymes' active sites, a hypothesis that is explored in light of their known crystal structures.

In the Ames Salmonella typhimurium reversion assay 1,6- and 1,8-dinitropyrenes (1,6- and 1,8-DNPs) are much more potent mutagens than 1-nitropyrene (1-NP). Genetic experiments established that certain differences in the metabolism of the DNPs, which in turn result in increased DNA adduction, play a role. It remained unclear, however, if the DNP adducts, N-(guanin-8-yl)-1-amino-6 (8)-nitropyrene (Gua-C8-1,6-ANP and Gua-C8-1,8-ANP), which contain a nitro group on the pyrene ring covalently linked to the guanine C8, are more mutagenic than the major 1-NP adduct, N-(guanin-8-yl)-1-aminopyrene (Gua-C8-AP). In order to address this, we have compared the mutation frequency of the three guanine C8 adducts, Gua-C8-AP, Gua-C8-1,6-ANP, and GuaC8-1,8-ANP in a CGCG*CG sequence. Single-stranded M13mp7L2 vectors containing these adducts and a control were constructed and replicated in Escherichia coli. A remarkable difference in the induced CpG deletion frequency between these adducts was noted. In repair-competent cells the 1-NP adduct induced 1.7% CpG deletions without SOS, whereas the 1,6- and 1,8-DNP adducts induced 6.8 and 10.0% two-base deletions, respectively. With SOS, CpG deletions increased up to 1.9, 11.1, and 15.1% by 1-NP, 1,6-, and 1,8-DNP adducts, respectively. This result unequivocally established that DNP adducts are more mutagenic than the 1-NP adduct in the repetitive CpG sequence. In each case the mutation frequency was significantly increased in a mutS strain, which is impaired in methyl-directed mismatch repair, and a dnaQ strain, which carries a defect in proofreading activity of the DNA polymerase III. Modeling studies showed that the nitro group on the pyrene ring at the 8-position can provide additional stabilization to the two-nucleotide extrahelical loop in the promutagenic slipped frameshift intermediate through its added hydrogen-bonding capability. This could account for the increase in CpG deletions in the M13 vector with the nitro-containing adducts compared with the Gua-C8-AP adduct itself.

High-resolution nuclear magnetic resonance (NMR) and crystallographic data have been taken to refine the force field used in the torsion angle space nucleic acids molecular mechanics program DUPLEX. The population balance deduced from NMR studies of two carcinogen-modified DNA conformers in equilibrium was used to fine tune a sigmoidal, distance-dependent dielectric function so that reasonable relative energies could be obtained. In addition, the base-pair and backbone geometry from high-resolution crystal structures of the Dickerson-Drew dodecamer was used to re-evaluate the deoxyribose pseudorotation profile and the Lennard-Jones nonbonded energy terms. With a modified dielectric function that assumes a very steep distance-dependent form, a deoxyribose pseudorotation profile with reduced energy barriers between C2'- and C3'-endo minima, and a shift of the Lennard-Jones potential energy minimum to a distance approximately 0.4 A greater than the sum of the van der Waals' radii, the sequence-dependent conformational features of the Dickerson-Drew dodecamer in both the solid state and the aqueous liquid crystalline phase are well reproduced. The robust performance of the revised force field, in conjunction with its efficiency through implicit treatment of solvent and counterions, provides a valuable tool for elucidating conformations and structure-function relationships of DNA, including those of molecules modified by carcinogens and other ligands.

The carcinogenic heterocyclic amine (HA) 2-amino-1-methyl-6-phenylimidazo[4,5- b ]pyridine (PhIP) is formed during the cooking of various meats. To enable structure/activity studies aimed at understanding how DNA damaged by a member of the HA class of compounds can ultimately lead to cancer, we have determined the first solution structure of an 11-mer duplex containing the C8-dG adduct formed by reaction with N -acetoxy-PhIP. A slow conformational exchange is observed in which the PhIP ligand either intercalates into the DNA helix by denaturing and displacing the modified base pair (main form) or is located outside the helix in a minimally perturbed B-DNA duplex (minor form). In the main base-displaced intercalation structure, the minor groove is widened, and the major groove is compressed at the lesion site because of the location of the bulky PhIP- N -methyl and phenyl ring in the minor groove; this distortion causes significant bending of the helix. The PhIP phenyl ring interacts with the phosphodiester-sugar ring backbone of the complementary strand and its fast rotation with respect to the intercalated imidazopyridine ring causes substantial distortions at this site, such as unwinding and bulging-out of the strand. The glycosidic torsion angle of the [PhIP]dG residue is syn , and the displaced guanine base is directed toward the 3′ end of the modified strand. This study contributes, to our knowledge, the first structural information on the biologically relevant HA class to a growing body of knowledge about how conformational similarities and differences for a variety of types of lesions can influence protein interactions and ultimately biological outcome.

We report below on the solution structures of stereoisomeric “fjord” region trans-anti-benzo[c]phenanthrene-N2-guanine (designated (BPh)G) adducts positioned opposite cytosine within the (C-(BPh)G-C)·(G-C-G) sequence context. We observe intercalation of the phenanthrenyl ring with stereoisomer-dependent directionality, without disruption of the modified (BPh)G·C base-pair. Intercalation occurs to the 5′ side of the modified strand for the 1S stereoisomeric adduct and to the 3′ side for the 1R stereoisomeric adduct, with the S and R-trans-isomers related to one another by inversion in a mirror plane at all four chiral carbon atoms on the benzylic ring. Intercalation of the fjord region BPh ring into the helix without disruption of the modified base-pair is achieved through buckling of the (BPh)G·C base-pair, displacement of the linkage bond from the plane of the (BPh)G base, adaptation of a chair pucker by the BPh benzylic ring and the propeller-like deviation from planarity of the BPh phenanthrenyl ring. It is noteworthy that intercalation without base-pair disruption occurs from the minor groove side for S and R-trans-anti BPh-N2-guanine adducts opposite C, in contrast to our previous demonstration of intercalation without modified base-pair disruption from the major groove side for S and R-trans-anti BPh-N6-adenine adducts opposite T. Further, these results on fjord region 1S and 1R-trans-anti (BPh)G adducts positioned opposite C are in striking contrast to earlier research with “bay” region benzo[a]pyrene-N2-guanine (designated (BP)G) adducts positioned opposite cytosine, where 10S and 10R-trans-anti stereoisomers were positioned with opposite directionality in the minor groove without modified base-pair disruption. They also are in contrast to the 10S and 10R-cis-anti stereoisomers of (BP)G adducts opposite C, where the pyrenyl ring is intercalated into the helix with directionality, but the modified base and its partner on the opposite strand are displaced out of the helix. These results are especially significant given the known greater tumorigenic potential of fjord region compared to bay region polycyclic aromatic hydrocarbons. The tumorigenic potential has been linked to repair efficiency such that bay region adducts can be readily repaired while their fjord region counterparts are refractory to repair. Our structural results propose a link between DNA adduct conformation and repair-dependent mutagenic activity, which could ultimately translate into structure-dependent differences in tumorigenic activities. We propose that the fjord region minor groove-linked BPh-N2-guanine and major groove-linked BPh-N6-adenine adducts are refractory to repair based on our observations that the phenanthrenyl ring intercalates into the helix without modified base-pair disruption. The helix is therefore minimally perturbed and the phenanthrenyl ring is not available for recognition by the repair machinery. By contrast, the bay region BP-N2-G adducts are susceptible to repair, since the repair machinery can recognize either the pyrenyl ring positioned in the minor groove for the trans-anti groove-aligned stereoisomers, or the disrupted modified base-pair for the cis-anti base-displaced intercalated stereoisomers.

The non-steroidal anti-estrogen tamoxifen [TAM] has been in clinical use over the last two decades as a potent adjunct chemotherapeutic agent for treatment of breast cancer. It has also been given prophylactically to women with a strong family history of breast cancer. However, tamoxifen treatment has also been associated with increased endometrial cancer, possibly resulting from the reaction of metabolically activated tamoxifen derivatives with cellular DNA. Such DNA adducts can be mutagenic and the activities of isomeric adducts may be conformation-dependent. We therefore investigated the high resolution NMR solution conformation of one covalent adduct (cis-isomer, S-epimer of [TAM]G) formed from the reaction of tamoxifen [TAM] to N2-of guanine in the d(C-[TAM]G-C)·d(G-C-G) sequence context at the 11-mer oligonucleotide duplex level. Our NMR results establish that the S-cis [TAM]G lesion is accomodated within a widened minor groove without disruption of the Watson-Crick [TAM]G·C and flanking Watson-Crick G·C base-pairs. The helix axis of the bound DNA oligomer is bent by about 30° and is directed away from the minor groove adduct site. The presence of such a bulky [TAM]G adduct with components of the TAM residue on both the 5′- and the 3′-side of the modified base could compromise the fidelity of the minor groove polymerase scanning machinery.

A solution structural study has been undertaken on the aminofluorene-C8-dG ([AF]dG) adduct located at a single-strand-double-strand d(A1-A2-C3-[AF]G4-C5-T6-A7-C8-C9-A10-T11-C12-C13). d(G14-G15-A16-T17-G18-G19-T20- A21-G22-N23) 13/10-mer junction (N = C or A) using proton-proton distance restraints derived from NMR data in combination with intensity-based relaxation matrix refinement computations. This single-strand-double-strand junction models one arm of a replication fork composed of a 13-mer template strand which contains the [AF]dG modification site and a 10-mer primer strand which has been elongated up to the modified guanine with either its complementary dC partner or a dA mismatch. The solution structures establish that the duplex segment retains a minimally perturbed B-DNA conformation with Watson-Crick hydrogen-bonding retained up to the dC5.dG22 base pair. The guanine ring of the [AF]dG4 adduct adopts a syn glycosidic torsion angle and is displaced into the major groove when positioned opposite dC or dA residues. This base displacement of the modified guanine is accompanied by stacking of one face of the aminofluorene ring of [AF]dG4 with the dC5.dG22 base pair, while the other face of the aminofluorene ring is stacked with the purine ring of the nonadjacent dA2 residue. By contrast, the dC and dA residues opposite the junctional [AF]dG4 adduct site adopt distinctly different alignments. The dC23 residue positioned opposite the adduct site is looped out into the minor groove by the aminofluorene ring. The syn displaced orientation of the modified dG with stacking of the aminofluorene and the looped out position of the partner dC could be envisioned to cause polymerase stalling associated with subsequent misalignment leading to frameshift mutations in appropriate sequences. The dA23 residue positioned opposite the adduct site is positioned in the major groove with its purine ring aligned face down over the van der Waals surface of the major groove and its amino group directed toward the T6.A21 base pair. The Hoogsteen edge of the modified guanine of [AF]dG4 and the Watson-Crick edge of dA23 positioned opposite it are approximately coplanar and directed toward each other but are separated by twice the hydrogen-bonding distance required for pairing. This structure of [AF]dG opposite dA at a model template-primer junctional site can be compared with a previous structure of [AF]dG opposite dA within a fully paired duplex [Norman, D., Abuaf, P., Hingerty, B. E., Live, D. , Grunberger, D., Broyde, S., and Patel, D. J. (1989) Biochemistry 28, 7462-7476]. The alignment of the Hoogsteen edge of [AF]dG (syn) positioned opposite the Watson-Crick edge of dA (anti) has been observed for both systems with the separation greater in the case of the junctional alignment in the model template-primer system. However, the aminofluorene ring is positioned in the minor groove in the fully paired duplex while it stacks over the junctional base pair in the template-primer system. This suggests that the syn [AF]dG opposite dA junctional alignment can be readily incorporated within a duplex by a translation of this entity toward the minor groove.

While there is now available a great deal of information on double stranded DNA from X-ray crytallography, high resolution NMR and computer modeling, very little is known about structures that are representative of the DNA core of replication intermediates. DNA replication occurs at a single strand/double strand junction and bulged out intermediates near the junction can lead to frameshift mutations. The single stranded domains are particularly challenging. Our interest is focused on strategies for modeling the DNA of these types of replication intermediates. Modeling such structures presents special problems in addressing the multiple minimum problem and in treating the electrostatic component of the force field. We are testing a number of search strategies for locating low energy structures of these types and we are also investigating two different distance dependent dielectric functions in the coulombic term of the force field. We are studying both unmodified DNA and DNA damaged by aromatic amines, carcinogens present in the environment in tobacco smoke, barbecued meats and automobile exhaust. The nature of the structure adopted by the carcinogen modified DNA at the replication fork plays a key role in determining whether the carcinogen will cause a mutation during replication that can initiate the carcinogenic process. In the present work results are presented for unmodified DNA in the NarI mutational hotspot sequence. Interestingly, we find that a slipped bulged mutagenic intermediate is much more stable relative to its normally extended counterpart when the guanine mutational hotspot is syn instead of anti. This suggests that carcinogens such as 2-(acetyl)aminofluorene which induce the syn conformation enhance mutagenesis at the hotspot by stabilizing a slipped intermediate in the syn orientation.

While the one-ring amine aniline (AN) has only slight genetic activity, the polycyclic aromatic amines 2-aminofluorene (AF) and 1-aminopyrene (AP) are significant mutagens and carcinogens. Moreover, the bulkier AP is more mutagenic per adduct than AF in the tetracycline-resistance gene of plasmid pBR322 [Melchior et al. (1994) Carcinogenesis 15, 889]. To elucidate possible conformational origins of the differing mutagenic effects of these three adducts, which may stem from their differing ring sizes, we have examined their conformations in two mutation-susceptible sequences from the above gene: TTGAG*GCCG (sequence I) and GAATG*GTGC (sequence II), where G* = C8-modified guanine. No experimental high-resolution NMR data are yet available for the aniline adduct in a DNA duplex. Minimized potential energy calculations were carried out, using the molecular mechanics program DUPLEX to explore the conformation space of these adducts. In the case of AN, a relatively unperturbed B-DNA helix with the amine in the major groove was strongly favored in both sequences. In the case of AF- and AP-modified DNA, however, several differing conformations were competitive in energy. They included major groove structures, as well as conformations with syn-modified guanine and the polycyclic amine in the minor groove, or the amine rings intercalated into the helix with displacement of the modified guanine, in overall harmony with high-resolution NMR solution structures. Thus, aniline distorts DNA structure to a lesser extent than larger aromatic amine ring systems, since a number of different conformations are energetically feasible and have been observed for the larger systems. This result may be relevant to their enhanced mutagenicity and their repair propensity, in contrast to aniline's low mutagenic effect.

A molecular dynamics simulation has been carried out with DNA polymerase beta (beta pol) complexed with a DNA primer-template. The templating guanine at the polymerase active site was covalently modified by the carcinogenic metabolite of benzo[a]pyrene, (+)-anti-benzo[a]pyrene diol epoxide, to form the major (+)-trans-anti-benzo[a]pyrene diol epoxide covalent adduct. Thus, the benzo[a]pyrenyl moiety (BP) is situated in the single-stranded template at the junction between double- and single-stranded DNA. The starting structure was based on the X-ray crystal structure of the rat beta pol primer-template and ddCTP complex [Pelletier, H., Sawaya, M. R., Kumar, A., Wilson, S. H., and Kraut, J. (1994) Science 264, 1891-1903]. During the simulation, the BP and its attached templating guanine rearrange to form a structure in which the BP is closer to parallel with the adjacent base pair. In addition, the templating attached guanine is displaced toward the major groove side and access to its Watson-Crick edge is partly obstructed. This structure is stabilized, in part, by new hydrogen bonds between the BP and beta pol Asn279 and Arg283. These residues are within hydrogen bonding distance to the incoming ddCTP and templating guanine, respectively, in the crystal structure of the beta pol ternary complex. Site-directed mutagenesis has confirmed their role in dNTP binding, discrimination, and catalytic efficiency [Beard, W. A., Osheroff, W. P., Prasad, R., Sawaya, M. R., Jaju, M., Wood, T. G., Kraut, J., Kunkel, T. A., and Wilson, S. H. (1996) J. Biol. Chem. 271, 12141-12144]. The predominant biological effect of the BP is DNA polymerase blockage. Consistent with this biological effect, the computed structure suggests the possibility that the BP's main deleterious impact on DNA synthesis might result at least in part from its specific interactions with key polymerase side chains. Moreover, relatively modest movement of BP and its attached guanine, with some concomitant enzyme motion, is necessary to relieve the obstruction and permit the observed rare incorporation of a dATP opposite the guanine lesion.

The Escherichia coli NarI restriction enzyme recognition site 5'G1G2C3G4C5C63' is a mutational hotspot for -2 deletions in E. coli plasmid pBR322, resulting in the sequence 5'GGCC3' when G4 is modified by the aromatic amine N-2-(acetyl)aminofluorene (AAF) [Burnouf, D., Koehl, P., and Fuchs, R. P. P. (1995) Proc. Natl. Acad. Sci. U.S.A. 86, 4147-4151] even though each G shows similar reactivity [Fuchs, R. P. P. (1984) J. Mol. Biol. 177, 173-180]. Modification at G4 by the related aromatic amine 2-aminofluorene (AF), which lacks the acetyl group of AAF, can also cause -2 deletions, but at a lower frequency [Bichara, M., and Fuchs, R. P. P. (1985) J. Mol. Biol. 183, 341-351]. A specific mechanism has been proposed to explain the double-base frameshifts in the NarI sequence in which the GC deletion results from a slipped mutagenic intermediate formed during replication [Schaaper, B. M., Koffel-Schwartz, N., and Fuchs, R. P. P. (1990) Carcinogenesis 11, 1087-1095]. We address the following key questions in this study. Why does AAF modification dramatically increase the mutagenicity at the NarI G4 position, and why does AAF enhance the mutagenicity more than AF? We studied two intermediates which model replication at one arm of a fork, using a fragment of DNA modified by AF or AAF at G4 in the NarI sequence: Intermediate I can be converted into intermediate II by misalignment. Elongation of intermediate I leads to error-free translesion synthesis, while elongation of intermediate II leads to a -2 frameshift mutation. Minimized potential energy calculations were carried out using the molecular mechanics program DUPLEX to investigate the conformations of the AF and AAF adducts at G4 in these two intermediates. We find that the slipped mutagenic intermediate is quite stable relative to its normally extended counterpart in the presence of AF and AAF in an abnormal syn orientation of the damaged base. An enhanced probability of elongation from a stable slipped structure rather than a properly aligned one would favor increased -2 frameshift mutations. Furthermore, AAF-modified DNA has a greater tendency to adopt the syn orientation than AF because of its greater bulk, which could explain its greater propensity to cause -2 deletions in the NarI sequence.

This paper addresses structural issues related to the capacity of aminofluorene [AF] for frameshift mutations of the -2 type on C-8 covalent adduct formation at the G(3) site in the d(C-G(1)-G(2)-C-G(3)-C-C) NarI hot spot sequence. This problem has been approached from a combined NMR and relaxation matrix analysis computational structural study of the [AF]dG adduct in the d(C-G-G-C-[AF]G-C-C).d(G-G-C-C-G) sequence context at the 12/10-mer adduct level (designated [AF]dG.del(-2) 12/10-mer). The proton spectra of this system are of exceptional quality and are consistent with the formation of an AF-intercalated conformer with the modified guanine in a syn alignment displaced along with the 5'-flanking cytosine residue into the major groove. The solution structure has been determined by initially incorporating intramolecular and intermolecular proton-proton distances defined by lower and upper bound deduced from NOESY spectra as restraints in molecular mechanics computations in torsion angle space and subsequently refined through restrainted molecular dynamics calculations based on a NOE distance and intensity refinement protocol. Strikingly, the [AF]dG.del(-2) 12/10-mer duplex adopts only one of two potential AF-intercalation alignments for the [AF]dG adduct opposite the -2 deletion site in the NarI sequence context with the extrusion of the dC-[AF]dG step favored completely over extrusion of the [AF]dG-dC step at the lesion site. This polarity establishes that the structural perturbation extends 5' rather than 3' to the [AF]dG lesion site in the adduct duplex. This structure of the [AF]dG adduct opposite a -2 deletion site shows distinct differences with conclusions reported on the alignment of the related acetylaminofluorene [AAF]dG adduct opposite a -2 deletion site in the identical NarI sequence context [Milhe, C., Fuchs, R. P. P., and Lefevre, J. F. (1996) Eur. J. Biochem. 235, 120-127]. In that study, qualitative NMR data without computational analysis were employed to conclude that the extrusion at the lesion site occurs at the [AAF]dG-dC step for the AAF-intercalated conformer of the adduct duplex. The structure of the [AF]dG adduct opposite a -2 deletion site determined in our group provides molecular insights into the architecture of extended slipped mutagenic intermediates involving aromatic amine intercalation and base-displaced syn modified guanines in AF and, by analogy, AAF-induced mutagenesis in the NarI hot spot sequence context.