X-Ray analysis of d(CGCGAATTXGCG)(2) containing a 2' or minute-deoxy-N(4)-methoxycytosine residue at X: a characteristic pattern of sugar puckers in the crystalline state of the Dickerson-Drew type DNA dodecamers

In a series of structural studies on damaged DNA, a modified Dickerson-Drew dodecamer with the sequence d(CGCGAATTmo(4)CGCG), where mo(4)C is 2'-deoxy-N(4)-methoxycytidine, was synthesized and its structure in a new crystal form has been determined by the X-ray diffraction method. The two dodecamers form a B-form duplex, in which the two mo(4)C residues, respectively, form a wobble pair and a Watson-Crick type pair with the guanine residues of the opposite strand. A comparison of the sugar conformations with those of the other related Dickerson-Drew dodecamers indicates a common feature of their puckering patterns. The sugar pucker of the third residue always adopts an intermediate state (C4'-exo-O4'-endo) between the A-form and B-form. This deviation is ascribed to the stacking interaction of the ribose ring at the third residue with the guanine base at the 12th residue, which is brought about by an extra G12:G2 interaction between two duplexes related by a crystallographic 2(1) symmetry.     

X-ray analyses of two DNA dodecamers containing N4-methoxy-cytosine paired with adenine or guanine.

To investigate mismatch of base-pairings in relation to mutagenesis by oxyamines, crystal structures of two DNA dodecamers with the sequence d(CGCZAA TTmo4CGCG) (Z = A or G), containing N4-methoxy-cytosine (mo4C), have been determined by X-ray analysis. These dodecamers essentially form right-handed B-form duplexes, respectively. In the dodecamer with Z = A, the two mo4C residues are adapted in imino form with the anti methoxyl group to form pairs with A on the opposite strand in a manner of Watson-Crick fashion. While in the dodecamer with Z = G, one mo4C in amino form with the anti methoxyl group forms a normal Watson-Crick pair with G, but the other one in imino form with syn methoxyl conformation wobbles with G. Based on these results, possible mutation mechanism has been proposed.     

Influence of an exocyclic guanine adduct on the thermal stability, conformation, and melting thermodynamics of a DNA duplex.

As part of an overall program to characterize the impact of mutagenic lesions on the physiochemical properties of DNA, we report here the results of a comparative spectroscopic study on pairs of DNA duplexes both with and without an exocyclic guanine lesion. Specifically, we have studied a family of four 13-mer duplexes of the form d(CGCATGYGTACGC).d(GCGTACZCATGCG) in which Y is either the normal deoxyguanosine residue (G) or the exocyclic guanine adduct 1,N2-propanodeoxyguanosine (X), while Z is either deoxycytosine (C) or deoxyadenosine (A). Thus, the four duplexes studied, which can be designated by the identity of their central Y.Z base pair, are a Watson-Crick duplex (GC), a duplex with a central mismatch (GA), and two duplexes with exocyclic guanine lesions (X), that differ only by the base opposite the lesion (XC and XA). The data derived from our spectroscopic measurements on these four duplexes have allowed us to evaluate the influence of the exocyclic guanine lesion, as well as the base opposite the lesion, on the conformation, thermal stability, and melting energetics of the host DNA duplex. To be specific, our circular dichroism (CD) spectra show that the exocyclic guanine lesion induces alterations in the duplex structure, while our temperature-dependent optical measurements reveal that these lesion-induced structural alterations reduce the thermal stability, the transition enthalpy, and the transition free energy of the duplex.(ABSTRACT TRUNCATED AT 250 WORDS)     

A functional analysis of the DNA glycosylase activity of mouse MUTYH protein excising 2-hydroxyadenine opposite guanine in DNA

