Base sequence context effects on nucleotide excision repair

Moderation of Structural DNA Properties by Coupled Dinucleotide Contents in Eukaryotes

Dinucleotides are known as determinants for various structural and physiochemical properties of DNA and for binding affinities of proteins to DNA. These properties (e.g., stiffness) and bound proteins (e.g., transcription factors) are known to influence important biological functions, such as transcription regulation and 3D chromatin organization. Accordingly, the question arises of how the considerable variations in dinucleotide contents of eukaryotic chromosomes could still provide consistent DNA properties resulting in similar functions and 3D conformations. In this work, we investigate the hypothesis that coupled dinucleotide contents influence DNA properties in opposite directions to moderate each other's influences. Analyzing all 2478 chromosomes of 155 eukaryotic species, considering bias from coding sequences and enhancers, we found sets of correlated and anti-correlated dinucleotide contents. Using computational models, we estimated changes of DNA properties resulting from this coupling. We found that especially pure A/T dinucleotides (AA, TT, AT, TA), known to influence histone positioning and AC/GT contents, are relevant moderators and that, e.g., the Roll property, which is known to influence histone affinity of DNA, is preferably moderated. We conclude that dinucleotide contents might indirectly influence transcription and chromatin 3D conformation, via regulation of histone occupancy and/or other mechanisms.

Individual Genetic Heterogeneity

Genetic variation has been widely covered in literature, however, not from the perspective of an individual in any species. Here, a synthesis of genetic concepts and variations relevant for individual genetic constitution is provided. All the different levels of genetic information and variation are covered, ranging from whether an organism is unmixed or hybrid, has variations in genome, chromosomes, and more locally in DNA regions, to epigenetic variants or alterations in selfish genetic elements. Genetic constitution and heterogeneity of microbiota are highly relevant for health and wellbeing of an individual. Mutation rates vary widely for variation types, e.g., due to the sequence context. Genetic information guides numerous aspects in organisms. Types of inheritance, whether Mendelian or non-Mendelian, zygosity, sexual reproduction, and sex determination are covered. Functions of DNA and functional effects of variations are introduced, along with mechanism that reduce and modulate functional effects, including TARAR countermeasures and intraindividual genetic conflict. TARAR countermeasures for tolerance, avoidance, repair, attenuation, and resistance are essential for life, integrity of genetic information, and gene expression. The genetic composition, effects of variations, and their expression are considered also in diseases and personalized medicine. The text synthesizes knowledge and insight on individual genetic heterogeneity and organizes and systematizes the central concepts.

DNA Repair in Haploid Context

DNA repair is a well-covered topic as alteration of genetic integrity underlies many pathological conditions and important transgenerational consequences. Surprisingly, the ploidy status is rarely considered although the presence of homologous chromosomes dramatically impacts the repair capacities of cells. This is especially important for the haploid gametes as they must transfer genetic information to the offspring. An understanding of the different mechanisms monitoring genetic integrity in this context is, therefore, essential as differences in repair pathways exist that differentiate the gamete’s role in transgenerational inheritance. Hence, the oocyte must have the most reliable repair capacity while sperm, produced in large numbers and from many differentiation steps, are expected to carry de novo variations. This review describes the main DNA repair pathways with a special emphasis on ploidy. Differences between Saccharomyces cerevisiae and Schizosaccharomyces pombe are especially useful to this aim as they can maintain a diploid and haploid life cycle respectively.

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.

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.

Nucleotide Excision Repair and Impact of Site-Specific 5',8-Cyclopurine and Bulky DNA Lesions on the Physical Properties of Nucleosomes

