A hypothesis for what conformation of the major adduct of (+)-anti-B[a]PDE (N2-dG) causes G-->T versus G-->A mutations based upon a correlation between mutagenesis and molecular modeling results

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.

Structural and energetic characterization of the major DNA adduct formed from the food mutagen ochratoxin A in the NarI hotspot sequence: influence of adduct ionization on the conformational preferences and implications for the NER propensity

The nephrotoxic food mutagen ochratoxin A (OTA) produces DNA adducts in rat kidneys, the major lesion being the C8-linked-2′-deoxyguanosine adduct (OTB-dG). Although research on other adducts stresses the importance of understanding the structure of the associated adducted DNA, site-specific incorporation of OTB-dG into DNA has yet to be attempted. The present work uses a robust computational approach to determine the conformational preferences of OTB-dG in three ionization states at three guanine positions in the NarI recognition sequence opposite cytosine. Representative adducted DNA helices were derived from over 2160 ns of simulation and ranked via free energies. For the first time, a close energetic separation between three distinct conformations is highlighted, which indicates OTA-adducted DNA likely adopts a mixture of conformations regardless of the sequence context. Nevertheless, the preferred conformation depends on the flanking bases and ionization state due to deviations in discrete local interactions at the lesion site. The structural characteristics of the lesion thus discerned have profound implications regarding its repair propensity and mutagenic outcomes, and support recent experiments suggesting the induction of double-strand breaks and deletion mutations upon OTA exposure. This combined structural and energetic characterization of the OTB-dG lesion in DNA will encourage future biochemical experiments on this potentially genotoxic lesion.

Biological properties of single chemical-DNA adducts: a twenty year perspective

The genome and its nucleotide precursor pool are under sustained attack by radiation, reactive oxygen and nitrogen species, chemical carcinogens, hydrolytic reactions, and certain drugs. As a result, a large and heterogeneous population of damaged nucleotides forms in all cells. Some of the lesions are repaired, but for those that remain, there can be serious biological consequences. For example, lesions that form in DNA can lead to altered gene expression, mutation, and death. This perspective examines systems developed over the past 20 years to study the biological properties of single DNA lesions.

Integrating S-phase checkpoint signaling with trans-lesion synthesis of bulky DNA adducts

Bulky adducts are DNA lesions generated in response to environmental agents including benzo[a]pyrene (a combustion product) and solar ultraviolet radiation. Error-prone replication of adducted DNA can cause mutations, which may result in cancer. To minimize the detrimental effects of bulky adducts and other DNA lesions, S-phase checkpoint mechanisms sense DNA damage and integrate DNA repair with ongoing DNA replication. The essential protein kinase Chk1 mediates the S-phase checkpoint, inhibiting initiation of new DNA synthesis and promoting stabilization and recovery of stalled replication forks. Here we review the mechanisms by which Chk1 is activated in response to bulky adducts and potential mechanisms by which Chk1 signaling inhibits the initiation stage of DNA synthesis. Additionally, we discuss mechanisms by which Chk1 signaling facilitates bypass of bulky lesions by specialized Y-family DNA polymerases, thereby attenuating checkpoint signaling and allowing resumption of normal cell cycle progression.

Mutagenesis studies of the major benzo[a]pyrene N2-dG adduct in a 5'-TG versus a 5'-UG sequence: removal of the methyl group causes a modest decrease in the [G->T/G->A] mutational ratio

The potent mutagen/carcinogen benzo[a]pyrene (B[a]P) is metabolically activated to (+)-anti-B[a]PDE, which induces a full spectrum of mutations primarily at the G:C base pairs (e.g. GC-->TA, GC-->AT, etc.). Each of these mutations can be induced by its major adduct [+ta]-B[a]P-N(2)-dG, where DNA sequence context appears to influence both the quantitative and qualitative pattern of mutagenesis. We noted previously that 5'-TG sequences tend to have a higher fraction of G-->T mutations for both [+ta]-B[a]P-N(2)-dG and (+)-anti-B[a]PDE in comparison with 5'-CG, 5'-GG or 5'-AG sequences. To investigate a possible structural element for this trend, the role (if any) of the methyl group on the 5'-T is considered. Using adduct site-specific means, the [G-->T/G-->A] mutational ratio for [+ta]-B[a]P-N(2)-dG is determined to be approximately 1.08 in a 5'-TGT sequence, and approximately 0.60 in a 5'-UGT sequence. (G-->C mutations are minor.) Although this modest approximately 1.8-fold decrease in [G-->T/G-->A] ratio is statistically significant (P = 0.03), it suggests that the methyl group on the 5'-T is not the main reason why a 5'-T tends to enhance G-->T mutations. This study was prompted by an adduct conformational hypothesis, which predicted that the removal of the methyl group in a 5'-TG sequence would lower the fraction of G-->T mutations; however, the approximately 1.8-fold decrease is too small to do additional experiments to assess whether this conformational hypothesis, or other hypotheses, are the true cause of the decrease, which is discussed in this paper.

