Mutagenic nucleotide incorporation and hindered translocation by a food carcinogen C8-dG adduct in Sulfolobus solfataricus P2 DNA polymerase IV (Dpo4): modeling and dynamics studies.

Bulky carcinogen-DNA adducts commonly cause replicative polymerases to stall, leading to a switch to bypass polymerases. We have investigated nucleotide incorporation opposite the major adduct of 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine (PhIP) in the DinB family polymerase, Dpo4, using molecular modeling and molecular dynamics (MD) simulations. PhIP, the most prevalent heterocyclic aromatic amine formed by cooking of proteinaceous food, is mutagenic in mammalian cells and is implicated in mammary and colon tumors. Our results show that the dG-C8-PhIP adduct can be accommodated in the spacious major groove Dpo4 open pocket, with Dpo4 capable of incorporating dCTP, dTTP or dATP opposite the adduct reasonably well. However, the PhIP ring system on the minor groove side would seriously disturb the active site, regardless of the presence and identity of dNTP. Furthermore, the simulations indicate that dATP and dTTP are better incorporated in the damaged system than in their respective mismatched but unmodified controls, suggesting that the PhIP adduct enhances incorporation of these mismatches. Finally, bulky C8-dG adducts, situated in the major groove, are likely to impede translocation in this polymerase (Rechkoblit et al. (2006), PLoS Biol., 4, e11). However, N2-dG adducts, which can reside on the minor groove side, appear to cause less hindrance when in this position.

INTRODUCTION

DNA polymerases of the Y-family have in recent years been shown to play a predominant role in synthesis past DNA bulky lesions, such as those derived from polycyclic aromatic chemicals present in tobacco smoke, automobile exhaust, and broiled meat and fish (1–5). High fidelity replicative DNA polymerases are usually impeded by such damage, leading to a switch to one or more low fidelity bypass polymerases for translesion synthesis (6–11).

The model Y-family polymerase, DNA polymerase IV (Dpo4), is from the archaeon bacterium Sulfolobus solfataricus P2. It is a member of the DinB family (12), of which human pol κ is also a member (13). It has been extensively investigated by X-ray crystallography in binary complexes with primer/template DNA and in ternary complexes in the presence of dNTP, both with and without DNA damage (14–23). These structures have revealed a spacious, water accessible active site that is capable of accommodating two templating bases, in contrast to high fidelity replicative polymerases whose ternary complexes show tight fit of the nascent base pair with exclusion of solvent (24). Dpo4 has an open pocket on the major groove side of the template, as well as a smaller open space on the minor groove side; this is in contrast to replicative polymerases, which only have an open pocket on the major groove side, while the minor groove side is packed with protein–DNA interactions critical for polymerase fidelity (25–29). In addition like other Y-family polymerases, Dpo4 has a unique little finger domain, also called wrist or polymerase associated domain (PAD) at the C-terminus (15,30,31). The flexibility of this little finger domain is believed to play an important role in accommodating specific types of DNA lesions (14,15,21). Crystal structures of Dpo4 binary and ternary complexes also reveal that the little finger domain plays a key role in translocation (23).

Recently, it has been suggested that the DinB family polymerases may be specifically suited for bypass of N2-dG adducts (32). The presence of a conserved ‘steric gate’, usually phenylalanine or tyrosine, is shown to be crucial in bypass of N2-dG minor groove adducts (32). However, the DinB polymerases may be less well suited for bypassing C8-dG bulky adducts. Such adducts would normally be expected, from their position of substitution on guanine, to reside in the major groove of double-stranded DNA. However, rotation of the modified guanine from anti to syn would place the adduct in the minor groove area, in a position roughly similar to an N2-substituted guanine. Despite this formal possibility, experimental primer extension data for two C8-dG adducts derived from 2-acetylaminofluorene (AAF), namely N-(deoxyguanosin-8-yl)-2-acetylaminofluorene (dG-C8-AAF) with Dpo4 (12) and N-(deoxyguanosin-8-yl)-2-aminofluorene (dG-C8-AF) with Escherichia coli DinB and human pol κΔc (33), indicates that these adducts cause polymerase stalling or blockage, with only small amounts of primer extension beyond the lesion site. A recent molecular dynamics (MD) study from our group has provided structural rationale for the case of dG-C8-AAF in Dpo4 (34). However, this adduct contains an acetyl group, which adds to the steric hindrance in C8 adducts.

Here, we investigate a C8-dG adduct (Figure 1a) derived from the most prevalent heterocyclic aromatic amine formed by cooking proteinaceous food, 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine (PhIP) (35–43). This substance causes predominantly G to T transversions in mammalian cells, with some G to A transitions, and a few deletions (44–52). In a previous molecular modeling and MD study of this adduct in the Pol α family replicative polymerase RB69, we found considerable active site distortion caused by this lesion at the templating base or in the first double-stranded extension position, when partnered by C or A (53). These results suggested that the replicative polymerase would be stalled by the lesion, affording opportunity for switch to one or more bypass polymerases. We now use this dG-C8-PhIP adduct, N2-(2′-deoxyguanosin-8-yl)-PhIP, as a model to study the structural feasibility of accommodating a C8-dG lesion lacking the acetyl group in Dpo4, employing Dpo4 type I (16) and type II (15) crystal structures as initial models (Figure 1).