2-Hydroxy-2-deoxyadenosine triphosphate (2-OH-dATP), generated by the oxidation of dATP, can be misincorporated by DNA polymerases opposite guanine in template DNA during DNA replication, thus causing spontaneous mutagenesis. We demonstrated that mouse MUTYH (mMUTYH) has a DNA glycosylase activity excising not only adenine opposite 8-oxoguanine (8-oxoG) but also 2-hydroxyadenine (2-OH-A) opposite guanine, using purified recombinant thioredoxin-mMUTYH fusion protein. mMUTYH formed a stable complex with duplex oligonucleotides containing an adenine:8-oxoG pair, but the binding of mMUTYH to oligonucleotides containing a 2-OH-A:guanine pair was barely detectable, thus suggesting that mMUTYH recognizes and interacts with these two substrates in a different manner which may reflect the difference in the base excision repair process for each substrate. Mutant mMUTYH with G365D amino acid substitution, corresponding to a G382D germline mutation of human MUTYH found in familial adenomatous polyposis patients, almost completely retained its DNA glycosylase activity excising adenine opposite 8-oxoG; however, it possessed 1.5% of the wild-type activity excising 2-OH-A opposite guanine. Our results imply that the reduced repair capacity of the mutant hMUTYH(G382D), which inefficiently excises 2-OH-A opposite guanine, results in an increased occurrence of somatic G:C to T:A transversion mutations in the APC gene as well as tumorigenesis in the colon.     

Crystallographic studies on damaged DNAs. I. An N(6)-methoxyadenine residue forms a Watson-Crick pair with a cytosine residue in a B-DNA duplex

Oxyamines such as hydroxylamine and methoxylamine disturb DNA replication and act as potent mutagens, causing nucleotide transition from one purine to another or one pyrimidine to another. In order to investigate mismatch base-pairing in DNA damaged with oxyamines, a dodecamer with the sequence d(CGCGmo(6)AATCCGCG), where mo(6) A is 2'-deoxy-N(6)-methoxyadenosine, was synthesized and its crystal structure determined. No significant conformational changes are found between the present dodecamer and the original undamaged B-form dodecamer. Electron density maps clearly show that the mo(6)A residue forms a base-pair with a 2'-deoxycytidine residue through hydrogen bonds similar to a Watson-Crick G.C base-pair. For these hydrogen bonds to be made, N(6)-methoxyadenine must chemically take the imino form. The methoxylation thus enables the adenine base to mimic a guanine base. As a result, misincorporation of 2'-deoxycytidine instead of thymidine, or 2'-deoxyadenosine instead of 2'-deoxyguanosine, can occur in DNA replication.     

Efficient and error-free replication past a minor-groove DNA adduct by the sequential action of human DNA polymerases iota and kappa

DNA polymerase iota (Poliota) is a member of the Y family of DNA polymerases, which promote replication through DNA lesions. The role of Poliota in lesion bypass, however, has remained unclear. Poliota is highly unusual in that it incorporates nucleotides opposite different template bases with very different efficiencies and fidelities. Since interactions of DNA polymerases with the DNA minor groove provide for the nearly equivalent efficiencies and fidelities of nucleotide incorporation opposite each of the four template bases, we considered the possibility that Poliota differs from other DNA polymerases in not being as sensitive to distortions of the minor groove at the site of the incipient base pair and that this enables it to incorporate nucleotides opposite highly distorting minor-groove DNA adducts. To check the validity of this idea, we examined whether Poliota could incorporate nucleotides opposite the gamma-HOPdG adduct, which is formed from an initial reaction of acrolein with the N(2) of guanine. We show here that Poliota incorporates a C opposite this adduct with nearly the same efficiency as it does opposite a nonadducted template G residue. The subsequent extension step, however, is performed by Polkappa, which efficiently extends from the C incorporated opposite the adduct. Based upon these observations, we suggest that an important biological role of Poliota and Polkappa is to act sequentially to carry out the efficient and accurate bypass of highly distorting minor-groove DNA adducts of the purine bases.     

Solution structures of aminofluorene [AF]-stacked conformers of the syn [AF]-C8-dG adduct positioned opposite dC or dA at a template-primer junction.

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.     

Solution structure by NMR and molecular dynamics of a duplex containing a guanine opposite a N-(2-deoxy-beta-D-erythro-pentofuranosyl)formamide lesion