Supercoiled DNA promotes formation of intercalated cis-N2-deoxyguanine adducts and base-stacked trans-N2-deoxyguanine adducts by (+)-7R,8S-dihydrodiol-9S,10R-epoxy-7,8,9,10-tetra- hydrobenzo[a]pyrene

The highly reactive and mutagenic benzo[a]pyrene metabolite, (+)-7R,8S-dihydroxy-9S,10R-epoxy-7,8,9,10-tetrahydrobenzo[a]pyrene (BPDE), forms predominantly N2-deoxyguanine DNA adducts in two stereoisomeric configurations (cis and trans). In previous in vitro assays using oligonucleotide substrates site specifically modified with cis- and trans-BPDE adducts, the nucleotide excision repair (NER) systems of eukaryotes and prokaryotes incise cis-BPDE adducts more efficiently than trans-BPDE adducts [Hess, et al. (1997) Mol. Cell Biol 17, 7069; Zou, et al. (2001) Biochemistry 40, 2923). We investigated the influence of DNA secondary structure on stereospecificity of BPDE adduct formation, and incision of BPDE adducts by the prokaryotic UvrABC NER endonuclease was examined. BPDE adducts formed at low density on supercoiled plasmids were incised 6-7-fold better by the thermoresistant Bacillus caldotenaxUvrABC than were BPDE adducts formed on linear DNA. Linearizing supercoiled plasmid DNAs after BPDE adduct formation did not diminish incision efficiency. These results suggested that configuration and/or conformation of adducts formed on linear and supercoiled DNAs differed. This hypothesis was confirmed by low temperature fluorescence spectroscopy of adducted supercoiled and linear DNAs. Spectroscopic results indicated that intercalated cis-BPDE adducts as well as base-stacked trans-BPDE adducts formed more abundantly in supercoiled DNA than in linear DNA. A higher cis to trans adduct ratio in supercoiled DNA was confirmed by high resolution [32P]postlabeling analyses. These results demonstrate that DNA secondary structure influences both configuration and conformation of BPDE adducts formed at low density (approximately 1 adduct/kbp) and suggests that the ratio of cis- to trans-BPDE adducts and amount of base-stacked trans adducts formed under physiological exposure conditions may be higher than inferred from high dose experiments.

Molecular modeling of the major benzo[a]pyrene N2-dG adduct in cases where mutagenesis results are known in double stranded DNA

The potent mutagen/carcinogen benzo[a]pyrene (B[a]P) is metabolically activated to (+)-anti-B[a]PDE, which induces a full spectrum of mutations (e.g. GC-->TA, GC-->AT, etc.). One hypothesis for this complexity is that different mutations are induced by different conformations of its major adduct [+ta]-B[a]P-N2-dG when bypassed during DNA replication (probably by different DNA polymerases). Previous molecular modeling studies suggested that B[a]P-N2-dG adducts can in principle adopt at least 16 potential conformational classes in ds-DNA. Herein we report on molecular modeling studies with the eight conformations most likely to be relevant to base substitution mutagenesis in 10 cases where mutagenesis has been studied in ds-DNA plasmids in E. coli with B[a]P-N2-dG adducts of differing stereoisomers and DNA sequence contexts, as well as in five cases where the conformation is known by NMR. Of the approximately 11,000 structures generated in this study, the computed lowest energy structures are reported for 120 cases (i.e. eight conformations and 15 examples), and their conformations compared. Of the eight conformations, four are virtually always computed to be high in energy. The remaining four lower energy conformations include two with the BP moiety in the minor groove (designated: BPmi5 and BPmi3), and two base-displaced conformations, one with the dG moiety in the major groove (designated: Gma5) and one with the dG in the minor groove (designated: Gmi3). Interestingly, these four are the only conformations that have been observed for B[a]P-N2-dG adducts in NMR studies. Independent of sequence contexts and adduct stereochemistry, BPmi5 structures tend to look reasonably similar, as do BPmi3 structures, while the base-displaced structures Gma5 and BPmi3 tend to show greater variability in structure. A correlation was sought between modeling and mutagenesis results in the case of the low energy conformations BPmi5, BPmi3, Gma5 and Gma3. Plots of log[(G-->T)/(G-->A)] versus energy[(conformation X)-(conformation Y)] were constructed for all six pairwise combinations of these four conformations, and the only plot giving a straight line involved Gma5 and Gmi3. While this finding is striking, its significance is unclear (as discussed).