Figure 1

(a) Structure and numbering scheme of the dG-C8-PhIP adduct. Torsion angle definitions are as follows: χ, O4′-C1′-N9-C4; α′, N9-C8-[PhIP]N2-[PhIP]C2; β′, C8-[PhIP]N2-[PhIP]C2-[PhIP]N1; γ′, [PhIP]C5-[PhIP]C6-[PhIP]C1′-[PhIP]C2′. (b) dG-C8-PhIP-modified DNA sequence used in G·dNTP. Models are derived from a Type I crystal structure of Dpo4 (PDBID: 1S0M) (16). (c) dG-C8-PhIP-modified DNA sequence used in −1 deletion studies. Models are derived from a Type II crystal structure of Dpo4 (PDBID: 1JXL) (15).

Our results show that this dG-C8 adduct, like dG-C8-AAF, if accommodated on the Dpo4 minor groove side pocket, would cause serious distortions to the active site region. However, the PhIP ring system can be accommodated in the spacious major groove pocket of this polymerase, and Dpo4 appears capable of incorporating dCTP, dTTP or dATP reasonably well. In addition, our results show that the PhIP-modified lesion produces less distorted structures for anti-G*·anti-dTTP and anti-G*·syn-dATP than for their respective unmodified controls (anti-G·anti-dTTP and anti-G·syn-dATP). These findings may be relevant to observed mutagenic behavior of PhIP in inducing G to T transversions and G to A transitions in mammalian systems (44–52), if the human DinB polymerase pol κ is involved in the error-prone incorporation opposite the lesion, and its structural properties prove to be similar to those of its prokaryotic homolog Dpo4. Finally, based on Dpo4 binary and ternary complex structures (23), we suggest that translocation may be seriously inhibited by bulky dG-C8 adducts positioned on the major groove side of the DNA duplex region, while N2-dG adducts residing on the minor groove side (54) would be less hindering.

COMPUTATIONAL METHODS

Molecular modeling

G·dNTP initial models. The crystal structure of Dpo4 polymerase with 10R(+)-cis-anti-benzo[a]pyrene(BP)-N6-dA located on the major groove side of the modified DNA (BP-2 complex) (16) was used to obtain initial structures for the G·dNTP models (PDB ID: 1S0M) (55). Of the many type I structures of Dpo4, the BP-2 complex has the most reaction-ready active site, since it contains two metal ions in the active site and the 3′ hydroxyl group at the primer end. This structure was further remodeled, as described in previous work from our group (34), to achieve ideal Mg2+ coordination (Supplementary Table S1) and an O3′-Pα distance of 3.1 Å, in the reaction-ready range.

The DNA sequence was then adjusted to match the sequence in the adenomatous polyposis coli (Apc) gene mutational hotspot for PhIP, codon 635 (Figure 1b), with the guanine selected for PhIP modification situated at the templating position opposite the incoming dNTP. The initial models for dynamics were then constructed by locating structures with minimal crowding between the PhIP moiety and the polymerase, by rotating α′ and β′ at 10° intervals, in combination (Figure 1a). γ′ was initiated at 26° as in the NMR solution structure (56). Both anti and syn conformations of the glycosidic torsion χ were investigated, with χ adjusted to achieve optimal hydrogen bonding and stacking in the nascent base pair for each χ domain. The torsion angles in the initial models are summarized in Supplementary Table S2.

For the G*·dATP mismatch, three starting models were obtained. Two of them featured the PhIP-modified guanine (G*) anti and the PhIP rings on the major groove side of the template, with the incoming nucleotide dATP either anti or syn. A syn-dATP was investigated since it has been observed opposite a lesion in a Dpo4 crystal structure (17). A high resolution NMR solution structure containing a dG-C8-PhIP in an 11mer duplex shows a base-displaced intercalation model (56). In this conformation, the modified guanine is syn and displaced into the major groove, while the PhIP rings intercalate within the DNA duplex. However, the PhIP rings in such a base-displaced intercalation conformation would occupy the position of the incoming nucleotide in Dpo4. A PhIP-modified syn-guanine model could be built with the PhIP rings on the minor groove side of the dNTP, which does allow the accommodation of the dNTP at the active site. With a syn-G* as the template, only anti-dATP is favorable, since the bases in the syn-G*·syn-dATP pair are too far apart for hydrogen bonding. Hydrogen bonding in the nascent base pair is an important consideration in building the initial models, since Dpo4 as well as other DinB family polymerases, such as Dbh and pol κ rely more on hydrogen bonding for catalytic efficiency than the high fidelity polymerases (57,58).

In the G*·dCTP/dTTP/dGTP series, the G*s were all constructed in the same anti conformation as in the G*·dATP mismatch. A syn-G* was eliminated based on results for the syn-G*·anti-dATP simulation; these show that a syn-G* generated much more disturbance to the polymerase active site than an anti one, as described in Results, and these disturbances are essentially independent of the dNTP. In the G*·dCTP model, the incoming dCTP is anti, in order to form Watson–Crick hydrogen bonds with the anti-G*. In the G*·dTTP model, the dTTP is anti to achieve a wobble pair with the anti-G*. A wobble paired T·dGTP has been observed in a Dpo4 ternary complex (18). A syn-dTTP can form no hydrogen bonds with anti-G* and syn pyrimidines have scarcely been observed (59). In the G*·dGTP model, a syn-dGTP was employed, since an anti-dGTP could not form hydrogen bonds with anti-G* in the polymerase, and would collide with either the anti-G* or the minor groove side protein residues.

