Differential nucleotide excision repair susceptibility of bulky DNA adducts in different sequence contexts: hierarchies of recognition signals

Utilizing biological experimental data and molecular dynamics for the classification of mutational hotspots through machine learning

Motivation

Benzo[a]pyrene, a notorious DNA-damaging carcinogen, belongs to the family of polycyclic aromatic hydrocarbons commonly found in tobacco smoke. Surprisingly, nucleotide excision repair (NER) machinery exhibits inefficiency in recognizing specific bulky DNA adducts including Benzo[a]pyrene Diol-Epoxide (BPDE), a Benzo[a]pyrene metabolite. While sequence context is emerging as the leading factor linking the inadequate NER response to BPDE adducts, the precise structural attributes governing these disparities remain inadequately understood. We therefore combined the domains of molecular dynamics and machine learning to conduct a comprehensive assessment of helical distortion caused by BPDE-Guanine adducts in multiple gene contexts. Specifically, we implemented a dual approach involving a random forest classification-based analysis and subsequent feature selection to identify precise topological features that may distinguish adduct sites of variable repair capacity. Our models were trained using helical data extracted from duplexes representing both BPDE hotspot and nonhotspot sites within the TP53 gene, then applied to sites within TP53, cII, and lacZ genes.

Results

We show our optimized model consistently achieved exceptional performance, with accuracy, precision, and f1 scores exceeding 91%. Our feature selection approach uncovered that discernible variance in regional base pair rotation played a pivotal role in informing the decisions of our model. Notably, these disparities were highly conserved among TP53 and lacZ duplexes and appeared to be influenced by the regional GC content. As such, our findings suggest that there are indeed conserved topological features distinguishing hotspots and nonhotpot sites, highlighting regional GC content as a potential biomarker for mutation.

Availability and implementation

Code for comparing machine learning classifiers and evaluating their performance is available at https://github.com/jdavies24/ML-Classifier-Comparison, and code for analysing DNA structure with Curves+ and Canal using Random Forest is available at https://github.com/jdavies24/ML-classification-of-DNA-trajectories.

5-Formylcytosine mediated DNA-peptide cross-link induces predominantly semi-targeted mutations in both Escherichia coli and human cells

Histone proteins can become trapped on DNA in the presence of 5-formylcytosine (5fC) to form toxic DNA-protein conjugates. Their repair may involve proteolytic digestion resulting in DNA-peptide cross-links (DpCs). Here, we have investigated replication of a model DpC comprised of an 11-mer peptide (NH2-GGGKGLGK∗GGA) containing an oxy-lysine residue (K∗) conjugated to 5fC in DNA. Both CXG and CXT (where X = 5fC-DpC) sequence contexts were examined. Replication of both constructs gave low viability (<10%) in Escherichia coli, whereas TLS efficiency was high (72%) in HEK 293T cells. In E. coli, the DpC was bypassed largely error-free, inducing only 2 to 3% mutations, which increased to 4 to 5% with SOS. For both sequences, semi-targeted mutations were dominant, and for CXG, the predominant mutations were G→T and G→C at the 3'-base to the 5fC-DpC. In HEK 293T cells, 7 to 9% mutations occurred, and the dominant mutations were the semi-targeted G → T for CXG and T → G for CXT. These mutations were reduced drastically in cells deficient in hPol η, hPol ι or hPol ζ, suggesting a role of these TLS polymerases in mutagenic TLS. Steady-state kinetics studies using hPol η confirmed that this polymerase induces G → T and T → G transversions at the base immediately 3' to the DpC. This study reveals a unique replication pattern of 5fC-conjugated DpCs, which are bypassed largely error-free in both E. coli and human cells and induce mostly semi-targeted mutations at the 3' position to the lesion.

Effect of Guanine Adduct Size, Shape, and Linker Type on the Conformation of Adducted DNA: A DFT and Molecular Dynamics Study