Evading the proofreading machinery of a replicative DNA polymerase: induction of a mutation by an environmental carcinogen

DNA replication fidelity is dictated by DNA polymerase enzymes and associated proteins. When the template DNA is damaged by a carcinogen, the fidelity of DNA replication is sometimes compromized, allowing mispaired bases to persist and be incorporated into the DNA, resulting in a mutation. A key question in chemical carcinogenesis by metabolically activated polycyclic aromatic hydrocarbons (PAHs) is the nature of the interactions between the carcinogen-damaged DNA and the replicating polymerase protein that permits the mutagenic misincorporation to occur. PAHs are environmental carcinogens that, upon metabolic activation, can react with DNA to form bulky covalently linked combination molecules known as carcinogen-DNA adducts. Benzo[a]pyrene (BP) is a common PAH found in a wide range of material ingested by humans, including cigarette smoke, car exhaust, broiled meats and fish, and as a contaminant in other foods. BP is metabolically activated into several highly reactive intermediates, including the highly tumorigenic (+)-anti-benzo[a]pyrene diol epoxide (BPDE). The primary product of the reaction of (+)-anti-BPDE with DNA, the (+)-trans-anti-benzo[a]pyrene diol epoxide-N(2)-dG ((+)-ta-[BP]G) adduct, is the most mutagenic BP adduct in mammalian systems and primarily causes G-to-T transversion mutations, resulting from the mismatch of adenine with BP-damaged guanine during replication. In order to elucidate the structural characteristics and interactions between the DNA polymerase and carcinogen-damaged DNA that allow a misincorporation opposite a DNA lesion, we have modeled a (+)-ta-[BP]G adduct at a primer-template junction within the replicative phage T7 DNA polymerase containing an incoming dATP, the nucleotide most commonly mismatched with the (+)-ta-[BP]G adduct during replication. A one nanosecond molecular dynamics simulation, using AMBER 5.0, has been carried out, and the resultant trajectory analyzed. The modeling and simulation have revealed that a (+)-ta-[BP]G:A mismatch can be accommodated stably in the active site so that the fidelity mechanisms of the polymerase are evaded and the polymerase accepts the incoming mutagenic base. In this structure, the modified guanine base is in the syn conformation, with the BP moiety positioned in the major groove, without interfering with the normal protein-DNA interactions required for faithful polymerase function. This structure is stabilized by a hydrogen bond between the modified guanine base and dATP partner, hydrophobic interactions between the BP moiety and the polymerase, a hydrogen bond between the modified guanine base and the polymerase, and several hydrogen bonds between the BP moiety and polymerase side-chains. Moreover, the G:A mismatch in this system closely resembles the size and shape of a normal Watson-Crick pair. These features reveal how the polymerase proofreading machinery may be evaded in the presence of a mutagenic carcinogen-damaged DNA, so that a mismatch can be accommodated readily, allowing bypass of the adduct by the replicative T7 DNA polymerase.

Factors that influence the mutagenic patterns of DNA adducts from chemical carcinogens

