Modeling and molecular mechanical studies of the cis-thymine glycol radiation damage lesion in DNA

Computer graphics and energy minimization techniques were used to construct a model of DNA containing cis-thymine glycol, an oxidation product of thymine formed in DNA by ionizing radiation. The model simulated an experimental DNA substrate used to study the effects of this lesion on DNA synthesis in vitro. The results derived from the model indicate that cis-thymine glycol lesions introduce localized perturbations of DNA structure. Specifically the model shows that interactions with the neighboring base pair on the 5' side are significantly destabilized by thymine glycol whereas interactions with the 3' base pair are stabilized by the lesion. The magnitude of these effects is modulated by the nucleotide sequence around the lesion, particularly by the nature of the base on the 3' side. The base pair formed between adenine and thymine glycol is energetically stable and shows minimal distortion, suggesting that this lesion retains the ability to direct the insertion of the correct nucleotide during DNA synthesis.     

Sequence-dependent distortions induced in DNA by monofunctional platinum(II) binding.

The effects on thermal stability and conformation of DNA produced by the monofunctional adducts of chlorodiethylenetriamineplatinum(II) chloride ([Pt(dien)Cl]Cl) have been investigated. Oligodeoxyribonucleotide duplexes of varying lengths (9-20 base pairs) and of varying central trinucleotide sequences were prepared and characterized that contained site-specific and unique N(7)-guanine adducts. Included are adducts at the sequences of d(AGC), d(AGT), d(CGA), d(TGA), d(TGC), and d(TGT). All these monofunctional adducts decrease the melting temperature (Tm) of the duplexes. This destabilization effect exhibits a sequence-dependent variability. The highest lowering of Tm is observed for the modified duplexes containing the central sequence of pyrimidine-guanine-pyrimidine. The destabilization effect is reduced with decreasing concentrations of Na+. Polarography, circular dichroism, phenanthroline-copper, and chemical probes reveal conformational distortions spreading over several base pairs around the adduct. The effects of monofunctional platinum(II) adducts on conformational distortions in DNA exhibit a sequence-dependent variability similar to those on thermal stability of DNA. The influence of the monofunctional adduct formed by cis-diamminemonoaquamonochloroplatinum(II) on the stability of the oligonucleotide duplex has been also studied. This lesion decreases thermal stability of DNA in the same way as does the adduct of [Pt(dien)Cl]Cl.     

DNA binding and mutation spectra of the carcinogen N-2-aminofluorene in Escherichia coli. A correlation between the conformation of the premutagenic lesion and the mutation specificity

When the chemical carcinogen N-2-acetylaminofluorene binds to DNA in vivo, two major adducts are formed, both at position C-8 of the guanine residue. One of these (the acetylaminofluorene adduct) retains the acetyl group, while the other (the aminofluorene adduct) is the corresponding deacetylated form. Unlike -AAF adducts, which trigger important structural changes of the DNA secondary structure (either the insertion-denaturation model or the induction of a Z-DNA structure, depending upon the local nucleotide sequence), -AF adducts bind to the C-8 of guanine residues without causing any major conformational change of the B-DNA structure. Well-defined adducts (either -AF or -AAF) can be formed in vitro by reacting DNA with either N-hydroxy-N-2-aminofluorene or N-acetoxy-N-2-acetylaminofluorene. Specific cleavage of the phosphodiester backbone at -AF adducts can be achieved by treating -AF-modified DNA in 1 M-piperidine at 90 degrees C. This observation led us to construct the spectrum for -AF binding to a defined DNA restriction fragment. It is found that only guanine residues react to form alkali-labile lesions and that the reactivity among the different guanines is similar. In a forward mutation assay, namely the inactivation of the tetracycline resistance gene, we found previously that more than 90% of mutations induced by -AAF adducts are frameshift mutations. Using the same assay, we show here that -AF adducts induce primarily base substitution mutations (85%), mainly of the G to T transversion type. There is therefore a strong correlation between the nature of the carcinogen-induced conformational change of the DNA structure and the corresponding mutation specificity. The -AF-induced base substitution mutations depend upon the umuC gene function(s). The data obtained in our forward mutation assay are compared to the data previously obtained in the histidine reversion assay (Ames test).     

DNA oligonucleotides with A, T, G or C opposite an abasic site: structure and dynamics

