The spacious active site of a Y-family DNA polymerase facilitates promiscuous nucleotide incorporation opposite a bulky carcinogen-DNA adduct: elucidating the structure-function relationship through experimental and computational approaches

Effects of N7-Alkylguanine Conformation and Metal Cofactors on the Translesion Synthesis by Human DNA Polymerase η

Non-enzymatic alkylation on DNA often generates N7-alkyl-2'-deoxyguanosine (N7alkylG) adducts as major lesions. N7alkylG adducts significantly block replicative DNA polymerases and can be bypassed by translesion synthesis (TLS) polymerases such as polymerase η (polη). To gain insights into the bypass of N7alkylG by TLS polymerases, we conducted kinetic and structural studies of polη catalyzing across N7BnG, a genotoxic lesion generated by the carcinogenic N-nitrosobenzylmethylamine. The presence of templating N7BnG in the polη catalytic site decreased the replication fidelity by ∼9-fold, highlighting the promutagenicity of N7BnG. The catalytic efficiency for dCTP incorporation opposite N7BnG decreased ∼22-fold and ∼7-fold compared to the incorporation opposite undamaged guanine in the presence of Mg²⁺ and Mn²⁺, respectively. A crystal structure of the complexes grown with polη, templating N7BnG, incoming dCTP, and Mg²⁺ ions showed the lack of the incoming nucleotide and metal cofactors in the polη catalytic site. Interestingly, the templating N7BnG adopted a syn conformation, which has not been observed in the published N7alkylG structures. The preferential formation of syn-N7BnG conformation at the templating site may deter the binding of an incoming dCTP, causing the inefficient bypass by polη. In contrast, the use of Mn²⁺ in place of Mg²⁺ in co-crystallization yielded a ternary complex displaying an anti-N7BnG:dCTP base pair and catalytic metal ions, which would be a close mimic of a catalytically competent state. We conclude that certain bulky N7-alkylG lesions can slow TLS polymerase-mediated bypass by adopting a catalytically unfavorable syn conformation in the replicating base pair site.

Mutagenic Replication of N2-Deoxyguanosine Benzo[a]pyrene Adducts by Escherichia coli DNA Polymerase I and Sulfolobus solfataricus DNA Polymerase IV

Benzo[a]pyrene, a potent human carcinogen, is metabolized in vivo to a diol epoxide that reacts with the N²-position of guanine to produce N²-BP-dG adducts. These adducts are mutagenic causing G to T transversions. These adducts block replicative polymerases but can be bypassed by the Y-family translesion synthesis polymerases. The mechanisms by which mutagenic bypass occurs is not well-known. We have evaluated base pairing structures using atomic substitution of the dNTP with two stereoisomers, 2'-deoxy-N-[(7R,8S,9R,10S)-7,8,9,10-tetrahydro-7,8,9-trihydroxybenzo[a]pyren-10-yl]guanosine and 2'-deoxy-N-[(7S,8R,9S,10R)-7,8,9,10-tetrahydro-7,8,9-trihydroxybenzo[a]pyren-10-yl]guanosine. We have examined the kinetics of incorporation of 1-deaza-dATP, 7-deaza-dATP, 2'-deoxyinosine triphosphate, and 7-deaza-dGTP, analogues of dATP and dGTP in which single atoms are changed. Changes in rate will occur if that atom provided a critical interaction in the transition state of the reaction. We examined two polymerases, Escherichia coli DNA polymerase I (Kf) and Sulfolobus solfataricus DNA polymerase IV (Dpo4), as models of a high fidelity and TLS polymerase, respectively. We found that with Kf, substitution of the nitrogens on the Watson-Crick face of the dNTPs resulted in decreased rate of reactions. This result is consistent with a Hoogsteen base pair in which the template N²-BP-dG flipped from the anti to syn conformation. With Dpo4, while the substitution did not affect the rate of reaction, the amplitude of the reaction decreased with all substitutions. This result suggests that Dpo4 bypasses N²-BP-dG via Hoogsteen base pairs but that the flipped nucleotide can be either the dNTP or the template.

Computational Evaluation of Nucleotide Insertion Opposite Expanded and Widened DNA by the Translesion Synthesis Polymerase Dpo4

Expanded (x) and widened (y) deoxyribose nucleic acids (DNA) have an extra benzene ring incorporated either horizontally (xDNA) or vertically (yDNA) between a natural pyrimidine base and the deoxyribose, or between the 5- and 6-membered rings of a natural purine. Far-reaching applications for (x,y)DNA include nucleic acid probes and extending the natural genetic code. Since modified nucleobases must encode information that can be passed to the next generation in order to be a useful extension of the genetic code, the ability of translesion (bypass) polymerases to replicate modified bases is an active area of research. The common model bypass polymerase DNA polymerase IV (Dpo4) has been previously shown to successfully replicate and extend past a single modified nucleobase on a template DNA strand. In the current study, molecular dynamics (MD) simulations are used to evaluate the accommodation of expanded/widened nucleobases in the Dpo4 active site, providing the first structural information on the replication of (x,y)DNA. Our results indicate that the Dpo4 catalytic (palm) domain is not significantly impacted by the (x,y)DNA bases. Instead, the template strand is displaced to accommodate the increased C1’–C1’ base-pair distance. The structural insights unveiled in the present work not only increase our fundamental understanding of Dpo4 replication, but also reveal the process by which Dpo4 replicates (x,y)DNA, and thereby will contribute to the optimization of high fidelity and efficient polymerases for the replication of modified nucleobases.