Carcinogens are generally mutagens, which is understandable given that tumor cells grow uncontrollably because they have mutations in critical genes involved in growth control. Carcinogens often induce a complex pattern of mutations (e.g., GC-->TA, GC-->AT, etc.). These mutations are thought to be initiated when a DNA polymerase encounters a carcinogen-DNA adduct during replication. In principle, mutational complexity could be due to either a collection of different adducts each inducing a single kind of mutation (Hypothesis 1a), or a single adduct inducing different kinds of mutations (Hypothesis 1b). Examples of each are discussed. Regarding Hypothesis 1b, structural factors (e.g., DNA sequence context) and biological factors (e.g., differing DNA polymerases) that can affect the pattern of adduct mutagenesis are discussed. This raises the question: how do structural and biological factors influence the pattern of adduct mutagenesis. For structural factors, three possibilities are considered: (Hypothesis 2a) a single conformation of an adduct giving rise to multiple mutations -- dNTP insertion by DNA polymerase being influenced by (e.g.) the surrounding DNA sequence context; (Hypothesis 2b) a variation on this ("dislocation mutagenesis"); or (Hypothesis 2c) a single adduct adopting multiple conformations, each capable of giving a different pattern of mutations. Hypotheses 2a, 2b and 2c can each in principle rationalize many mutational results, including how the pattern of adduct mutagenesis might be influenced by factors, such as DNA sequence context. Five lines of evidence are discussed suggesting that Hypothesis 2c can be correct for base substitution mutagenesis. For example, previous work from our laboratory was interpreted to indicate that [+ta]-B[a]P-N(2)-dG in a 5'-CGG sequence context (G115) could be trapped in a conformation giving predominantly G-->T mutations, but heating caused the adduct to equilibrate to its thermodynamic mixture of conformations, leading to a decrease in the fraction of G-->T mutations. New work is described suggesting that [+ta]-B[a]P-N(2)-dG at G115 can also be trapped predominantly in the G-->A mutational conformation, from which equilibration can also occur, leading to an increase in the fraction of G-->T mutations. Evidence is also presented that the fraction of G-->T mutations is higher when [+ta]-B[a]P-N(2)-dG at G115 is in ss-DNA ( approximately 89%) vs. ds-DNA ( approximately 66%), a finding that can be rationalized if the mixture of adduct conformations is different in ss- and ds-DNA. In summary, the factors affecting adduct mutagenesis are reviewed and five lines of evidence that support one hypothesis (2c: adduct conformational complexity can cause adduct mutational complexity) are discussed.

The processing of a Benzo(a)pyrene adduct into a frameshift or a base substitution mutation requires a different set of genes in Escherichia coli

Replication through a single DNA lesion may give rise to a panel of translesion synthesis (TLS) events, which comprise error-free TLS, base substitutions and frameshift mutations. In order to determine the genetic control of the various TLS events induced by a single lesion, we have chosen the major N2-dG adduct of (+)-anti-Benzo(a)pyrene diol epoxide [(+)-anti-BPDE] adduct located within a short run of guanines as a model lesion. Within this sequence context, in addition to the major event, i.e. error-free TLS, the adduct also induces base substitutions (mostly G --> T transversions) and -1 frameshift mutations. The pathway leading to G --> T base substitution mutagenesis appears to be SOS independent, suggesting that TLS is most probably performed by the replicative Pol III holoenzyme itself. In contrast, both error-free and frameshift TLS pathways are dependent upon SOS-encoded functions that belong to the pool of inducible DNA polymerases specialized in TLS (translesional DNA polymerases), namely umuDC (Pol V) and dinB (Pol IV). It is likely that, given the diversity of conformations that can be adopted by lesion-containing replication intermediates, cells use one or several translesional DNA polymerases to achieve TLS.

Sequence-dependent conformational perturbation in DNA duplexes containing an epsilonA.T mismatch using molecular dynamics simulation

Previous experiments from this laboratory showed that 1, N:(6)-ethenoadenine (epsilonA) in 15mer DNA oligonucleotide duplexes with GGepsilonAGG and CCepsilonACC central sequences is repaired 3-5-fold more efficiently than in duplexes containing AAepsilonAAA and TTepsilonATT central sequences. This sequence dependence in repair rates appeared to correlate with the observed thermodynamic stability of these duplexes [Hang et al. (1998) J. Biol. Chem., 273, 33406-33413]. In the present work, unrestrained molecular dynamics was used to evaluate the sequence-dependent structural features of these duplexes. Explicit solvent and the particle mesh Ewald method were applied for the accurate representation of the electrostatic interactions. The differences observed in the axis- and intra-base pair parameters were primarily localized at the epsilonA*T mismatch in all sequences and indicate conformational diversity between the structures. However, all four structures remained in the B-conformational family. In the tip, tilt and propeller twist parameters for the five central base pairs, larger perturbations were found for the two duplexes with epsilonA flanked by A or T bases than for duplexes with epsilonA flanked by G or C bases. As a result of these perturbations, the average global curvature of the AAepsilonAAA and TTvarepsilonATT DNA duplexes was larger by approximately 12 degrees than that of the duplexes with the GGepsilonAGG and CCepsilonACC central sequences. The observed conformational differences between the duplexes containing A or T and G or C neighbors of epsilonA may contribute to the observed differential enzymatic repair of the same sequences.