Abasic sites are common DNA lesions resulting from spontaneous depurination and excision of damaged nucleobases by DNA repair enzymes. However, the influence of the local sequence context on the structure of the abasic site and ultimately, its recognition and repair, remains elusive. In the present study, duplex DNAs with three different bases (G, C or T) opposite an abasic site have been synthesized in the same sequence context (5′-CCA AAG₆ XA₈C CGG G-3′, where X denotes the abasic site) and characterized by 2D NMR spectroscopy. Studies on a duplex DNA with an A opposite the abasic site in the same sequence has recently been reported [Chen,J., Dupradeau,F.-Y., Case,D.A., Turner,C.J. and Stubbe,J. (2007) Nuclear magnetic resonance structural studies and molecular modeling of duplex DNA containing normal and 4′-oxidized abasic sites. Biochemistry, 46, 3096–3107]. Molecular modeling based on NMR-derived distance and dihedral angle restraints and molecular dynamics calculations have been applied to determine structural models and conformational flexibility of each duplex. The results indicate that all four duplexes adopt an overall B-form conformation with each unpaired base stacked between adjacent bases intrahelically. The conformation around the abasic site is more perturbed when the base opposite to the lesion is a pyrimidine (C or T) than a purine (G or A). In both the former cases, the neighboring base pairs (G6-C21 and A8-T19) are closer to each other than those in B-form DNA. Molecular dynamics simulations reveal that transient H-bond interactions between the unpaired pyrimidine (C20 or T20) and the base 3′ to the abasic site play an important role in perturbing the local conformation. These results provide structural insight into the dynamics of abasic sites that are intrinsically modulated by the bases opposite the abasic site.     

Lesion search and recognition by thymine DNA glycosylase revealed by single molecule imaging

The ability of DNA glycosylases to rapidly and efficiently detect lesions among a vast excess of nondamaged DNA bases is vitally important in base excision repair (BER). Here, we use single molecule imaging by atomic force microscopy (AFM) supported by a 2-aminopurine fluorescence base flipping assay to study damage search by human thymine DNA glycosylase (hTDG), which initiates BER of mutagenic and cytotoxic G:T and G:U mispairs in DNA. Our data reveal an equilibrium between two conformational states of hTDG–DNA complexes, assigned as search complex (SC) and interrogation complex (IC), both at target lesions and undamaged DNA sites. Notably, for both hTDG and a second glycosylase, hOGG1, which recognizes structurally different 8-oxoguanine lesions, the conformation of the DNA in the SC mirrors innate structural properties of their respective target sites. In the IC, the DNA is sharply bent, as seen in crystal structures of hTDG lesion recognition complexes, which likely supports the base flipping required for lesion identification. Our results support a potentially general concept of sculpting of glycosylases to their targets, allowing them to exploit the energetic cost of DNA bending for initial lesion sensing, coupled with continuous (extrahelical) base interrogation during lesion search by DNA glycosylases.     

The p-benzoquinone DNA adducts derived from benzene are highly mutagenic

Benzene is a human leukemia carcinogen, resulting from its cellular metabolism. A major benzene metabolite is p-benzoquinone (pBQ), which can damage DNA by forming the exocyclic base adducts pBQ-dC, pBQ-dA, and pBQ-dG in vitro. To gain insights into the role of pBQ in benzene genotoxicity, we examined in vitro translesion synthesis and in vivo mutagenesis of these pBQ adducts. Purified REV1 and Polkappa were essentially incapable of translesion synthesis in response to the pBQ adducts. Opposite pBQ-dA and pBQ-dC, purified human Poliota was capable of error-prone nucleotide insertion, but was unable to perform extension synthesis. Error-prone translesion synthesis was observed with Poleta. However, DNA synthesis largely stopped opposite the lesion. Consistent with in vitro results, replication of site-specifically damaged plasmids was strongly inhibited by pBQ adducts in yeast cells, which depended on both Polzeta and Poleta. In wild-type cells, the majority of translesion products were deletions at the site of damage, accounting for 91%, 90%, and 76% for pBQ-dA, pBQ-dG, and pBQ-dC, respectively. These results show that the pBQ-dC, pBQ-dA, and pBQ-dG adducts are strong blocking lesions, and are highly mutagenic by predominantly inducing deletion mutations. These results are consistent with the lesion structures predicted by molecular dynamics simulation. Our results led to the following model. Translesion synthesis normally occurs by directly copying the lesion site through base insertion and extension synthesis. When the lesion becomes incompatible in accommodating a base opposite the lesion in DNA, translesion synthesis occurs by a less efficient lesion loop-out mechanism, resulting in avoiding copying the damaged base and leading to deletion.     