DNA is damaged through various exogenous sources (e.g., automobile exhaust, tobacco smoke, and processed foods), which can yield diverse C8-dG bulky aryl adducts. Adducts are known to induce structural changes to DNA that can lead to various biological outcomes, ranging from cell death to diseases such as cancer. Unfortunately, the relationship between the chemical composition of the damaged product, the adducted DNA structure, and the biological consequences is not well understood, which limits the development of disease detection and prevention strategies. The present study uses density functional theory (DFT) calculations and quintuplicate 1 μs molecular dynamics (MD) simulations to characterize the structure of DNA containing 21 model C8-dG adducts that systematically differ in size (phenyl to pyrenyl), shape (α (2,3), β (3,4) fusion, or ring substitution), and nucleobase-aryl group linkage (N, O, and C-linked). DFT calculations reveal that the inherent structural features of the G nucleobase adducts are impacted by linker type and bulky moiety shape, but not size, with the conformational flexibility reducing with α-ring fusion and linker composition as N > O > C. These structural properties are maintained in nucleoside models, which also reveal an increased propensity for anti-to-syn rotation about the glycosidic bond with N < O < C linker type. Although these diverse chemical features do not influence the global structure of adducted DNA, the adducts differentially impact the conformation local to the adducted site, including the relative populations of structures with the bulky moiety in the major groove (B conformer) and intercalated (stacked) into the helix (S conformer). Specifically, while the smallest phenyl adducts favor the B conformation and the largest pyrenyl-derived adducts stabilize the S conformation, the B/S ratio decreases with an increase in ring size and N > O > C linker composition. The shape and size (length) of the adduct can further finetune the B/S ratio, with β-fused naphthyl or α-fused phenanthryl N-linked adducts and O or C-linked adducts containing ring substitution increasing the prevalence of the S adducted DNA conformation. Overall, this work uncovers the significant effect of bulky moiety size and linker type, as well as the lesser impact of aryl group shape, on adducted DNA structure, which suggests differential replication and repair outcomes, and thereby represents an important step toward rationalizing connections between the structure and biological consequences of diverse DNA adducts.

Establishing Linkages Among DNA Damage, Mutagenesis, and Genetic Diseases

DNA damage by chemicals, radiation, or oxidative stress leads to a mutational spectrum, which is complex because it is determined in part by lesion structure, the DNA sequence context of the lesion, lesion repair kinetics, and the type of cells in which the lesion is replicated. Accumulation of mutations may give rise to genetic diseases such as cancer and therefore understanding the process underlying mutagenesis is of immense importance to preserve human health. Chemical or physical agents that cause cancer often leave their mutational fingerprints, which can be used to back-calculate the molecular events that led to disease. To make a clear link between DNA lesion structure and the mutations a given lesion induces, the field of single-lesion mutagenesis was developed. In the last three decades this area of research has seen much growth in several directions, which we attempt to describe in this Perspective.

"Flexible hinge" dynamics in mismatched DNA revealed by fluorescence correlation spectroscopy

Altered unwinding/bending fluctuations at DNA lesion sites are implicated as plausible mechanisms for damage sensing by DNA-repair proteins. These dynamics are expected to occur on similar timescales as one-dimensional (1D) diffusion of proteins on DNA if effective in stalling these proteins as they scan DNA. We examined the flexibility and dynamics of DNA oligomers containing 3 base pair (bp) mismatched sites specifically recognized in vitro by nucleotide excision repair protein Rad4 (yeast ortholog of mammalian XPC). A previous Forster resonance energy transfer (FRET) study mapped DNA conformational distributions with cytosine analog FRET pair primarily sensitive to DNA twisting/unwinding deformations (Chakraborty et al. Nucleic Acids Res. 46: 1240-1255 (2018)). These studies revealed B-DNA conformations for nonspecific (matched) constructs but significant unwinding for mismatched constructs specifically recognized by Rad4, even in the absence of Rad4. The timescales of these unwinding fluctuations, however, remained elusive. Here, we labeled DNA with Atto550/Atto647N FRET dyes suitable for fluorescence correlation spectroscopy (FCS). With these probes, we detected higher FRET in specific, mismatched DNA compared with matched DNA, reaffirming unwinding/bending deformations in mismatched DNA. FCS unveiled the dynamics of these spontaneous deformations at ~ 300 µs with no fluctuations detected for matched DNA within the ~ 600 ns-10 ms FCS time window. These studies are the first to visualize anomalous unwinding/bending fluctuations in mismatched DNA on timescales that overlap with the < 500 µs "stepping" times of repair proteins on DNA. Such "flexible hinge" dynamics at lesion sites could arrest a diffusing protein to facilitate damage interrogation and recognition.

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.

Sequence context effects of replication of Fapy•dG in three mutational hot spot sequences of the p53 gene in human cells