One- and two-dimensional NMR spectroscopy has been used combined with molecular dynamics to determine the fine structure of the DNA duplex 5'-d(AGGAGCCACG).d(CGTGGFTCCT) where F is the N-(2-deoxy-beta-D-erythro-pentofuranosyl)formamide residue which is a ring fragmentation product of thymine. The formamide deoxyribose exists as two isomers with respect to the orientation about the peptide bond. The two isomers (trans and cis) are observed in a ratio 3:2 in solution. For both species, the oligonucleotide adopts a globally B form structure although conformational changes are observed around the mismatch site. The formamide residue, whatever the isomer, is intrahelical and can pair with the guanine on the opposite strand with one hydrogen bond. For the cis isomer, the residue adopts a syn orientation and is able to form a second hydrogen bond with the guanine on the 5' side on the same strand. Off-resonance ROESY experiments have been used to investigate the chemical exchange observed at low temperature of the duplex. Conformational exchange has only been found for the oligonucleotide with the formamide residue in the trans conformation.     

Structure of an oligonucleotide containing a N-(2-deoxy-beta-D-erythro-pentofuranosyl)formamide residue facing a guanine

Formamide residue is a major oxidative DNA damage product from ionizing radiation on thymine residues in DNA. We report NMR and molecular modeling studies on a DNA duplex structure which contains guanine opposite formamide residue. Formamide residue exists as either the cis and trans isomer. For the trans and the cis isomers, we find that guanine and formamide are stacked inside the helix and are hydrogen bonded. The oligonucleotide adopts globally a B form structure for the two isomers. Conformational changes are observed between the two isomers.     

Interactions between an anthracycline antibiotic and DNA: molecular structure of daunomycin complexed to d(CpGpTpApCpG) at 1.2-A resolution

The crystal structure of a daunomycin-d(CGTACG) complex has been solved by X-ray diffraction analysis and refined to a final R factor of 0.175 at 1.2-A resolution. The crystals are in a tetragonal crystal system with space group P4(1)2(1)2 and cell dimensions of a = b = 27.86 A and c = 52.72 A. The self-complementary DNA forms a six base pair right-handed double helix with two daunomycin molecules intercalated in the d(CpG) sequences at either end of the helix. Daunomycin in the complex has a conformation different from that of daunomycin alone. The daunomycin aglycon chromophore is oriented at right angles to the long dimension of the DNA base pairs, and the cyclohexene ring A rests in the minor groove of the double helix. Substituents on this ring have hydrogen-bonding interactions to the base pairs above and below the intercalation site. O9 hydroxyl group of the daunomycin forms two hydrogen bonds with N3 and N2 of an adjacent guanine base. Two bridging water molecules between the drug and DNA stabilize the complex in the minor groove. In the major groove, a hydrated sodium ion is coordinated to N7 of the terminal guanine and the O4 and O5 of daunomycin with a distorted octahedral geometry. The amino sugar lies in the minor groove without bonding to the DNA. The DNA double helix is distorted with an asymmetrical rearrangement of the backbone conformation surrounding the intercalator drug. The sugar puckers are C1,C2'-endo, G2,C1'-endo, C11,C1'-endo, and G12,C3'-exo. Only the C1 residue has a normal anti-glycosyl torsion angle (chi = -154 degrees), while the other three residues are all in the high anti range (average chi = -86 degrees). This structure allows us to identify three principal functional components of anthracycline antibiotics: the intercalator (rings B-D), the anchoring functions associated with ring A, and the amino sugar. The structure-function relationships of daunomycin binding to DNA as well as other related anticancer drugs are discussed.     

Incorporation of oxidized guanine nucleoside 5'-triphosphates in DNA with DNA polymerases and preparation of single-lesion carrying DNA

