Nucleotide excision repair efficiencies of bulky carcinogen-DNA adducts are governed by a balance between stabilizing and destabilizing interactions

Diet derived polycyclic aromatic hydrocarbons and its pathogenic roles in colorectal carcinogenesis

Polycyclic aromatic hydrocarbon (PAHs) are molecules that contaminate meat products during the high-temperature cooking of meat. This study reviewed the pathogenic roles of meat derived polycyclic aromatic hydrocarbons in the carcinogenesis of colorectal cancer (CRC). Ingested PAHs undergo xenobiotic metabolism resulting in the activation of genotoxic metabolites that can induce DNA damage in the colorectum. Genetic polymorphisms in PAH xenobiotic enzymes are linked to the risk of CRC and suggest a role for PAH-meat ingestion in carcinogenesis of colorectal malignancies. Furthermore, PAH specific DNA adducts have been identified in colorectal cancer tissue and linked to high meat intake. DNA adduct resolution is mediated by the nucleotide excision repair, and polymorphisms within genes of this repair pathway and high meat intake are associated with increased CRC risk. In the literature, there is evidence from metabolic enzyme gene variants, DNA repair genes, PAH metabolites, and epidemiological studies suggesting PAH involvement in CRC.

Impact of DNA sequences on DNA 'opening' by the Rad4/XPC nucleotide excision repair complex

Rad4/XPC recognizes diverse DNA lesions to initiate nucleotide excision repair (NER). However, NER propensities among lesions vary widely and repair-resistant lesions are persistent and thus highly mutagenic. Rad4 recognizes repair-proficient lesions by unwinding ('opening') the damaged DNA site. Such 'opening' is also observed on a normal DNA sequence containing consecutive C/G's (CCC/GGG) when tethered to Rad4 to prevent protein diffusion. However, it was unknown if such tethering-facilitated DNA 'opening' could occur on any DNA or if certain structures/sequences would resist being 'opened'. Here, we report that DNA containing alternating C/G's (CGC/GCG) failed to be opened even when tethered; instead, Rad4 bound in a 180°-reversed manner, capping the DNA end. Fluorescence lifetime studies of DNA conformations in solution showed that CCC/GGG exhibits local pre-melting that is absent in CGC/GCG. In MD simulations, CGC/GCG failed to engage Rad4 to promote 'opening' contrary to CCC/GGG. Altogether, our study illustrates how local sequences can impact DNA recognition by Rad4/XPC and how certain DNA sites resist being 'opened' even with Rad4 held at that site indefinitely. The contrast between CCC/GGG and CGC/GCG sequences in Rad4-DNA recognition may help decipher a lesion's mutagenicity in various genomic sequence contexts to explain lesion-determined mutational hot and cold spots.

Carcinogen-induced DNA structural distortion differences in the RAS gene isoforms; the importance of local sequence

BACKGROUND

Local sequence context is known to have an impact on the mutational pattern seen in cancer. The RAS genes and a smoking carcinogen, Benzo[a]pyrene diol epoxide (BPDE), have been utilised to explore these context effects. BPDE is known to form an adduct at the guanines in a number of RAS gene sites, KRAS codons 12, 13 and 14, NRAS codon 12, and HRAS codons 12 and 14.

RESULTS

Molecular modelling techniques, along with multivariate analysis, have been utilised to determine the sequence influenced differences between BPDE-adducted RAS gene sequences as well as the local distortion caused by the adducts.

CONCLUSIONS

We conclude that G:C > T:A mutations at KRAS codon 12 in the tumours of lung cancer patients (who smoke), proposed to be predominantly caused by BPDE, are due to the effect of the interaction methyl group at the C5 position of the thymine base in the KRAS sequence with the BPDE carcinogen investigated causing increased distortion. We further suggest methylated cytosine would have a similar effect, showing the importance of methylation in cancer development.

Regularization of least squares problems in CHARMM parameter optimization by truncated singular value decompositions

We examine the use of the truncated singular value decomposition and Tikhonov regularization in standard form to address ill-posed least squares problems Ax = b that frequently arise in molecular mechanics force field parameter optimization. We illustrate these approaches by applying them to dihedral parameter optimization of genotoxic polycyclic aromatic hydrocarbon-DNA adducts that are of interest in the study of chemical carcinogenesis. Utilizing the discrete Picard condition and/or a well-defined gap in the singular value spectrum when A has a well-determined numerical rank, we are able to systematically determine truncation and in turn regularization parameters that are correspondingly used to produce truncated and regularized solutions to the ill-posed least squares problem at hand. These solutions in turn result in optimized force field dihedral terms that accurately parameterize the torsional energy landscape. As the solutions produced by this approach are unique, it has the advantage of avoiding the multiple iterations and guess and check work often required to optimize molecular mechanics force field parameters.

Inflammation response, oxidative stress and DNA damage caused by urban air pollution exposure increase in the lack of DNA repair XPC protein

Air pollution represents a considerable threat to health worldwide. The São Paulo Metropolitan area, in Brazil, has a unique composition of atmospheric pollutants with a population of nearly 20 million people and 9 million passenger cars. It is long known that exposure to particulate matter less than 2.5 µm (PM_2.5) can cause various health effects such as DNA damage. One of the most versatile defense mechanisms against the accumulation of DNA damage is the nucleotide excision repair (NER), which includes XPC protein. However, the mechanisms by which NER protects against adverse health effects related to air pollution are largely unknown. We hypothesized that reduction of XPC activity may contribute to inflammation response, oxidative stress and DNA damage after PM_2.5 exposure. To address these important questions, XPC knockout and wild type mice were exposed to PM_2.5 using the Harvard Ambient Particle concentrator. Results from one-single exposure have shown a significant increase in the levels of anti-ICAM, IL-1β, and TNF-α in the polluted group when compared to the filtered air group. Continued chronic PM_2.5 exposure increased levels of carbonylated proteins, especially in the lung of XPC mice, probably as a consequence of oxidative stress. As a response to DNA damage, XPC mice lungs exhibit increased γ-H2AX, followed by severe atypical hyperplasia. Emissions from vehicles are composed of hazardous substances, with polycyclic aromatic hydrocarbons (PAHs) and metals being most frequently cited as the major contributors to negative health impacts. This analysis showed that benzo[b]fluoranthene, 2-nitrofluorene and 9,10-anthraquinone were the most abundant PAHs and derivatives. Taken together, these findings demonstrate the participation of XPC protein, and NER pathway, in the protection of mice against the carcinogenic potential of air pollution. This implicates that DNA is damaged directly (forming adducts) or indirectly (Reactive Oxygen Species) by the various compounds detected in urban PM_2.5.