Biochemical characterization of eight genetic variants of human DNA polymerase κ involved in error-free bypass across bulky N(2)-guanyl DNA adducts

DNA polymerase (pol) κ, one of the Y-family polymerases, has been shown to function in error-free translesion DNA synthesis (TLS) opposite the bulky N(2)-guanyl DNA lesions induced by many carcinogens such as polycyclic aromatic hydrocarbons. We analyzed the biochemical properties of eight reported human pol κ variants positioned in the polymerase core domain, using the recombinant pol κ (residues 1-526) protein and the DNA template containing an N(2)-CH2(9-anthracenyl)G (N(2)-AnthG). The truncation R219X was devoid of polymerase activity, and the E419G and Y432S variants showed much lower polymerase activity than wild-type pol κ. In steady-state kinetic analyses, E419G and Y432S displayed 20- to 34-fold decreases in kcat/Km for dCTP insertion opposite G and N(2)-AnthG compared to that of wild-type pol κ. The L21F, I39T, and D189G variants, as well as E419G and Y432S, displayed 6- to 22-fold decreases in kcat/Km for next-base extension from C paired with N(2)-AnthG, compared to that of wild-type pol κ. The defective Y432S variant had 4- to 5-fold lower DNA-binding affinity than wild-type, while a slightly more efficient S423R variant possessed 2- to 3-fold higher DNA-binding affinity. These results suggest that R219X abolishes and the E419G, Y432S, L21F, I39T, and D189G variations substantially impair the TLS ability of pol κ opposite bulky N(2)-G lesions in the insertion step opposite the lesion and/or the subsequent extension step, raising the possibility that certain nonsynonymous pol κ genetic variations translate into individual differences in susceptibility to genotoxic carcinogens.

Variants of mouse DNA polymerase κ reveal a mechanism of efficient and accurate translesion synthesis past a benzo[a]pyrene dG adduct

DNA polymerase κ (Polκ) is the only known Y-family DNA polymerase that bypasses the 10S (+)-trans-anti-benzo[a]pyrene diol epoxide (BPDE)-N(2)-deoxyguanine adducts efficiently and accurately. The unique features of Polκ, a large structure gap between the catalytic core and little finger domain and a 90-residue addition at the N terminus known as the N-clasp, may give rise to its special translesion capability. We designed and constructed two mouse Polκ variants, which have reduced gap size on both sides [Polκ Gap Mutant (PGM) 1] or one side flanking the template base (PGM2). These Polκ variants are nearly as efficient as WT in normal DNA synthesis, albeit with reduced accuracy. However, PGM1 is strongly blocked by the 10S (+)-trans-anti-BPDE-N(2)-dG lesion. Steady-state kinetic measurements reveal a significant reduction in efficiency of dCTP incorporation opposite the lesion by PGM1 and a moderate reduction by PGM2. Consistently, Polκ-deficient cells stably complemented with PGM1 GFP-Polκ remained hypersensitive to BPDE treatment, and complementation with WT or PGM2 GFP-Polκ restored BPDE resistance. Furthermore, deletion of the first 51 residues of the N-clasp in mouse Polκ (mPolκ(52-516)) leads to reduced polymerization activity, and the mutant PGM2(52-516) but not PGM1(52-516) can partially compensate the N-terminal deletion and restore the catalytic activity on normal DNA. However, neither WT nor PGM2 mPolκ(52-516) retains BPDE bypass activity. We conclude that the structural gap physically accommodates the bulky aromatic adduct and the N-clasp is essential for the structural integrity and flexibility of Polκ during translesion synthesis.

Sequence context effect on strand slippage in natural DNA primer-templates

Strand slippage has been found to occur in primer-templates containing a templating thymine, cytosine, and guanine, leading to the formation of misaligned structures with a single-nucleotide bulge. If remained in the active site of low-fidelity polymerases during DNA replication, these misaligned structures can ultimately bring about deletion mutations. In this study, we performed NMR investigations on primer-template models containing a templating adenine. Similar to our previous results on guanine, adenine templates are also less prone to strand slippage than pyrimidine templates. Misalignment occurs only in primer-templates that form a terminal C·G or G·C base pair. Together with our previous findings on thymine, cytosine, and guanine templates, the present study reveals strand slippage can occur in any kind of natural templating bases during DNA replication, providing insights into the origin of mutation hotspots in natural DNA sequences. In addition to the type of incoming base upon misincorporation, the propensity of strand slippage in primer-templates depends also on the type of templating base, its upstream and downstream bases.

NMR investigation of DNA primer-template models: guanine templates are less prone to strand slippage upon misincorporation

Misaligned structures can result from strand slippage during DNA replication and, if not repaired, would lead to mutations. Previously, we showed that strand slippage can occur upon misincorporation of a dNTP opposite thymine and cytosine templates, resulting in a misaligned structure with a T- or C-bulge. The formation propensity for misaligned structures was found to depend on the type of terminal base pair. In this study, we performed NMR investigations on primer-template models containing a guanine template. Our results reveal guanine templates are less prone to strand slippage than pyrimidine templates. Misalignment was found to occur only in 5'-CG templates with a downstream purine. In addition to the significance of terminal base pair and upstream nucleotide, the present study reveals the importance of the templating base and its downstream nucleotide, which also determine the propensity of strand slippage in primer-templates.