The nonbulky 5',8-cyclopurine DNA lesions (cP) and the bulky, benzo[ a]pyrene diol epoxide-derived stereoisomeric cis- and trans- N2-guanine adducts (BPDE-dG) are good substrates of the human nucleotide excision repair (NER) mechanism. These DNA lesions were embedded at the In or Out rotational settings near the dyad axis in nucleosome core particles reconstituted either with native histones extracted from HeLa cells (HeLa-NCP) or with recombinant histones (Rec-NCP). The cP lesions are completely resistant to NER in human HeLa cell extracts. The BPDE-dG adducts are also NER-resistant in Rec-NCPs but are good substrates of NER in HeLa-NCPs. The four BPDE-dG adduct samples are excised with different efficiencies in free DNA, but in HeLa-NCPs, the efficiencies are reduced by a common factor of 2.2 ± 0.2 relative to the NER efficiencies in free DNA. The NER response of the BPDE-dG adducts in HeLa-NCPs is not directly correlated with the observed differences in the thermodynamic destabilization of HeLa-NCPs, the Förster resonance energy transfer values, or hydroxyl radical footprint patterns and is weakly dependent on the rotational settings. These and other observations suggest that NER is initiated by the binding of the DNA damage-sensing NER factor XPC-RAD23B to a transiently opened BPDE-modified DNA sequence that corresponds to the known footprint of XPC-DNA-RAD23B complexes (≥30 bp). These observations are consistent with the hypothesis that post-translational modifications and the dimensions and properties of the DNA lesions are the major factors that have an impact on the dynamics and initiation of NER in nucleosomes.

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.

Molecular basis for damage recognition and verification by XPC-RAD23B and TFIIH in nucleotide excision repair

Global genome nucleotide excision repair (GG-NER) is the main pathway for the removal of bulky lesions from DNA and is characterized by an extraordinarily wide substrate specificity. Remarkably, the efficiency of lesion removal varies dramatically and certain lesions escape repair altogether and are therefore associated with high levels of mutagenicity. Central to the multistep mechanism of damage recognition in NER is the sensing of lesion-induced thermodynamic and structural alterations of DNA by the XPC-RAD23B protein and the verification of the damage by the transcription/repair factor TFIIH. Additional factors contribute to the process: UV-DDB, for the recognition of certain UV-induced lesions in particular in the context of chromatin, while the XPA protein is believed to have a role in damage verification and NER complex assembly. Here we consider the molecular mechanisms that determine repair efficiency in GG-NER based on recent structural, computational, biochemical, cellular and single molecule studies of XPC-RAD23B and its yeast ortholog Rad4. We discuss how the actions of XPC-RAD23B are integrated with those of other NER proteins and, based on recent high-resolution structures of TFIIH, present a structural model of how XPC-RAD23B and TFIIH cooperate in damage recognition and verification.

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.

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.

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.

Resistance to Nucleotide Excision Repair of Bulky Guanine Adducts Opposite Abasic Sites in DNA Duplexes and Relationships between Structure and Function

The nucleotide excision repair of certain bulky DNA lesions is abrogated in some specific non-canonical DNA base sequence contexts, while the removal of the same lesions by the nucleotide excision repair mechanism is efficient in duplexes in which all base pairs are complementary. Here we show that the nucleotide excision repair activity in human cell extracts is moderate-to-high in the case of two stereoisomeric DNA lesions derived from the pro-carcinogen benzo[a]pyrene (cis- and trans-B[a]P-N2-dG adducts) in a normal DNA duplex. By contrast, the nucleotide excision repair activity is completely abrogated when the canonical cytosine base opposite the B[a]P-dG adducts is replaced by an abasic site in duplex DNA. However, base excision repair of the abasic site persists. In order to understand the structural origins of these striking phenomena, we used NMR and molecular spectroscopy techniques to evaluate the conformational features of 11mer DNA duplexes containing these B[a]P-dG lesions opposite abasic sites. Our results show that in these duplexes containing the clustered lesions, both B[a]P-dG adducts adopt base-displaced intercalated conformations, with the B[a]P aromatic rings intercalated into the DNA helix. To explain the persistence of base excision repair in the face of the opposed bulky B[a]P ring system, molecular modeling results suggest how the APE1 base excision repair endonuclease, that excises abasic lesions, can bind productively even with the trans-B[a]P-dG positioned opposite the abasic site. We hypothesize that the nucleotide excision repair resistance is fostered by local B[a]P residue-DNA base stacking interactions at the abasic sites, that are facilitated by the absence of the cytosine partner base in the complementary strand. More broadly, this study sets the stage for elucidating the interplay between base excision and nucleotide excision repair in processing different types of clustered DNA lesions that are substrates of nucleotide excision repair or base excision repair mechanisms.

DNA sequence context greatly affects the accuracy of bypass across an ultraviolet light 6-4 photoproduct in mammalian cells