Cellular levels and molecular dynamics simulations of estragole DNA adducts point at inefficient repair resulting from limited distortion of the double-stranded DNA helix

Estragole, naturally occurring in a variety of herbs and spices, can form DNA adducts after bioactivation. Estragole DNA adduct formation and repair was studied in in vitro liver cell models, and a molecular dynamics simulation was used to investigate the conformation dependent (in)efficiency of N²-(trans-isoestragol-3'-yl)-2'-deoxyguanosine (E-3'-N²-dG) DNA adduct repair. HepG2, HepaRG cells, primary rat hepatocytes and CHO cells (including CHO wild-type and three NER-deficient mutants) were exposed to 50 μM estragole or 1'-hydroxyestragole and DNA adduct formation was quantified by LC-MS immediately following exposure and after a period of repair. Results obtained from CHO cell lines indicated that NER plays a role in repair of E-3'-N²-dG adducts, however, with limited efficiency since in the CHO wt cells 80% DNA adducts remained upon 24 h repair. Inefficiency of DNA repair was also found in HepaRG cells and primary rat hepatocytes. Changes in DNA structure resulting from E-3'-N²-dG adduct formation were investigated by molecular dynamics simulations. Results from molecular dynamics simulations revealed that conformational changes in double-stranded DNA by E-3'-N²-dG adduct formation are small, providing a possible explanation for the restrained repair, which may require larger distortions in the DNA structure. NER-mediated enzymatic repair of E-3'-N²-dG DNA adducts upon exposure to estragole will be limited, providing opportunities for accumulation of damage upon repeated daily exposure. The inability of this enzymatic repair is likely due to a limited distortion of the DNA double-stranded helix resulting in inefficient activation of nucleotide excision repair.

Enhanced spontaneous DNA twisting/bending fluctuations unveiled by fluorescence lifetime distributions promote mismatch recognition by the Rad4 nucleotide excision repair complex

Rad4/XPC recognizes diverse DNA lesions including ultraviolet-photolesions and carcinogen-DNA adducts, initiating nucleotide excision repair. Studies have suggested that Rad4/XPC senses lesion-induced helix-destabilization to flip out nucleotides from damaged DNA sites. However, characterizing how DNA deformability and/or distortions impact recognition has been challenging. Here, using fluorescence lifetime measurements empowered by a maximum entropy algorithm, we mapped the conformational heterogeneities of artificially destabilized mismatched DNA substrates of varying Rad4-binding specificities. The conformational distributions, as probed by FRET between a cytosine-analog pair exquisitely sensitive to DNA twisting/bending, reveal a direct connection between intrinsic DNA deformability and Rad4 recognition. High-specificity CCC/CCC mismatch, free in solution, sampled a strikingly broad range of conformations from B-DNA-like to highly distorted conformations that resembled those observed with Rad4 bound; the extent of these distortions increased with bound Rad4 and with temperature. Conversely, the non-specific TAT/TAT mismatch had a homogeneous, B-DNA-like conformation. Molecular dynamics simulations also revealed a wide distribution of conformations for CCC/CCC, complementing experimental findings. We propose that intrinsic deformability promotes Rad4 damage recognition, perhaps by stalling a diffusing protein and/or facilitating ‘conformational capture’ of pre-distorted damaged sites. Surprisingly, even mismatched DNA specifically bound to Rad4 remains highly dynamic, a feature that may reflect the versatility of Rad4/XPC to recognize many structurally dissimilar lesions.

Lesion Sensing during Initial Binding by Yeast XPC/Rad4: Toward Predicting Resistance to Nucleotide Excision Repair

Nucleotide excision repair (NER) excises a variety of environmentally derived DNA lesions. However, NER efficiencies for structurally different DNA lesions can vary by orders of magnitude; yet the origin of this variance is poorly understood. Our goal is to develop computational strategies that predict and identify the most hazardous, repair-resistant lesions from the plethora of such adducts. In the present work, we are focusing on lesion recognition by the xeroderma pigmentosum C protein complex (XPC), the first and required step for the subsequent assembly of factors needed to produce successful NER. We have performed molecular dynamics simulations to characterize the initial binding of Rad4, the yeast orthologue of human XPC, to a library of 10 different lesion-containing DNA duplexes derived from environmental carcinogens. These vary in lesion chemical structures and conformations in duplex DNA and exhibit a wide range of relative NER efficiencies from repair resistant to highly susceptible. We have determined a promising set of structural descriptors that characterize initial binding of Rad4 to lesions that are resistant to NER. Key initial binding requirements for successful recognition are absent in the repair-resistant cases: There is little or no duplex unwinding, very limited interaction between the β-hairpin domain 2 of Rad4 and the minor groove of the lesion-containing duplex, and no conformational capture of a base on the lesion partner strand. By contrast, these key binding features are present to different degrees in NER susceptible lesions and correlate to their relative NER efficiencies. Furthermore, we have gained molecular understanding of Rad4 initial binding as determined by the lesion structures in duplex DNA and how the initial binding relates to the repair efficiencies. The development of a computational strategy for identifying NER-resistant lesions is grounded in this molecular understanding of the lesion recognition mechanism.

Impact of DNA lesion repair, replication and formation on the mutational spectra of environmental carcinogens: Aflatoxin B1 as a case study

In a multicellular organism, somatic mutations represent a permanent record of the past chemical and biochemical perturbations experienced by a cell in its local microenvironment. Akin to a perpetual recording device, with every replication, genomic DNA accumulates mutations in patterns that reflect: i) the sequence context-dependent formation of DNA damage, due to environmental or endogenous reactive species, including spontaneous processes; ii) the activity of DNA repair pathways, which, depending on the type of lesion, can erase, ignore or exacerbate the mutagenic consequences of that DNA damage; and iii) the choice of replication machinery that synthesizes the nascent genomic copy. These three factors result in a richly contoured sequence context-dependent mutational spectrum that, from appearances, is distinct for most individual forms of DNA damage. Such a mutagenic legacy, if appropriately decoded, can reveal the local history of genome-altering events such as chemical or pathogen exposures, metabolic stress, and inflammation, which in turn can provide an indication of the underlying causes and mechanisms of genetic disease. Modern tools have positioned us to develop a deep mechanistic understanding of the cellular factors and pathways that modulate a mutational process and, in turn, provide opportunities for better diagnostic and prognostic biomarkers, better exposure risk assessment and even actionable therapeutic targets. The goal of this Perspective is to present a bottom-up, lesion-centric framework of mutagenesis that integrates the contributions of lesion replication, lesion repair and lesion formation to explain the complex mutational spectra that emerge in the genome following exposure to mutagens. The mutational spectra of the well-studied hepatocarcinogen aflatoxin B₁ are showcased here as specific examples, but the implications are meant to be generalizable.