As an unmodified control for all the simulations, an unmodified anti-G opposite a Watson–Crick paired anti-dCTP was used. Additional controls containing an unmodified anti-G opposite anti-dATP for the G*·dATP mismatch, anti-dTTP for the G*·dTTP mismatch and syn-dGTP for the G*·dGTP mismatch were also obtained. Stereoviews of the active sites of the initial models are shown in Supplementary Figure S1.

−1 deletion initial models. The crystal structure of the Dpo4 type II ternary complex (PDB ID: 1JXL) (15) was used to obtain the starting models for the −1 deletion structures (Figure 1). In this crystal structure, the coordinates of the first base at the 5′ end of the template could not be resolved, probably due to its flexible conformation outside the polymerase. We therefore modeled it into the structure, using the conformation of an analogous terminal base in the Dpo4 type I structure (PDB ID: 1JX4). Hydrogen atoms absent in the crystal structures were added by the AMBER suite. The dideoxy group at the 3′ end of the primer was replaced by a hydroxyl group. The Ca2+ ion, residing in the position of the nucleotide binding metal ion, was replaced by a Mg2+ ion, and repositioned for proper octahedral coordination with water molecules and amino acid residues (Supplementary Table S1). The DNA sequence was also remodeled to match the sequence in the Apc gene mutational hotspot codon 635 (Figure 1c). The PhIP moiety was linked to the unpartnered guanine in the active site, while the incoming dCTP paired with the guanine on the 5′ side of the adduct. Starting models were obtained using the same approach as for G*·dNTP models, by rotating torsion angles α′ and β′ at 10° intervals, in combination. Structures were selected for anti and syn domains of χ based on minimal steric close contacts and optimal stacking. The torsion angles in these starting models are summarized in Supplementary Table S2 and stereoviews of the active sites of the initial models are shown in Supplementary Figure S1.

Force field parameterization

Parameters for the anti-dG-C8-PhIP adduct and the incoming anti-dCTP, anti-dATP, syn-dATP and anti-dTTP were the same as in earlier work (53,54). The parameters for the syn-dG-C8-PhIP adduct and the syn-dGTP were obtained using the same method described previously (53). Supplementary Tables S3 and S4 shows the AMBER atom type, connection type and partial charge assignment for these cases.

MD simulation and data analyses

Simulations were carried out using the SANDER module of the AMBER 6.0 MD software package (University of California, San Francisco), the Cornell et al. force field (60) and the parm99 parameter set (61). Electrostatic interactions were approximated by the particle mesh Ewald method (62), and a 10 Å cutoff was applied to Lennard–Jones interactions. All bonds involving hydrogen atoms were constrained by the SHAKE algorithm (63) with a tolerance of 0.0005 Å, and a 2 fs time-step was used in the dynamics simulation. Periodic boundary conditions were applied, and all MD simulations were carried out under constant temperature, 300 K, with a temperature coupling parameter of 4.0 ps and at constant atmospheric pressure with a 1.0 ps coupling parameter. The translational motion of the center of mass was removed every 1 ps. No obvious overall rotation of the system was observed during the simulation; thus, energy leakage from internal motion to global rotation through the ‘flying ice cube effect’ did not happen here (64).

Fifteen Na+ ions were added to each system for neutralization using the LEaP module of AMBER, followed by 600 steps of steepest descent (SD) and 400 steps of conjugate gradient (CG) to relax the added Na+ ions and crystal waters. Then ∼15 000 TIP3P water molecules were added to solvate each system creating a rectangular periodic box containing a total of ∼50 000 atoms. The systems were further minimized and heated up to 300 K. The equilibration protocol was conducted in the same manner as in earlier work (53) and the details are given in Supplementary Data. Production MD simulation was carried out for 3 ns.

The PTRAJ, ANAL and CARNAL modules of AMBER were employed to analyze the trajectories. The root-mean-square deviations (RMSDs) of each system were calculated relative to the starting structures as shown in Supplementary Figure S2. The overall structure and the active site region of the protein–DNA complexes became reasonably stable after ∼500 ps production MD simulation in each system. Therefore, the last 2.5 ns of MD simulation were used for structural analyses. Stereoviews of unmodified and modified systems after the total 3 ns of production MD simulations are shown in Supplementary Figure S3, and the active sites of these systems are shown in Figures 2 and 3, and Supplementary Figure S4.

Figure 2

The active sites after 3 ns MD in the G·dNTP models. The Dpo4 DNA polymerase, shown in ribbon representation, is colored as follows: thumb is in dark orange, palm in lime, finger in yellow, and little finger in blue; the DNA ternary complex is shown in stick representation. The primer and template DNA strands are colored in gray, the dCTP is colored in cyan, dATP in magenta, dTTP in marine, dGTP in pink, the template guanine is colored in green and the PhIP moiety is in red. Mg2+ ions A and B, shown in spheres representation, are colored in orange.

Figure 3