Translesion DNA synthesis (TLS) is a DNA damage tolerance mechanism carried out by low-fidelity DNA polymerases that bypass DNA lesions, which overcomes replication stalling. Despite the miscoding nature of most common DNA lesions, several of them are bypassed in mammalian cells in a relatively accurate manner, which plays a key role maintaining a low mutation load. Whereas it is generally agreed that TLS across the major UV and sunlight induced DNA lesion, the cyclobutane pyrimidine dimer (CPD), is accurate, there were conflicting reports on whether the same is true for the thymine-thymine pyrimidine-pyrimidone(6-4) ultraviolet light photoproduct (TT6-4PP), which represents the second most common class of UV lesions. Using a TLS assay system based on gapped plasmids carrying site-specific TT6-4PP lesions in defined sequence contexts we show that the DNA sequence context markedly affected both the extent and accuracy of TLS. The sequence exhibiting higher TLS exhibited also higher error-frequency, caused primarily by semi-targeted mutations, at the nearest nucleotides flanking the lesion. Our results resolve the discrepancy reported on TLS across TT6-4PP, and suggest that TLS is more accurate in human cells than in mouse cells.

Genomic landscape of carcinogen-induced and genetically induced mouse skin squamous cell carcinoma

Mouse models of cancers are routinely used to study cancer biology. However, it remains unclear whether carcinogenesis in mice is driven by the same spectrum of genomic alterations found in humans. Here we conducted a comprehensive genomic analysis of 9,10-dimethyl-1,2-benzanthracene (DMBA)-induced skin cancer, the most commonly used skin cancer model, which appears as benign papillomas that progress into squamous cell carcinomas (SCCs). We also studied genetically induced SCCs that expressed G12D mutant Kras (Kras G12D) but were deficient for p53. Using whole-exome sequencing, we uncovered a characteristic mutational signature of DMBA-induced SCCs. We found that the vast majority of DMBA-induced SCCs presented recurrent mutations in Hras, Kras or Rras2 and mutations in several additional putative oncogenes and tumor-suppressor genes. Similar genes were recurrently mutated in mouse and human SCCs that were from different organs or had been exposed to different carcinogens. Invasive SCCs, but not papillomas, presented substantial chromosomal aberrations, especially in DMBA-induced and genetically induced Trp53-mutated SCCs. Metastasis occurred through sequential spreading, with relatively few additional genetic events. This study provides a framework for future functional cancer genomic studies in mice.

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.

Unusual sequence effects on nucleotide excision repair of arylamine lesions: DNA bending/distortion as a primary recognition factor

The environmental arylamine mutagens are implicated in the etiology of various sporadic human cancers. Arylamine-modified dG lesions were studied in two fully paired 11-mer duplexes with a -G*CN- sequence context, in which G* is a C8-substituted dG adduct derived from fluorinated analogs of 4-aminobiphenyl (FABP), 2-aminofluorene (FAF) or 2-acetylaminofluorene (FAAF), and N is either dA or dT. The FABP and FAF lesions exist in a simple mixture of 'stacked' (S) and 'B-type' (B) conformers, whereas the N-acetylated FAAF also samples a 'wedge' (W) conformer. FAAF is repaired three to four times more efficiently than FABP and FAF. A simple A- to -T polarity swap in the G*CA/G*CT transition produced a dramatic increase in syn-conformation and resulted in 2- to 3-fold lower nucleotide excision repair (NER) efficiencies in Escherichia coli. These results indicate that lesion-induced DNA bending/thermodynamic destabilization is an important DNA damage recognition factor, more so than the local S/B-conformational heterogeneity that was observed previously for FAF and FAAF in certain sequence contexts. This work represents a novel 3'-next flanking sequence effect as a unique NER factor for bulky arylamine lesions in E. coli.

The efficiencies of damage recognition and excision correlate with duplex destabilization induced by acetylaminofluorene adducts in human nucleotide excision repair

Nucleotide excision repair (NER) removes lesions caused by environmental mutagens or UV light from DNA. A hallmark of NER is the extraordinarily wide substrate specificity, raising the question of how one set of proteins is able to recognize structurally diverse lesions. Two key features of good NER substrates are that they are bulky and thermodynamically destabilize DNA duplexes. To understand what the limiting step in damage recognition in NER is, we set out to test the hypothesis that there is a correlation of the degree of thermodynamic destabilization induced by a lesion, binding affinity to the damage recognition protein XPC-RAD23B, and overall NER efficiency. We chose to use acetylaminofluorene (AAF) and aminofluorene (AF) adducts at the C8 position of guanine in different positions within the NarI (GGCGCC) sequence, as it is known that the structures of the duplexes depend on the position of the lesion in this context. We found that the efficiency of NER and the binding affinity of the damage recognition factor XPC-RAD23B correlated with the thermodynamic destabilization induced by the lesion. Our study is the first systematic analysis correlating these three parameters and supports the idea that initial damage recognition by XPC-RAD23B is a key rate-limiting step in NER.

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.