Repair-Resistant DNA Lesions

The eukaryotic global genomic nucleotide excision repair (GG-NER) pathway is the major mechanism that removes most bulky and some nonbulky lesions from cellular DNA. There is growing evidence that certain DNA lesions are repaired slowly or are entirely resistant to repair in cells, tissues, and in cell extract model assay systems. It is well established that the eukaryotic DNA lesion-sensing proteins do not detect the damaged nucleotide, but recognize the distortions/destabilizations in the native DNA structure caused by the damaged nucleotides. In this article, the nature of the structural features of certain bulky DNA lesions that render them resistant to NER, or cause them to be repaired slowly, is compared to that of those that are good-to-excellent NER substrates. Understanding the structural features that distinguish NER-resistant DNA lesions from good NER substrates may be useful for interpreting the biological significance of biomarkers of exposure of human populations to genotoxic environmental chemicals. NER-resistant lesions can survive to replication and cause mutations that can initiate cancer and other diseases. Furthermore, NER diminishes the efficacy of certain chemotherapeutic drugs, and the design of more potent pharmaceuticals that resist repair can be advanced through a better understanding of the structural properties of DNA lesions that engender repair-resistance.

Nucleotide Excision Repair Lesion-Recognition Protein Rad4 Captures a Pre-Flipped Partner Base in a Benzo[a]pyrene-Derived DNA Lesion: How Structure Impacts the Binding Pathway

The xeroderma pigmentosum C protein complex (XPC) recognizes a variety of environmentally induced DNA lesions and is the key in initiating their repair by the nucleotide excision repair (NER) pathway. When bound to a lesion, XPC flips two nucleotide pairs that include the lesion out of the DNA duplex, yielding a productively bound complex that can lead to successful lesion excision. Interestingly, the efficiencies of NER vary greatly among different lesions, influencing their toxicity and mutagenicity in cells. Though differences in XPC binding may influence NER efficiency, it is not understood whether XPC utilizes different mechanisms to achieve productive binding with different lesions. Here, we investigated the well-repaired 10R-(+)-cis-anti-benzo[a]pyrene-N²-dG (cis-B[a]P-dG) DNA adduct in a duplex containing normal partner C opposite the lesion. This adduct is derived from the environmental pro-carcinogen benzo[a]pyrene and is likely to be encountered by NER in the cell. We have extensively investigated its binding to the yeast XPC orthologue, Rad4, using umbrella sampling with restrained molecular dynamics simulations and free energy calculations. The NMR solution structure of this lesion in duplex DNA has shown that the dC complementary to the adducted dG is flipped out of the DNA duplex in the absence of XPC. However, it is not known whether the "pre-flipped" base would play a role in its recognition by XPC. Our results show that Rad4 first captures the displaced dC, which is followed by a tightly coupled lesion-extruding pathway for productive binding. This binding path differs significantly from the one deduced for the small cis-syn cyclobutane pyrimidine dimer lesion opposite mismatched thymines [ Mu , H. , ( 2015 ) Biochemistry , 54 ( 34 ), 5263 - 7 ]. The possibility of multiple paths that lead to productive binding to XPC is consistent with the versatile lesion recognition by XPC that is required for successful NER.

Repair Rate of Clustered Abasic DNA Lesions by Human Endonuclease: Molecular Bases of Sequence Specificity

In the present contribution, the interaction between damaged DNA and repair enzymes is examined by means of molecular dynamics simulations. More specifically, we consider clustered abasic DNA lesions processed by the primary human apurinic/apyrimidinic (AP) endonuclease, APE1. Our results show that, in stark contrast with the corresponding bacterial endonucleases, human APE1 imposes strong geometrical constraints on the DNA duplex. As a consequence, the level of recognition and, hence, the repair rate is higher. Important features that guide the DNA/protein interactions are the presence of an extended positively charged region and of a molecular tweezers that strongly constrains DNA. Our results are on very good agreement with the experimentally determined repair rate of clustered abasic lesions. The lack of repair for one particular arrangement of the two abasic sites is also explained considering the peculiar destabilizing interaction between the recognition region and the second lesion, resulting in a partial opening of the molecular tweezers and, thus, a less stable complex. This contribution cogently establishes the molecular bases for the recognition and repair of clustered DNA lesions by means of human endonucleases.