We investigated the incorporation of oxidatively modified guanine residues in DNA using three DNA polymerases, Escherichia coli Kf exo+, Kf exo-, and Taq DNA polymerase. We prepared nucleoside 5'-triphosphates with modified bases (dN (ox)TP) including imidazolone associated with oxazolone (dIzTP/dZTP), dehydroguanidinohydantoin (dOGhTP), and oxaluric acid (dOxaTP). We showed that the single-nucleotide incorporation of these dN (ox)TP at the 3'-end of a primer DNA strand was possible opposite C or G for dIzTP/dZTP, opposite C for dOGhTP using the Klenow fragment, and opposite C for dOxaTP using Taq. The efficiency of these misincorporations was compared to that of the nucleoside 5'-triphosphate modified with the mutagenic guanine lesion 8-oxo-G opposite A or C as well as to that of the natural dNTPs. The reaction was found not competitive. However, the ability of Kf exo- to further copy the whole template DNA strand from the primer carrying one modified residue at the 3'-end proved to be easy and rapid. The two-step polymerization process consisting of the single-nucleotide extension followed by the full extension of a primer afforded a method for the preparation of tailored double-stranded DNA oligonucleotides carrying a single modified base at a precise site on any sequence. This very rapid method allowed the incorporation of unique residues in DNA that were not available before due to their unstable character.     

A unique unnatural base pair between a C analogue, pseudoisocytosine, and an A analogue, 6-methoxypurine, in replication

Pseudoisocytidine, a C-nucleoside analogue of cytosine, has two possible isomers of the H1- and H3-forms. Enzymatic incorporation experiments confirmed the existence of the two isomers in solution, and the 2'-deoxyribonucleoside triphosphate of pseudoisocytosine (PIC) was incorporated into DNA opposite both guanine and 6-methoxypurine (M) by the Klenow fragment of Escherichia coli DNA polymerase I. In addition to the PIC*M pairing in replication, M also functioned as an A analogue and T was efficiently incorporated opposite M. Thus, the PIC*M pair is regarded as a base pair between a C analogue and an A analogue, and can mediate the interconversion between the G*C and A*T base pairs. The combination of PIC and M could be used as a G*C<-->A*T transition mutagen.     

Crystallographic studies on damaged DNAs. II. N(6)-methoxyadenine can present two alternate faces for Watson-Crick base-pairing, leading to pyrimidine transition mutagenesis

In a previous paper, 2'-deoxy-N(6)-methoxyadenosine (mo(6)A) was shown to form a mismatch base-pair with 2'-deoxycytidine with a Watson-Crick-type geometry. To fully understand the structural basis of genetic mutations with damaged DNA, it is necessary to examine whether the methoxylated adenine residue still has the ability to form the regular Watson-Crick pairing with a thymine residue. Therefore, a DNA dodecamer with the sequence d(CGCGmo(6)AATTCGCG) has been synthesized and its crystal structure determined. The methoxylation has no significant effect on the overall DNA conformation, which is that of a standard B-form duplex. The methoxylated adenine moieties adopt the amino tautomer with an anti conformation around the C(6)-N(6) bond to the N(1) atom, and they form a Watson-Crick base-pair with thymine residues on the opposite strand, similar to an unmodified adenine residue. It is concluded that methoxylated adenine can present two alternate faces for base-pairing, thanks to the amino<-->imino tautomerism allowed by methoxylation. Based on this property, two gene transition routes are proposed.     

Propeller-twisted adenine.thymine and guanine.cytosine base pairs tend to buckle and stagger in opposite directions.

Base pairs are propeller-twisted, buckled and staggered in DNA fragment crystals. These deformations were analyzed with isolated Watson-Crick base pairs using empirical potentials and buckle was found to almost linearly correlate with propeller. Interestingly, the thymine.adenine pair favours negative buckling for propellers mostly observed in DNA crystals while positive buckling is preferred by the cytosine.guanine pair. The propeller also induces opposite staggers in the adenine.thymine and guanine.cytosine base pairs.     

Mechanisms of dCMP transferase reactions catalyzed by mouse Rev1 protein.