Toward an understanding of the role of DNA adduct conformation in defining mutagenic mechanism based on studies of the major adduct (formed at N(2)-dG) of the potent environmental carcinogen, benzo[a]pyrene

The process of carcinogenesis is initiated by mutagenesis, which often involves replication past damaged DNA. One question - what exactly is a DNA polymerase seeing when it incorrectly copies a damaged DNA base (e.g., inserting dATP opposite a dG adduct)? - has not been answered in any case. Herein, we reflect on this question, principally by considering the mutagenicity of one activated form of benzo[a]pyrene, (+)-anti-B[a]PDE, and its major adduct [+ta]-B[a]P-N(2)-dG. In previous work, [+ta]-B[a]P-N(2)-dG was shown to be capable of inducing>95% G-->T mutations in one sequence context (5'-TGC), and approximately 95% G-->A mutations in another (5'-AGA). This raises the question - how can a single chemical entity induce different mutations depending upon DNA sequence context? Our current working hypothesis is that adduct conformational complexity causes adduct mutational complexity, where DNA sequence context can affect the former, thereby influencing the latter. Evidence supporting this hypothesis was discussed recently (Seo et al., Mutation Res. [in press]). Assuming this hypothesis is correct (at least in some cases), one goal is to consider what these mutagenic conformations might be. Based on molecular modeling studies, 16 possible conformations for [+ta]-B[a]P-N(2)-dG are proposed. A correlation between molecular modeling and mutagenesis work suggests a hypothesis (Hypothesis 3): a base displaced conformation with the dG moiety of the adduct in the major vs. minor groove gives G-->T vs. G-->A mutations, respectively. (Hypothesis 4, which is a generalized version of Hypothesis 3, is also proposed, and can potentially rationalize aspects of both [+ta]-B[a]P-N(2)-dG and AP-site mutagenesis, as well as the so-called "A-rule".) Finally, there is a discussion of how conformational complexity might explain some unusual mutagenesis results that suggest [+ta]-B[a]P-N(2)-dG can become trapped in different conformations, and why we think it makes sense to interpret adduct mutagenesis results by modeling ds-DNA (at least in some cases), even though the mutagenic event must occur at a ss/ds-DNA junction in the presence of a DNA polymerase.

Molecular modeling of the major adduct of (+)-anti-B[a]PDE (N2-dG) in the eight conformations and the five DNA sequences most relevant to base substitution mutagenesis

The potent mutagen/carcinogen 7R,8S-dihydroxy-9S, 10R-epoxy-7,8,9,10-tetrahydrobenzo[a]pyrene [(+)-anti-B[a]PDE], which is the activated form of benzo[a]pyrene (B[a]P), is able to induce different kinds of mutations (G-->T, G-->A, etc.). One hypothesis for this is that different mutations are induced depending upon the conformation of its major adduct ([+ta]-B[a]P-N2-dG) when bypassed during DNA replication. Based on molecular modeling, there appear to be at least 16 potential conformations that the major adduct [+ta]-B[a]P-N2-dG can adopt in dsDNA. Regarding base substitution mutagenesis, eight conformations are most likely to be relevant. In two conformations the dG moiety of the adduct is base paired with its complementary dC and the B[a]P moiety is in the minor groove. In two others the dG moiety of the adduct is in the Hoogsteen orientation and the B[a]P moiety is in the major groove. There are four base displaced structures in which the B[a]P moiety of the adduct is stacked with the surrounding base pairs, two with dG in the major groove and two with dG in the minor groove. Using a simulated annealing protocol, these eight conformations were evaluated in five different DNA sequence contexts (5'-TGC-3', 5'-CGT-3', 5'-AGA-3', 5'-CGG-3' and 5'-GGG-3'); the latter were chosen because they may be particularly revealing about mutagenic mechanism based on studies with [+ta]-B[a]P-N2-dG and (+)-anti-B[a]PDE. For each conformation and each sequence context, 25 simulated annealing runs were conducted by systematically varying several parameters (such as the initial annealing temperature) based on a protocol established recently. The goal of this work was to exclude conformations that are clearly inferior. Three conformations are virtually always high in energy, including the two Hoogsteen oriented species and one of the base displaced species with dG in the major groove. Remarkably, the remaining five conformations are often quite close in energy and are deemed most likely to be relevant to mutagenesis (see accompanying paper).