The active sites after 3 ns MD in −1 deletion models. The protein and two Mg2+ ions are colored the same as in Figure 2. The template guanine is colored in purple; the dCTP is in cyan; the template guanine without partner is in green, and the PhIP moiety is in red.

Molecular modeling was carried out with Insight II 97.0 (Accelrys, Inc., a subsidiary of Pharmacopeia, Inc.). Figures of structures were prepared with PyMOL (DeLano Scientific, LLC.).

RESULTS

Initial models

Our initial models (see Computational Methods) investigated both anti and syn conformations of the adduct in the active site of Dpo4, partnered with dATP or unpaired in the −1 deletion simulations. In addition, anti-G* opposite anti-dCTP, anti-dTTP or syn-dGTP were investigated using results from dATP to guide selection of the starting structures. Production MD simulations were carried out for 3 ns and the ensembles of structures derived from the last 2.5 ns trajectory were analyzed.

Structural analyses

A number of structural criteria were evaluated in each modified system and compared to the undamaged control simulation. These include: (i) the number and occupancy of the hydrogen bonds in the nascent base pair (Supplementary Figure S5 and Supplementary Tables S5 and S10); (ii) the distance between C1′ of the template and C1′ of the incoming nucleotide (dNTP), normally near 10.8 Å in a Watson–Crick pair (65) (Supplementary Table S7); (iii) stacking interactions between the nascent base pair and the primer-terminus base pair (Figures 2 and 3, and Supplementary Tables S8 and S9); (iv) the frequency of sampling a near reaction-ready distance (3.1 to 3.5 Å) between O3′ of the primer-terminus and Pα of the dNTP (Supplementary Table S7); (v) angle O3′ (primer 3′ end)-Pα(dNTP)-O3α(dNTP), ideally 180° for in-line attack of O3′ on Pα (66) (Supplementary Table S7); (vi) the chelation of Mg2+ ions (Figures 2 and 3); (vii) the distance between the two Mg2+ ions (Supplementary Table S7) and (viii) hydrogen bonds between the nascent base pair and the polymerase. These are summarized in Tables 1 and 2, using a distortion scoring function to evaluate the lesion impact on the polymerase active site (34). While we analyze the structural features individually, we realize there is interdependence among them; how they affect each other remains unclear but their composite quality sheds light on a given model's structural feasibility. While scoring has subjective criteria elements, it reflects our current stage for evaluating polymerase active site distortions. We hypothesize that the more distorted the active site is, the less efficient the nucleotidyl transfer reaction will be.

Table 1

Composite evaluation of structural distortion in G·dNTP models

	Control	anti-G*·dCTP	anti-G·syn-dATP	anti-G*·anti-dATP	anti-G*·syn-dATP	syn-G*·anti-dATP	anti-G·anti-dTTP	anti-G*·anti-dTTP	anti-G·syn-dGTP	anti-G*·syn-dGTP
Pα-O3′a	√	√	×	√	√	×	×	√	√	×
HB geometryb	√	√	××	××	×	×××	××	×	××	××
C1′-C1′c	√	√	√	×	√	×	√	√	×	×
Key amino acids-DNA interactiond	√	××	×××	×××	×	×××	×	×	××	××
Base pair stackinge	√	√	√	√	√	×	√	√	√	√
In-line attackf	√	√	√	√	√	√	√	√	√	√
Mg2+ ions distanceg	√	√	√	√	√	√	√	√	√	√
Mg2+ ions octahedral coordinationh	√	√	√	√	√	√	√	√	√	√
Composite scorei	0	−2	−6	−6	−2	−9	−4	−2	−5	−6

aOne cross denotes a decrease in occupancy compared to the control by more than 15%, two crosses denote a decrease by >30%, and three denote a decrease by >45%. See Supplementary Table S7 for details.

bEach cross denotes a decrease in occupancy compared to the control by >10%, or a loss of one hydrogen bond. See Supplementary Table S6 for details.

cOne cross denotes an increase in C1′–C1′ distance compared to the control (10.9 Å) by more than 1 Å. See Supplementary Table S7 for details.

dEach cross denotes a decrease in occupancy compared to the control by >10%, or a loss of one stacking interaction. See Supplementary Table S6 for details.

eOne cross denotes a disruption of stacking interaction between the nascent and the primer-terminus base pairs, indicated by diminished coplanarity of the base pairs determined visually, and a weakened van der Waals interaction energy between them compared to the control (−15.1 kcal/mol) by ≥5.0 kcal/mol. See Figure 2 and Supplementary Table S8 for details.

fOne cross denotes occupancy ≤50%. See Supplementary Table S7 for details.

gOne cross denotes occupancy ≤90%. See Supplementary Table S7 for details.

hNo disruptions of the chelation interactions involving the two Mg2 ions are observed (the distance between chelating atom and Mg2 ≥2.5 Å is considered to be a disruption). See Figure 2 for details.

iScore is based on the unmodified control, which was assigned a base score of 0. Each cross penalizes the base score by 1.

Table 2

Composite evaluation of structural distortion in the −1 deletion models

	Control	anti-G* −1 deletion	syn-G* −1 deletion
Pα-O3′a	√	××	×××
HB geometryb	√	√	√
C1′-C1′c	√	√	√
Key amino acids-DNA interactiond	√	×××	××××
Base stackinge	√	√	×
In-line attackf	√	√	×
Mg2 ions distanceg	√	√	√
Mg2 ions octahedral coordinationh	√	√	√
Composite scorei	0	−5	−9

