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

p53 triggers mitochondrial apoptosis following DNA damage-dependent replication stress by the hepatotoxin methyleugenol

Liver cancer is one of the most frequent tumor entities worldwide, which is causally linked to viral infection, fatty liver disease, life-style factors and food-borne carcinogens, particularly aflatoxins. Moreover, genotoxic plant toxins including phenylpropenes are suspected human liver carcinogens. The phenylpropene methyleugenol (ME) is a constituent of essential oils in many plants and occurs in herbal medicines, food, and cosmetics. Following its uptake, ME undergoes Cytochrome P450 (CYP) and sulfotransferase 1A1 (SULT1A1)-dependent metabolic activation, giving rise to DNA damage. However, little is known about the cellular response to the induced DNA adducts. Here, we made use of different SULT1A1-proficient cell models including primary hepatocytes that were treated with 1'-hydroxymethyleugenol (OH-ME) as main phase I metabolite. Firstly, mass spectrometry showed a concentration-dependent formation of N²-MIE-dG as major DNA adduct, strongly correlating with SULT1A1 expression as attested in cells with and without human SULT1A1. ME-derived DNA damage activated mainly the ATR-mediated DNA damage response as shown by phosphorylation of CHK1 and histone 2AX, followed by p53 accumulation and CHK2 phosphorylation. Consistent with these findings, the DNA adducts decreased replication speed and caused replication fork stalling. OH-ME treatment reduced viability particularly in cell lines with wild-type p53 and triggered apoptotic cell death, which was rescued by pan-caspase-inhibition. Further experiments demonstrated mitochondrial apoptosis as major cell death pathway. ME-derived DNA damage caused upregulation of the p53-responsive genes NOXA and PUMA, Bax activation, and cytochrome c release followed by caspase-9 and caspase-3 cleavage. We finally demonstrated the crucial role of p53 for OH-ME triggered cell death as evidenced by reduced pro-apoptotic gene expression, strongly attenuated Bax activation and cell death inhibition upon genetic knockdown or pharmacological inhibition of p53. Taken together, our study demonstrates for the first time that ME-derived DNA damage causes replication stress and triggers mitochondrial apoptosis via the p53-Bax pathway.

A Unique B-Family DNA Polymerase Facilitating Error-Prone DNA Damage Tolerance in Crenarchaeota

Sulfolobus islandicus codes for four DNA polymerases: three are of the B-family (Dpo1, Dpo2, and Dpo3), and one is of the Y-family (Dpo4). Western analysis revealed that among the four polymerases, only Dpo2 exhibited DNA damage-inducible expression. To investigate how these DNA polymerases could contribute to DNA damage tolerance in S. islandicus, we conducted genetic analysis of their encoding genes in this archaeon. Plasmid-borne gene expression revealed that Dpo2 increases cell survival upon DNA damage at the expense of mutagenesis. Gene deletion studies showed although dpo1 is essential, the remaining three genes are dispensable. Furthermore, although Dpo4 functions in housekeeping translesion DNA synthesis (TLS), Dpo2, a B-family DNA polymerase once predicted to be inactive, functions as a damage-inducible TLS enzyme solely responsible for targeted mutagenesis, facilitating GC to AT/TA conversions in the process. Together, our data indicate that Dpo2 is the main DNA polymerase responsible for DNA damage tolerance and is the primary source of targeted mutagenesis. Given that crenarchaea encoding a Dpo2 also have a low-GC composition genome, the Dpo2-dependent DNA repair pathway may be conserved in this archaeal lineage.

Molecular Insights into the Translesion Synthesis of Benzyl-Guanine from Molecular Dynamics Simulations: Structural Evidence of Mutagenic and Nonmutagenic Replication