Extent of helix perturbation associated with DNA modification by the o-acetyl derivative of the carcinogen 4-hydroxyaminoquinoline-1-oxide

Duplex unwinding associated with DNA modification by 4-acetoxyaminoquinoline-1-oxide, a model ultimate carcinogen of 4-nitroquinoline-1-oxide, has been determined by the agarose gel electrophoresis band-shift method. An average unwinding angle per stable adduct of -15.1 degrees +/- 1.5 degrees for negatively supercoiled topoisomers and -6.5 degrees +/- 1.4 degrees for positively supercoiled topoisomers was obtained. Because of the different proportion of stable adducts (dGuo-N2-AQO, dGuo-C8-AQO, dAdo-N6-AQO) between negatively (8:1.5:0.5) and positively (5:2.5:1) supercoiled topoisomers, the difference in unwinding angles is suggestive of a diverse contribution of the various adducts to the overall conformational change. Since the largest unwinding angle was coupled with the highest proportion of dGuo-N2-AQO adduct, it is likely that this adduct is the most distortive lesion. A contribution of sites of base loss to DNA unwinding was also observed.     

Energy Minimized Structures of Carcinogen-DNA. Adducts: 2-Acetylaminofluorene and 2-Aminofluorene

Abstract

Energy minimized structures of DNA modified by the aromatic amines 2-acetylaminofluorene (AAF) and 2-aminofuorene (AF), for which no experimental atomic resolution data exist, are presented. These have been computed with a new molecular mechanics program specifically designed to define distortions imposed by such adducts, and employing a rational strategy for searching the conformation space of a DNA molecule with covalently linked carcinogen. In alternating G-C sequences, the AAF adduct prefers to reside at the exterior of an undeformed Z-helix. It can also induce base displacement with attendant denaturation and helix bending in sequences that disfavor the Z form, but undeformed B helices are excluded. The AF adduct, by contrast, prefers the major groove of an unperturbed B-helix, but can also induce carcinogen-base stacking in single stranded regions of the DNA, such as at the replication fork. The different biological properties of these two adducts may be related to their distinct conformational features.     

Theoretical analysis of competitive conformational transitions in torsionally stressed DNA

This paper presents a theoretical analysis of the conformation of a torsionally deformed segment of DNA containing two sites susceptible to stress-induced transitions in secondary structure. A mechanical analysis of the ensuing competitive behavior is developed and applied to several phenomena of possible biological relevance. First, a molecular lesion which disrupts base pairing without strand breakage (such as a pyrimidine dimer) is shown to provide an effective nucleation site for further stress-induced denaturation. A competition is established between strand separation at this lesion site and at the A + T-richest portion of the molecule. The relative importance of these two forms of melting is shown to depend upon the A + T-content of the sites involved, segment length, local environmental conditions and the magnitude of the imposed torsional deformation. A possible alternative mode of behavior of a stressed segment of DNA involves transitions from B-form to Z-form. The second application of this theory analyzes the interplay between B → Z transitions and local denaturation in torsionally stressed DNA. Finally, local melting is shown to be energetically preferred over transitions to A-form under physiologically reasonable conditions in vitro, due primarily to the greater degree of unwinding involved in melting.The mechanical theory presented here makes several simplifying assumptions regarding the nature of the transitions and the sequences involved. First, the theory is developed explicitly for the competition between two sites of possible transition, with no further consideration given to conformational degeneracy or sequence effects. These sites are regarded as being uniform in composition. A multistate, heteropolymeric statistical mechanical transition theory is required to account rigorously for degeneracy and the influence of base sequence. A preliminary formulation of such a theory is used to analyze the denaturation of a segment containing a site of disrupted base pairing.     

Recognition of a cytosine base lesion by a human damage-specific DNA binding protein.

Sodium bisulfite reacts with cytosine and 5-methylcytosine, forming the 5,6-dihydrosulfonate adducts which deaminate to the uracil and thymine adducts, respectively. At alkaline pH, the sulfonate groups are then released, generating uracil and thymine. In DNA, the resulting G:U and G:T base mismatches generated are potential sites of mutagenesis. Using a human damage-specific DNA binding protein as a probe, we have found protein-recognizable lesions in bisulfite-treated DNA and poly d(I-C), but not in treated poly d(A-T) or poly d(A-U). Although this suggests that the lesion recognized is cytosine-derived, there was no correlation between the number of uracils induced and the number of binding sites, suggesting that the protein-bound damage is not a uracil-containing mismatch. Modification of the treatment protocol to reduce elimination of the bisulfite from the base adducts increased the level of binding, suggesting that the protein recognizes a base-sulfonate adduct.     