DNA Sequence Modulates Geometrical Isomerism of the trans-8,9- Dihydro-8-(2,6-diamino-4-oxo-3,4-dihydropyrimid-5-yl-formamido)- 9-hydroxy Aflatoxin B1 Adduct

Aflatoxin B₁ (AFB₁), a mycotoxin produced by Aspergillus flavus, is oxidized by cytochrome P450 enzymes to aflatoxin B₁-8,9-epoxide, which alkylates DNA at N7-dG. Under basic conditions, this N7-dG adduct rearranges to yield the trans-8,9-dihydro-8-(2,6-diamino-4-oxo-3,4-dihydropyrimid-5-yl-formamido)-9-hydroxy aflatoxin B₁ (AFB₁–FAPY) adduct. The AFB₁–FAPY adduct exhibits geometrical isomerism involving the formamide moiety. NMR analyses of duplex oligodeoxynucleotides containing the 5′-XA-3′, 5′-XC-3′, 5′-XT-3′, and 5′-XY-3′ sequences (X = AFB₁–FAPY; Y = 7-deaza-dG) demonstrate that the equilibrium between E and Z isomers is controlled by major groove hydrogen bonding interactions. Structural analysis of the adduct in the 5′-XA-3′ sequence indicates the preference of the E isomer of the formamide group, attributed to formation of a hydrogen bond between the formyl oxygen and the N⁶ exocyclic amino group of the 3′-neighbor adenine. While the 5′-XA-3′ sequence exhibits the E isomer, the 5′-XC-3′ sequence exhibits a 7:3 E:Z ratio at equilibrium at 283 K. The E isomer is favored by a hydrogen bond between the formyl oxygen and the N⁴-dC exocyclic amino group of the 3′-neighbor cytosine. The 5′-XT-3′ and 5′-XY-3′ sequences cannot form such a hydrogen bond between the formyl oxygen and the 3′-neighbor T or Y, respectively, and in these sequence contexts the Z isomer is favored. Additional equilibria between α and β anomers and the potential to exhibit atropisomers about the C5–N⁵ bond do not depend upon sequence. In each of the four DNA sequences, the AFB₁–FAPY adduct maintains the β deoxyribose configuration. Each of these four sequences feature the atropisomer of the AFB₁ moiety that is intercalated above the 5′-face of the damaged guanine. This enforces the R_a axial conformation for the C5–N⁵ bond.

Impact of thymine glycol damage on DNA duplex energetics: Correlations with lesion-induced biochemical and structural consequences

The magnitude and nature of lesion-induced energetic perturbations empirically correlate with mutagenicity/cytotoxicity profiles and can be predictive of lesion outcomes during polymerase-mediated replication in vitro. In this study, we assess the sequence and counterbase-dependent energetic impact of the Thymine glycol (Tg) lesion on a family of deoxyoligonucleotide duplexes. Tg damage arises from thymine and methyl-cytosine exposure to oxidizing agents or radiation-generated free-radicals. The Tg lesion blocks polymerase-mediated DNA replication in vitro and the unrepaired site elicits cytotoxic lethal consequences in vivo. Our combined calorimetric and spectroscopic characterization correlates Tg -induced energetic perturbations with biological and structural properties. Specifically, we incorporate a 5R-Tg isomer centered within the tridecanucleotide sequence 5'-GCGTACXCATGCG-3' (X = Tg or T) which is hybridized with the corresponding complementary sequence 5'-CGCATGNGTACGC-3' (N = A, G, T, C) to generate families of Tg -damaged (Tg ·N) and lesion-free (T·N) duplexes. We demonstrate that the magnitude and nature of the Tg destabilizing impact is dependent on counterbase identity (i.e., A ∼ G < T < C). The observation that a Tg lesion is less destabilizing when positioned opposite purines suggests that favorable counterbase stacking interactions may partially compensate lesion-induced perturbations. Moreover, the destabilizing energies of Tg ·N duplexes parallel their respective lesion-free T·N mismatch counterparts (i.e., G < T < C). Elucidation of Tg-induced destabilization relative to the corresponding undamaged mismatch energetics allows resolution of lesion-specific and sequence-dependent impacts. The Tg-induced energetic perturbations are consistent with its replication blocking properties and may serve as differential recognition elements for discrimination by the cellular repair machinery.

Poor response to platinum-based chemotherapy is associated with KRAS mutation and concomitant low expression of BRAC1 and TYMS in NSCLC

Objective

To evaluate treatment response, survival, and the associations between KRAS mutation status and tumour expression levels of BRCA1, TYMS and SRC retrospectively in a cohort of patients with non-small cell lung cancer (NSCLC), treated exclusively with conjunctive platinum-based doublet chemotherapy.

Methods

KRAS mutation status was determined via amplification refractory mutation and multiple quantitative polymerase chain reaction (PCR) analysis. Tumour expression levels of BRCA1, TYMS and SRC were determined via real time quantitative PCR.

Results

Patients with KRAS mutations (n = 3) had significantly shorter survival duration than patients with wild type KRAS (n = 42). Tumour expression levels of BRCA1 and TYMS, but not SRC, were significantly lower in patients with, than in those without, KRAS mutations. Tumour expression level of BRCA1 was positively correlated with survival duration.

Conclusions

KRAS mutation status and BRCA1 tumour expression are potential biomarkers for tailoring chemotherapy and predicting clinical outcome.

Nucleotide Excision Repair and Transcription-coupled DNA Repair Abrogate the Impact of DNA Damage on Transcription