Lack of recognition by global-genome nucleotide excision repair accounts for the high mutagenicity and persistence of aristolactam-DNA adducts

Exposure to aristolochic acid (AA), a component of Aristolochia plants used in herbal remedies, is associated with chronic kidney disease and urothelial carcinomas of the upper urinary tract. Following metabolic activation, AA reacts with dA and dG residues in DNA to form aristolactam (AL)-DNA adducts. These mutagenic lesions generate a unique TP53 mutation spectrum, dominated by A:T to T:A transversions with mutations at dA residues located almost exclusively on the non-transcribed strand. We determined the level of AL-dA adducts in human fibroblasts treated with AA to determine if this marked strand bias could be accounted for by selective resistance to global-genome nucleotide excision repair (GG-NER). AL-dA adduct levels were elevated in cells deficient in GG-NER and transcription-coupled NER, but not in XPC cell lines lacking GG-NER only. In vitro, plasmids containing a single AL-dA adduct were resistant to the early recognition and incision steps of NER. Additionally, the NER damage sensor, XPC-RAD23B, failed to specifically bind to AL-DNA adducts. However, placing AL-dA in mismatched sequences promotes XPC-RAD23B binding and renders this adduct susceptible to NER, suggesting that specific structural features of this adduct prevent processing by NER. We conclude that AL-dA adducts are not recognized by GG-NER, explaining their high mutagenicity and persistence in target tissues.

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

The nucleotide excision repair (NER) machinery, the primary defense against cancer-causing bulky DNA lesions, is surprisingly inefficient in recognizing certain mutagenic DNA adducts and other forms of DNA damage. However, the biochemical basis of resistance to repair remains poorly understood. To address this problem, we have investigated a series of intercalated DNA-adenine lesions derived from carcinogenic polycyclic aromatic hydrocarbon (PAH) diol epoxide metabolites that differ in their response to the mammalian NER apparatus. These stereoisomeric PAH-derived adenine lesions represent ideal model systems for elucidating the effects of structural, dynamic, and thermodynamic properties that determine the recognition of these bulky DNA lesions by NER factors. The objective of this work was to gain a systematic understanding of the relation between aromatic ring topology and adduct stereochemistry with existing experimental NER efficiencies and known thermodynamic stabilities of the damaged DNA duplexes. For this purpose, we performed 100 ns molecular dynamics studies of the lesions embedded in identical double-stranded 11-mer sequences. Our studies show that, depending on topology and stereochemistry, stabilizing PAH-DNA base van der Waals stacking interactions can compensate for destabilizing distortions caused by these lesions that can, in turn, cause resistance to NER. The results suggest that the balance between helix stabilizing and destabilizing interactions between the adduct and nearby DNA residues can account for the variability of NER efficiencies observed in this class of PAH-DNA lesions.

Structural, energetic and dynamic properties of guanine(C8)-thymine(N3) cross-links in DNA provide insights on susceptibility to nucleotide excision repair

The one-electron oxidation of guanine in DNA by carbonate radical anions, a decomposition product of peroxynitrosocarbonate which is associated with the inflammatory response, can lead to the formation of intrastrand cross-links between guanine and thymine bases [Crean et al. (Oxidation of single-stranded oligonucleotides by carbonate radical anions: generating intrastrand cross-links between guanine and thymine bases separated by cytosines. Nucleic Acids Res. 2008; 36: 742-755.)]. These involve covalent bonds between the C8 positions of guanine (G*) and N3 of thymine (T*) in 5'-d(…G*pT*…) and 5'-d(…G*pCpT*…) sequence contexts. We have performed nucleotide excision repair (NER) experiments in human HeLa cell extracts which show that the G*CT* intrastrand cross-link is excised with approximately four times greater efficiency than the G*T* cross-link embedded in 135-mer DNA duplexes. In addition, thermal melting studies reveal that both lesions significantly destabilize duplex DNA, and that the destabilization induced by the G*CT* cross-link is considerably greater. Consistent with this difference in NER, our computations show that both lesions dynamically distort and destabilize duplex DNA. They disturb Watson-Crick base-pairing and base-stacking interactions, and cause untwisting and minor groove opening. These structural perturbations are much more pronounced in the G*CT* than in the G*T* cross-link. Our combined experimental and computational studies provide structural and thermodynamic understanding of the features of the damaged duplexes that produce the most robust NER response.