DNA can be damaged by many compounds in our environment, and the resulting damaged DNA is commonly replicated by translesion synthesis (TLS) polymerases. Because the mechanism and efficiency of TLS are affected by the type of DNA damage, obtaining information for a variety of DNA adducts is critical. However, there is no structural information for the insertion of a dNTP opposite an O6-dG adduct, which is a particularly harmful class of DNA lesions. We used molecular dynamics (MD) simulations to investigate structural and energetic parameters that dictate preferred dNTP insertion opposite O6-benzyl-guanine (Bz-dG) by DNA polymerase IV, a prototypical TLS polymerase. Specifically, MD simulations were completed on all possible ternary insertion complexes and ternary -1 base deletion complexes with different Bz-dG conformations. Our data suggests that the purines are unlikely to be inserted opposite anti- or syn-Bz-dG, and dTTP is unlikely to be inserted opposite syn-Bz-dG, because of changes in the active site conformation, including critical hydrogen-bonding interactions and/or reaction-ready parameters compared to natural dG replication. In contrast, a preserved active site conformation suggests that dCTP can be inserted opposite either anti- or syn-Bz-dG and dTTP can be inserted opposite anti-Bz-dG. This is the first structural explanation for the experimentally observed preferential insertion of dCTP and misincorporation of dTTP opposite Bz-dG. Furthermore, we provide atomic level insight into why Bz-dG replication does not lead to deletion mutations, which is in contrast with the replication outcomes of other adducts. These findings provide a basis for understanding the replication of related O6-dG adducts.

DNA damage response curtails detrimental replication stress and chromosomal instability induced by the dietary carcinogen PhIP

PhIP is an abundant heterocyclic aromatic amine (HCA) and important dietary carcinogen. Following metabolic activation, PhIP causes bulky DNA lesions at the C8-position of guanine. Although C8-PhIP-dG adducts are mutagenic, their interference with the DNA replication machinery and the elicited DNA damage response (DDR) have not yet been studied. Here, we analyzed PhIP-triggered replicative stress and elucidated the role of the apical DDR kinases ATR, ATM and DNA-PKcs in the cellular defense response. First, we demonstrate that PhIP induced C8-PhIP-dG adducts and DNA strand breaks. This stimulated ATR-CHK1 signaling, phosphorylation of histone 2AX and the formation of RPA foci. In proliferating cells, PhIP treatment increased the frequency of stalled replication forks and reduced fork speed. Inhibition of ATR in the presence of PhIP-induced DNA damage strongly promoted the formation of DNA double-strand breaks, activation of the ATM-CHK2 pathway and hyperphosphorylation of RPA. The abrogation of ATR signaling potentiated the cell death response and enhanced chromosomal aberrations after PhIP treatment, while ATM and DNA-PK inhibition had only marginal effects. These results strongly support the notion that ATR plays a key role in the defense against cancer formation induced by PhIP and related HCAs.

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.

Tolerance of base pair size and shape in postlesion DNA synthesis

The influence of base pair size and shape on the fidelity of DNA polymerase-mediated extension past lesion-containing mispairs was examined. Primer extension analysis was performed with synthetic nucleosides paired opposite the pro-mutagenic DNA lesion O(6)-benzylguanine (O(6)-BnG). These data indicate that the error-prone DNA polymerase IV (Dpo4) inefficiently extended past the larger Peri:O(6)-BnG base pair, and in contrast, error-free extension was observed for the smaller BIM:O(6)-BnG base pair. Steady-state kinetic analysis revealed that Dpo4 catalytic efficiency was strongly influenced by the primer:template base pair. Compared to the C:G pair, a 1.9- and 79,000-fold reduction in Dpo4 efficiency was observed for terminal C:O(6)-BnG and BIM:G base pairs respectively. These results demonstrate the impact of geometrical size and shape on polymerase-mediated mispair extension.

Preferred WMSA catalytic mechanism of the nucleotidyl transfer reaction in human DNA polymerase κ elucidates error-free bypass of a bulky DNA lesion

Human DNA Pol κ is a polymerase enzyme, specialized for near error-free bypass of certain bulky chemical lesions to DNA that are derived from environmental carcinogens present in tobacco smoke, automobile exhaust and cooked food. By employing ab initio QM/MM-MD (Quantum Mechanics/Molecular Mechanics-Molecular Dynamics) simulations with umbrella sampling, we have determined the entire free energy profile of the nucleotidyl transfer reaction catalyzed by Pol κ and provided detailed mechanistic insights. Our results show that a variant of the Water Mediated and Substrate Assisted (WMSA) mechanism that we previously deduced for Dpo4 and T7 DNA polymerases is preferred for Pol κ as well, suggesting its broad applicability. The hydrogen on the 3'-OH primer terminus is transferred through crystal and solvent waters to the γ-phosphate of the dNTP, followed by the associative nucleotidyl transfer reaction; this is facilitated by a proton transfer from the γ-phosphate to the α,β-bridging oxygen as pyrophosphate leaves, to neutralize the evolving negative charge. MD simulations show that the near error-free incorporation of dCTP opposite the major benzo[a]pyrene-derived dG lesion is compatible with the WMSA mechanism, allowing for an essentially undisturbed pentacovalent phosphorane transition state, and explaining the bypass of this lesion with little mutation by Pol κ.

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.