DNA adducts derived from carcinogenic polycyclic aromatic hydrocarbons like benzo[a]pyrene (B[a]P) and benzo[c]phenanthrene (B[c]Ph) impede replication and transcription, resulting in aberrant cell division and gene expression. Global nucleotide excision repair (NER) and transcription-coupled DNA repair (TCR) are among the DNA repair pathways that evolved to maintain genome integrity by removing DNA damage. The interplay between global NER and TCR in repairing the polycyclic aromatic hydrocarbon-derived DNA adducts (+)-trans-anti-B[a]P-N(6)-dA, which is subject to NER and blocks transcription in vitro, and (+)-trans-anti-B[c]Ph-N(6)-dA, which is a poor substrate for NER but also blocks transcription in vitro, was tested. The results show that both adducts inhibit transcription in human cells that lack both NER and TCR. The (+)-trans-anti-B[a]P-N(6)-dA lesion exhibited no detectable effect on transcription in cells proficient in NER but lacking TCR, indicating that NER can remove the lesion in the absence of TCR, which is consistent with in vitro data. In primary human cells lacking NER, (+)-trans-anti-B[a]P-N(6)-dA exhibited a deleterious effect on transcription that was less severe than in cells lacking both pathways, suggesting that TCR can repair the adduct but not as effectively as global NER. In contrast, (+)-trans-anti-B[c]Ph-N(6)-dA dramatically reduces transcript production in cells proficient in global NER but lacking TCR, indicating that TCR is necessary for the removal of this adduct, which is consistent with in vitro data showing that it is a poor substrate for NER. Hence, both global NER and TCR enhance the recovery of gene expression following DNA damage, and TCR plays an important role in removing DNA damage that is refractory to NER.

Effect of base sequence context on the conformational heterogeneity of aristolactam-I adducted DNA: structural and energetic insights into sequence-dependent repair and mutagenicity

Aristolochic acids (AAs) are nephrotoxic and potentially carcinogenic plant mutagens that form bulky DNA adducts at the exocyclic amino groups of the purines. The present work utilizes classical molecular dynamics simulations and free energy calculations to investigate the role of lesion site sequence context in dictating the conformational outcomes of DNA containing ALI-N6-dA, the most persistent and mutagenic adduct arising from the AAs. Our calculations reveal that the anti base-displaced intercalated conformer is the lowest energy conformer of damaged DNA in all sequence contexts considered (CXC, CXG, GXC and GXG). However, the experimentally-observed greater mutagenicity of the adduct in the CXG sequence context does not correlate with the relative thermodynamic stability of the adduct in different sequences. Instead, AL-N6-dA adducted DNA is least distorted in the CXG sequence context, which points toward a possible differential repair propensity of the lesion in different sequences. Nevertheless, the structural deviations between adducted DNA with different lesion site sequences are small, and therefore other factors (such as interactions between the adducted DNA and lesion-bypass polymerases during replication) are likely more important for dictating the observed sequence-dependent mutagenicity of ALI-N6-dA.

Base damage, local sequence context and TP53 mutation hotspots: a molecular dynamics study of benzo[a]pyrene induced DNA distortion and mutability

The mutational pattern for the TP53 tumour suppressor gene in lung tumours differs to other cancer types by having a higher frequency of G:C>T:A transversions. The aetiology of this differing mutation pattern is still unknown. Benzo[a]pyrene,diol epoxide (BPDE) is a potent cigarette smoke carcinogen that forms guanine adducts at TP53 CpG mutation hotspot sites including codons 157, 158, 245, 248 and 273. We performed molecular modelling of BPDE-adducted TP53 duplex sequences to determine the degree of local distortion caused by adducts which could influence the ability of nucleotide excision repair. We show that BPDE adducted codon 157 has greater structural distortion than other TP53 G:C>T:A hotspot sites and that sequence context more distal to adjacent bases must influence local distortion. Using TP53 trinucleotide mutation signatures for lung cancer in smokers and non-smokers we further show that codons 157 and 273 have the highest mutation probability in smokers. Combining this information with adduct structural data we predict that G:C>T:A mutations at codon 157 in lung tumours of smokers are predominantly caused by BPDE. Our results provide insight into how different DNA sequence contexts show variability in DNA distortion at mutagen adduct sites that could compromise DNA repair at well characterized cancer related mutation hotspots.

Recognition of Damaged DNA for Nucleotide Excision Repair: A Correlated Motion Mechanism with a Mismatched cis-syn Thymine Dimer Lesion

Mammalian global genomic nucleotide excision repair requires lesion recognition by XPC, whose detailed binding mechanism remains to be elucidated. Here we have delineated the dynamic molecular pathway and energetics of lesion-specific and productive binding by the Rad4/yeast XPC lesion recognition factor, as it forms the open complex [Min, J. H., and Pavletich, N. P. (2007) Nature 449, 570–575; Chen, X., et al. (2015) Nat. Commun. 6, 5849] that is required for excision. We investigated extensively a cis-syn cyclobutane pyrimidine dimer in mismatched duplex DNA, using high-level computational approaches. Our results delineate a preferred correlated motion mechanism, which provides for the first time an atomistic description of the sequence of events as Rad4 productively binds to the damaged DNA.

Adenine versus guanine DNA adducts of aristolochic acids: role of the carcinogen-purine linkage in the differential global genomic repair propensity

Computational modeling is employed to provide a plausible structural explanation for the experimentally-observed differential global genome repair (GGR) propensity of the ALII-N(2)-dG and ALII-N(6)-dA DNA adducts of aristolochic acid II. Our modeling studies suggest that an intrinsic twist at the carcinogen-purine linkage of ALII-N(2)-dG induces lesion site structural perturbations and conformational heterogeneity of damaged DNA. These structural characteristics correlate with the relative repair propensities of AA-adducts, where GGR recognition occurs for ALII-N(2)-dG, but is evaded for intrinsically planar ALII-N(6)-dA that minimally distorts DNA and restricts the conformational flexibility of the damaged duplex. The present analysis on the ALII adduct model systems will inspire future experimental studies on these adducts, and thereby may extend the list of structural factors that directly correlate with the propensity for GGR recognition.

Theoretical considerations for thresholds in chemical carcinogenesis