Resistance of bulky DNA lesions to nucleotide excision repair can result from extensive aromatic lesion-base stacking interactions

The molecular basis of resistance to nucleotide excision repair (NER) of certain bulky DNA lesions is poorly understood. To address this issue, we have studied NER in human HeLa cell extracts of two topologically distinct lesions, one derived from benzo[a]pyrene (10R-(+)-cis-anti-B[a]P-N(2)-dG), and one from the food mutagen 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine (C8-dG-PhIP), embedded in either full or 'deletion' duplexes (the partner nucleotide opposite the lesion is missing). All lesions adopt base-displaced intercalated conformations. Both full duplexes are thermodynamically destabilized and are excellent substrates of NER. However, the identical 10R-(+)-cis-anti-B[a]P-N(2)-dG adduct in the deletion duplex dramatically enhances the thermal stability of this duplex, and is completely resistant to NER. Molecular dynamics simulations show that B[a]P lesion-induced distortion/destabilization is compensated by stabilizing aromatic ring system-base stacking interactions. In the C8-dG-PhIP-deletion duplex, the smaller size of the aromatic ring system and the mobile phenyl ring are less stabilizing and yield moderate NER efficiency. Thus, a partner nucleotide opposite the lesion is not an absolute requirement for the successful initiation of NER. Our observations are consistent with the hypothesis that carcinogen-base stacking interactions, which contribute to the local DNA stability, can prevent the successful insertion of an XPC β-hairpin into the duplex and the normal recruitment of other downstream NER factors.

Intercalative conformations of the 14R (+)- and 14S (-)-trans-anti-DB[a,l]P-N⁶-dA adducts: molecular modeling and MD simulations

Among the polycyclic aromatic hydrocarbon class of chemical carcinogens, dibenzo[a,l]pyrene (DB[a,l]P) is the most potent tumorigen that has been identified to date. Structurally, it is bulky with six aromatic rings, and it contains the nonplanar fjord-region. The conformational properties of DB[a,l]P-derived DNA adducts responsible for its extraordinary carcinogenicity are hence of great interest. We have carried out molecular modeling and MD simulations for the 14R (+)- and 14S (-)-trans-anti-DB[a,l]P-N⁶-dA adducts derived from the reactions of the DB[a,l]P diol epoxides with adenine in double-stranded DNA. The structures are based on the classically intercalated NMR solution structures of the analogous fjord-region benzo[c]phenanthrene-derived-N⁶-dA adducts. One objective was to gain insight on the impact of the more bulky DB[a,l]P ring system on the structural characteristics of the intercalative adduct conformations. A further objective was to elucidate the effect of the flexible twist associated with the sterically hindered aromatic ring in the fjord-region on the intercalated conformations, for comparison with the intercalated but planar bay-region benzo[a]pyrene-derived-N⁶-dA adducts. For the DB[a,l]P-N⁶-dA adducts, our results show that the 14R (+)-adduct is more favorably intercalated on the 5'-side of the modified adenine than the stereoisomeric 14S (-)-adduct, intercalated on its 3'-side. The 14R (+)-adduct manifests better van der Waals stacking interactions with flanking base pairs, less perturbed Watson-Crick hydrogen bonding, less local groove enlargement, less unwinding, and a lower solvent exposure than the 14S (-)-adduct. These structural findings are consistent with observed thermodynamic melting data, UV absorption properties, and fluorescence quenching studies. By contrast, the NMR solution structures for the analogous but less bulky B[c]Ph-derived adducts reveal no such stereoisomeric effect, while the planar bay-region benzo[a]pyrene-derived-N⁶-dA adducts do. Differences in nucleotide excision repair susceptibilities of the fjord and bay region adducts stem from distinctions in their intercalative conformations, produced by the intrinsic topological variations in their polycyclic aromatic ring systems.