aOne cross denotes a decrease in occupancy compared to the control by >15%, two crosses denote a decrease by >30%, and three denote a decrease by >45%. See Supplementary Table S7 for details.

bEach cross denotes a decrease in occupancy compared to the control by >10%, or a loss of one hydrogen bond. See Supplementary Table S11 for details.

cOne cross denotes an increase in C1′–C1′ distance compared to the control (10.8 Å) by >1 Å. See Supplementary Table S7 for details.

dEach cross denotes a decrease in occupancy compared to the control by >10%, or a loss of one stacking interaction. See Supplementary Table S11 for details.

eOne cross denotes a disruption of the stacking interactions between G4 and A5, indicated by diminished coplanarity of the bases determined visually and a weakened van der Waals interaction energy between them compared to the control (−5.3 kcal/mol) by ≥2.5 kcal/mol. See Figure 3 and Supplementary Table S9 for details.

fOne cross denotes occupancy ≤50%. See Supplementary Table S7 for details.

gOne cross denotes occupancy ≤90%. See Supplementary Table S7 for details.

hNo disruptions of the chelation interactions involving the two Mg2 ions are observed (the distance between chelating atom and Mg2 ≥2.5 Å is considered to be a disruption). See Figure 3 for details.

iScore is based on the unmodified control, which was assigned a base score of 0. Each cross penalizes the base score by 1.

G·dNTP models

Unmodified control of anti-G·anti-dCTP. In the unmodified control simulation, anti-G Watson–Crick pairs with an incoming anti-dCTP. Views of this complex and its active site after 3 ns production MD simulation are shown in Supplementary Figures S3a and S4a, and Figure 2a. The active site remains essentially undisturbed throughout the simulation. The Watson–Crick hydrogen bonds in the nascent base pair have occupancies >90% (Supplementary Table S5); the C1′–C1′ distance in the nascent base pair remains normal (10.9 ± 0.1 Å); the chelation of the two Mg2+ ions is preserved (Figure 2a); the distance between the two Mg2+ has an average value of 3.7 ± 0.1 Å; the Pα–O3′ distance samples the near reaction-ready range frequently, during 78.7% of the time; and the in-line attack angle is 167.1 ± 5.1°, close to the ideal value of 180°. These active site features are similar in polymerase ternary complex crystal structures (67,68). In the active site of Dpo4 polymerase, the flat face of the nascent base pair is topped by protein residues with unusually small and hydrophobic side chains: Val32, Ala42, Ala44, Ala57 and Gly58 (Figure 2a). These residues form the ‘ceiling’ of the dNTP binding pocket (25). The ‘floor’ of this pocket is considered to be the flat face of the base pair at the primer terminus (25). The geometry of the dNTP binding pocket remains normal here (Figure 2a). A highly conserved aromatic residue in the DinB family (Tyr12 in Dpo4, Phe12 in Dbh and Phe13 in DinB) is tightly packed on the sugar ring of the dNTP. This has the effect of placing the dNTP properly for catalysis and acting as a ‘steric gate’ to prevent binding of ribo-NTPs through steric exclusion of a 2′-OH (32,69). In our unmodified control simulation, the proper position of Tyr12 is preserved. Polymerase residues Tyr10, Phe11, Tyr12, Ser34, Thr45, Arg51, Lys159, Thr250 and Arg331 form hydrogen bonds with the nascent base pair in the unmodified control simulation (Supplementary Table S6).

Syn-G*·anti-dATP minor groove. The syn-G*·anti-dATP minor groove model is highly disturbed. With G* syn, the PhIP rings are on the minor groove side and inserted between the nascent base pair and the base pair at the primer terminus (Figure 2f). The location of the PhIP rings in the syn-G*·anti-dATP simulation disrupts the flat face of the primer-terminus base pair, and also causes poor stacking interactions between the nascent and the primer-terminus base pairs, as shown by energetic consideration: a partial energetic assessment of differential stacking interactions can be obtained from van der Waals interactions between bases for the same sequences. In the syn-G*·anti-dATP simulation, this interaction between the nascent and the primer-end base pairs is only −2.3 ± 0.5 kcal/mol, significantly weaker than −16.8 ± 0.9 kcal/mol in the anti-G*·syn-dATP, or −15.8 ± 0.9 kcal/mol in the anti-G*·anti-dATP case. The syn-G*·anti-dATP also manifests an enlarged C1′–C1′ distance (12.2 ± 0.4 Å), and a less frequently sampled near reaction-ready Pα-O3′ distance (59.8%). The hydrogen bond in the nascent base pair in the syn-G*·anti-dATP initial model, between O6 of the syn-G* and N6 of the anti-dATP, is ruptured during equilibration and remains broken throughout the simulation. To achieve the hydrogen bond in the initial model, the PhIP moiety was positioned in a crowded region near the ceiling of the dNTP binding pocket, which proved unfavorable (Supplementary Figure S1f). These disturbances, triggered by the syn-G* adduct, are essentially independent of the specific dNTP.

