{'fraction_non_alphanumeric': 0.03963480128893663, 'fraction_numerical': 0.032545649838882924, 'mean_word_length': 4.54924750409775, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Exocyclic amino groups of flanking guanines govern sequence-dependent adduct conformations and local structural distortions for minor groove-aligned benzo[a]pyrenyl-guanine lesions in a GG mutation hotspot context

The environmental carcinogen benzo[a]pyrene (BP) is metabolized to reactive diol epoxides that bind to cellular DNA by predominantly forming N2-guanine adducts (G*). Mutation hotspots for these adducts are frequently found in 5′- ··· GG ··· dinucleotide sequences, but their origins are poorly understood. Here we used high resolution NMR and molecular dynamics simulations to investigate differences in G* adduct conformations in 5′- ··· CG*GC ··· and 5′- ··· CGG* C··· sequence contexts in otherwise identical 12-mer duplexes. The BP rings are positioned 5′ along the modified strand in the minor groove in both cases. However, subtle orientational differences cause strong distinctions in structural distortions of the DNA duplexes, because the exocyclic amino groups of flanking guanines on both strands compete for space with the BP rings in the minor groove, acting as guideposts for placement of the BP. In the 5′- ··· CGG* C ··· case, the 5′-flanking G · C base pair is severely untwisted, concomitant with a bend deduced from electrophoretic mobility. In the 5′- ··· CG*GC ··· context, there is no untwisting, but there is significant destabilization of the 5′-flanking Watson–Crick base pair. The minor groove width opens near the lesion in both cases, but more for 5′- ··· CGG*C···. Differential sequence-dependent removal rates of this lesion result and may contribute to the mutation hotspot phenomenon.


INTRODUCTION
Polycyclic aromatic hydrocarbons (PAHs) are products of combustion of organic matter and are therefore ubiquitous in the environment, and have long been associated with induced human cancers (1,2). Benzo [a]pyrene (BP) is the most widely studied and representative compound of this class of chemical carcinogens (3). While BP is relatively inert, it is metabolically activated to highly reactive and genotoxic diol epoxide derivatives (4) that bind chemically to DNA and form mutagenic adducts (5). These DNA lesions have been correlated with the initiation of animal and human cancers (6). Among the most active and tumorigenic metabolites of BP is the (þ)-7R,8S,9S,10R enantiomer of 7,8-dihydroxy-9, 10-epoxy-7,8,9,10-tetrahydrobenzo[a]pyrene ((þ)-anti-BPDE); this metabolite binds predominantly to the exocyclic amino group of guanine in DNA to form the 10S (þ)-trans-anti- [BP]-N 2 -dG adduct shown in Figure 1 (7). If DNA adducts are not excised by cellular repair mechanisms, they can persist until replication occurs and cause mutations. Multiple DNA mutations in critically important genes such as ras and p53 constitute genetic alterations that play key roles in the regulation of cell cycle control and cancer (1).
The diol epoxide (þ)-anti-BPDE is highly mutagenic (8) and tumorigenic (9,10). Mutation hotspots associated with adduct formation derived from the reactions of BPDE with DNA are frequently found in runs of guanines, especially in GG sequence contexts (11)(12)(13)(14). Examples of GG mutation hotspot phenomena include mutations in plasmids containing human c-Ha-ras1 sequences, reacted with anti-BPDE and transfected into mouse embryonic NIH 3T3 cells (15) in the coding region of the BPDE-modified HPRT gene in mammalian V-79 cells (14), and in the SupF gene of an E. coli plasmid (13). The mutational specificity and the mutation frequencies depend not only on base sequence contexts (5), but also on the nature of the host or the polymerases bypassing the lesion (16). Furthermore, in cellular environments, nucleotide excision repair of these bulky lesions may also depend on base sequence contexts and thus contribute to the mutation hotspot phenomena, by efficiently excising lesions in some sequence contexts but not in others. However, the structural basis of base sequence effects on mutagenesis and DNA repair is poorly understood. Although the conformations of 10S (þ)-trans-anti- [BP]-N 2 -dG adducts in double-stranded oligonucleotides have been studied by NMR methods (17,18), previous emphasis has been on understanding stereochemical rather than base sequence effects. In this work, we have used high resolution NMR methods to investigate the conformational properties of the mutagenic 10S (þ)-trans-anti- [BP]-N 2 -dG adduct ( Figure 1A) at either of the two guanines in a GG sequence context in a double-stranded oligonucleotide ( Figure 1B). The central 6-mer sequence of this duplex contains the SupF mutation hotspot (13) 5 0 -ÁÁÁ GCGGCC ÁÁÁ with the BP lesion positioned either at the first (5 0 -side) or the second G (3 0 -side). In both sequence contexts, the BP rings are positioned 5 0 -directed along the modified strand in the minor groove (17,19), but with subtle differences in orientations which produce markedly different local duplex distortions. These distinctions arise from differential steric hindrance from the exocyclic amino groups of the neighboring undamaged guanine bases that compete with the bulky aromatic BP ring systems for space in the minor groove. This leads to significant differences in local helical untwisting, minor groove width opening near the lesion site, and weakening of local Watson-Crick base pairs. The unusual flexible bend observed in duplexes with the 10S (þ)-trans-anti- [BP]-N 2 -dG adduct positioned only on the 3 0 -G in the 5 0 -ÁÁÁ GG ÁÁÁ sequence context (20) is consistent with severe local untwisting found only for this sequence. The pronounced differences in the local structural distortions affect the processing of the lesions in GG mutation hotspots by cellular enzymes.
were recorded at mixing times of 50, 100, 150 and 200 ms using a Bruker 800 MHz (for the G6*G7 sample) and a 600 MHz (for G6G7* sample) NMR spectrometers, equipped with cryoprobes, at the New York Structural Biology Center (NYSBC). The NOESY experiments in H 2 O at mixing times of 90, 175 and 250 ms at 58C to visualize the imino protons were conducted using a Bruker 500 MHz NMR instrument at New York University. Full details are provided in Supplementary Data.

Molecular dynamics computations and intensity refinement
Starting models. Initial models were created from a high resolution NMR solution structure for the 10S (þ)-trans-anti- [BP]-N 2 -dG adduct (17). We remodeled the base sequence to the two sequences listed in Figure 1B, using INSIGHT II 97.0 (Accelrys, Inc.).
Force field. MD simulations and intensity refinement against the NOE distances and volumes were carried out using SANDER in the AMBER 8.0 simulation package (University of San Francisco), the Cornell et al. force field (24), with the parm99.dat parameter set (25). Partial charges and all parameters for the 10S (þ)-trans-anti- [BP]-N 2 -dG adduct employed in this work were computed previously (26).
Molecular dynamics computation protocols and best representative structures. Details of the molecular dynamics protocols and the computation of the best representative structure, using the cluster analyses option in MOIL-View (27) are provided in Supplementary Data. NMR intensity refinement and structural analyses. The refinement started with the best representative structure of the final 1.5 ns unrestrained MD simulation. This structure was subjected to distance-restrained, and subsequently, intensity-restrained molecular dynamics computations in order to elucidate the solution structures of the 10S (þ)-trans-anti- [BP]-N 2 -dG adduct in the G6*G7 and G6G7* sequence contexts. The SANDER module in AMBER 8.0 was employed. Details of the structural analyses are provided in Supplementary Data.

Exchangeable proton spectra
The G6*G7 duplex. The exchangeable 1D proton NMR spectrum (12.0-14.0 ppm) of the G6*G7 duplex (sequence context shown in Figure 1B) in H 2 O buffer solution (pH ¼ 6.8) at 58C is plotted in Figure 2A. The imino protons were assigned via their NOE connectivities to the amino protons of the hydrogen bond partner and the flanking base pairs, as described in earlier publications (17,(28)(29)(30)(31). The observed distinguishable imino proton resonances indicate a single adduct conformation. An expanded contour plot of the NOESY spectrum in H 2 O solution determined at 250 ms mixing time at 58C is shown in Figure 2B. NOEs between deoxyguanosine imino and cytosine amino protons are characteristic of a G Á C Watson-Crick basepair (12.4-13.1 ppm), while NOEs between thymidine imino and adenine H2 protons are characteristic of an A Á T Watson-Crick base pair (13.4-13.9 ppm). The observed NOE pattern establishes Watson-Crick pairing for G Á C pairs (peaks labeled A to C) and A Á T pairs (peaks labeled E and F), including the G6* Á C19 base pair, in the 12-mer duplex. However, somewhat weaker NOE cross-peaks are observed at the lesioned G6* Á C19 base pair and the 5 0 -flanking C5* Á G20 base pair. Imino-imino proton sequential connectivities can be traced throughout the duplex except at G20(NH1)-G6*(NH1), G7(NH1)-G17(NH1) and G17(NH1)-G16(NH1). Some of these connectivities cannot be observed, presumably because they are too close to the diagonal, while imino protons at the terminal base pairs are difficult to observe due to fraying effects. None of the imino proton resonances are upfield shifted (Figure 2A), which indicates that the aromatic BP ring system is not intercalated between adjacent base pairs in the duplex (28,29) ( Figure 2C). An observable BP(H8)-G6*(H1) cross-peak is also observed in this region of the spectrum (labeled D in Figure 2B).
The G6G7* duplex. The 1D imino proton NMR spectrum (12.0-14.0 ppm) in H 2 O buffer of the G6G7* duplex, together with the assignments, is depicted in Figure 2D.
The imino proton resonances are distinguishable and establish the presence of a single conformation in this duplex. An expanded contour plot of the NOESY spectrum (250 ms mixing time) in H 2 O buffer at 58C is shown in Figure 2E. Characteristic interstrand Watson-Crick NOE cross-peaks are observed between guanine imino and cytosine amino protons of G Á C base pairs (peaks labeled A to C), and thymine imino and adenine H2 protons (peaks labeled E to H). Although Watson-Crick hydrogen bonding is observed at the modified base pair (peak A in Figure 2E) and the two flanking base pairs G6 Á C19 (peak B) and C8 Á G17 (peak C), the NOE at the lesion site (G7* Á C18), is relatively weak (peak A) when compared to the NOE connectivities between other base pairs. This implies that hydrogen bonding at the lesion site is weakened. A similar conclusion can be drawn for the 5 0 -flanking G6 Á C19 base pair where hydrogen bonding is also diminished. On the other hand, the amino-imino NOE cross-peak for the C8 Á G17 pair is quite strong (peak C), most likely reflecting a unidirectional orientation of the BP ring system. Imino-imino proton sequential assignments can be traced for the entire duplex except at the G6(NH1)-G7*(NH1), G7*(NH1)-G17(NH1) and the terminal base pairs. None of the imino proton resonances are upfield shifted ( Figure 2D) suggesting that the aromatic BP ring system is not intercalated between adjacent base pairs in this duplex.
Temperature dependence of imino proton resonances in G6*G7 and G6G7* duplexes Imino proton chemical shifts and linewidth broadening of the central three base pair segments around the lesion as a function of temperature show interesting differences between the two adducts ( Figure S1, Supplementary Data). Increases in linewidths of the imino protons reflect enhanced solvent exposure, and thus increased exchange rates with solvent water molecules. These phenomena are therefore indicators of the relative stability of the local duplex environment (32). As the temperature is increased from 5 to 208C, there is a pronounced broadening of the imino proton resonances of G6*, G7 and G20; the latter two bases belong to base pairs flanking G6* Á C19 in the G6*G7 duplex. However, in the case of the G6G7* duplex, the changes in the linewidths of imino protons of G7*, and neighboring bases G6 and G17 are less pronounced ( Figure S1A, Supplementary Data). Upfield chemical shift changes of imino protons are attributed to a premelting conformational transition associated with increased duplex opening rates (32). We have examined the chemical shift values of the imino protons at the lesion site and the two adjacent G Á C base pairs as a function of temperature ( Figure S1B, Supplementary Data). In the case of the G6G7* duplex, rather small and monotonous upfield chemical shifts of the imino proton resonances are observed at the site of the G7* lesion and the two flanking base pairs as the temperature increases from 5 to 458C ( Figure S1B, Supplementary Data). However, in the G6*G7 duplex, pronounced upfield chemical shift changes occur at the 5 0 -side of G6*, i.e. at the C5 Á G20 base pair and, to a lesser extent, at the G7 Á C18 base pair flanking the lesion on both sides in the G6*G7 duplex ( Figure S1B, Supplementary Data). All three base pairs, C5 Á G20, G6* Á C19, and G7 Á C18, appear to premelt at temperatures of $20-308C, well below the global melting point, T m , of either of the two duplexes determined from measurements of the UV absorbance as a function of temperature ( Figure S1C, Supplementary Data). When the absorbance is monitored at the DNA absorption maximum of 260 nm, the melting curves are cooperative and the T m values are 57 AE 1 and 55 AE 18C for the G6G7* and G6*G7 duplexes, respectively. When the temperaturedependent absorbance is monitored at 346 nm, corresponding to the absorption maximum of the aromatic BP rings, the local premelting of the duplexes in the immediate vicinity of the lesions can be monitored. The local T m is 50 AE 18C in both duplexes, which is significantly lower than the global melting point for the entire duplex measured at 260 nm. It is evident that the region of the duplex in the immediate vicinity of the lesions is destabilized by the 10S (þ)-trans-anti- [BP]-N 2 -dG adducts, causing local premelting of the doublestranded DNA. The pronounced upfield shift changes of the G20(H1), G6*(H1) and G7(H1) imino protons in the G6*G7 duplex suggest that significant local changes occur at temperatures much lower than the global T m , i.e. between 20 and 308C. Changes in local duplex structure associated with increased opening rates occur at these relatively low temperatures. The lack of similar changes in linewidths and chemical shifts of the analogous base pairs in the G6G7* duplex, indicates that the G6 Á C19, G7* Á C18, and C8 Á G17 are more stable in the same 20-308C temperature interval ( Figure S1A and B, Supplementary Data).

Non-exchangeable proton spectra
The G6*G7 duplex. A portion of the 2D NOESY contour plot (mixing time 250 ms) of the G6*G7 duplex in 2 H 2 O buffer solution (pH ¼ 6.8) at 158C is plotted in Figure 3A. The characteristic NOE patterns between base protons and their own and 5 0 -flanking base sugar H1 0 protons for a right-handed B-DNA duplex (33) are traced for the C5-G6*-G7-C8 segment of the modified strand (solid lines) and for the G17-C18-C19-G20-C21-A22 segment of the non-modified strand (broken lines). Internucleotide base to sugar NOEs could not be identified unequivocally in the G20-C21 segment of the complementary strand because of severe overlap of the C19(H6), G20(H8) and T23(H6) base protons. We tentatively assigned the G20(H1 0 ) and G20(H3 0 ) sugar protons at 4.18 and 3.52 ppm, respectively. These protons exhibit NOE cross-peaks to the G20(H2 0 ) and G20(H2 00 ) sugar protons. Large upfield shifts of the base (6.77 ppm) and sugar H1 0 (4.55 ppm) proton of C21, as well as similar upfield shifts for all other sugar protons of G20 and C21 ( Figure 3B), indicate that the aromatic BP ring system of G6* is in contact with the G20-C21 sugar rings, and is therefore positioned in the minor groove. These results also indicate that the aromatic BP ring system is oriented toward the 5 0 -end of the modified strand in the G6*G7 duplex. The aromatic BP protons, assigned on the basis of comparisons of NOESY and TOCSY patterns for these protons and their chemical shifts, are shown in Table 1.
The G6G7* duplex. Representative 2D NOESY spectra of the G6G7* duplex (mixing time 250 ms) in 2 H 2 O buffer solution (pH ¼ 6.8) at 158C are shown in Figure 4A. The NOE connectivities between base protons and their own and 5 0 -flanking H1 0 sugar protons are traced for the C5-G6-G7*-C8 segment of the modified strand (solid lines), and for the G17-C18-C19-G20-C21-A22 segment of the unmodified strand (dashed lines). All connectivities are present, confirming that the G6G7* duplex is consistent with a B-DNA helix conformation. The upfield shift observed for the C19 base proton (6.71 ppm) and its sugar H1 0 proton (3.87 ppm) indicate that the aromatic BP ring stacks mainly over the sugar ring of the C19 residue of the G6 Á C19 base pair flanking the lesion G7* on the 5 0 -side of the modified strand. The BP aromatic ring system extends less toward the G20 sugar protons, since smaller upfield chemical shifts are observed for this base than for C19 ( Figure 4B). The aromatic BP protons were assigned on the basis of comparisons of NOESY and TOCSY patterns for the indicated protons, and their chemical shift values are shown in Table 1.
Comparisons of nucleic acid chemical shifts in the G6*G7 and G6G7* duplexes In order to compare the effects of the aromatic BP rings on chemical shifts of DNA proton resonances, plots of chemical shifts, relative to those of the unmodified control duplexes, are presented in Figures 3B and 4B. The relative chemical shifts for sugar protons H1 0 , H2 0 , H2 00 and H3 0 , and H6/H8 base protons are generally unperturbed except in the vicinity of the lesion (C5-G7), and up to three bases on the 5 0 -side of the lesion in the complementary strand. These observations are consistent with a 5 0 -orientation of the BP moiety in both the G6*G7 and G6G7* duplexes, with the aromatic ring system contacting mostly sugar protons of the complementary strand in the minor groove (17). However, there are significant differences. In the case of the G6G7* duplex, the largest upfield shifts, caused by aromatic ring current effects associated with the BP rings, are observed at the sugar protons of the 5 0 -adjacent base C19 ( Figure 4B). However, in the case of G6*G7, the effects of the BP rings extend to protons on the 5 0 -flanking G20 and C21 connectivities are shown between a base (purine H8 and pyrimidine H6) proton and their own and 5 0 -flanking sugar H1 0 protons from dC5 to dC8 on the modified strand (solid arrows) and from dG17 to dA22 on the non-modified strand (broken arrows). The NOEs between the 2.45 Å fixed distance H6 and H5 of the cytosine are designated by crosses. Internucleotide base to sugar NOEs, at the non-modified strand, were not identified clearly for G20 ! C21 as a result of the overlap of the C19-H6, G20-H8 and T23-H6 protons. There is a tentative assignment of the G20(H1 0 ) and G20(H3 0 ) sugar protons at 4.18 and 3.52 ppm, respectively through NOEs between theses protons and the G20(H2 0 ) and G20(H2 00 ) sugar protons. in the complementary strand, up to two bases from the modification site (G6*) ( Figure 3B). These results suggest that the long axes of the BP aromatic ring systems in the minor grooves are oriented differently in the G6*G7 and G6G7* duplexes.
BP-DNA proton NOE cross-peaks. We have identified and assigned all BP aliphatic and aromatic protons in both 12-mer duplexes. A total of 12 NOE cross-peaks between the BP rings and DNA protons were identified for the G6*G7 sample. Several of the NOEs between minor groove sugar protons and the BP rings are shown in Figure 3A and C for the G6*G7 duplex. The intermolecular NOE connectivities involve protons on both the modified and the unmodified strands. In the case of the G6G7* duplex, a total of 14 intermolecular NOE crosspeaks between the BP rings and DNA protons were identified. Several of these intermolecular NOEs are shown in Figure 4A and C, and involve sugar and BP aromatic ring protons. A list of all of the BP-DNA intermolecular NOE cross-peaks for the G6*G7 and the G6G7* duplexes are summarized in Table 1. These intermolecular connectivities provide further evidence that the BP aromatic ring systems are positioned in the minor groove in both duplexes and are directed toward the 5 0 -side of the modified strand. The results are fully consistent with the chemical shift data ( Figures 3B  and 4B) showing that the BP aromatic ring system stacks predominantly over the G20 and C21 sugar rings in the G6*G7 duplex, and over the C19 sugar ring in the case of the G6G7* duplex.
Determination of the conformations of the 10S (þ)-trans-anti- [BP]-N 2 -dG adducts in the G6*G7 and G6G7* duplexes NMR distance and intensity-restrained molecular dynamics computations were employed to define the solution structures of the 10S (þ)-trans-anti- [BP]-N 2 -dG adducts in the G6*G7 and G6G7* sequence contexts at the 12-mer duplex level. Details concerning the initial models and molecular dynamics computation protocols for unrestrained and distance and intensity-restrained MD simulations are given in the Materials and Methods section. The restraints and refinement statistics are listed in Table 2.

Solution structures
Stick views of five superpositioned intensity-refined structures of the 10S (þ)-trans-anti- [BP]-N 2 -dG adducts in G6*G7 and G6G7* sequence contexts are shown in Figure 5. The views shown are into the minor grooves of the G6*G7 ( Figure 5A) and the G6G7* duplexes ( Figure 5B). In both cases, the BP moiety is in the B-DNA minor groove, 5 0 -directed along the modified strand. However, there are subtle differences in the orientations of the aromatic BP rings within the minor groove of the DNA. Figure 6 shows the origin of this phenomenon. In each case, the BP rings are positioned to avoid collisions with adjacent exocyclic amino groups of the flanking guanine bases that are protruding into the minor groove. In the case of the G6*G7 duplex, the relevant amino groups are those belonging to G7 on the 3 0 -side of the lesion G6* and positioned on the same  strand, and G20 on the complementary strand. These two amino groups are in contact with the aliphatic ring and the aromatic ring system of BP, respectively ( Figure 6A). In G6G7*, the relevant amino groups are those belonging to G6 on the 5 0 -side to the lesion G7* and positioned on the same strand, and G17 on the complementary strand. These two amino groups are also in contact with the aromatic ring and the aliphatic ring system of BP, respectively. The difference between the two cases is that the amino group of G6 flanking the modified guanine G7* is on the same strand and on the 5 0 -side of G7* ( Figure 6B), while in the G6*G7 duplex the 5 0 -flanking amino group is part of G20 on the complementary strand ( Figure 6A). These variations in positions of the amino groups relative to the modified guanines produce steric hindrance differences; subtle orientational differences result, manifested by distinctions in the distribution of the torsion angles at the carcinogen-DNA linkage site a 0 , b 0 and the glycosidic torsion angles w ( Figure 7A and Table 3). As a result, the G6*G7 and G6G7* duplexes have distinct structural properties. Most notably, there are significant differences in the Twist angles flanking the modified guanines in the G6*G7 and G6G7* duplexes ( Figure 7B). In the G6*G7 duplex, the Twist angle varies within the narrow range of 29-328 from base pair A2 Á T23 to base pair T10 Á A15. However, in the case of the G6G7* duplex, there is a significantly larger Twist (408) at the modified base pair G7* Á C18, and a strikingly lower Twist angle of only 138 at the 5 0 -flanking base pair G6 Á C19 ( Figure 7B). These results indicate that there is a significant impact of guanine G6 flanking the modified guanine G7* on the 5 0 -side, since the base pair G7* Á C18 is over-twisted while the 5 0 -flanking unmodified G6 Á C19 base pair is severely undertwisted. Undertwisting is accompanied by increased Roll (Figure S2, Supplementary Data). In addition, there are differences in the minor groove dimensions near the lesion ( Figure 7C), with average values of 10.6 Å for P7-P22 (defined in Figure 7C) in the G6*G7 sequence context, and 12.3 Å for P8-P21 in the G6G7* duplex.
The amino groups governing the position of the BP moiety in the minor groove facilitate extra hydrogen bonding between the hydroxyl groups of the BP benzylic ring and O2 atoms of cytosines in the immediate vicinity of the lesion; this is seen in both the G6*G7 and the G6G7* sequence contexts in the computed structures ( Figure S3, Supplementary Data). These cytosines are either partner bases to the modified guanines, or their 3 0 side-flanking neighbors. In the case of G6*G7, there is one hydrogen bond between the HO7 and O2 of C19, the partner to the modified G6*. (Figure S3A, Supplementary Data). In the G6G7* case, there are two hydrogen bonds: one between the HO7 of the benzylic ring and O2 of C18, the partner base to the modified guanine, G7*; the other is between the HO8 and O2 of C8, on the 3 0 -side of G7* ( Figure S3B, Supplementary Data). These intermolecular hydrogen bonds between the benzylic ring and cytosines further anchor the BP moiety in the position governed by the guanine amino groups. Thus the  Figure 5. View into the minor groove of superpositioned intensityrefined structures of the (A) G6*G7 and (B) G6G7* sequence contexts at the 12-mer duplex level. The structures shown are the five best representative conformations for each sequence context from the final 1 ps of unrestrained MD simulation after intensity refinement. The BP rings are in yellow, the modified base guanine and its partner cytosine are in cyan, and the rest of the DNA is in white except for the phosphorus atoms, which are colored red.
differential orientation of the BP in the different sequence contexts is dictated by the competing amino groups in the minor groove, with possible additional stability derived from hydrogen-bonding interactions of the hydroxyl groups on the BP benzylic ring with keto groups of neighboring cytosines.

Steric competition in minor groove produces sequence-governed conformational differences
The present study provides detailed insights into the structures of two BP-DNA adducts that have guanine bases flanking the modified G*, on either the 3 0 -side or the 5 0 -side. The exocyclic amino groups on the modified and complementary strands in the immediate vicinity of the adducts protrude into the minor groove and strongly affect the structural properties of the 10S (þ)-trans-anti- [BP]-N 2 -dG adduct conformations in the minor groove. Like the structures in the CG*C sequence context studied previously by Cosman et al. (17) and in a TG*C sequence investigated by Fountain and Krugh (19), the G6*G7 and G6G7* duplexes exhibit 5 0 -directed BP minor groove conformations. However, the orientations of the BP rings are subtly different ( Figure 5). In both the G6*G7 and G6G7* duplexes, the exocyclic amino groups on the 3 0 -side contact the aliphatic ring of BP, while the amino group on the 5 0 -side contacts the aromatic ring system. Yet, different steric hindrance effects by the unmodified guanine amino groups combine to produce distinct positioning of the BP moiety in each case. This is evident from the intermolecular NOE patterns ( Table 1). The large upfield chemical shifts due to stacking of the BP rings over the sugar protons of the two bases G20 and C21 in the G6*G7 duplex ( Figure 3B), and predominantly over one complementary base C19 only, in the G6G7* duplex ( Figure 4B), is evidence for the different orientations of the BP rings in these duplexes. These differences are illustrated schematically in Figure S4, Supplementary Data. The chemical shifts of the aliphatic and aromatic BP ring protons are not sensitive to their interactions with the different sugar protons in the G6*G7 and G6G7* duplexes, which is also consistent with minor groove conformations (17,19) ( Figure S5, Supplementary Data). Intercalative structures are characterized by significant upfield shifts of the BP aromatic ring protons due to base stacking interactions (36). We compared the conformations of the BP residues in the G6*G7 and G6G7* duplexes with the structure of the same 10S (þ)-trans-anti- [BP]-N 2 -dG adduct in a CG*C double-stranded sequence context, previously studied by NMR (17), using an MD-derived ensemble which reproduced well the NOE-derived distance restraints (26) and employed the same force field as in the current work. There are subtle differences in the orientation of the BP rings in the minor groove in all three sequence contexts G6*G7, G6G7* and CG*C, with attendant small differences in w, a 0 and b 0 torsion angles (defined in Figure 1A) ( Figure 7A and Table 3). The combination of subtle differences in these three torsion angles together reflects the subtly different positioning of the pyrenyl ring system in the three sequences ( Figure 6), because the steric hindrance effects differ for each case. G6G7* is most distinct because it uniquely contains a 5 0 -flanking amino group on the same strand. Specifically, in the G6G7* duplex, the amino group of G6 on the 5 0 -side of G7* sterically forces the BP rings towards the complementary strand ( Figure 6B). In addition, the bulky amino group of G17, the partner base of C8, is in van der Waals contact with the benzylic ring of BP and thus also affects the positioning of the BP aromatic ring system in the minor groove. In the case of the G6*G7 duplex, the amino group of G7 which is 3 0 to G6*, as well as the amino group of base G20 which is the Watson-Crick partner of C5, must be avoided by the bulky benzylic ring ( Figure 6A). In the case of the CG*C duplex, the pyrenyl and benzylic rings are again positioned to avoid collisions with the amino groups of G16 and G18 in the complementary strand, which are the Watson-Crick partners of C5 and C7, respectively ( Figure 6C). In all cases, the BP rings are 5 0 -directed along the modified strand, the torsion angles a 0 b 0 remaining in the low energy domains (a 0 ¼ 180 AE 408, b 0 ¼ 270 AE 408) computed for this 10S (þ)-trans-anti- [BP]-N 2 -dG adduct (37) and observed in all NMR solution structures in which the cytosine partner to the lesion is present (18), and also in a crystal structure within a replicative polymerase (38). Therefore, modest flexibility in the a 0 , b 0 torsion angles within the allowed domains (37) permits the subtle differential accommodation of the bulky BP rings in the minor groove.

Sequence-governed duplex distortions
It has been shown that oligonucleotide duplexes with 10S (þ)-trans-anti- [BP]-N 2 -dG adducts in the G6G7* sequence context exhibit a bend at the site of the lesion that are characterized by remarkably diminished electrophoretic mobilities relative to the same duplexes but with a G6*G7 sequence context (20). The most striking difference in the G6*G7 and G6G7* duplexes is the remarkable untwisting in the latter case ( Figure 7B). There is a striking decrease in the Twist angle of the G6 Á C19 base pair to 138 relative to that of the adjacent modified G7* Á C18 base pair whose Twist angle is somewhat above normal (408). The Twist angles are about 308 in the same region in the G6*G7 duplex ( Figure 7B). This large untwisting effect in the G6G7* duplex is consistent with a bend at the lesion site, which markedly decreases the electrophoretic mobilities of such duplexes with 10S (þ)-trans-anti- [BP]-N 2 -dG adducts in 5 0 -ÁÁÁ GG*ÁÁÁ sequence contexts (20). Untwisting has been shown to correlate with increased Roll (seen also in our G6G7* case, Figure S2, Supplementary Data) and hence DNA bending (39)(40)(41). In contrast, the electrophoretic mobilities are normal in the 5 0 -ÁÁÁ GG*ÁÁÁ sequence context (20). The bend in the G6G7* duplex is attributed to the steric hindrance by the exocyclic amino group in the minor groove at G6, on the same strand, which causes the untwisting needed to accommodate BP, with the attendant bend. In support of this interpretation, The definitions of the minor groove widths for the G6G7 12-mer duplexes is illustrated above the figure. The minor groove widths are the distance between phosphate groups P5 and P24, P6 and P23, P7 and P22, etc. The values shown were obtained after subtracting 5.8 Å from the pairwise distances to account for the van der Waals radius of the P atoms (34). The minor groove widths of the unmodified control were obtained from an unrestrained 3.0 ns simulation of the unmodified sequence (Cai et al., unpublished data), using the last 1.5 ns of production. (C) Computed ensemble average of Twist angles for the G6G7* and the G6G7* duplexes. The numbering scheme for the nucleotide base pair steps is that C1 Á G24 to A2 Á T23 is step 1, A2 Á T23 to T3 Á A22 is step 2, . . . and so on. The minor groove widths and Twist angles were calculated using MD Toolchest (35). Ensemble average values based on the last 1.5 ns of an MD simulation (26), employing the NMR solution structure of Cosman et al. (17) as initial structure and utilizing the same AMBER force field as in the current study. This MD simulation reproduced well the interproton distances given in (17).
it was found that the replacement of the 5 0 -flanking G in GG* duplexes by inosine, lacking the amino group, abolishes the unusual electrophoretic mobilities (Ruan,Q. and Geacintov,N.E., unpublished data). However, in the CG*C sequence context, no significant effect on electrophoretic mobility (20) and untwisting (26) is observed. Another difference we observe between the G6*G7 and G6G7* duplexes is the marked local helix destabilization in only the former case, which exhibits low-temperature (20-308C) perturbation of Watson-Crick base pairing. While these differences are important in the base pair opening local dynamics around the lesion sites, the effect on the overall cooperative melting of the two duplexes is small ( Figure S1D, Supplementary Data). Another structural difference we find is widening of the minor groove centered about the lesion site, which is greater in the G6G7* case ( Figure 7C). We find that an amino group of a guanine base flanking the 10S (þ)-transanti- [BP]-N 2 -dG adduct on the 5 0 -side causes a greater distortion of the minor groove and prevents the bulky BP rings from aligning as well within the normal groove contours. However, the strain caused by the greater groove opening does not strongly affect the thermal stabilities at the modified base pair nor the flanking base pairs in G6G7* duplexes. The space filling views in Figure S6, Supplementary Data, show the overall structures of the G6*G7 and G6G7* duplexes and the differences in groove widths. In sum, the steric effects of the amino groups differ depending on their exact position, and produce different structural distortions. In the G6G7* case, the steric hindrance produces prominent untwisting with attendant bending. In the G6*G7 case, the necessary accommodation of the BP rings to the amino groups causes instead local helix destabilization. Local minor groove width opening centered about the lesion site occurs in both cases, but is somewhat greater for G6G7*. In all cases, these amino groups act as guideposts for the placement of the BP in the minor groove.

Biological implications
The flanking guanine bases strongly affect the structural characteristics of the 10S (þ)-trans-anti- [BP]-N 2 -dG lesions in double-stranded DNA. The G6G7* duplex is characterized by a wider minor groove centered about the lesion site than the G6*G7 duplex, but the base pairing at the lesion site and flanking base pairs is more stable in the former than the latter. Local untwisting with concomitant bending occurs only in the G6G7* duplex. Differences in distortions caused by the exocyclic amino groups of guanines surrounding these lesions on either side may be distinguished by DNA repair enzymes that are known to probe the structural integrity of the damaged DNA (42). Indeed, recent experiments in our laboratory have shown that human nucleotide excision repair enzymes in cell-free extracts exhibit distinct differences in excising sequences with the 10S (þ)-trans-anti- [BP]-N 2 -dG adduct positioned at either one or the other G in GG sequence contexts. The dual incisions are more efficient in the G6*G7 than the G6G7* sequence context by a factor of $2 (manuscript in preparation). While this difference is not large, the lesions that are more slowly repaired in cellular environments and in tissues are more likely to survive until DNA replication occurs. Finally, sequence effects can also manifest themselves in DNA replication. The relevance of the minor groove orientation of the 10S (þ)-trans-anti- [BP]-N 2 -dG lesion to polymerase structure and function has recently been shown by crystallographic analyses of a replicative DNA polymerase containing the 10S (þ)trans-anti- [BP]-N 2 -dG adduct (38).