Dynamics of a benzo[a]pyrene-derived guanine DNA lesion in TGT and CGC sequence contexts: enhanced mobility in TGT explains conformational heterogeneity, flexible bending, and greater susceptibility to nucleotide excision repair

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.

"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.

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.

Molecular Dynamics Study of One-Base Deletion Duplexes Containing the Major DNA Adduct Formed by Ochratoxin A: Effects of Sequence Context and Adduct Ionization State on Lesion Site Structure and Mutagenicity

Ochratoxin A (OTA) is a ubiquitous food toxin associated with chronic nephropathy in humans and renal carcinogenicity in rodents. The mutational spectra of cells exposed to OTA reveal that one-base deletions comprise the largest percentage (73%) of the total mutations that occur upon OTA exposure. To contribute toward understanding the prevalence of OTA-induced one-base deletion mutations, the present work uses molecular dynamics (MD) simulations to analyze the conformational preferences of one-base deletion duplexes containing OT-G, the major OTA adduct (addition product) at the C⁸-site of guanine. Specifically, the influence of OT-G in four possible ionization states and three sequence contexts (G¹, G² and G³ in the NarI (5'-G¹G²CG³CC-3'), a prokaryotic mutational hotspot sequence) on the structure of the adducted DNA is investigated. Our data reveal that the damaged helices are stable in two (B-type (B) and stacked (S)) conformations that are structurally similar to those adopted by common N-linked C⁸-guanine lesions. However, the adduct ionization state and sequence context affect the degree of helical distortion and the B/S conformational heterogeneity, which will impact the lesion repair and replication outcomes. This finding correlates with the experimentally reported tissue-specific mutagenicity of OTA exposure. Furthermore, regardless of the adduct conformation, ionization state, or sequence context, more stable lesion-site interactions and lack of disruption of the flanking base pairs in the one-base deletion duplexes compared to the corresponding two-base deletion helices rationalize the greater abundance of OTA induced one-base deletions. Overall, our work provides valuable structural insights that help explain the experimentally observed mutagenicity associated with OTA.

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.

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.

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.

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.

BPDE-induced genotoxicity: relationship between DNA adducts, mutagenicity in the in vitro PIG-A assay, and the transcriptional response to DNA damage in TK6 cells

Benzo[a]pyrene is a known human carcinogen. As underlying mechanism, the induction of stable DNA adducts and mutations have been repeatedly demonstrated. Also, the activation of cellular stress response on the transcriptional level has been described. Nevertheless, the interrelationship between these different events is less well understood, especially at low, for human exposure relevant concentrations. Within the present study, we applied the reactive metabolite benzo[a]pyrene diolepoxide (BPDE) in the nanomolar, non-cytotoxic concentration range in human TK6 cells and quantified the induction and repair of stable DNA adducts at the N²-position of guanine by HPLC with fluorescence detection. Significant levels of DNA lesions were detected even at the lowest concentration of 10 nM BPDE, with a linear increase up to 50 nM. Relative repair was similar at all damage levels, reaching about 30% after 8 h and 60% after 24 h. Mutation frequencies were quantified as GPI-deficient cells by the recently established in vitro PIG-A mutagenicity assay. Again, a linear dose–response-relationship in the before-mentioned concentration range was observed, also when plotting the number of GPI-deficient cells against the number of DNA adducts. Furthermore, we explored the time- and concentration-dependent DNA damage response on the transcriptional level via a high-throughput RT-qPCR technique by quantifying the impact of BPDE on the transcription of 95 genes comprising DNA damage response, DNA repair factors, oxidative stress response, cell cycle arrest, cell proliferation, and apoptosis. As expected, BPDE activated DNA damage signaling, p53 and AP-1 dependent signaling, oxidative stress response, and apoptosis. However, in contrast to DNA adducts and mutations, the onset of the transcriptional DNA damage response was restricted to higher concentrations, indicating that its respective activations require a certain level of DNA lesions. Altogether, the results indicate that in case of BPDE, DNA lesions and mutations were correlated at all concentrations, suggesting that repair is not complete even at low levels of DNA damage. Considering the ongoing discussion on potential thresholds also for genotoxic carcinogens, the results are of major relevance, both with respect to basic research as well as to risk assessment of chemical carcinogens.

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.

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

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

Molecular dynamics of dibenz[a,h]anthracene and its metabolite interacting with lung surfactant phospholipid bilayers

Dibenz[a,h]anthracene and its metabolite may form aggregates, which have implications in the clearance process of the lung surfactant phospholipid bilayers.

Ribonucleotides as nucleotide excision repair substrates

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

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.

Electrochemical study on the effects of epigenetic cytosine methylation on anti-benzo[a]pyrene diol epoxide damage at TP53 oligomers