Visualizing sequence-governed nucleotide selectivities and mutagenic consequences through a replicative cycle: processing of a bulky carcinogen N2-dG lesion in a Y-family DNA polymerase

Understanding how DNA polymerases process lesions remains fundamental to determining the molecular origins of mutagenic translesion bypass. We have investigated how a benzo[a]pyrene-derived N(2)-dG adduct, 10S-(+)-trans-anti-[BP]-N(2)-dG ([BP]G*), is processed in Dpo4, the well-characterized Y-family bypass DNA polymerase. This polymerase has a slippage-prone spacious active site region. Experimental results in a 5'-C[BP]G*G-3' sequence context reveal significant selectivity for dGTP insertion that predominantly yields -1 deletion extension products. A less pronounced error-prone nonslippage pathway that leads to full extension products with insertion of A > C > G opposite the lesion is also observed. Molecular modeling and dynamics simulations follow the bypass of [BP]G* through an entire replication cycle for the first time in Dpo4, providing structural interpretations for the experimental observations. The preference for dGTP insertion is explained by a 5'-slippage pattern in which the unmodified G rather than G* is skipped, the incoming dGTP pairs with the C on the 5'-side of G*, and the -1 deletion is produced upon further primer extension which is more facile than nucleotide insertion. In addition, the simulations suggest that the [BP]G* may undergo an anti/syn conformational rearrangement during the stages of the replication cycle. In the minor nonslippage pathway, the nucleotide insertion preferences opposite the lesion are explained by relative distortions to the active site region. These structural insights, provided by the modeling and dynamics studies, augment kinetic and limited available crystallographic investigations with bulky lesions, by providing molecular explanations for lesion bypass activities over an entire replication cycle.

Genomic analysis of cancer tissue reveals that somatic mutations commonly occur in a specific motif

Somatic mutations are hallmarks of cancer progression. We sequenced 26 matched human prostate tumor and constitutional DNA samples for somatic alterations in the SRD5A2, HPRT, and HSD3B2 genes, and identified 71 nucleotide substitutions. Of these substitutions, 79% (56/71) occur within a WKVnRRRnVWK sequence (a novel motif we call THEMIS [from the ancient Greek goddess of prophecy]: W=A/T, K=G/T, V=G/A/C, R=purine (A/G), and n=any nucleotide), with one mismatch allowed. Literature searches identified this motif with one mismatch allowed in 66% (37/56) of the somatic prostate cancer mutations and in 74% (90/122) of the somatic breast cancer mutations found in all human genes analyzed. We also found the THEMIS motif with one allowed mismatch in 88% (23/26) of the ras1 gene somatic mutations formed in the sensitive to skin carcinogenesis (SENCAR) mouse model, after induction of error-prone DNA repair following mutagenic treatment. The high prevalence of the motif in each of the above mentioned cases cannot be explained by chance (P<0.046). We further identified 27 somatic mutations in the error-prone DNA polymerase genes pol eta, pol kappa, and pol beta in these prostate cancer patients. The data suggest that most somatic nucleotide substitutions in human cancer may occur in sites that conform to the THEMIS motif. These mutations may be caused by "mutator" mutations in error-prone DNA polymerase genes.

Versatility of Y-family Sulfolobus solfataricus DNA polymerase Dpo4 in translesion synthesis past bulky N2-alkylguanine adducts