The Rev1 protein, a member of a large family of translesion DNA polymerases, catalyzes a dCMP transfer reaction. Recombinant mouse Rev1 protein was found to insert a dCMP residue opposite guanine, adenine, thymine, cytosine, uracil, and an apurinic/apyrimidinic site and to have weak ability for transfer to a mismatched terminus. The mismatch-extension ability was strongly enhanced by a guanine residue on the template near the mismatched terminus; this was not the case with an apurinic/apyrimidinic site and the other template nucleotides. Kinetic analysis of the dCMP transferase reaction provided evidence for high affinity for dCTP with template G but not the other templates, whereas the template nucleotide did not much affect the V(max) value. Furthermore, it could be established that the mouse Rev1 protein inserts dGMP and dTMP residues opposite template guanine at a V(max) similar to that for dCMP.     

Mutagenicity associated with O6-methylguanine-DNA damage and mechanism of nucleotide flipping by AGT during repair

Methylated guanine damage at O6 position (i.e. O6MG) is dangerous due to its mutagenic and carcinogenic character that often gives rise to G:C-A:T mutation. However, the reason for this mutagenicity is not known precisely and has been a matter of controversy. Further, although it is known that O6-alkylguanine-DNA alkyltransferase (AGT) repairs O6MG paired with cytosine in DNA, the complete mechanism of target recognition and repair is not known completely. All these aspects of DNA damage and repair have been addressed here by employing high level density functional theory in gas phase and aqueous medium. It is found that the actual cause of O6MG mediated mutation may arise due to the fact that DNA polymerases incorporate thymine opposite to O6MG, misreading the resulting O6MG:T complex as an A:T base pair due to their analogous binding energies and structural alignments. It is further revealed that AGT mediated nucleotide flipping occurs in two successive steps. The intercalation of the finger residue Arg128 into the DNA double helix and its interaction with the O6MG:C base pair followed by rotation of the O6MG nucleotide are found to be crucial for the damage recognition and nucleotide flipping.     

Observing translesion synthesis of an aromatic amine DNA adduct by a high-fidelity DNA polymerase

Aromatic amines have been studied for more than a half-century as model carcinogens representing a class of chemicals that form bulky adducts to the C8 position of guanine in DNA. Among these guanine adducts, the N-(2'-deoxyguanosin-8-yl)-aminofluorene (G-AF) and N-2-(2'-deoxyguanosin-8-yl)-acetylaminofluorene (G-AAF) derivatives are the best studied. Although G-AF and G-AAF differ by only an acetyl group, they exert different effects on DNA replication by replicative and high-fidelity DNA polymerases. Translesion synthesis of G-AF is achieved with high-fidelity polymerases, whereas replication of G-AAF requires specialized bypass polymerases. Here we have presented structures of G-AF as it undergoes one round of accurate replication by a high-fidelity DNA polymerase. Nucleotide incorporation opposite G-AF is achieved in solution and in the crystal, revealing how the polymerase accommodates and replicates past G-AF, but not G-AAF. Like an unmodified guanine, G-AF adopts a conformation that allows it to form Watson-Crick hydrogen bonds with an opposing cytosine that results in protrusion of the bulky fluorene moiety into the major groove. Although incorporation opposite G-AF is observed, the C:G-AF base pair induces distortions to the polymerase active site that slow translesion synthesis.     

Human mitochondrial DNA polymerase γ exhibits potential for bypass and mutagenesis at UV-induced cyclobutane thymine dimers