Anti-benzo[a]pyrene-r-7,t-8-dihydrodiol-t-9,10-epoxide (anti-BPDE) is a known carcinogen that damages DNA, and this damage is influenced by the DNA sequence and epigenetic factors. The influence of epigenetic cytosine methylation on the reaction with anti-BPDE at a known hotspot DNA damage site was studied electrochemically. Gold electrodes were modified with thiolated DNA oligomers spanning codons 270-276 of the TP53 gene. The oligomers exhibited 5-carbon cytosine methylation at the codon 273 location on the bound probe, the acquired complementary target, or both. Redox active diviologen compounds of the form C(12)H(25)V(2+)C(6)H(12)V(2+)C(12)H(25) (V(2+) = 4,4'-bipyridyl or viologen, C12-Viologen) were employed to detect anti-BPDE damage to DNA. DNA was exposed to racemic (±)- or enantiomerically pure (+)-anti-BPDE solutions followed by electrochemical interrogation in the presence of C12-Viologen. Background subtracted square wave voltammograms (SWV) showed the appearance of two peaks at approximately -0.38 V and -0.55 V vs Ag/AgCl upon anti-BPDE exposure. The acquired voltammetry is consistent with singly reduced C12-Viologen dimers bound at two different DNA environments, which arise from BPDE damage and are influenced by cytosine methylation and BPDE stereochemical considerations. UV spectroscopic and mass spectrometric methods employed to validate the electrochemical responses showed that (+)-anti-BPDE primarily adopts a minor groove bound orientation within the oligomers while selectively targeting the nontranscribed ssDNA sequence within the duplexes.

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.

Frontiers in molecular dynamics simulations of DNA

It has been known for decades that DNA is extremely flexible and polymorphic, but our knowledge of its accessible conformational space remains limited. Structural data, primarily from X-ray diffraction studies, is sparse in comparison to the manifold configurations possible, and direct experimental examinations of DNA's flexibility still suffer from many limitations. In the face of these shortcomings, molecular dynamics (MD) is now an essential tool in the study of DNA. It affords detailed structural and dynamical insights, which explains its recent transition from a small number of highly specialized laboratories to a large variety of groups dealing with challenging biological problems. MD is now making an irreversible journey to the mainstream of research in biology, with the attendant opportunities and challenges. But given the speed with which MD studies of DNA have spread, the roots remain somewhat shallow: in many cases, there is a lack of deep knowledge about the foundations, strengths, and limits of the technique. In this Account, we discuss how MD has become the most important source of structural and flexibility data on DNA, focusing on advances since 2007 of atomistic MD in the description of DNA under near-physiological conditions and highlighting the possibilities and shortcomings of the technique. The evolution in the field over the past four years is a prelude to the ongoing revolution. The technique has gained in robustness and predictive power, which when coupled with the spectacular improvements in software and hardware has enabled the tackling of systems of increasing complexity. Simulation times of microseconds have now been achieved, with even longer times when specialized hardware is used. As a result, we have seen the first real-time simulation of large conformational transitions, including folding and unfolding of short DNA duplexes. Noteworthy advances have also been made in the study of DNA-ligand interactions, and we predict that a global thermodynamic and kinetic picture of the binding landscape of DNA will become available in a few years. MD will become a crucial tool in areas such as biomolecular engineering and synthetic biology. MD has also been shown to be an excellent source of parameters for mesoscopic models of DNA flexibility. Such models can be refined through atomistic MD simulations on small duplexes and then applied to the study of entire chromosomes. Recent evidence suggests that MD-derived elastic models can successfully predict the position of regulatory regions in DNA and can help advance our understanding of nucleosome positioning and chromatin plasticity. If these results are confirmed, MD simulations can become the ultimate tool to decipher a physical code that can contribute to gene regulation. We are entering the golden age of MD simulations of DNA. Undoubtedly, the expectations are high, but the challenges are also enormous. These include the need for more accurate potential energy functionals and for longer and more complex simulations in more realistic systems. The joint research effort of several groups will be crucial for adapting the technique to the requirements of the coming decade.

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.

Conformational flexibility of C8-phenoxylguanine adducts in deoxydinucleoside monophosphates

M06-2X/6-31G(d,p) is used to calculate the structure of all natural deoxydinucleoside monophosphates with G in the 5' or 3' position, the anti or syn conformation, and each natural (A, C, G, T) base in the corresponding flanking position. When the ortho or para C8-phenoxyl-2'-deoxyguanosine (C8-phenoxyl-dG) adduct replaces G in each model, there is little change in the relative base-base orientation or backbone conformation. However, the orientation of the C8-phenoxyl group can be characterized according to the position (5' versus 3'), conformation (anti versus syn), and isomer (ortho versus para) of damage. Although the degree of coplanarity between the phenoxyl ring and G base in the ortho adduct is highly affected by the sequence since the hydroxyl group can interact with neighboring bases, the para adduct generally does not exhibit discrete interactions with flanking bases. For both adducts, steric clashes between the phenoxyl group and the backbone or flanking base destabilize the anti conformation preferred by the natural nucleotide and thereby result in a clear preference for the syn conformation regardless of the sequence or position. This contrasts the conclusions drawn from smaller (nucleoside, nucleotide) models previously used in the literature, which stresses the importance of using models that address the steric constraints present due to the surrounding environment. Since replication errors for other C8-dG bulky adducts have been linked to a preference for the syn conformation, our findings provide insight into the possible mutagenicity of phenolic adducts.

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.

Conformational and thermodynamic impact of bulky aminofluorene adduction on simulated translesion DNA synthesis

We report a systematic spectroscopic investigation on the conformational evolution during primer extension of a bulky fluoroaminofluorene-modified dG adduct (FAF-dG) in chemically simulated translesion synthesis. FAF-dG was paired either with dC or dA (dC-match and dA-mismatch series, respectively). Dynamic (19)F NMR/CD results showed that the FAF-adduct exists in a syn/anti equilibrium and that its conformational characteristics are modulated by the identity of an inserted nucleotide at the lesion site and the extent of primer elongation. At the pre-insertion site, the adduct adopted preferentially a syn conformation where FAF stacked with preceding bases. Insertion of the correct nucleotide dC at the lesion site and subsequent elongation resulted in a gradual transition to the anti conformation. By contrast, the syn conformer was persistent along with primer extension in the dA-mismatch series. In the dC-match series, FAF-induced thermal (T(m)) and thermodynamic (-ΔG°(37 °C)) stabilities were significantly reduced relative to those of the controls. However, the corresponding T(m) and -ΔG°(37 °C) values were increased in the FAF-modified mismatched dA series. The lesion impact persisted up to three 5'-nucleotides from the lesion. Occupation of the minor groove of the W-conformer with the bulky carcinogenic fluorene moiety not only would limit the DNA mobility but also would impose a serious difficulty for the active site of a polymerase throughout the replication process. Our spectroscopic results are consistent with reported data on AF, which showed dramatic (~10(4)-fold) differences in the nucleotide insertion rates between the dC-match and dA-mismatch series. The results emphasize the importance of adduct-induced steric constraints for determining the replication fidelity of a polymerase.

Modeling the interactions of the nucleotide excision repair UvrA(2) dimer with DNA

The UvrA protein initiates the DNA damage recognition process by the bacterial nucleotide excision repair (NER) system. Recently, crystallographic structures of holo-UvrA(2) dimers from two different microorganisms have been released (Protein Data Bank entries 2r6f , 2vf7 , and 2vf8 ). However, the details of the DNA binding by UvrA(2) and other peculiarities involved in the damage recognition process remain unknown. We have undertaken a molecular modeling approach to appraise the possible modes of DNA-UvrA(2) interaction using molecular docking and short-scale guided molecular dynamics [continuum field, constrained, and/or unrestricted simulated annealing (SA)], taking into account the three-dimensional location of a series of mutation-identified UvrA residues implicated in DNA binding. The molecular docking was based on the assumptions that the UvrA(2) dimer is preformed prior to DNA binding and that no major protein conformational rearrangements, except moderate domain reorientations, are required for binding of undamaged DNA. As a first approximation, DNA was treated as a rigid ligand. From the electrostatic relief of the ventral surface of UvrA(2), we initially identified three, noncollinear DNA binding paths. Each of the three resulting nucleoprotein complexes (C1, C2, and C3) was analyzed separately, including calculation of binding energies, the number and type of interaction residues (including mutated ones), and the predominant mode of translational and rotational motion of specific protein domains after SA to ensure improved DNA binding. The UvrA(2) dimer can accommodate DNA in all three orientations, albeit with different binding strengths. One of the UvrA(2)-DNA complexes (C1) fulfilled most of the requirements (high interaction energy, proximity of DNA to mutated residues, etc.) expected for a natural, high-affinity DNA binding site. This nucleoprotein presents a structural organization that is designed to clamp and bend double-stranded DNA. We examined the binding site in more detail by docking DNAs of significantly different (AT- vs CG-enriched) sequences and by submitting the complexes to DNA-unrestricted SA. It was found that in a manner independent of the DNA sequence and applied MD protocols, UvrA(2) favors binding of a bent and unwound undamaged DNA, with a kink positioned in the proximity of the Zn3 hairpins, anticollinearly aligned at the bottom of the ventral protein surface. It is further hypothesized that the Zn3 modules play an essential role in the damage recognition process and that the apparent existence of a family of DNA binding sites might be biologically relevant. Our data should prove to be useful in rational (structure-based) mutation studies.

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.

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.

Enthalpy-entropy contribution to carcinogen-induced DNA conformational heterogeneity

DNA damage by adduct formation is a critical step for the initiation of carcinogenesis. Aromatic amines are strong inducers of environmental carcinogenesis. Their DNA adducts are known to exist in an equilibrium between the major groove (B) and base-displaced stacked (S) conformations. However, the factors governing such heterogeneity remain unclear. Here we conducted extensive calorimetry/NMR/CD studies on the model DNA lesions caused by fluorinated 2-aminfluorene (FAF) and 4-aminobiphenyl (FABP) in order to gain thermodynamic and kinetic insights into the S/B conformational equilibrium. We demonstrate that there are large differences in enthalpy-entropy compensations for FABP and FAF. The small and flexible FABP exclusively adopts the less perturbed B conformer with small enthalpy (DeltaDeltaH-2.7 kcal/mol)/entropy (DeltaDeltaS-0.7 eu) change. In contrast, FAF stacks better and exists as a mixture of B and S conformers with large enthalpy (DeltaDeltaH-13.4 kcal/mol)/entropy (DeltaDeltaS-34.2 eu) compensation. van't Hoff analysis of dynamic (19)F NMR data indicated DeltaH(B<-->S) = 4.1 kcal/mol in favor of the B conformer and DeltaS(B<-->S) = 15.6 cal mol(-1) K(-1) in favor of the intercalated S conformer. These findings demonstrate that the favorable entropy of the S conformer over B conformer determines the S/B population ratios at physiological temperatures.

Influence of flanking sequence context on the conformational flexibility of aminofluorene-modified dG adduct in dA mismatch DNA duplexes

When positioned opposite to a dA in a DNA duplex, the prototype arylamine-DNA adduct [N-(2'-deoxyguanosin-yl)-7-fluoro-2-aminofluorene (FAF)] adopts the so-called 'wedge' (W) conformation, in which the carcinogen resides in the minor groove of the duplex. All 16 FAF-modified 12-mer NG*N/NAN dA mismatch duplexes (G* = FAF, N = G, A, C, T) exhibited strongly positive induced circular dichroism in the 290-360 nm range (ICD(290-360 nm)), which supports the W conformation. The ICD(290-360 nm) intensities were the greatest for duplexes with a 3'-flanking T. The AG*N duplex series showed little adduct-induced destabilization. An exception was the AG*T duplex, which displayed two well-resolved signals in the (19)F NMR spectra. This was presumably due to a strong lesion-destabilizing effect of the 3'-T. The flanking T effect was substantiated further by findings with the TG*T duplex, which exhibited greater lesion flexibility and nucleotide excision repair recognition. Adduct conformational heterogeneity decreased in order of TG*T > AG*T > CG*T > AG*A > AG*G > AG*C. The dramatic flanking T effect on W-conformeric duplexes is consistent with the strong dependence of the ICD(290-360) on both temperature and salt concentration and could be extended to the arylamine food mutagens that are biologically relevant in humans.

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.

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

The structural origin underlying differential nucleotide excision repair (NER) susceptibilities of bulky DNA lesions remains a challenging problem. We investigated the 10S (+)-trans-anti-[BP]-N(2)-2'-deoxyguanosine (G*) adduct in double-stranded DNA. This adduct arises from the reaction, in vitro and in vivo, of a major genotoxic metabolite of benzo[a]pyrene (BP), (+)-(7R,8S,9S,10R)-7,8-dihydroxy-9,10-epoxy-7,8,9,10-tetrahydrobenzo[a]pyrene, with the exocyclic amino group of guanine. Removal of this lesion by the NER apparatus in cell-free extracts has been found to depend on the base sequence context in which the lesion is embedded, providing an excellent opportunity for elucidating the properties of the damaged DNA duplexes that favor NER. While the BP ring system is in the B-DNA minor groove, 5' directed along the modified strand, there are orientational distinctions that are sequence dependent and are governed by flanking amino groups [Nucleic Acids Res.35 (2007), 1555-1568]. To elucidate sequence-governed NER susceptibility, we conducted molecular dynamics simulations for the 5'-...CG*GC..., 5'-...CGG*C..., and 5'-...TCG*CT... adduct-containing duplexes. We also investigated the 5'-...CG*IC... and 5'-...CIG*C... sequences, which contain "I" (2'-deoxyinosine), with hydrogen replacing the amino group in 2'-deoxyguanosine, to further characterize the structural and dynamic roles of the flanking amino groups in the damaged duplexes. Our results pinpoint explicit roles for the amino groups in tandem GG sequences on the efficiency of NER and suggest a hierarchy of destabilizing structural features that differentially facilitate NER of the BP lesion in the sequence contexts investigated. Furthermore, combinations of several locally destabilizing features in the hierarchy, consistent with a multipartite model, may provide a relatively strong recognition signal.