Fapy•dG and 8-OxodGuo are formed in DNA from a common N7-dG radical intermediate by reaction with hydroxyl radical. Although cellular levels of Fapy•dG are often greater, its effects on replication are less well understood than those of 8-OxodGuo. In this study plasmid DNA containing Fapy•dG in three mutational hotspots of human cancers, codons 248, 249, and 273 of the p53 tumor suppressor gene, was replicated in HEK 293T cells. TLS efficiencies for the Fapy•dG containing plasmids varied from 72 to 89%, and were further reduced in polymerase-deficient cells. The mutation frequency (MF) of Fapy•dG ranged from 7.3 to 11.6%, with G→T and G→A as major mutations in codons 248 and 249 compared to primarily G→T in codon 273. Increased MF in hPol ι-, hPol κ-, and hPol ζ-deficient cells suggested that these polymerases more frequently insert the correct nucleotide dC opposite Fapy•dG, whereas decreased G→A in codons 248 and 249 and reduction of all mutations in codon 273 in hPol λ-deficient cells indicated hPol λ's involvement in Fapy•dG mutagenesis. In vitro kinetic analysis using isolated translesion synthesis polymerases and hPol λ incompletely corroborated the mutagenesis experiments, indicating codependence on other proteins in the cellular milieu. In conclusion, Fapy•dG mutagenesis is dependent on the DNA sequence context, but its bypass by the TLS polymerases is largely error-free.

Impact of DNA Adduct Size, Number, and Relative Position on the Toxicity of Aromatic Amines: A Molecular Dynamics Case Study of ANdG- and APdG-Containing DNA Duplexes