There is increasing evidence for non-linear relationships for gene mutations, chromosomal aberrations and even tumor incidences in response to low doses of genotoxic carcinogens. To attain the biological relevance of such non-linear responses, there is a need to identify the underlying defense mechanisms that allow tolerance to low doses of genotoxicants. This communication discusses presumptive cancer prevention mechanisms that may contribute to thresholds, i.e. points of departure, for each endpoint, from initial DNA lesion to tumor formation. We discuss a sequential order of genome protection during carcinogenesis where genotoxicant scavenging, cellular efflux, DNA repair, elimination of damaged cells by apoptosis, autophagy, silencing by DNA damage-triggered replicative senescence, and finally, elimination of transformed (premalignant) cells by the immune system are thought to be responsible for a threshold in tumor formation. We highlight DNA repair, for which experimental evidence has been recently provided to dictate a role in PoDs. In conclusion, from a theoretical perspective it is reasonable to posit that tolerance to low dose levels exists for each requisite step of tumor formation and these tolerance mechanisms are critical in determining thresholds in chemical carcinogenesis.

(DNS)C: a fluorescent, environmentally sensitive cytidine derivative for the direct detection of GGG triad sequences

With the goal of developing a fluorescent nucleoside sensitive to its environment, in this study we synthesized (DNS)C, a novel modified 2'-deoxycytidine bearing a 5-(dimethylamino)naphthalene-1-sulfonyl (dansyl) moiety at the N4 position, and tested its properties in monomeric and oligomeric states. (DNS)C undergoes intramolecular photoinduced electron transfer between its dansyl and cytosine units, resulting in remarkable changes in fluorescence that depend on the choice of solvent. In addition, the fluorescence behavior and thermal stability of oligonucleotides containing (DNS)C are dependent on the nature of the flanking and neighboring bases. Notably, (DNS)C exhibits fluorescence enhancement only in fully matched duplex DNA containing a GGG triad sequence. The environmental sensitivity of (DNS)C can be exploited as a fluorescence tool for monitoring the interactions of DNA with other biomolecules, including DNA, RNA, and proteins.

Nuclear magnetic resonance studies of an N2-guanine adduct derived from the tumorigen dibenzo[a,l]pyrene in DNA: impact of adduct stereochemistry, size, and local DNA sequence on solution conformations

The dimensions and arrangements of aromatic rings (topology) in adducts derived from the reactions of polycyclic aromatic hydrocarbon (PAH) diol epoxide metabolites with DNA influence the distortions and stabilities of double-stranded DNA, and hence their recognition and processing by the human nucleotide excision repair (NER) system. Dibenzo[a,l]pyrene (DB[a,l]P) is a highly tumorigenic six-ring PAH, which contains a nonplanar and aromatic fjord region that is absent in the structurally related bay region five-ring PAH benzo[a]pyrene (B[a]P). The PAH diol epoxide–DNA adducts formed include the stereoisomeric 14S and 14Rtrans-anti-DB[a,l]P-N²-dG and the stereochemically analogous 10S- and 10R-B[a]P-N²-dG (B[a]P-dG) guanine adducts. However, nuclear magnetic resonance (NMR) solution studies of the 14S-DB[a,l]P-N²-dG adduct in DNA have not yet been presented. Here we have investigated the 14S-DB[a,l]P-N²-dG adduct in two different sequence contexts using NMR methods with distance-restrained molecular dynamics simulations. In duplexes with dC opposite the adduct deleted, a well-resolved base-displaced intercalative adduct conformation can be observed. In full duplexes, in contrast to the intercalated 14R stereoisomeric adduct, the bulky DB[a,l]P residue in the 14S adduct is positioned in a greatly widened and distorted minor groove, with significant disruptions and distortions of base pairing at the lesion site and two 5′-side adjacent base pairs. These unique structural features are significantly different from those of the stereochemically analogous but smaller B[a]P-dG adduct. The greater size and different topology of the DB[a,l]P aromatic ring system lead to greater structurally destabilizing DNA distortions that are partially compensated by stabilizing DB[a,l]P-DNA van der Waals interactions, whose combined effects impact the NER response to the adduct. These structural results broaden our understanding of the structure–function relationship in NER.

Structural basis for the recognition of diastereomeric 5',8-cyclo-2'-deoxypurine lesions by the human nucleotide excision repair system

The hydroxyl radical is a powerful oxidant that generates DNA lesions including the stereoisomeric R and S 5′,8-cyclo-2′-deoxyadenosine (cdA) and 5′,8-cyclo-2′-deoxyguanosine (cdG) pairs that have been detected in cellular DNA. Unlike some other oxidatively generated DNA lesions, cdG and cdA are repaired by the human nucleotide excision repair (NER) apparatus. The relative NER efficiencies of all four cyclopurines were measured and compared in identical human HeLa cell extracts for the first time under identical conditions, using identical sequence contexts. The cdA and cdG lesions were excised with similar efficiencies, but the efficiencies for both 5′R cyclopurines were greater by a factor of ∼2 than for the 5′S lesions. Molecular modeling and dynamics simulations have revealed structural and energetic origins of this difference in NER-incision efficiencies. These lesions cause greater DNA backbone distortions and dynamics relative to unmodified DNA in 5′R than in 5′S stereoisomers, producing greater impairment in van der Waals stacking interaction energies in the 5′R cases. The locally impaired stacking interaction energies correlate with relative NER incision efficiencies, and explain these results on a structural basis in terms of differences in dynamic perturbations of the DNA backbone imposed by the R and S covalent 5′,8 bonds.

Base-displaced intercalation of the 2-amino-3-methylimidazo[4,5-f]quinolone N2-dG adduct in the NarI DNA recognition sequence