In contrast to replicative DNA polymerases, Sulfolobus solfataricus Dpo4 showed a limited decrease in catalytic efficiency (k(cat)/Km) for insertion of dCTP opposite a series of N2-alkylguanine templates of increasing size from (methyl (Me) to (9-anthracenyl)-Me (Anth)). Fidelity was maintained with increasing size up to (2-naphthyl)-Me (Naph). The catalytic efficiency increased slightly going from the N2-NaphG to the N2-AnthG substrate, at the cost of fidelity. Pre-steady-state kinetic bursts were observed for dCTP incorporation throughout the series (N2-MeG to N2-AnthG), with a decrease in the burst amplitude and k(pol), the rate of single-turnover incorporation. The pre-steady-state kinetic courses with G and all of the six N2-alkyl G adducts could be fit to a general DNA polymerase scheme to which was added an inactive complex in equilibrium with the active ternary Dpo4.DNA.dNTP complex, and only the rates of equilibrium with the inactive complex and phosphodiester bond formation were altered. Two crystal structures of Dpo4 with a template N2-NaphG (in a post-insertion register opposite a 3'-terminal C in the primer) were solved. One showed N2-NaphG in a syn conformation, with the naphthyl group located between the template and the Dpo4 "little finger" domain. The Hoogsteen face was within hydrogen bonding distance of the N4 atoms of the cytosine opposite N2-NaphG and the cytosine at the -2 position. The second structure showed N2-Naph G in an anti conformation with the primer terminus largely disordered. Collectively these results explain the versatility of Dpo4 in bypassing bulky G lesions.

The N-clasp of human DNA polymerase kappa promotes blockage or error-free bypass of adenine- or guanine-benzo[a]pyrenyl lesions

DNA bypass polymerases are utilized to transit bulky DNA lesions during replication, but the process frequently causes mutations. The structural origins of mutagenic versus high fidelity replication in lesion bypass is therefore of fundamental interest. As model systems, we investigated the molecular basis of the experimentally observed essentially faithful bypass of the guanine 10S-(+)-trans-anti-benzo[a]pyrene-N(2)-dG adduct by the Y-family human DNA polymerase kappa, and the observed blockage of pol kappa produced by the adenine 10S-(+)-trans-anti-benzo[a]pyrene-N(2)-dA adduct. These lesions are derived from the most tumorigenic metabolite of the ubiquitous cancer-causing pollutant, benzo[a]pyrene. We compare our results for the dG adduct with our earlier studies for the pol kappa archaeal homolog Dpo4, which processes the same lesion in an error-prone manner. Molecular modeling, molecular mechanics calculations and molecular dynamics simulations were utilized. Our results show that the pol kappa N-clasp is a key structural feature that accounts for the dA adduct blockage and the near-error-free bypass of the dG lesion. Absence of the N-clasp in Dpo4 explains the error-prone processing of the same lesion by this enzyme. Thus, our studies elucidate structure-function relationships in the fidelity of lesion bypass.

Quantum mechanics/molecular mechanics investigation of the chemical reaction in Dpo4 reveals water-dependent pathways and requirements for active site reorganization

The nucleotidyl-transfer reaction coupled with the conformational transitions in DNA polymerases is critical for maintaining the fidelity and efficiency of DNA synthesis. We examine here the possible reaction pathways of a Y-family DNA polymerase, Sulfolobus solfataricus DNA polymerase IV (Dpo4), for the correct insertion of dCTP opposite 8-oxoguanine using the quantum mechanics/molecular mechanics (QM/MM) approach, both from a chemistry-competent state and a crystal closed state. The latter examination is important for understanding pre-chemistry barriers to interpret the entire enzyme mechanism, since the crystal closed state is not an ideal state for initiating the chemical reaction. The most favorable reaction path involves initial deprotonation of O3'H via two bridging water molecules to O1A, overcoming an overall potential energy barrier of approximately 20.0 kcal/mol. The proton on O1A-P(alpha) then migrates to the gamma-phosphate oxygen of the incoming nucleotide as O3' attacks P(alpha), and the P(alpha)-O3A bond breaks. The other possible pathway in which the O3'H proton is transferred directly to O1A on P(alpha) has an overall energy barrier of 25.0 kcal/mol. In both reaction paths, the rate-limiting step is the initial deprotonation, and the trigonal-bipyramidal configuration for P(alpha) occurs during the concerted bond formation (O3'-P(alpha)) and breaking (P(alpha)-O3A), indicating the associative nature of the chemical reaction. In contrast, the Dpo4/DNA complex with an imperfect active-site geometry corresponding to the crystal state must overcome a much higher activation energy barrier (29.0 kcal/mol) to achieve a tightly organized site due to hindered O3'H deprotonation stemming from larger distances and distorted conformation of the proton acceptors. This significant difference demonstrates that the pre-chemistry reorganization in Dpo4 costs approximately 4.0 to 9.0 kcal/mol depending on the primer terminus environment. Compared to the higher fidelity DNA polymerase beta from the X-family, Dpo4 has a higher chemical reaction barrier (20.0 vs 15.0 kcal/mol) due to the more solvent-exposed active site.

Lesion processing: high-fidelity versus lesion-bypass DNA polymerases

When a high-fidelity DNA polymerase encounters certain DNA-damage sites, its progress can be stalled and one or more lesion-bypass polymerases are recruited to transit the lesion. Here, we consider two representative types of lesions: (i) 7,8-dihydro-8-oxoguanine (8-oxoG), a small, highly prevalent lesion caused by oxidative damage; and (ii) bulky lesions derived from the environmental pre-carcinogen benzo[a]pyrene, in the high-fidelity DNA polymerase Bacillus fragment (BF) from Bacillus stearothermophilus and in the lesion-bypass DNA polymerase IV (Dpo4) from Sulfolobus solfataricus. The tight fit of the BF polymerase around the nascent base pair contrasts with the more spacious, solvent-exposed active site of Dpo4, and these differences in architecture result in distinctions in their respective functions: one-step versus stepwise polymerase translocation, mutagenic versus accurate bypass of 8-oxoG, and polymerase stalling versus mutagenic bypass at bulky benzo[a]pyrene-derived lesions.

Nucleotide selectivity opposite a benzo[a]pyrene-derived N2-dG adduct in a Y-family DNA polymerase: a 5'-slippage mechanism

The Y-family DNA polymerase Dpo4, from the archaeon bacterium Sulfolobus solfataricus, is a member of the DinB family, which also contains human Pol kappa. It has a spacious active site that can accommodate two templating bases simultaneously, with one of them skipped by the incoming dNTP. Assays of single dNTP insertion opposite a benzo[ a]pyrene-derived N (2)-dG adduct, 10 S(+)- trans- anti-[BP]- N (2)-dG ([BP]G*), reveal that an incoming dATP is significantly preferred over the other three dNTPs in the TG 1*G 2 sequence context. Molecular modeling and dynamics simulations were carried out to interpret this experimental observation on a molecular level. Modeling studies suggest that the significant preference for dATP insertion observed experimentally can result from two possible dATP incorporation modes. The dATP can be inserted opposite the T on the 5' side of the adduct G 1*, using an unusual 5'-slippage pattern, in which the unadducted G 2, rather than G 1*, is skipped, to produce a -1 deletion. In addition, the dATP can be misincorporated opposite the adduct. The 5'-slippage pattern may be generally facilitated in cases where the base 3' to the lesion is the same as the adducted base.