Anti-G*·anti-dATP and anti-G*·syn-dATP major groove. The anti-G*·anti-dATP and anti-G*·syn-dATP major groove models are more disturbed than the unmodified control, but less than the syn-G*·anti-dATP minor groove model. In the anti-G*·syn-dATP complex, the phenyl ring of the PhIP is directed towards the 3′ end of the template, along the template backbone (Figure 2e), while in the anti-G*·anti-dATP complex, it is oriented toward the 5′ end of the template (Figure 2d). A ∼180° difference in α′ is responsible for this difference in orientation (Supplementary Figure S6). With an incoming syn-dATP, one full and one bifurcated hydrogen bond are formed in the nascent base pair (Supplementary Figure S5 and Supplementary Table S5), with a close to normal C1′–C1′ distance (11.1 ± 0.2 Å). On the other hand, with the incoming nucleotide anti, the C1′–C1′ distance is significantly enlarged (12.7 ± 0.2 Å), and the full hydrogen bond between O6 of G* and N6 of anti-dATP has an occupancy 20% lower than for syn-dATP (Supplementary Table S5). A bifurcated hydrogen bond is formed between the little finger residue Arg332 and the PhIP rings only in the anti-G*·syn-dATP simulation (Supplementary Table S6). The frequencies of sampling the near reaction-ready Pα-O3′ distance and the geometries of the binding pocket are not dramatically different in these two major groove complexes (Supplementary Table S7 and Figure 2).

Anti-G*·anti-dTTP and anti-G*·syn-dGTP major groove. The anti-G*·anti-dTTP major groove model is little distorted, while the anti-G*·syn-dGTP major groove structure is very distorted. At this stage, the syn-G* was eliminated due to poor accommodation of the minor groove positioned PhIP regardless of dNTP, as revealed in the syn-G*·anti-dATP simulation (Figure 2f and Table 1). Therefore, only anti-G* was employed to build the models for G*·dTTP and G*·dGTP. With the incoming dTTP anti, the anti-G* forms a wobble pair with the anti-dTTP (Supplementary Figure S5). During the anti-G*·anti-dTTP simulation, the two hydrogen bonds in the wobble nascent base pair are highly occupied (>94% of the time); the C1′–C1′ distance has a normal average value of 10.8 ± 0.2 Å; the Pα-O3′ distance remains in the near reaction-ready range almost all the time (93.0%), and the nascent base pair is well stacked with the base pair at the primer-terminus. In the G*·G mismatch model, the incoming dGTP is syn and forms two hydrogen bonds with the anti-G* (Supplementary Figure S5 and Supplementary Table S5). However, the C1′–C1′ distance is enlarged to 11.9 ± 0.3 Å, and a somewhat less frequently sampled near reaction-ready Pα-O3′ distance is noted (68.1%).

Anti-G*·anti-dCTP major groove. With anti-G* opposite a Watson–Crick paired dCTP in the active site, the structural features are comparable to those of the unmodified control. In this anti-G*·anti-dCTP simulation, the Watson–Crick hydrogen bonds in the nascent base pair are highly occupied (>95% of the time); the C1′–C1′ distance has a normal average value of 10.9 ± 0.9 Å; the nascent base pair stacks well with the base pair at the primer-terminus, and the near reaction-ready Pα-O3′ distance is frequently sampled (82.4% of the time). However, two protein–DNA interactions have lower hydrogen bond occupancies in comparison to the unmodified control (Supplementary Table S6).

Controls of anti-G·syn-dATP, anti-G·anti-dTTP and anti-G·syn-dGTP. These control models containing G·dATP, G·dTTP and G·dGTP mismatches without PhIP modification show distortions at the active sites in comparison to the Watson–Crick G·dCTP paired unmodified control (Table 1). In the control of anti-G·syn-dATP, although the C1′–C1′ distance has a close to normal value of 11.2 ± 0.2 Å (Supplementary Table S7), the full and the bifurcated hydrogen bonds in the nascent base pair have lower occupancies (<90% of the time, Supplementary Table S5), and the near reaction-ready Pα-O3′ distance is less frequently sampled (61.4% of the time, Supplementary Table S7). In the control of anti-G·anti-dTTP, the C1′–C1′ distance has a normal value of 10.8 ± 0.2 Å (Supplementary Table S7). However, one of the two hydrogen bonds in the nascent pair is less frequently occupied (80.2% of the time, Supplementary Table S5), and the Pα-O3′ distance sampled the near reaction-ready range during only 59.7% of the time (Supplementary Table S5). In the control of anti-G·syn-dGTP, the C1′–C1′ distance is enlarged to 11.9 ± 0.3 Å and one of the two hydrogen bonds in the nascent base pair is poorly occupied (only 46.6% of the time, Supplementary Table S5). However, the near reaction-ready Pα-O3′ distance is frequently sampled (76.7% of the time, Supplementary Table S5).

−1 Deletion models