A theoretical model of structural chromosome rearrangement induction]

A theoretical model to explain the occurrence of chromosome rearrangements at the DNA molecular level is proposed. It involves: (10) breakage of one strand of each DNA molecule; (20) enzymatic digestion by an exonuclease of a short sequence of bases adjacent to the lesion, hence the presence of a short monocatenar segment on the complementary strand; (30) association of two monocatenar sequences if they are sufficiently complementary hence the formation of an X-shaped heteroduplex; (40) repair synthesis of DNA and, during this repair, an exhange of strands identical to that occuring in sister-chromatid exchanges, hence the rearrangement. This model has the advantage of necessitating one single breakpoint per DNA molecule, and explains the efficiency of radiations with low linear transfer of energy, such as X- and gamma-rays. It also explains how rearrangements such as translocations can be totally balanced since they would occur base to base within identical sequences.     

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

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

DNA Sequence Determinants Controlling Affinity,Stability and Shape of DNA Complexes Bound by the Nucleoid Protein Fis

The abundant Fis nucleoid protein selectively binds poorly related DNA sequences with high affinities to regulate diverse DNA reactions. Fis binds DNA primarily through DNA backbone contacts and selects target sites by reading conformational properties of DNA sequences, most prominently intrinsic minor groove widths. High-affinity binding requires Fis-stabilized DNA conformational changes that vary depending on DNA sequence. In order to better understand the molecular basis for high affinity site recognition, we analyzed the effects of DNA sequence within and flanking the core Fis binding site on binding affinity and DNA structure. X-ray crystal structures of Fis-DNA complexes containing variable sequences in the noncontacted center of the binding site or variations within the major groove interfaces show that the DNA can adapt to the Fis dimer surface asymmetrically. We show that the presence and position of pyrimidine-purine base steps within the major groove interfaces affect both local DNA bending and minor groove compression to modulate affinities and lifetimes of Fis-DNA complexes. Sequences flanking the core binding site also modulate complex affinities, lifetimes, and the degree of local and global Fis-induced DNA bending. In particular, a G immediately upstream of the 15 bp core sequence inhibits binding and bending, and A-tracts within the flanking base pairs increase both complex lifetimes and global DNA curvatures. Taken together, our observations support a revised DNA motif specifying high-affinity Fis binding and highlight the range of conformations that Fis-bound DNA can adopt. The affinities and DNA conformations of individual Fis-DNA complexes are likely to be tailored to their context-specific biological functions.     

Base pair conformation-dependent excision of benzo[a]pyrene diol epoxide-guanine adducts by human nucleotide excision repair enzymes.

Human nucleotide excision repair processes carcinogen-DNA adducts at highly variable rates, even at adjacent sites along individual genes. Here, we identify conformational determinants of fast or slow repair by testing excision of N2-guanine adducts formed by benzo[a]pyrene diol epoxide (BPDE), a potent and ubiquitous mutagen that induces mainly G x C-->T x A transversions and frameshift deletions. We found that human nucleotide excision repair processes the predominant (+)-trans-BPDE-N2-dG adduct 15 times less efficiently than a standard acetylaminofluorene-C8-dG lesion in the same sequence. No difference was observed between (+)-trans- and (-)-trans-BPDE-N2-dG, but excision was enhanced about 10-fold by changing the adduct configurations to either (+)-cis- or (-)-cis-BPDE-N2-dG. Conversely, excision of (+)-cis- and (-)-cis- but not (+)-trans-BPDE-N2-dG was reduced about 10-fold when the complementary cytosine was replaced by adenine, and excision of these BPDE lesions was essentially abolished when the complementary deoxyribonucleotide was missing. Thus, a set of chemically identical BPDE adducts yielded a greater-than-100-fold range of repair rates, demonstrating that nucleotide excision repair activity is entirely dictated by local DNA conformation. In particular, this unique comparison between structurally highly defined substrates shows that fast excision of BPDE-N2-dG lesions is correlated with displacement of both the modified guanine and its partner base in the complementary strand from their normal intrahelical positions. The very slow excision of carcinogen-DNA adducts located opposite deletion sites reveals a cellular strategy that minimizes the fixation of frameshifts after mutagenic translesion synthesis.     

Conformation and dynamics of abasic sites in DNA investigated by time-resolved fluorescence of 2-aminopurine