DNA adduct structure-function relationships: comparing solution with polymerase structures

It has now been nearly two decades since the first solution structures of DNA duplexes covalently damaged by metabolically activated polycyclic aromatic hydrocarbons and amines were determined by NMR. Dozens of such high-resolution structures are now available, and some broad structural themes have been uncovered. It has been hypothesized that the solution structures are relevant to the biochemical processing of the adducts. The structural features of the adducts are considered to determine their mutational properties in DNA polymerases and their repair susceptibilities. In recent years, a number of crystal structures of DNA adducts of interest to our work have been determined in DNA polymerases. Accordingly, it is now timely to consider how NMR solution structures relate to structures within DNA polymerases. The NMR solution structural themes for polycyclic aromatic adducts are often observed in polymerase crystal structures. While the polymerase interactions can on occasion override the solution preferences, intrinsic adduct conformations favored in solution are often manifested within polymerases and likely play a significant role in lesion processing.

Optimization of non-natural nucleotides for selective incorporation opposite damaged DNA

The promutagenic process known as translesion DNA synthesis reflects the ability of a DNA polymerase to misinsert a nucleotide opposite a damaged DNA template. To study the underlying mechanism of nucleotide selection during this process, we quantified the incorporation of various non-natural nucleotide analogs opposite an abasic site, a non-templating DNA lesion. Our kinetic studies using the bacteriophage T4 DNA polymerase reveal that the pi-electron surface area of the incoming nucleotide substantially contributes to the efficiency of incorporation opposite an abasic site. A remaining question is whether the selective insertion of these non-hydrogen-bonding analogs can be achieved through optimization of shape and pi-electron density. In this report, we describe the synthesis and kinetic characterization of four novel nucleotide analogs, 5-cyanoindolyl-2'-deoxyriboside 5'-triphosphate (5-CyITP), 5-ethyleneindolyl-2'-deoxyriboside 5'-triphosphate (5-EyITP), 5-methylindolyl-2'-deoxyriboside 5'-triphosphate (5-MeITP), and 5-ethylindolyl-2'-deoxyriboside 5'-triphosphate (5-EtITP). Kinetic analyses indicate that the overall catalytic efficiencies of all four nucleotides are related to their base-stacking properties. In fact, the catalytic efficiency for nucleotide incorporation opposite an abasic site displays a parabolic trend in the overall pi-electron surface area of the non-natural nucleotide. In addition, each non-natural nucleotide is incorporated opposite templating DNA approximately 100-fold worse than opposite an abasic site. These data indicate that selectivity for incorporation opposite damaged DNA can be achieved through optimization of the base-stacking properties of the incoming nucleotide.

What a difference a decade makes: insights into translesion DNA synthesis

Living organisms are continually under attack from a vast array of DNA-damaging agents that imperils their genomic integrity. As a consequence, cells possess an army of enzymes to repair their damaged chromosomes. However, DNA lesions often persist and pose a considerable threat to survival, because they can block the cell's replicase and its ability to complete genome duplication. It has been clear for many years that cells must possess a mechanism whereby the DNA lesion could be tolerated and physically bypassed. Yet it was only within the past decade that specialized DNA polymerases for "translesion DNA synthesis" or "TLS" were identified and characterized. Many of the TLS enzymes belong to the recently described "Y-family" of DNA polymerases. By possessing a spacious preformed active site, these enzymes can physically accommodate a variety of DNA lesions and facilitate their bypass. Flexible DNA-binding domains and a variable binding pocket for the replicating base pair further allow these TLS polymerases to select specific lesions to bypass and favor distinct non-Watson-Crick base pairs. Consequently, TLS polymerases tend to exhibit much lower fidelity than the cell's replicase when copying normal DNA, which results in a dramatic increase in mutagenesis. Occasionally this can be beneficial, but it often speeds the onset of cancer in humans. Cells use both transcriptional and posttranslational regulation to keep these low-fidelity polymerases under strict control and limit their access to a replication fork. Our perspective focuses on the mechanistic insights into TLS by the Y-family polymerases, how they are regulated, and their effects on genomic (in)stability that have been described in the past decade.

A structural gap in Dpo4 supports mutagenic bypass of a major benzo[a]pyrene dG adduct in DNA through template misalignment

Erroneous replication of lesions in DNA by DNA polymerases leads to elevated mutagenesis. To understand the molecular basis of DNA damage-induced mutagenesis, we have determined the x-ray structures of the Y-family polymerase, Dpo4, in complex with a DNA substrate containing a bulky DNA lesion and incoming nucleotides. The DNA lesion is derived from an environmentally widespread carcinogenic polycyclic aromatic hydrocarbon, benzo[a]pyrene (BP). The potent carcinogen BP is metabolized to diol epoxides that form covalent adducts with cellular DNA. In the present study, the major BP diol epoxide adduct in DNA, BP-N(2)-deoxyguanosine (BP-dG), was placed at a template-primer junction. Three ternary complexes reveal replication blockage, extension past a mismatched lesion, and a -1 frameshift mutation. In the productive structures, the bulky adduct is flipped/looped out of the DNA helix into a structural gap between the little finger and core domains. Sequestering of the hydrophobic BP adduct in this new substrate-binding site permits the DNA to exhibit normal geometry for primer extension. Extrusion of the lesion by template misalignment allows the base 5' to the adduct to serve as the template, resulting in a -1 frameshift. Subsequent strand realignment produces a mismatched base opposite the lesion. These structural observations, in combination with replication and mutagenesis data, suggest a model in which the additional substrate-binding site stabilizes the extrahelical nucleotide for lesion bypass and generation of base substitutions and -1 frameshift mutations.