Anti-G*, syn-G* and unmodified control. In these models, the damaged template guanine has been skipped and the incoming dCTP pairs with the guanine on the 5′ side of the adduct (Figure 1c). In comparison to the unmodified −1 control system, the two adduct systems, the anti-G* −1 deletion and syn-G* −1 deletion are significantly distorted. In the anti-G* −1 deletion model, the PhIP moiety is positioned on the major groove side, near-perpendicular to the template strand (Figure 3b); in the syn-G* −1 deletion model, the PhIP moiety is positioned on the minor groove side (Figure 3c). Views of these ternary complexes and their active sites after 3 ns production MD simulations are shown in Supplementary Figures S3 and S4, and Figure 3. For the syn-G* −1 deletion system, the PhIP rings are reoriented through rotation of α′ and β′ after about 320 ps of MD (Supplementary Figure S7) to avoid crowding between the distal phenyl ring and Lys78. This rearrangement results in a favorable environment around the PhIP ring system. The phenyl ring is in van der Waals contact with hydrophobic residue Ile104; the imidazo ring N3 and N4 atoms have favorable electrostatic interactions with the amino group of Lys78 (actual hydrogen bonds are formed ∼9% of the time). As shown in Supplementary Figure S8, the PhIP moiety is pocketed by the protein residues, the nascent base pair and the primer-terminus base pair.

Watson–Crick hydrogen bonds in the nascent base pair are preserved in all three simulations (Supplementary Table S10). The C1′–C1′ distances are also normal in all cases (Supplementary Table S7), as is the chelation of the two Mg2+ ions (Figure 3). Compared to the unmodified control, the primer end is repositioned towards the 5′ direction and away from the active site in the anti-G* and syn-G*adduct simulations (Supplementary Figure S9); this is more prominent in the syn-G* −1 deletion simulation, where the terminal primer base is shifted by about 3 Å. The backbone of the primer strand is concomitantly shifted, accompanied by small rotation of the thumb domain, which maintains contacts with the relocated backbone phosphate groups of the primer strand. In order to retain the B-DNA geometry, the backbone of the template strand is also adjusted, followed by small rotation of the little finger domain in order to maintain contacts with the same phosphate groups of the template strand. A similar relocation of the primer terminal base was observed in a Dpo4 crystal structure containing the 10R(+)-cis-anti-[BP]-N6-dA adduct, with BP aromatic rings intercalated between the nascent and the primer terminal base pair (16). Here, the primer-terminus was positioned >10 Å away from the incoming dNTP. In both the anti-G* and syn-G* −1 deletion simulations, the dCTP loses its proper stacking contact with the ‘steric gate’ Tyr12 from the finger domain, and its hydrogen bonding with the finger residue Arg51. The separation of the dCTP from the finger domain could be a result of hydrogen bonds formed between the dCTP and the adduct (Supplementary Table S11). The Pα-O3′ distance is also affected by this primer side relocation. In the unmodified control system for the −1 deletion, the near reaction-ready Pα-O3′ distance is sampled during the whole simulation (51.4%, Supplementary Figure S10k). In the anti-G* −1 deletion model, this distance is only achieved in the 380–480 ps time frame, and comes close to 4 Å occasionally after this time range (Supplementary Figure S10l). However, in the syn-G* −1 deletion case, the distance has an average value of 7.5 ± 0.7 Å throughout the simulation, and does not approach closer than 5.5 Å (Supplementary Figure S10m). In the syn-G* −1 deletion system, the stacking between the PhIP-modified syn-G and its 3′-neighbour is also disrupted (Figure 3c). Thus, the −1 deletion simulations suggest various distortions likely to hinder nucleotide incorporation.

DISCUSSION

Our simulations suggest that regardless of the incoming dNTP or whether the damaged base has a partner, the dG-C8-PhIP adduct is not likely to be accommodated in the minor groove side pocket of the Dpo4 DinB family polymerase (Tables 1 and 2). When situated in this position, with the G* in the syn conformation, the bulky aromatic rings disrupt the geometry of the active site, particularly on the primer side (Figure 2f and Figure 3c). With the G* anti, however, the aromatic rings fit well on the major groove side. With dCTP opposite the anti-G*, Watson–Crick pairing is preserved and simulations show only modest distortions compared to the unmodified control (Supplementary Table S5 and Table 1). An incoming dTTP affords a well-formed wobble pair with only small distortions. In the case of the dATP, greater distortions are observed, with syn-dATP providing the more favorable structure, which contains hydrogen bonds employing the Hoogsteen edge of the dATP. Incoming dGTP was most distorted because even the syn conformation for dGTP caused enlargement of the nascent base pair. Furthermore, in the case of structures with no partner opposite the lesion, the −1 deletion models, the PhIP aromatic rings can again reside in the large major groove side open space; however, the active site region is notably distorted, and the approach of the primer-terminus to the α-phosphate of dNTP is inhibited by the PhIP moiety, particularly by its protruding methyl group (Supplementary Figure S9 and Supplementary Table S7).

In contrast, an earlier study of the 10S(+)-trans-anti-[BP]-N2-dG showed that this adduct can be accommodated reasonably well on the minor groove side of the Dpo4, irrespective of incoming dNTP, with only modest enzyme perturbation including opening of the little finger and some small rearrangement of active site region residues (54). In this case, the modified dG adopts the anti conformation. The major groove position with the N2-dG adduct in syn conformation also provides a good accommodation for the BP ring system. Experimental kinetic studies revealed promiscuous nucleotide incorporation for this case. Running-start primer extension experiments indicated that the damage can be bypassed to a significant extent (54).