Human exposure to aromatic amines (AAs) can result in carcinogenic DNA adducts. To complement previous work geared toward understanding the mutagenicity of AA-derived adducts, which has almost exclusively studied (monoadducted) DNA containing a single lesion, the present work provides the first in-depth comparison of the structure of monoadducted and diadducted DNA duplexes. Specifically, molecular dynamics (MD) simulations were initially performed on DNA containing the nonmutagenic single-ringed N-(deoxyguanosin-8-yl)-aniline (^ANdG) or the mutagenic four-ringed N-(deoxyguanosin-8-yl)-1-aminopyrene (^APdG) lesion at G₁, G₂, or G₃ in the AA deletion hotspot (5'-G₁G₂CG₃CC) in the anti or syn glycosidic orientation (B/S duplex conformation). Subsequently, diadducted strands were assessed that span each combination of damaged sites (G₁G₂ (nearest neighbors), G₂G₃ (next-nearest neighbors), and G₁G₃ (two intervening nucleotides)) and anti/syn lesion glycosidic orientations. Despite other N-linked C8-dG adducts exhibiting sequence dependence conformational heterogeneity, a single ^ANdG or ^APdG lesion induces helical conformational homogeneity that is exclusively controlled by aryl moiety size. However, the preferred damaged DNA conformation can change upon the addition of a second adduct depending on lesion separation, with neighboring lesions stabilizing a nonmutagenic conformation and next-nearest damaged sites stabilizing a promutagenic conformation regardless of adduct size. As a result, diadducted DNA is found to adopt conformations that are unfavored for the corresponding monoadducted system, pointing to differential replication and repair outcomes for diadducted DNA compared to those for monoadducted DNA. Thus, although the toxicity of monoadducted DNA is most significantly dictated by lesion size, the toxicity can increase or decrease upon a second damaging event depending on lesion size and relative position. Overall, our work adds the number of lesions and their spatial separation to the growing list of factors that determine the structure and biological outcomes of adducted DNA.

Structure and mechanism of pyrimidine-pyrimidone (6-4) photoproduct recognition by the Rad4/XPC nucleotide excision repair complex

Failure in repairing ultraviolet radiation-induced DNA damage can lead to mutations and cancer. Among UV-lesions, the pyrimidine–pyrimidone (6-4) photoproduct (6-4PP) is removed from the genome much faster than the cyclobutane pyrimidine dimer (CPD), owing to the more efficient recognition of 6-4PP by XPC-RAD23B, a key initiator of global-genome nucleotide excision repair (NER). Here, we report a crystal structure of a Rad4–Rad23 (yeast XPC-Rad23B ortholog) bound to 6-4PP-containing DNA and 4-μs molecular dynamics (MD) simulations examining the initial binding of Rad4 to 6-4PP or CPD. This first structure of Rad4/XPC bound to a physiological substrate with matched DNA sequence shows that Rad4 flips out both 6-4PP-containing nucleotide pairs, forming an ‘open’ conformation. The MD trajectories detail how Rad4/XPC initiates ‘opening’ 6-4PP: Rad4 initially engages BHD2 to bend/untwist DNA from the minor groove, leading to unstacking and extrusion of the 6-4PP:AA nucleotide pairs towards the major groove. The 5′ partner adenine first flips out and is captured by a BHD2/3 groove, while the 3′ adenine extrudes episodically, facilitating ensuing insertion of the BHD3 β-hairpin to open DNA as in the crystal structure. However, CPD resists such Rad4-induced structural distortions. Untwisting/bending from the minor groove may be a common way to interrogate DNA in NER.

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.

Adaptation of an amplicon-based human cancer next-generation sequencing panel assay for murine tumors

Unlike humans, inbred genetically engineered mice have minimal inter-individual variation and, consequently, offer substantially increased statistical power for robust definition of recurrent cooperating cancer mutations. While technically feasible, whole exome sequencing is expensive and extremely data-intensive. Somatic mutation analysis using panels of 25-75 genes now provides detailed insight into the biology of human tumors. Here we report an adaptation for mouse tumors of a human PCR amplicon-based panel (Ion Torrent Cancer Hotspot Panel v2) allowing analysis of 18 cancer genes, including Kras, Nras, Hras, Pten, Pik3ca and Smad4, and encompassing regions homologous to more than 2000 known human cancer mutations.

Effects of Black Raspberry Extract and Berry Compounds on Repair of DNA Damage and Mutagenesis Induced by Chemical and Physical Agents in Human Oral Leukoplakia and Rat Oral Fibroblasts

Black raspberries (BRB) have been shown to inhibit carcinogenesis in a number of systems, with most studies focusing on progression. Previously we reported that an anthocyanin-enriched black raspberry extract (BE) enhanced repair of dibenzo-[a,l]-pyrene dihydrodiol (DBP-diol)-induced DNA adducts and inhibited DBP-diol and DBP-diolepoxide (DBPDE)-induced mutagenesis in a lacI rat oral fibroblast cell line, suggesting a role for BRB in the inhibition of initiation of carcinogenesis. Here we extend this work to protection by BE against DNA adduct formation induced by dibenzo-[a,l]-pyrene (DBP) in a human oral leukoplakia cell line (MSK) and to a second carcinogen, UV light. Treatment of MSK cells with DBP and DBPDE led to a dose-dependent increase in DBP-DNA adducts. Treatment of MSK cells with BE after addition of DBP reduced levels of adducts relative to cells treated with DBP alone, and treatment of rat oral fibroblasts with BE after addition of DBPDE inhibited mutagenesis. These observations showed that BE affected repair of DNA adducts and not metabolism of DBP. As a proof of principle we also tested aglycones of two anthocyanins commonly found in berries, delphinidin chloride and pelargonidin chloride. Delphinidin chloride reduced DBP-DNA adduct levels in MSK cells, while PGA did not. These results suggested that certain anthocyanins can enhance repair of bulky DNA adducts. As DBP and its metabolites induced formation of bulky DNA adducts, we investigated the effects of BE on genotoxic effects of a second carcinogen that induces bulky DNA damage, UV light. UV irradiation produced a dose-dependent increase in cyclobutanepyrimidine dimer levels in MSK cells, and post-UV treatment with BE resulted in lower cyclobutanepyrimidine dimer levels. Post-UV treatment of the rat lacI cells with BE reduced UV-induced mutagenesis. Taken together, the results demonstrate that BE extract reduces bulky DNA damage and mutagenesis and support a role for BRB in the inhibition of initiation of carcinogenesis.

Human genome-wide repair map of DNA damage caused by the cigarette smoke carcinogen benzo[a]pyrene

Benzo[a]pyrene (BaP), a polycyclic aromatic hydrocarbon, is the major cause of lung cancer. BaP forms covalent DNA adducts after metabolic activation and induces mutations. We have developed a method for capturing oligonucleotides carrying bulky base adducts, including UV-induced cyclobutane pyrimidine dimers (CPDs) and BaP diol epoxide-deoxyguanosine (BPDE-dG), which are removed from the genome by nucleotide excision repair. The isolated oligonucleotides are ligated to adaptors, and after damage-specific immunoprecipitation, the adaptor-ligated oligonucleotides are converted to dsDNA with an appropriate translesion DNA synthesis (TLS) polymerase, followed by PCR amplification and next-generation sequencing (NGS) to generate genome-wide repair maps. We have termed this method translesion excision repair-sequencing (tXR-seq). In contrast to our previously described XR-seq method, tXR-seq does not depend on repair/removal of the damage in the excised oligonucleotides, and thus it is applicable to essentially all DNA damages processed by nucleotide excision repair. Here we present the excision repair maps for CPDs and BPDE-dG adducts generated by tXR-Seq for the human genome. In addition, we report the sequence specificity of BPDE-dG excision repair using tXR-seq.

Removal of oxidatively generated DNA damage by overlapping repair pathways

It is generally believed that the mammalian nucleotide excision repair pathway removes DNA helix-distorting bulky DNA lesions, while small non-bulky lesions are repaired by base excision repair (BER). However, recent work demonstrates that the oxidativly generated guanine oxidation products, spiroimininodihydantoin (Sp), 5-guanidinohydantoin (Gh), and certain intrastrand cross-linked lesions, are good substrates of NER and BER pathways that compete with one another in human cell extracts. The oxidation of guanine by peroxynitrite is known to generate 5-guanidino-4-nitroimidazole (NIm) which is structurally similar to Gh, except that the 4-nitro group in NIm is replaced by a keto group in Gh. However, unlike Gh, NIm is an excellent substrate of BER, but not of NER. These and other related results are reviewed and discussed in this article.

Single-stranded DNA Binding by the Helix-Hairpin-Helix Domain of XPF Protein Contributes to the Substrate Specificity of the ERCC1-XPF Protein Complex

The nucleotide excision repair protein complex ERCC1-XPF is required for incision of DNA upstream of DNA damage. Functional studies have provided insights into the binding of ERCC1-XPF to various DNA substrates. However, because no structure for the ERCC1-XPF-DNA complex has been determined, the mechanism of substrate recognition remains elusive. Here we biochemically characterize the substrate preferences of the helix-hairpin-helix (HhH) domains of XPF and ERCC-XPF and show that the binding to single-stranded DNA (ssDNA)/dsDNA junctions is dependent on joint binding to the DNA binding domain of ERCC1 and XPF. We reveal that the homodimeric XPF is able to bind various ssDNA sequences but with a clear preference for guanine-containing substrates. NMR titration experiments and in vitro DNA binding assays also show that, within the heterodimeric ERCC1-XPF complex, XPF specifically recognizes ssDNA. On the other hand, the HhH domain of ERCC1 preferentially binds dsDNA through the hairpin region. The two separate non-overlapping DNA binding domains in the ERCC1-XPF heterodimer jointly bind to an ssDNA/dsDNA substrate and, thereby, at least partially dictate the incision position during damage removal. Based on structural models, NMR titrations, DNA-binding studies, site-directed mutagenesis, charge distribution, and sequence conservation, we propose that the HhH domain of ERCC1 binds to dsDNA upstream of the damage, and XPF binds to the non-damaged strand within a repair bubble.

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.

The somatic autosomal mutation matrix in cancer genomes

DNA damage in somatic cells originates from both environmental and endogenous sources, giving rise to mutations through multiple mechanisms. When these mutations affect the function of critical genes, cancer may ensue. Although identifying genomic subsets of mutated genes may inform therapeutic options, a systematic survey of tumor mutational spectra is required to improve our understanding of the underlying mechanisms of mutagenesis involved in cancer etiology. Recent studies have presented genome-wide sets of somatic mutations as a 96-element vector, a procedure that only captures the immediate neighbors of the mutated nucleotide. Herein, we present a 32 × 12 mutation matrix that captures the nucleotide pattern two nucleotides upstream and downstream of the mutation. A somatic autosomal mutation matrix (SAMM) was constructed from tumor-specific mutations derived from each of 909 individual cancer genomes harboring a total of 10,681,843 single-base substitutions. In addition, mechanistic template mutation matrices (MTMMs) representing oxidative DNA damage, ultraviolet-induced DNA damage, (5m)CpG deamination, and APOBEC-mediated cytosine mutation, are presented. MTMMs were mapped to the individual tumor SAMMs to determine the maximum contribution of each mutational mechanism to the overall mutation pattern. A Manhattan distance across all SAMM elements between any two tumor genomes was used to determine their relative distance. Employing this metric, 89.5% of all tumor genomes were found to have a nearest neighbor from the same tissue of origin. When a distance-dependent 6-nearest neighbor classifier was used, 10.4% of the SAMMs had an Undetermined tissue of origin, and 92.2% of the remaining SAMMs were assigned to the correct tissue of origin. [corrected]. Thus, although tumors from different tissues may have similar mutation patterns, their SAMMs often display signatures that are characteristic of specific tissues.

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.

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.

Prokaryotic nucleotide excision repair

Nucleotide excision repair (NER) has allowed bacteria to flourish in many different niches around the globe that inflict harsh environmental damage to their genetic material. NER is remarkable because of its diverse substrate repertoire, which differs greatly in chemical composition and structure. Recent advances in structural biology and single-molecule studies have given great insight into the structure and function of NER components. This ensemble of proteins orchestrates faithful removal of toxic DNA lesions through a multistep process. The damaged nucleotide is recognized by dynamic probing of the DNA structure that is then verified and marked for dual incisions followed by excision of the damage and surrounding nucleotides. The opposite DNA strand serves as a template for repair, which is completed after resynthesis and ligation.

C8-linked bulky guanosine DNA adducts: experimental and computational insights into adduct conformational preferences and resulting mutagenicity

Bulky DNA adducts are formed through the covalent attachment of aryl groups to the DNA nucleobases. Many of these adducts are known to possess conformational heterogeneity, which is responsible for the variety of mutagenic outcomes associated with these lesions. The present contribution reviews several conformational and mutagenic themes that are prevalent among the DNA adducts formed at the C8-site of the guanine nucleobase. The most important conclusions obtained (to date) from experiments are summarized including the anti/syn conformational preference of the adducts, their potential to inflict DNA mutations and mismatch stabilization, and their interactions with DNA polymerases and repair enzymes. Additionally, the unique role that computer calculations can play in understanding the structural properties of these adducts are highlighted.

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.

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.

Chemistry and structural biology of DNA damage and biological consequences

The formation of adducts by the reaction of chemicals with DNA is a critical step for the initiation of carcinogenesis. The structural analysis of various DNA adducts reveals that conformational and chemical rearrangements and interconversions are a common theme. Conformational changes are modulated both by the nature of adduct and the base sequences neighboring the lesion sites. Equilibria between conformational states may modulate both DNA repair and error-prone replication past these adducts. Likewise, chemical rearrangements of initially formed DNA adducts are also modulated both by the nature of adducts and the base sequences neighboring the lesion sites. In this review, we focus on DNA damage caused by a number of environmental and endogenous agents, and biological consequences.

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.

Probing for DNA damage with β-hairpins: similarities in incision efficiencies of bulky DNA adducts by prokaryotic and human nucleotide excision repair systems in vitro

Nucleotide excision repair (NER) is an important prokaryotic and eukaryotic defense mechanism that removes a large variety of structurally distinct lesions in cellular DNA. While the proteins involved are completely different, the mode of action of these two repair systems is similar, involving a cut-and-patch mechanism in which an oligonucleotide sequence containing the lesion is excised. The prokaryotic and eukaryotic NER damage-recognition factors have common structural features of β-hairpin intrusion between the two DNA strands at the site of the lesion. In the present study, we explored the hypothesis that this common β-hairpin intrusion motif is mirrored in parallel NER incision efficiencies in the two systems. We have utilized human HeLa cell extracts and the prokaryotic UvrABC proteins to determine their relative NER incision efficiencies. We report here comparisons of relative NER efficiencies with a set of stereoisomeric DNA lesions derived from metabolites of benzo[a]pyrene and equine estrogens in different sequence contexts, utilizing 21 samples. We found a general qualitative trend toward similar relative NER incision efficiencies for ∼65% of these substrates; the other cases deviate mostly by ∼30% or less from a perfect correlation, although several more distant outliers are also evident. This resemblance is consistent with the hypothesis that lesion recognition through β-hairpin insertion, a common feature of the two systems, is facilitated by local thermodynamic destabilization induced by the lesions in both cases. In the case of the UvrABC system, varying the nature of the UvrC endonuclease, while maintaining the same UvrA/B proteins, can markedly affect the relative incision efficiencies. These observations suggest that, in addition to recognition involving the initial modified duplexes, downstream events involving UvrC can also play a role in distinguishing and processing different lesions in prokaryotic NER.

Base sequence context effects on nucleotide excision repair

Nucleotide excision repair (NER) plays a critical role in maintaining the integrity of the genome when damaged by bulky DNA lesions, since inefficient repair can cause mutations and human diseases notably cancer. The structural properties of DNA lesions that determine their relative susceptibilities to NER are therefore of great interest. As a model system, we have investigated the major mutagenic lesion derived from the environmental carcinogen benzo[a]pyrene (B[a]P), 10S (+)-trans-anti-B[a]P-N²-dG in six different sequence contexts that differ in how the lesion is positioned in relation to nearby guanine amino groups. We have obtained molecular structural data by NMR and MD simulations, bending properties from gel electrophoresis studies, and NER data obtained from human HeLa cell extracts for our six investigated sequence contexts. This model system suggests that disturbed Watson-Crick base pairing is a better recognition signal than a flexible bend, and that these can act in concert to provide an enhanced signal. Steric hinderance between the minor groove-aligned lesion and nearby guanine amino groups determines the exact nature of the disturbances. Both nearest neighbor and more distant neighbor sequence contexts have an impact. Regardless of the exact distortions, we hypothesize that they provide a local thermodynamic destabilization signal for repair.

DNA-Destabilizing Agents as an Alternative Approach for Targeting DNA: Mechanisms of Action and Cellular Consequences

DNA targeting drugs represent a large proportion of the actual anticancer drug pharmacopeia, both in terms of drug brands and prescription volumes. Small DNA-interacting molecules share the ability of certain proteins to change the DNA helix's overall organization and geometrical orientation via tilt, roll, twist, slip, and flip effects. In this ocean of DNA-interacting compounds, most stabilize both DNA strands and very few display helix-destabilizing properties. These types of DNA-destabilizing effect are observed with certain mono- or bis-intercalators and DNA alkylating agents (some of which have been or are being developed as cancer drugs). The formation of locally destabilized DNA portions could interfere with protein/DNA recognition and potentially affect several crucial cellular processes, such as DNA repair, replication, and transcription. The present paper describes the molecular basis of DNA destabilization, the cellular impact on protein recognition, and DNA repair processes and the latter's relationships with antitumour efficacy.

Distant neighbor base sequence context effects in human nucleotide excision repair of a benzo[a]pyrene-derived DNA lesion

The effects of non-nearest base sequences, beyond the nucleotides flanking a DNA lesion on either side, on nucleotide excision repair (NER) in extracts from human cells were investigated. We constructed two duplexes containing the same minor groove-aligned 10S (+)-trans-anti-B[a]P-N(2)-dG (G*) DNA adduct, derived from the environmental carcinogen benzo[a]pyrene (B[a]P): 5'-C-C-A-T-C-G*-C-T-A-C-C-3' (CG*C-I), and 5'-C-A-C3-A4-C5-G*-C-A-C-A-C-3' (CG*C-II). We used polyacrylamide gel electrophoresis to compare the extent of DNA bending, and molecular dynamics simulations to analyze the structural characteristics of these two DNA duplexes. The NER efficiencies are 1.6(+/-0.2)-fold greater in the case of the CG*C-II than the CG*C-I sequence context in 135-mer duplexes. Gel electrophoresis and self-ligation circularization experiments revealed that the CG*C-II duplex is more bent than the CG*C-I duplex, while molecular dynamics simulations showed that the unique -C3-A4-C5- segment in the CG*C-II duplex plays a key role. The presence of a minor groove-positioned guanine amino group, the Watson-Crick partner to C3, acts as a wedge; facilitated by a highly deformable local -C3-A4- base step, this amino group allows the B[a]P ring system to produce a more enlarged minor groove in CG*C-II than in CG*C-I, as well as a local untwisting and enlarged and flexible Roll only in the CG*C-II sequence. These structural properties fit well with our earlier findings that in the case of the family of minor groove 10S (+)-trans-anti-B[a]P-N(2)-dG lesions, flexible bends and enlarged minor groove widths constitute NER recognition signals, and extend our understanding of sequence context effects on NER to the neighbors that are distant to the lesion.

DNA interstrand crosslink repair in mammalian cells: step by step

Interstrand DNA crosslinks (ICLs) are formed by natural products of metabolism and by chemotherapeutic reagents. Work in E. coli identified a two cycle repair scheme involving incisions on one strand on either side of the ICL (unhooking) producing a gapped intermediate with the incised oligonucleotide attached to the intact strand. The gap is filled by recombinational repair or lesion bypass synthesis. The remaining monoadduct is then removed by nucleotide excision repair (NER). Despite considerable effort, our understanding of each step in mammalian cells is still quite limited. In part this reflects the variety of crosslinking compounds, each with distinct structural features, used by different investigators. Also, multiple repair pathways are involved, variably operative during the cell cycle. G(1) phase repair requires functions from NER, although the mechanism of recognition has not been determined. Repair can be initiated by encounters with the transcriptional apparatus, or a replication fork. In the case of the latter, the reconstruction of a replication fork, stalled or broken by collision with an ICL, adds to the complexity of the repair process. The enzymology of unhooking, the identity of the lesion bypass polymerases required to fill the first repair gap, and the functions involved in the second repair cycle are all subjects of active inquiry. Here we will review current understanding of each step in ICL repair in mammalian cells.

NMR and computational studies of stereoisomeric equine estrogen-derived DNA cytidine adducts in oligonucleotide duplexes: opposite orientations of diastereomeric forms

The equine estrogens equilin (EQ) and equilenin (EN) are the active components in the widely prescribed hormone replacement therapy formulation Premarin. Metabolic activation of EQ and EN generates the catechol 4-hydroxyequilenin (4-OHEN) that autoxidizes to the reactive o-quinone form in aerated aqueous solutions. The o-quinones react predominantly with C, and to a lesser extent with A and G, to form premutagenic cyclic covalent DNA adducts in vitro and in vivo. To obtain insights into the structural properties of these biologically important DNA lesions, we have synthesized site-specifically modified oligonucleotides containing the stereoisomeric 1'S,2'R,3'R-4-OHEN-C3 and 1'R,2'S,3'S-4-OHEN-C4 adducts derived from the reaction of 4-OHEN with the C in the oligonucleotide 5'-GGTAGCGATGG in aqueous solution. A combined NMR and computational approach was utilized to determine the conformational characteristics of the two major 4-OHEN-C3 and 4-OHEN-C4 stereoisomeric adducts formed in this oligonucleotide hybridized with its complementary strand. In both cases, the modified C adopts an anti glycosidic bond conformation; the equilenin distal ring protrudes into the minor groove while its two proximal hydroxyl groups are exposed on the major groove side of the DNA duplex. The bulky 4-OHEN-C adduct distorts the duplex within the central GC*G portion, but Watson-Crick pairing is maintained adjacent to C* in both stereoisomeric adducts. For the 4-OHEN-C3 adduct, the equilenin rings are oriented toward the 5'-end of the modified strand, while in 4-OHEN-C4 the equilenin is 3'-directed. Correspondingly, the distortions of the double-helical structures are more pronounced on the 5'- or the 3'-side of the lesion, respectively. These differences in stereoisomeric adduct conformations may play a role in the processing of these lesions in cellular environments.

The sequence dependence of human nucleotide excision repair efficiencies of benzo[a]pyrene-derived DNA lesions: insights into the structural factors that favor dual incisions

Nucleotide excision repair (NER) is a vital cellular defense system against carcinogen-DNA adducts, which, if not repaired, can initiate cancer development. The structural features of bulky DNA lesions that account for differences in NER efficiencies in mammalian cells are not well understood. In vivo, the predominant DNA adduct derived from metabolically activated benzo[a]pyrene (BP), a prominent environmental carcinogen, is the 10S (+)-trans-anti-[BP]-N(2)-dG adduct (G*), which resides in the B-DNA minor groove 5'-oriented along the modified strand. We have compared the structural distortions in double-stranded DNA, imposed by this adduct, in the different sequence contexts 5'-...CGG*C..., 5'-...CG*GC..., 5'-...CIG*C... (I is 2'-deoxyinosine), and 5'-...CG*C.... On the basis of electrophoretic mobilities, all duplexes manifest moderate bends, except the 5'-...CGG*C...duplex, which exhibits an anomalous, slow mobility attributed to a pronounced flexible kink at the site of the lesion. This kink, resulting from steric hindrance between the 5'-flanking guanine amino group and the BP aromatic rings, both positioned in the minor groove, is abolished in the 5'-...CIG*C...duplex (the 2'-deoxyinosine group, I, lacks this amino group). In contrast, the sequence-isomeric 5'-...CG*GC...duplex exhibits only a moderate bend, but displays a remarkably increased opening rate at the 5'-flanking base pair of G*, indicating a significant destabilization of Watson-Crick hydrogen bonding. The NER dual incision product yields were compared for these different sequences embedded in otherwise identical 135-mer duplexes in cell-free human HeLa extracts. The yields of excision products varied by a factor of as much as approximately 4 in the order 5'-...CG*GC...>5'...CGG*C...>or=5'...CIG*C...>or=5'-...CG*C.... Overall, destabilized Watson-Crick hydrogen bonding, manifested in the 5'-...CG*GC...duplex, elicits the most significant NER response, while the flexible kink displayed in the sequence-isomeric 5'-...CGG*C...duplex represents a less significant signal in this series of substrates. These results demonstrate that the identical lesion can be repaired with markedly variable efficiency in different local sequence contexts that differentially alter the structural features of the DNA duplex around the lesion site.