2-Amino-3-methylimidazo[4,5-f]quinolone (IQ), a heterocyclic amine found in cooked meats, undergoes bioactivation to a nitrenium ion, which alkylates guanines at both the C8-dG and N2-dG positions. The conformation of a site-specific N2-dG-IQ adduct in an oligodeoxynucleotide duplex containing the iterated CG repeat restriction site of the NarI endonuclease has been determined. The IQ moiety intercalates, with the IQ H4a and CH3 protons facing the minor groove, and the IQ H7a, H8a and H9a protons facing the major groove. The adducted dG maintains the anti-conformation about the glycosyl bond. The complementary dC is extruded into the major groove. The duplex maintains its thermal stability, which is attributed to stacking between the IQ moiety and the 5'- and 3'-neighboring base pairs. This conformation is compared to that of the C8-dG-IQ adduct in the same sequence, which also formed a 'base-displaced intercalated' conformation. However, the C8-dG-IQ adopted the syn conformation placing the Watson-Crick edge of the modified dG into the major groove. In addition, the C8-dG-IQ adduct was oriented with the IQ CH3 group and H4a and H5a facing the major groove. These differences may lead to differential processing during DNA repair and replication.

Ribonucleotides as nucleotide excision repair substrates

The incorporation of ribonucleotides in DNA has attracted considerable notice in recent years, since the pool of ribonucleotides can exceed that of the deoxyribonucleotides by at least 10-20-fold, and single ribonucleotide incorporation by DNA polymerases appears to be a common event. Moreover ribonucleotides are potentially mutagenic and lead to genome instability. As a consequence, errantly incorporated ribonucleotides are rapidly repaired in a process dependent upon RNase H enzymes. On the other hand, global genomic nucleotide excision repair (NER) in prokaryotes and eukaryotes removes damage caused by covalent modifications that typically distort and destabilize DNA through the production of lesions derived from bulky chemical carcinogens, such as polycyclic aromatic hydrocarbon metabolites, or via crosslinking. However, a recent study challenges this lesion-recognition paradigm. The work of Vaisman et al. (2013) [34] reveals that even a single ribonucleotide embedded in a deoxyribonucleotide duplex is recognized by the bacterial NER machinery in vitro. In their report, the authors show that spontaneous mutagenesis promoted by a steric-gate pol V mutant increases in uvrA, uvrB, or uvrC strains lacking rnhB (encoding RNase HII) and to a greater extent in an NER-deficient strain lacking both RNase HI and RNase HII. Using purified UvrA, UvrB, and UvrC proteins in in vitro assays they show that despite causing little distortion, a single ribonucleotide embedded in a DNA duplex is recognized and doubly-incised by the NER complex. We present the hypothesis to explain the recognition and/or verification of this small lesion, that the critical 2'-OH of the ribonucleotide - with its unique electrostatic and hydrogen bonding properties - may act as a signal through interactions with amino acid residues of the prokaryotic NER complex that are not possible with DNA. Such a mechanism might also be relevant if it were demonstrated that the eukaryotic NER machinery likewise incises an embedded ribonucleotide in DNA.

Role of structural and energetic factors in regulating repair of a bulky DNA lesion with different opposite partner bases

Extensive molecular modeling with molecular dynamics simulations and van der Waals energy analyses were used to elucidate the striking finding that a mutagenic benzo[a]pyrene-derived DNA lesion, the base-displaced intercalated 10R-(+)-cis-anti-B[a]P-N(2)-dG (G*), manifests large differences in nucleotide excision repair (NER) efficiencies in DNA duplexes, which depend on the identities of the partner base opposite G*. The nature of the partner base causes marked differences in the extent of its major groove extrusion and dynamics, as well as energetic stability of the intercalation pocket that parallels the relative NER efficiencies.

Free energy profiles of base flipping in intercalative polycyclic aromatic hydrocarbon-damaged DNA duplexes: energetic and structural relationships to nucleotide excision repair susceptibility

The crystal structure of Rad4/Rad23, the yeast homolog of the human nucleotide excision repair (NER) lesion recognition factor XPC-RAD23B ( Min , J. H. and Pavletich , N. P. ( 2007 ) Nature 449 , 570 - 575 ) reveals that the lesion-partner base is flipped out of the helix and binds to amino acids of the protein. This suggests the hypothesis that the flipping of this partner base must overcome a free energy barrier, which constitutes one element contributing to changes in the thermodynamic properties induced by the DNA damage and sensed by the recognition protein. We explored this hypothesis by computing complete flipping free energy profiles for two lesions derived from the procarcinogenic polycyclic aromatic hydrocarbons (PAHs), dibenzo[a,l]pyrene (DB[a,l]P) and benzo[a]pyrene (B[a]P), R-trans-anti-DB[a,l]P-N(6)-dA (R-DB[a,l]P-dA) and R-trans-anti-B[a]P-N(6)-dA (R-B[a]P-dA), and the corresponding unmodified duplex. The DB[a,l]P and B[a]P adducts differ in number and organization of their aromatic rings. We integrate these results with prior profiles for the R-trans-anti-DB[a,l]P-dG adduct ( Zheng , H. et al. ( 2010 ) Chem. Res. Toxicol. 23 , 1868 - 1870 ). All adopt conformational themes involving intercalation of the PAH aromatic ring system into the DNA duplex; however, R-DB[a,l]P-dA and R-B[a]P-dA intercalate from the major groove, while R-DB[a,l]P-dG intercalates from the minor groove. These structural differences produce different computed van der Waals stacking interaction energies between the flipping partner base with the lesion aromatic ring system and adjacent bases; we find that the better the stacking, the higher the relative flipping free energy barrier and hence lower flipping probability. The better relative NER susceptibilities correlate with greater ease of flipping in these three differently intercalated lesions. In addition to partner base flipping, the Rad4/Rad23 crystal structure shows that a protein-β-hairpin, BHD3, intrudes from the major groove side between the DNA strands at the lesion site. We present a molecular modeling study for the R-DB[a,l]P-dG lesion in Rad4/Rad23 showing BHD3 β-hairpin intrusion with lesion eviction, and we hypothesize that lesion steric effects play a role in the recognition of intercalated adducts.

Adenine-DNA adducts derived from the highly tumorigenic Dibenzo[a,l]pyrene are resistant to nucleotide excision repair while guanine adducts are not