Figure 4 illustrates the differential accommodation of the C8 and N2-dG adducts in the Dpo4 minor groove, showing why the C8 adduct is not well positioned there, while the N2 adduct is, regardless of dNTP (see Supplementary Figure S11 for stereoview). Specifically, this N2 adduct ring system is directed 5′ along the template strand in our models; however, the C8 adduct rings are oriented 3′ along the template strand, protrude to the primer side and thereby disrupt the active site directly where the nucleotidyl transfer reaction is to take place. Perhaps C8-dG bulky adducts, such as dG-C8-PhIP and dG-C8-AAF might less readily allow nucleotide incorporation by Dpo4, in part because of their poor accommodation in the enzyme minor groove side pocket; thus, only one site is available to harbor the lesion, while N2 adducts, such as 10S(+)-trans-anti-[BP]-N2-dG can reside in either major or minor groove side pockets.

Figure 4

Minor groove positioned adducts in Dpo4. The view is into the minor groove side with the major groove in back. (a) syn-dG-C8-PhIP·anti-dATP after 3 ns MD simulation and (b) anti-dG-N-BP·syn-dATP as in Figure 3b of (54). The carcinogen-modified guanine is green; the PhIP moiety is red; the BP moiety is brown; the incoming dATP is magenta, and the primer-terminus base pair is yellow. The protein is colored the same as in Figure 2.

Additionally, our structures suggest that translocation may be difficult for such major groove positioned multi-ringed adducts. Specifically, a recent crystal structure of Dpo4 has suggested that translocation starts with movement of the little finger during nucleotide binding, followed by thumb movement during the chemical reaction (23). The replication cycle involves a screw-like counterclockwise rotation/translocation of the polymerase along the DNA helix axis (viewed in the 5′ to 3′ direction of the template strand). This threading of the DNA through the polymerase places the next templating base in the active site, awaiting entry of its complementary dNTP. Prior study of the dG-C8-AAF adduct suggested that the preferred 3′-directed orientation of the fluorenyl ring system along the template strand, in the major groove, would impede the rotation of the little finger (34); this impediment is not due to the acetyl group. The current study of the dG-C8-PhIP adduct suggests a similar impact on translocation. As shown in Figure 5a and b, in the insertion position the adduct in the anti conformation is enveloped by the little finger without collision (see Supplementary Figure S12 for stereoview). However, once in the post-insertion site, the adduct would be in collision with the little finger, indicating that translocation to the post-insertion site would be difficult. Very large rearrangement of the PhIP adduct would be required to move it to another orientation that allows translocation. Additional bypass polymerases may be required to further extend beyond the lesion (4). Furthermore, it appears that N2 adducts, such as 10S(+)-trans-anti-[BP]-N2-dG, would be less likely to impede translocation by colliding with the little finger, when in the anti conformation and placed in the minor groove side polymerase pocket (Figure 5c and d). The number of rings in the adduct together with their orientation will determine the extent of these effects on translocation and this remains to be elucidated.

Figure 5

Major groove positioned anti-dG-C8-PhIP adducts and minor groove positioned anti-dG-N2-BP adducts in Dpo4. For (a) and (b), the view is into the major groove with the minor groove in back, while for (c) and (d), the view is into the minor groove with the major groove in back. (a) anti-dG-C8-PhIP at the insertion position. (b) Modeled anti-dG-C8-PhIP at the post-insertion position with linkage site torsion angles α′, β′ and γ′ same as in (a); collisions between the PhIP rings and the little finger domain are indicated. (c) anti-dG-N2-BP adduct at insertion position opposite anti-dCTP as in Figure 3a of (54). (d) Modeled anti-dG-N2-BP at the post-insertion position in same linkage conformation as (c). The color code is the same as in Figure 4.

Primer extension studies of dG-N2-AAF and dG-C8-AAF with hpol κ have shown that the frequency of bypass of dG-N2-AAF is at least 4 orders of magnitude higher than for dG-C8-AAF, and it was proposed that the observed difference may be due to the different linkage site (70). This is consistent with our suggested poorer accommodation on the minor groove side of C8 than N2 adducts. However, we do not know exactly how structurally similar pol κ and Dpo4 ternary complexes are, since to date, only a crystal structure of the apo-pol κ catalytic core is available (71).

Furthermore, our results show that the PhIP-modified lesion produces less distorted structures for G*·dTTP and G*·dATP (anti-G*·anti-dTTP and anti-G*·syn-dATP) than for their respective unmodified controls (anti-G·anti-dTTP and anti-G·syn-dATP) (Table 1). The structural origin of this interesting phenomenon is at least partly in the formation of hydrogen bonds between the PhIP rings and the little finger residue Arg332, which stabilize the anti conformation in the adduct (Figure 2 and Supplementary Table S6). These findings would be relevant to observed mutagenic behavior of PhIP in inducing G to T transversions and G to A transitions in mammalian systems (44–52), if error-prone incorporation opposite the lesion involved the human DinB polymerase pol κ, and its structural properties prove similar to those of its prokaryotic homolog Dpo4.

In conclusion, our modeling and MD simulations for dG-C8-PhIP suggest that the adduct would increase the infidelity of Dpo4 and hinder translocation by the enzyme. We hope that the hypotheses resulting from our modeling studies will provide useful suggestions for future experimental investigations.

SUPPLEMENTARY DATA

Supplementary Data are available at NAR online.