Translesion synthesis past the C8- and N2-deoxyguanosine adducts of the dietary mutagen 2-Amino-3-methylimidazo[4,5-f]quinoline in the NarI recognition sequence by prokaryotic DNA polymerases

2-Amino-3-methylimidazo[4,5-f]quinoline (IQ) is found in cooked meats and forms DNA adducts at the C8- and N2-positions of dGuo after appropriate activation. IQ is a potent inducer of frameshift mutations in bacteria and is carcinogenic in laboratory animals. We have incorporated both IQ-adducts into the G1- and G3-positions of the NarI recognition sequence (5'-G1G2CG3CC-3'), which is a hotspot for arylamine modification. The in vitro replication of the oligonucleotides was examined with Escherichia coli pol I Klenow fragment exo-, E. coli pol II exo-, and Sulfolobus solfataricus P2 DNA polymerase IV (Dpo4), and the extension products were sequenced by tandem mass spectrometry. Replication of the C8-adduct at the G3-position resulted in two-base deletions with all three polymerases, whereas error-free bypass and extension was observed at the G1-position. The N2-adduct was bypassed and extended by all three polymerases when positioned at the G1-position, and the error-free product was observed. The N2-adduct at the G3-position was more blocking and was bypassed and extended only by Dpo4 to produce an error-free product. These results indicate that the replication of the IQ-adducts of dGuo is strongly influenced by the local sequence and the regioisomer of the adduct. These results also suggest a possible role for pol II and IV in the error-prone bypass of the C8-IQ-adduct leading to frameshift mutations in reiterated sequences, whereas noniterated sequences result in error-free bypass.

Use of damaged DNA and dNTP substrates by the error-prone DNA polymerase X from African swine fever virus

The structural specificity that translesion DNA polymerases often show for a particular class of lesions suggests that the predominant criterion of selection during their evolution has been the capacity for lesion tolerance and that the error-proneness they display when copying undamaged templates may simply be a byproduct of this adaptation. Regardless of selection criteria/evolutionary history, at present both of these properties coexist in these enzymes, and both properties confer a fitness advantage. The repair polymerase, Pol X, encoded by the African swine fever virus (ASFV) is one of the most error-prone polymerases known, leading us to previously hypothesize that it may work in tandem with the exceptionally error-tolerant ASFV DNA ligase to effect viral mutagenesis. Here, for the first time, we test whether the error-proneness of Pol X is coupled with a capacity for lesion tolerance by examining its ability to utilize the types of damaged DNA and dNTP substrates that are expected to be relevant to ASFV. We (i) test Pol X's ability to both incorporate opposite to and extend from ubiquitous oxidative purine (7,8-dihydro-8-oxoguanine), oxidative pyrimidine (5,6-dihydroxy-5,6-dihydrothymine), and noncoding (AP site) lesions, in addition to 5,6-dihydrothymine, (ii) determine the catalytic efficiency and dNTP specificity of Pol X when catalyzing incorporation opposite to, and when extending from, 7,8-dihydro-8-oxoguanine in a template/primer context, and (iii) quantitate Pol X-catalyzed incorporation of the damaged nucleotide 8-oxo-dGTP opposite to undamaged templates in the context of both template/primer and a single-nucleotide gap. Our findings are discussed in light of ASFV biology and the mutagenic DNA repair hypothesis described above.

Increased flexibility enhances misincorporation: temperature effects on nucleotide incorporation opposite a bulky carcinogen-DNA adduct by a Y-family DNA polymerase

The Y-family DNA polymerase Dpo4, from the thermophilic crenarchaeon Sulfolobus solfataricus P2, offers a valuable opportunity to investigate the effect of conformational flexibility on the bypass of bulky lesions because of its ability to function efficiently at a wide range of temperatures. Combined molecular modeling and experimental kinetic studies have been carried out for 10S-(+)-trans-anti-[BP]-N2-dG ((+)-ta-[BP]G), a lesion derived from the covalent reaction of a benzo[a]pyrene metabolite with guanine in DNA, at 55 degrees C and results compared with an earlier study at 37 degrees C (Perlow-Poehnelt, R. A., Likhterov, I., Scicchitano, D. A., Geacintov, N. E., and Broyde, S. (2004) J. Biol. Chem. 279, 36951-36961). The experimental results show that there is more overall nucleotide insertion opposite (+)-ta-[BP]G due to particularly enhanced mismatch incorporation at 55 degrees C compared with 37 degrees C. The molecular dynamics simulations suggest that mismatched nucleotide insertion opposite (+)-ta-[BP]G is increased at 55 degrees C compared with 37 degrees C because the higher temperature shifts the preference of the damaged base from the anti to the syn conformation, with the carcinogen on the more open major groove side. The mismatched dNTP structures are less distorted when the damaged base is syn than when it is anti, at the higher temperature. However, with the normal partner dCTP, the anti conformation with close to Watson-Crick alignment remains more favorable. The molecular dynamics simulations are consistent with the kcat values for nucleotide incorporation opposite the lesion studied, providing structural interpretation of the experimental observations. The observed temperature effect suggests that conformational flexibility plays a role in nucleotide incorporation and bypass fidelity opposite (+)-ta-[BP]G by Dpo4.

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, N²-dG adducts, which can reside on the minor groove side, appear to cause less hindrance when in this position.