Abasic sites are highly mutagenic lesions in DNA that arise as intermediates in the excision repair of modified bases. These sites are generated by the action of damage-specific DNA glycosylases and are converted into downstream intermediates by the specific activity of apurinic/apyrimidinic endonucleases. Enzymes in both families have been observed in crystal structures to impose deformations on the abasic-site DNA, including DNA kinking and base flipping. On the basis of these apparent protein-induced deformations, we propose that altered conformation and dynamics of abasic sites may contribute to the specificity of these repair enzymes. Previously, measurements of the steady-state fluorescence of the adenine analogue 2-aminopurine (2AP) opposite an abasic site demonstrated that binding of divalent cations could induce a conformational change that increased the accessibility of 2AP to solute quenching [Stivers, J. T. (1998) Nucleic Acids Res. 26, 3837-44]. We have performed time-resolved fluorescence experiments to characterize the states involved in this conformational change. Interpretation of these studies is based on a recently developed model attributing the static and dynamic fluorescence quenching of 2AP in DNA to aromatic stacking and collisional interactions with neighboring bases, respectively (see the preceding paper in this issue). The time-resolved fluorescence results indicate that divalent cation binding shifts the equilibrium of the abasic site between two conformations: a "closed" state, characterized by short average fluorescence lifetime and complex decay kinetics, and an "open" state, characterized by monoexponential decay with lifetime approximately that of the free nucleoside. Because the lifetime and intensity decay kinetics of 2AP incorporated into DNA are sensitive primarily to collisional interactions with the neighboring bases, the absence of dynamic quenching in the open state strongly suggests that the fluorescent base is extrahelical in this conformation. Consistent with this interpretation, time-resolved quenching studies reveal that the open state is accessible to solute quenching by potassium iodide, but the closed state is not. Greater static quenching is observed in the abasic site when the fluorescent base is flanked by 5'- and 3'-thymines than in the context of 5'- and 3'-adenines, indicating that 2AP is more stacked with the neighboring bases in the former sequence. These results imply that the conformation of the abasic site varies in a sequence-dependent manner. Undamaged sequences in which the abasic site is replaced by thymine do not exhibit an open state and have different levels of both static and dynamic quenching than their damaged homologues. These differences in structure and dynamics may be significant determinants of the high specific affinity of repair enzymes for the abasic site.     

Non-B DNA-forming sequences and WRN deficiency independently increase the frequency of base substitution in human cells

Although alternative DNA secondary structures (non-B DNA) can induce genomic rearrangements, their associated mutational spectra remain largely unknown. The helicase activity of WRN, which is absent in the human progeroid Werner syndrome, is thought to counteract this genomic instability. We determined non-B DNA-induced mutation frequencies and spectra in human U2OS osteosarcoma cells and assessed the role of WRN in isogenic knockdown (WRN-KD) cells using a supF gene mutation reporter system flanked by triplex- or Z-DNA-forming sequences. Although both non-B DNA and WRN-KD served to increase the mutation frequency, the increase afforded by WRN-KD was independent of DNA structure despite the fact that purified WRN helicase was found to resolve these structures in vitro. In U2OS cells, ～70% of mutations comprised single-base substitutions, mostly at G·C base-pairs, with the remaining ～30% being microdeletions. The number of mutations at G·C base-pairs in the context of NGNN/NNCN sequences correlated well with predicted free energies of base stacking and ionization potentials, suggesting a possible origin via oxidation reactions involving electron loss and subsequent electron transfer (hole migration) between neighboring bases. A set of ～40,000 somatic mutations at G·C base pairs identified in a lung cancer genome exhibited similar correlations, implying that hole migration may also be involved. We conclude that alternative DNA conformations, WRN deficiency and lung tumorigenesis may all serve to increase the mutation rate by promoting, through diverse pathways, oxidation reactions that perturb the electron orbitals of neighboring bases. It follows that such "hole migration" is likely to play a much more widespread role in mutagenesis than previously anticipated.     

DNA polymerase activity at the single-molecule level

DNA polymerases are essential enzymes responsible for replication and repair of DNA in all organisms. To replicate DNA with high fidelity, DNA polymerases must select the correct incoming nucleotide substrate during each cycle of nucleotide incorporation, in accordance with the templating base. When an incorrect nucleotide is sometimes inserted, the polymerase uses a separate 3'→5' exonuclease to remove the misincorporated base (proofreading). Large conformational rearrangements of the polymerase-DNA complex occur during both the nucleotide incorporation and proofreading steps. Single-molecule fluorescence spectroscopy provides a unique tool for observation of these dynamic conformational changes in real-time, without the need to synchronize a population of DNA-protein complexes.