Effect of oxidatively damaged DNA on the active site preorganization during nucleotide incorporation in a high fidelity polymerase from Bacillus stearothermophilus

We study the effect of the oxidative lesion 8-oxoguanine (8oxoG) on the preorganization of the active site for DNA replication in the closed (active) state of the Bacillus fragment (BF), a Klenow analog from Bacillus stearothermophilus. Our molecular dynamics and free energy simulations of explicitly solvated model ternary complexes of BF bound to correct dCTP/incorrect dATP opposite guanine (G) and 8oxoG bases in DNA suggest that the lesion introduces structural and energetic changes at the catalytic site to favor dATP insertion. Despite the formation of a stable Watson-Crick pairing in the 8oxoG:dCTP system, the catalytic geometry is severely distorted to possibly slow down catalysis. Indeed, our calculated free energy landscapes associated with active site preorganization suggest additional barriers to assemble an efficient catalytic site, which need to be overcome during dCTP incorporation opposite 8oxoG relative to that opposite undamaged G. In contrast, the catalytic geometry for the Hoogsteen pairing in the 8oxoG:dATP system is highly organized and poised for efficient nucleotide incorporation via the "two-metal-ion" catalyzed phosphoryl transfer mechanism. However, the free energy calculations suggest that the catalytic geometry during dATP incorporation opposite 8oxoG is considerably less plastic than that during dCTP incorporation opposite G despite a very similar, well organized catalytic site for both systems. A correlation analysis of the dynamics trajectories suggests the presence of significant coupling between motions of the polymerase fingers and the primary distance for nucleophilic attack (i.e., between the terminal primer O3' and the dNTP P(alpha.) atoms) during correct dCTP incorporation opposite undamaged G. This coupling is shown to be disrupted during nucleotide incorporation by the polymerase with oxidatively damaged DNA/dNTP substrates. We also suggest that the lesion affects DNA interactions with key polymerase residues, thereby affecting the enzymes ability to discriminate against non-complementary DNA/dNTP substrates. Taken together, our results provide a unified structural, energetic, and dynamic platform to rationalize experimentally observed relative nucleotide incorporation rates for correct dCTP/incorrect dATP insertion opposite an undamaged/oxidatively damaged template G by BF.

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.

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.

Examination of the long-range effects of aminofluorene-induced conformational heterogeneity and its relevance to the mechanism of translesional DNA synthesis

Adduct-induced conformational heterogeneity complicates the understanding of how DNA adducts exert mutation. A case in point is the N-deacetylated AF lesion [N-(2'-deoxyguanosin-8-yl)-2-aminofluorene], the major adduct derived from the strong liver carcinogen N-acetyl-2-aminofluorene. Three conformational families have been previously characterized and are dependent on the positioning of the aminofluorene rings: B is in the "B-DNA" major groove, S is "stacked" into the helix with base-displacement, and W is "wedged" into the minor groove. Here, we conducted (19)F NMR, CD, T(m), and modeling experiments at various primer positions with respect to a template modified by a fluorine tagged AF-adduct (FAF). In the first set, the FAF-G was paired with C and in the second set it was paired with A. The FAF-G:C oligonucleotides were found to preferentially adopt the B or S-conformers while the FAF-G:A mismatch ones preferred the B and W-conformers. The conformational preferences of both series were dependent on temperature and complementary strand length; the largest differences in conformation were displayed at lower temperatures. The CD and T(m) results are in general agreement with the NMR data. Molecular modeling indicated that the aminofluorene moiety in the minor groove of the W-conformer would impose a steric clash with the tight-packing amino acid residues on the DNA binding area of the Bacillus fragment (BF), a replicative DNA polymerase. In the case of the B-type conformer, the carcinogenic moiety resides in the solvent-exposed major groove throughout the replication/translocation process. The present dynamic NMR results, combined with previous primer extension kinetic data by Miller & Grollman, support a model in which adduct-induced conformational heterogeneities at positions remote from the replication fork affect polymerase function through a long-range DNA-protein interaction.