Free-energy perturbation methods to study structure and energetics of DNA adducts: results for the major N2-dG adduct of benzo[a]pyrene in two conformations and different sequence contexts

The potent mutagen/carcinogen benzo[a]pyrene (B[a]P) is activated to (+)-anti-B[a]PDE, which induces a variety of mutations (e.g., G --> T, G --> A, etc.) via its major adduct [+ta]-B[a]P-N2-dG. One hypothesis is that adducts (such as [+ta]-B[a]P-N2-dG) induce different mutations via different conformations, probably when replicated by different lesion-bypass DNA polymerases (DNAPs). We showed that Escherichia coli DNAP V was responsible for G --> T mutations with [+ta]-B[a]P-N2-dG in a 5'-TGT sequence (Yin et al., (2004) DNA Repair 3, 323), so we wish to study conformations of this adduct/sequence context by molecular modeling. The development of a CHARMM-based molecular dynamics (MD) simulations protocol with free-energy calculations in the presence of solvent and counterions is described. A representative base-pairing and base-displaced conformation of [+ta]-B[a]P-N2-dG in the 5'-TGT sequence are used: (1) BPmi5, which has the B[a]P moiety in the minor groove pointing toward the base on the 5'-side of the adduct, and (2) Gma5, which has the B[a]P moiety stacked with the surrounding base pairs and the dG moiety displaced into the major groove. The MD output structures are reasonable when compared to known NMR structures. Changes in DNA sequence context dramatically affect the biological consequences (e.g., mutagenesis) of [+ta]-B[a]P-N2-dG. Consequently, we also developed a MD-based free-energy perturbation (FEP) protocol to study DNA sequence changes. FEP involves the gradual "fading-out" of atoms in a starting structure (A) and "fading-in" of atoms in a final structure (B), which allows a realistic assessment of the energetic and structural changes when two structures A and B are closely related. Two DNA sequence changes are described: (1) 5'-TGT --> 5'-TGG, which involves two steps [T:A --> T:C --> G:C], and (2) 5'-TGT --> 5'-TGC, which involves three steps [T:A --> T:2AP --> C:2AP --> C:G], where 2AP (2-aminopurine) is included, because T:2AP and C:2AP retain more-or-less normal pairing orientations between complementary bases. FEP is also used to evaluate the impact that a 5'-TGT to 5'-UGT sequence change might have on mutagenesis with [+ta]-B[a]P-N2-dG. In summary, we developed (1) a CHARMM-based molecular dynamics (MD) simulations protocol with free-energy calculations in the presence of solvent and counterions to study B[a]P-N2-dG adducts in DNA duplexes, and (2) a MD-based free-energy perturbation (FEP) protocol to study DNA sequence context changes around B[a]P-N2-dG adducts.

Mismatch-induced conformational distortions in polymerase beta support an induced-fit mechanism for fidelity

Molecular dynamics simulations of DNA polymerase (pol) beta complexed with different incorrect incoming nucleotides (G x G, G x T, and T x T template base x incoming nucleotide combinations) at the template-primer terminus are analyzed to delineate structure-function relationships for aberrant base pairs in a polymerase active site. Comparisons, made to pol beta structure and motions in the presence of a correct base pair, are designed to gain atomically detailed insights into the process of nucleotide selection and discrimination. In the presence of an incorrect incoming nucleotide, alpha-helix N of the thumb subdomain believed to be required for pol beta's catalytic cycling moves toward the open conformation rather than the closed conformation as observed for the correct base pair (G x C) before the chemical reaction. Correspondingly, active-site residues in the microenvironment of the incoming base are in intermediate conformations for non-Watson-Crick pairs. The incorrect incoming nucleotide and the corresponding template residue assume distorted conformations and do not form Watson-Crick bonds. Furthermore, the coordination number and the arrangement of ligands observed around the catalytic and nucleotide binding magnesium ions are mismatch specific. Significantly, the crucial nucleotidyl transferase reaction distance (P(alpha)-O3') for the mismatches between the incoming nucleotide and the primer terminus is not ideally compatible with the chemical reaction of primer extension that follows these conformational changes. Moreover, the extent of active-site distortion can be related to experimentally determined rates of nucleotide misincorporation and to the overall energy barrier associated with polymerase activity. Together, our studies provide structure-function insights into the DNA polymerase-induced constraints (i.e., alpha-helix N conformation, DNA base pair bonding, conformation of protein residues in the vicinity of dNTP, and magnesium ions coordination) during nucleotide discrimination and pol beta-nucleotide interactions specific to each mispair and how they may regulate fidelity. They also lend further support to our recent hypothesis that additional conformational energy barriers are involved following nucleotide binding but prior to the chemical reaction.