Cyclobutane thymine dimers (T-T) comprise the majority of DNA damage caused by short wavelength ultraviolet radiation. These lesions generally block replicative DNA polymerases and are repaired by nucleotide excision repair or bypassed by translesion polymerases in the nucleus. Mitochondria lack nucleotide excision repair, and therefore, it is important to understand how the sole mitochondrial DNA polymerase, pol γ, interacts with irreparable lesions such as T-T. We performed in vitro DNA polymerization assays to measure the kinetics of incorporation opposite the lesion and bypass of the lesion by pol γ with a dimer-containing template. Exonuclease-deficient pol γ bypassed thymine dimers with low relative efficiency; bypass was attenuated but still detectable when using exonuclease-proficient pol γ. When bypass did occur, pol γ misincorporated a guanine residue opposite the 3'-thymine of the dimer only 4-fold less efficiently than it incorporated an adenine. Surprisingly, the pol γ exonuclease-proficient enzyme excised the incorrectly incorporated guanine at similar rates irrespective of the nature of the thymines in the template. In the presence of all four dNTPs, pol γ extended the primer after incorporation of two adenines opposite the lesion with relatively higher efficiency compared with extension past either an adenine or a guanine incorporated opposite the 3'-thymine of the T-T. Our results suggest that T-T usually stalls mitochondrial DNA replication but also suggest a mechanism for the introduction of point mutations and deletions in the mitochondrial genomes of chronically UV-exposed cells.     

Toward understanding the mutagenicity of an environmental carcinogen: structural insights into nucleotide incorporation preferences

Bulky carcinogen-DNA adducts, including (+)-trans-anti-[BP]-N(2)-dG derived from the reaction of (+)-anti-benzo[a]pyrene diol epoxide with guanine, often block the progression of DNA polymerases. However, when rare bypass of the lesions does occur, they may be misreplicated. Experimental results have shown that nucleotides are inserted opposite the (+)-trans-anti-[BP]-N(2)-dG adduct by bacteriophage T7 DNA polymerase with the order of preference A>T>or=G>C. To gain structural insights into the effects of the bulky adduct on nucleotide incorporation within the polymerase active site, molecular modeling and molecular dynamics simulations were carried out using T7 DNA polymerase to permit the relation of function to structure. We modeled the (+)-trans-anti-[BP]-N(2)-dG adduct opposite incoming dGTP, dTTP and dCTP nucleotides, as well as unmodified guanine opposite its normal partner dCTP as a control, to compare with our previous simulation with dATP opposite the adduct. The modeling required that the (+)-trans-anti-[BP]-N(2)-dG adduct adopt the syn conformation in each case to avoid deranging essential protein-DNA interactions. While the dATP: (+)-trans-anti-[BP]-N(2)-dG pair was well accommodated within the active site of T7 DNA polymerase, dCTP fit poorly opposite the adduct, adopting an orientation perpendicular to the plane of the syn modified guanine during the simulation. Rotation about the glycosidic bond of the dCTP residue to this abnormal position was allowed because only one hydrogen bond between dCTP and the (+)-trans-anti-[BP]-N(2)-dG residue evolved during the simulation, and this hydrogen bond was directly across from the dCTP glycosidic bond. The dTTP and dGTP nucleotides, incorporated with an intermediate preference opposite (+)-trans-anti-[BP]-N(2)-dG, were accommodated reasonably well, but not as stably as the dATP nucleotide, due to a skewed primer-template alignment and more exposed BP moiety, respectively. In addition, the extent of stabilizing interactions between the nascent base-pair in each simulation was correlated positively with the incorporation preference of that particular nucleotide. The dATP nucleotide is accommodated most stably opposite the adduct, with protein-DNA hydrogen bonding interactions and an active-site pocket size that do not deviate significantly from those of the control simulation. The simulations of dTTP and dGTP opposite (+)-trans-anti-[BP]-N(2)-dG exhibited more instability in interactions between the protein and the nascent base-pair than the dATP system. However, the active-site pocket size of the dTTP and dGTP simulations remained stable. The dCTP: (+)-trans-anti-[BP]-N(2)-dG system had the least number of stabilizing interactions, and the active-site pocket of this system increased in size significantly compared to the control and other dNTPs opposite the adduct. These simulations elucidated why A is inserted opposite (+)-trans-anti-[BP]-N(2)-dG most frequently, while T and G are inserted opposite the adduct to an extent intermediate between A and C, and C is most rarely incorporated. Structural rationalization of the incorporation preference opposite (+)-trans-anti-[BP]-N(2)-dG by T7 DNA polymerase contributes to providing a molecular explanation for mutations caused by this carcinogen-DNA adduct in a model system.