The structural origins of differences in susceptibilities of various DNA lesions to nucleotide excision repair (NER) are poorly understood. Here we compared, in the same sequence context, the relative NER dual incision efficiencies elicited by two stereochemically distinct pairs of guanine (N(2)-dG) and adenine (N(6)-dA) DNA lesions, derived from enantiomeric genotoxic diol epoxides of the highly tumorigenic fjord region polycyclic aromatic hydrocarbon dibenzo[a,l]pyrene (DB[a,l]P). Remarkably, in cell-free HeLa cell extracts, the guanine adduct with R absolute chemistry at the N(2)-dG linkage site is ∼35 times more susceptible to NER dual incisions than the stereochemically identical N(6)-dA adduct. For the guanine and adenine adducts with S stereochemistry, a similar but somewhat smaller effect (factor of ∼15) is observed. The striking resistance of the bulky N(6)-dA in contrast to the modest to good susceptibilities of the N(2)-dG adducts to NER is interpreted in terms of the balance between lesion-induced DNA distorting and DNA stabilizing van der Waals interactions in their structures, that are partly reflected in the overall thermal stabilities of the modified duplexes. Our results are consistent with the hypothesis that the high genotoxic activity of DB[a,l]P is related to the formation of NER-resistant and persistent DB[a,l]P-derived adenine adducts in cellular DNA.

Biomonitoring of polycyclic aromatic hydrocarbons from coke oven emissions and reproductive toxicity in nonsmoking workers

The objective of the cross-sectional study was to assess whether exposure to polycyclic aromatic hydrocarbons (PAHs) from coke oven emissions contributed to alteration of semen quality and sperm DNA integrity in nonsmoking workers. Nonsmoking coke oven workers from a steel plant in Taiwan served as the exposure groups (topside-oven workers for the high exposure group and side-oven workers for the low exposure group), and administrators and security personnel in the plant served as the control. An exposure assessment was conducted to determine both particulate and gaseous phase of PAH levels and urinary 1-hydroxypyrene (1-OHP) levels. Semen quality was analyzed according to WHO guidelines. DNA fragmentation and bulky DNA adducts were measured to assess sperm DNA integrity. There was no significant difference in sperm concentrations, vitality, and DNA fragmentation between the exposed group and the control. The high exposure group experienced significantly lower percentages of normal morphology as compared with the control (p=0.0001). Bulky DNA adducts were detected in the exposed group that were significant higher than the control (p=0.04). Exposure to PAHs from coke-oven emissions could contribute to increased levels of bulky DNA adducts in sperm.

Nuclear magnetic resonance solution structure of an N(2)-guanine DNA adduct derived from the potent tumorigen dibenzo[a,l]pyrene: intercalation from the minor groove with ruptured Watson-Crick base pairing

The most potent tumorigen identified among the polycyclic aromatic hydrocarbons (PAH) is the nonplanar fjord region dibenzo[a,l]pyrene (DB[a,l]P). It is metabolically activated in vivo through the widely studied diol epoxide (DE) pathway to form covalent adducts with DNA bases, predominantly guanine and adenine. The (+)-11S,12R,13R,14S DE enantiomer forms adducts via its C14 position with the exocyclic amino group of guanine. Here, we present the first nuclear magnetic resonance solution structure of a DB[a,l]P-derived adduct, the 14R-(+)-trans-anti-DB[a,l]P-N(2)-dG (DB[a,l]P-dG) lesion in double-stranded DNA. In contrast to the stereochemically identical benzo[a]pyrene-derived N(2)-dG adduct (B[a]P-dG) in which the B[a]P rings reside in the B-DNA minor groove on the 3'-side of the modifed deoxyguanosine, in the DB[a,l]P-derived adduct the DB[a,l]P rings intercalate into the duplex on the 3'-side of the modified base from the sterically crowded minor groove. Watson-Crick base pairing of the modified guanine with the partner cytosine is broken, but these bases retain some stacking with the bulky DB[a,l]P ring system. This new theme in PAH DE-DNA adduct conformation differs from (1) the classical intercalation motif in which Watson-Crick base pairing is intact at the lesion site and (2) the base-displaced intercalation motif in which the damaged base and its partner are extruded from the helix. The structural considerations that lead to the intercalated conformation of the DB[a,l]P-dG lesion in contrast to the minor groove alignment of the B[a]P-dG adduct, and the implications of the DB[a,l]P-dG conformational motif for the recognition of such DNA lesions by the human nucleotide excision repair apparatus, are discussed.

Nucleotide excision repair of 2-acetylaminofluorene- and 2-aminofluorene-(C8)-guanine adducts: molecular dynamics simulations elucidate how lesion structure and base sequence context impact repair efficiencies

Nucleotide excision repair (NER) efficiencies of DNA lesions can vary by orders of magnitude, for reasons that remain unclear. An example is the pair of N-(2′-deoxyguanosin-8-yl)-2-aminofluorene (dG-C8-AF) and N-(2′-deoxyguanosin-8-yl)-2-acetylaminofluorene (dG-C8-AAF) adducts that differ by a single acetyl group. The NER efficiencies in human HeLa cell extracts of these lesions are significantly different when placed at G₁, G₂ or G₃ in the duplex sequence (5′-CTCG₁G₂CG₃CCATC-3′) containing the NarI mutational hot spot. Furthermore, the dG-C8-AAF adduct is a better substrate of NER than dG-C8-AF in all three NarI sequence contexts. The conformations of each of these adducts were investigated by Molecular dynamics (MD) simulation methods. In the base-displaced conformational family, the greater repair susceptibility of dG-C8-AAF in all sequences stems from steric hindrance effects of the acetyl group which significantly diminish the adduct-base stabilizing van der Waals stacking interactions relative to the dG-C8-AF case. Base sequence context effects for each adduct are caused by differences in helix untwisting and minor groove opening that are derived from the differences in stacking patterns. Overall, the greater NER efficiencies are correlated with greater extents of base sequence-dependent local untwisting and minor groove opening together with weaker stacking interactions.