Subtle but variable conformational rearrangements in the replication cycle of Sulfolobus solfataricus P2 DNA polymerase IV (Dpo4) may accommodate lesion bypass

The possible conformational changes of DNA polymerase IV (Dpo4) before and after the nucleotidyl-transfer reaction are investigated at the atomic level by dynamics simulations to gain insight into the mechanism of low-fidelity polymerases and identify slow and possibly critical steps. The absence of significant conformational changes in Dpo4 before chemistry when the incoming nucleotide is removed supports the notion that the "induced-fit" mechanism employed to interpret fidelity in some replicative and repair DNA polymerases does not exist in Dpo4. However, significant correlated movements in the little finger and finger domains, as well as DNA sliding and subtle catalytic-residue rearrangements, occur after the chemical reaction when both active-site metal ions are released. Subsequently, Dpo4's little finger grips the DNA through two arginine residues and pushes it forward. These metal ion correlated movements may define subtle, and possibly characteristic, conformational adjustments that operate in some Y-family polymerase members in lieu of the prominent subdomain motions required for catalytic cycling in other DNA polymerases like polymerase beta. Such subtle changes do not easily provide a tight fit for correct incoming substrates as in higher-fidelity polymerases, but introduce in low-fidelity polymerases different fidelity checks as well as the variable conformational-mobility potential required to bypass different lesions.

Accommodation of a 1S-(-)-benzo[c]phenanthrenyl-N6-dA adduct in the Y-family Dpo4 DNA polymerase active site: structural insights through molecular dynamics simulations

Molecular modeling and molecular dynamics simulations have been performed to elucidate feasible structures in the Y-family Dpo4 DNA polymerase for the 1S-(-)-trans-anti-B[c]Ph-N6-dA adduct, derived from the fjord region polycyclic aromatic hydrocarbon (PAH) benzo[c]phenanthrene. Three types of models were delineated as follows: an intercalation model, a model with the aromatic ring system in the polymerase major groove open pocket, and a -1 deletion major groove model. All four 2'-deoxyribonucleoside 5'-triphosphates (dNTPs) were considered in the first two cases, and a normal Watson-Crick partner positioned to have skipped the modified template was employed as the incoming dNTP in the -1 deletion case. The trajectories derived from the dynamics simulations were analyzed in detail to evaluate the extents of distortion for each system. Overall, our results suggest that the major groove model is the least distorted, followed by the -1 deletion model, while the intercalation model is perturbed the most. The syn-dGTP and syn-dATP mismatches opposite the lesion are well-accommodated in the major groove model, as is the normal Watson-Crick partner dTTP. The intercalation model appears most likely to impede the polymerase. More broadly, these models look reasonable for other PAH metabolite-derived adducts to adenine with similar 1S stereochemistry. Furthermore, these models suggest how error-prone translesion synthesis by Y-family polymerases might produce mutations that may play a role in the initiation of cancer.

Structure of a high fidelity DNA polymerase bound to a benzo[a]pyrene adduct that blocks replication

Of the carcinogens to which humans are most frequently exposed, the polycyclic aromatic hydrocarbon benzo[a]pyrene (BP) is one of the most ubiquitous. BP is a byproduct of grilled foods and tobacco and fuel combustion and has long been linked to various human cancers, particularly lung and skin. BP is metabolized to diol epoxides that covalently modify DNA bases to form bulky adducts that block DNA synthesis by replicative or high fidelity DNA polymerases. Here we present the structure of a high fidelity polymerase from a thermostable strain of Bacillus stearothermophilus (Bacillus fragment) bound to the most common BP-derived N2-guanine adduct base-paired with cytosine. The BP adduct adopts a conformation that places the polycyclic BP moiety in the nascent DNA minor groove and is the first structure of a minor groove adduct bound to a polymerase. Orientation of the BP moiety into the nascent DNA minor groove results in extensive disruption to the interactions between the adducted DNA duplex and the polymerase. The disruptions revealed by the structure of Bacillus fragment bound to a BP adduct provide a molecular basis for rationalizing the potent blocking effect on replication exerted by BP adducts.

Dynamic conformational heterogeneities of carcinogen-DNA adducts and their mutagenic relevance

Arylamines and polycyclic aromatic hydrocarbons (PAHs), which are known as "bulky" carcinogens, have been studied extensively and upon activation in vivo, react with cellular DNA to form DNA-adducts. The available structure data accumulated thus far has revealed that conformational heterogeneity is a common theme among duplex DNA modified with these carcinogens. Several conformationally diverse structures have been elucidated and found to be in equilibrium in certain cases. The dynamics of the heterogeneity appear to be modulated by the nature of the adduct structure and the base sequences neighboring the lesion site. These can be termed as "adduct- and sequence-induced conformational heterogeneities," respectively. Due to the small energy differences, the population levels of these conformers could readily be altered within the active sites of repair or replicate enzymes. Thus, the complex role of "enzyme-induced conformational heterogeneity" must also be taken into consideration for the establishment of a functional structure-mutation relationship. Ultimately, a major challenge in mutation structural biology is to carry out adduct- and site-specific experiments in a conformationally specific manner within biologically relevant environments. Results from such experiments should provide an accurate account of how a single chemically homogenous adduct gives rise to complex multiple mutations, the earliest step in the induction of cancer.