|
Advertisement | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 279, Issue 14, 13435-13446, April 2, 2004
S-(2-Chloroethyl)glutathione-generated p53 Mutation Spectra Are Influenced by Differential Repair Rates More than Sites of Initial DNA Damage*![]() ¶ ||
From the
Received for publication, November 11, 2003 , and in revised form, January 5, 2004.
Several steps occur between the reaction of a chemical with DNA and a mutation, and each may influence the resulting mutation spectrum, i.e. nucleotides at which the mutations occur. The half-mustard S-(2-bro-moethyl)glutathione is the reactive conjugate implicated in ethylene dibromide-induced mutagenesis attributed to the glutathione-dependent pathway. A human p53-driven Ade reporter system in yeast was used to study the factors involved in producing mutations. The synthetic analog S-(2-chloroethyl)glutathione was used to produce DNA damage; the damage to the p53 exons was analyzed using a new fluorescence-based modification of ligation-mediated polymerase chain reaction and an automated sequencer. The mutation spectrum was strongly dominated by the G to A transition mutations seen in other organisms with S-(2-chloroethyl)glutathione or ethylene dibromide. The mutation spectrum clearly differed from the spontaneous spectrum or that derived from N-ethyl,N-nitrosourea. Distinct differences were seen between patterns of modification of p53 DNA exposed to the mutagen in vitro versus in vivo. In the four p53 exons in which mutants were analyzed, the major sites of mutation matched the sites with long half-lives of repair much better than the sites of initial damage. However, not all slowly repaired sites yielded mutations in part because of the lack of effect of mutations on phenotype. We conclude that the rate of DNA repair at individual nucleotides is a major factor in influencing the mutation spectra in this system. The results are consistent with a role of N7-guanyl adducts in mutagenesis.
The somatic theory of cancer holds that carcinogens are mutagens, at least with regard to what are classically termed tumor initiators (13). In a general model, cells are initiated by damage resulting from a DNA-alkylating agent, yielding mutations that are fixed by subsequent rounds of replication. Mutation spectra are generally observed when chemical damage to the DNA occurs; i.e. the sites of the mutation are not random and are often characteristic of the DNA-damaging agent (4, 5). Sometimes mutant spectra can be measured in tumors (in experimental animals or humans), although the existence of any "hotspots" for mutations does not necessarily indicate that the particular mutation is the cause of the tumor. The biochemical basis for where mutations occur (i.e. the mutant spectrum) is generally not well understood. Mutation resulting from chemical damage to DNA is a multistep process, and the observed mutant spectrum may be the result of events occurring during five or possibly more individual processes: (i) reaction of the chemical with DNA, (ii) further non-enzymatic reactions of the DNA adduct (e.g. opening of the imidazole ring in the case of guanyl-N7 alkyl adducts), (iii) DNA repair (or, more specifically, resistance to repair), (iv) the action of DNA polymerases, and (v) the biological effect of a particular nucleotide/amino acid change in the phenotypic assay used. Ethylene dibromide (1,2-dibromoethane, BrCH2CH2Br) was extensively used as an insecticide, fumigant, and anti-knock agent in leaded gasoline prior to its classification as a carcinogen in the 1970s (68). The chemical is mutagenic in a number of bacterial and eukaryotic systems, producing primarily GC to AT transitions (912). Ethylene dibromide appears not to be DNA-reactive itself; enzymatic conjugation with GSH yields a half-mustard (GSCH2CH2Br) that generates an episulfonium ion that reacts with DNA (13, 14) and leads to mutations (15). Recent work has demonstrated that similar episul-fonium ion chemistry appears to be involved in the reaction of the active site cysteine of the DNA repair protein O6-alkylguanine-DNA alkyltransferase with ethylene dibromide (16). The GSH half-mustard (GSCH2CH2Cl) forms primarily S-(2-(N7-guanyl)ethyl)GSH adducts (17, 18) and also the minor adducts S-[2-(N2-guanyl)ethyl]GSH, S-[2-(O6-guanyl)ethyl]GSH, and S-[2-(N1-adenyl)ethyl]GSH (12, 19). This reactive half-mustard can be synthesized, under anhydrous conditions, and used as a model for GSCH2CH2Br, with similar results in terms of DNA alkylation and mutagenesis (15). This mustard was used in some previous work on sequence selectivity with the lacZ gene in a bacteriophage M13mp18/Salmonella typhimurium system. The patterns of both DNA alkylation and mutation were determined in the lacZ sequence. Some similarity between the patterns of highest alkylation and mutation was seen, although the overlap was not particularly strong (12). These results with the GSH half-mustard can be considered in context of the general question of the role of DNA adduct profiles and their influence on mutation spectra (20). Studies with other physical and chemical agents have been done, particularly with UV light and PAHs1 in p53, ras, and hprt models (2125). In the case of damage by UV light, a role for differential rates of DNA repair along the gene has been concluded to be a dominant factor in explaining the mutation spectra that appear in human tumors (2628). With PAHs, the initial sites of DNA modification seem more important and are explained at least partly by C-5 methylation patterns at CpG sites in mammalian systems (2931).
We considered some of the issues of the mutation spectra of ethylene dibromide using the half-mustard GSCH2CH2Cl in Saccharomyces cerevisiae with a plasmid containing the human tumor suppressor gene p53 and an Ade reporter responsive to the transactivation activity of p53 (3234). This yeast system has been used previously with several physical and chemical agents (3537). We developed a variation of an LMPCR method to identify and quantify sites of DNA damage. This mutagen has some advantages, in the context of general questions about sequence selectivity, in that DNA alkylation by GSCH2CH2Cl is largely (but not exclusively) restricted to guanine N7 atoms, in contrast to UV- and PAH-derived DNA modifications. We found differences in the damage patterns for isolated DNA and the same DNA harbored inside of yeast cells. The FLMPCR method also allowed the quantitation of piperidine-labile sites in DNA in cells as a function of time. Our results provide evidence that the mutation sites are highly localized to regions of persistent DNA adducts (i.e. regions of slow repair).
In Vivo Exposure to S-(2-Chloroethyl)glutathione or ENUThe yDM56p strain uses a double selection for mutant p53 transactivation activity. Both Can1 and Ade2 gene expression are regulated by a p53-responsive promoter. Cells with wild type p53 express the Can1 (arginine permease) gene, resulting in the uptake of canavanine, which is toxic to the cells and limits growth on selective media. Cells with mutant p53 that cannot drive expression of Ade2 accumulate an adenine metabolic intermediate and are red in color. S. cerevisiae yDM56p strain was grown to mid-log phase (A600 = 0.3; 3 x 106 cells/ml) in 250 ml of high adenine complete minimal medium as described (36). Cells were washed three times with 30 ml of PBS buffer and resuspended in 4 ml of PBS buffer (2 x 108 cells/ml). ENU (Sigma) was added to a final concentration of 100 mM as described (38). S-(2-Chloroethyl)GSH (GSCH2CH2Cl), synthesized as described (15), was freshly dissolved in chilled 0.10 M Tris-HCl buffer (pH 7.7, containing 1 mM EDTA) to a final concentration of 1 M and immediately added to yeast cells to a final concentration of 120 mM. In both cases, incubation was at 30 °C for 30 min with shaking. Cells were washed three times with 10 ml of PBS buffer, resuspended in 4 ml of PBS buffer, and diluted in the same buffer. Appropriate dilutions (100 µl, usually 104 to 106 dilutions) were spread onto high adenine complete minimal medium plates to determine total colony number and on selective plates (low adenine complete minimal medium, containing 0.73 mM canavanine (Ref. 36)) to determine the mutant fraction. Incubation was at 30 °C until red colonies appeared (usually 34 days). For in vivo adduct detection experiments (see below), yeast cells were exposed to 200 mM GSCH2CH2Cl and DNA was immediately extracted after washing the cells. p53 DNA PreparationTotal DNA from the yDM56p yeast strain was purified using a DNeasy Tissue kit (Qiagen, Valencia, CA). Purified DNA was employed as a template for p53 cDNA amplification (exons 310) using the primers P3 and P4 (Operon Technologies, Alameda, CA) and conditions reported by Flaman et al. (34). The resulting PCR product was purified using a High Pure PCR Purification System (Roche Molecular Biochemicals) and collected with two 30-µl washes of elution buffer. This PCR product was employed either as a template for detection of mutations by the yeast functional assay described by Moshinsky et al. (36), for in vivo adduct detection by FLMPCR (after insertion into yeast cells and subsequent modification by GSCH2CH2Cl), or for in vitro exposure to GSCH2CH2Cl, and further detection of adducts as described below. In Vitro Exposure of Isolated p53 cDNA to GSCH2CH2ClSamples containing 1.5 µg of purified wild-type p53 PCR product were exposed to GSCH2CH2Cl (0 or 2 mM, final concentrations) in a 100-µl reaction for 20 min at 37 °C. As a positive control, a parallel reaction done with dimethylsulfate (Aldrich) was run for 10 min at 25 °C. Piperidine (Fisher Scientific) was added to give a final concentration of 1 M, as well as 100 µl of Nujol mineral oil (PerkinElmer Life Sciences) and incubated for 30 min at 90 °C (3941). After piperidine treatment, DNA fragments (a collection of DNA molecules of variable size) were purified using the Roche kit (see above). These fragments served as templates for FLMPCR. FLMPCRA general scheme of the FLMPCR technique is shown in Fig. 1. Piperidine-treated DNA fragments were used as templates in a first primer extension assay. P3 and P4 primers were employed in separate reactions for adduct detection on the transcribed and on the non-transcribed strands of p53, respectively. This procedure generated a family of variably sized, blunt-ended duplex molecules, which were substrates for T4 DNA ligase-catalyzed addition of an asymmetric linker (42). This linker provided each DNA fragment with a common defined end. The bulk of the DNA was then denatured and used as template in a second primer extension reaction with the same oligonucleotides (P3 or P4), so that the extension products included the sequence complementary to the added linker. After each member of the bulk DNA pool had two defined ends, the mixture was used as templates in an exponential PCR reaction (30 cycles) using as primers the largest oligonucleotide of the common linker, and one of four different internal (LA) primers, 5' end-labeled with 6-FAM (Operon). Primers LA1 (5'-CCCAGAAAACCTACCAGGGCAGCTACGG) and LA2 (5'-CCGAGTGGAAGGAAATTTGCGTGTGGAG) were used to localize the adducts on the transcribed strand of p53; primers LA3 (5'-GCCGCCCATGCAGGAACTGTTACACATG) and LA4 (5'-TCATTCAGCTCTCGGAACATCACGAAGCG) were used for adduct detection on the non-transcribed strand. After every step of the FLMPCR assay, the DNA fragments generated were purified using the Roche kit (see above). Finally, an aliquot (1.0%) of the exponential PCR reactions was mixed with appropriate amounts of loading buffer, deionized formamide, and ROX-350-labeled size markers (PerkinElmer Life Sciences), loaded onto an 4% (w/v) acrylamide sequencing gel, and analyzed using an ABI Prism 377 DNA sequencer (Applied Biosystems, Foster City, CA) in the Vanderbilt facility.
Adduct QuantitationThe GeneScan software of the ABI Prism 377 DNA sequencer generated an image of the fluorescent bands in the acrylamide gel. The analysis of this image resulted in a profile (electropherogram) of peaks (the DNA fragments labeled with 6-FAM) and size markers (labeled with ROX-350), and the corresponding tabular display of data indicated the size (in base pairs) and the height (in arbitrary units) of each peak. The position of the peaks relative to the markers in these profiles revealed the nucleotide position where each GSCH2CH2Cl adduct was formed. Furthermore, by calculating the ratio between the heights of a given peak in the GSCH2CH2Cl-treated sample and its equivalent in the negative control, it was possible to establish the relative adduction level of every nucleotide along both p53 strands (exons 58).
Adduct t
Mutagenesis with Alkylating AgentsExperiments were done with S. cerevisiae yDM56p to determine the mutation frequency in the p53 gene induced by two alkylating agents, ENU and GSCH2CH2Cl. Cells were exposed to 100 mM ENU or 120 mM GSCH2CH2Cl as described under "Experimental Procedures." Survival in nonselective medium at these doses rendered similar results ( 17% survival with respect to non-exposed, control cells). Exposure of yDM56p yeast to ENU produced a 20-fold increase in the induced mutant fraction. However, the mutant fraction resulting from GSCH2CH2Cl treatment was only 2.3-fold higher than the spontaneous fraction (Table I). This increase was less than the 10-fold level previously observed in the M13mp18 bacteriophage study (12), where the bacterial cells were deficient in nucleotide excision repair. It has been reported that a combined treatment with a given cytotoxic agent and the canavanine present in the selective medium kills a larger proportion of cells than the cytotoxic agent alone (36). The present results, therefore, are probably an underestimate of induced mutant fractions.
Mutant CharacterizationTo determine the proportion of red colonies that actually contain mutations in p53, the DNA from selected red colonies was PCR-amplified and checked in the functional assay (36). PCR products that rendered red colonies were considered p53 mutants, whereas those that produced white colonies were presumed to contain mutations outside the p53 region (e.g. in the Ade2 reporter gene). All PCR products from red colonies were analyzed by agarose gel electrophoresis to detect the presence of large insertions and deletions. Two out of 10 (20%) spontaneous red colonies were found to contain large insertions (20 and 70 bp long, determined by further sequencing; Table II). In contrast, only 2 of 32 (6%) ENU-induced clones had large deletions (Table II), and only 1 of 45 (2%) GSCH2CH2Cl-induced red colonies contained a deletion (13 bp, determined by nucleotide sequencing, Table III).
Results of the p53 functional assay are presented in Table I. Of 23 spontaneous red colonies tested, 17 (74%) contained mutant p53, and the sequencing of exons 58 revealed that 10 of them (59%) harbored mutations in this region of the gene. Among ENU-induced colonies, 51 of 71 (72%) contained mutant p53 as judged by the functional assay, and 34 (67%) of these had detectable mutations in the DNA binding domain of the gene. Of the GSCH2CH2Cl-induced red colonies, 95 of 133 (71%) were mutants in p53, but only 46 of these 95 (48%) harbored mutations in the central domain of the gene. Spontaneous and ENU-induced mutations found in the sequenced region of the gene are listed in Table II and the GSCH2CH2Cl-induced mutations in Table III. In the spontaneous mutants (in addition to the two large insertions), there were two 2-bp deletions, three transitions, and three transversions. In the ENU- and GSCH2CH2Cl-induced mutants, the vast majority of p53 mutations (78 and 87%, respectively) were G to A transitions. Mutations were observed in a similar number of codons after exposure to ENU and GSCH2CH2Cl (25 and 24, respectively) in the region between the 3' end of exon 4 and the 5' end of exon 9. ENU and GSCH2CH2Cl both produced mutations preferentially in purine-rich sites (see Tables II and III). The determined consensus sequence around the mutated base for ENU was (5')G32G45G84A39G45; whereas that for GSCH2CH2Cl was (5') G35G70G100G63G/C33 (where the indicated value in each position is the percentage of cases where that nucleotide appears in the mutants and the mutated base is in the central position (in bold font)). The distribution of mutations induced by the alkylating agents tested was not random. Some codons were more frequently changed than others. Furthermore, each treatment produced a distinctive and statistically different fingerprint of hotspots (p values < 0.05 calculated with the computer program designed by Cariello and co-workers (Refs. 44 and 45)). After ENU-treatment, codons 141, 177, 194, 199, and 278 were found mutated more than once (Table II). Codons 152, 177, 244, and 279 were frequent targets for mutation by GSCH2CH2Cl (Table III). Taken together, the mutations detected were symmetrically distributed between the transcribed and non-transcribed strands ( 50% in each strand for both chemicals tested). However, when individual exons were analyzed, it was evident that GSCH2CH2Cl produced mutations almost exclusively in the transcribed strand of exon 5 and in the non-transcribed strands of exons 7 and 8 (Table III). ENU treatment rendered rather diverse results; all mutations in exon 7 were in the non-transcribed strand, and mutations in other exons were rather equally distributed (Table II). To our knowledge, the p53 mutations present in clones SPO2, SPO5, ENU1, ENU13, and ENU54 and in GSCH2CH2Cl clones 9, 43, 49, and 91 have not been reported before. They were not included in the list of 14,968 entries available on a web site collection (p53.curie.fr).
Analysis of GSCH2CH2Cl-produced DNA Adducts by LMPCR AssaysAmong the different types of DNA adducts produced by GSCH2CH2Cl, the most abundant ( The strand breaks and the free 5'-phosphate termini generated by the combined action of GSCH2CH2Cl and piperidine render treated DNA a suitable substrate for analysis using LMPCR. Modifications to this technique allowed the mapping of DNA damage induced by GSCH2CH2Cl in the central domain (exons 58) of human p53 at nucleotide resolution. LMPCR assays have been successfully used to detect the presence of adducts in several DNA sequences at nucleotide resolution (4648). The traditional assay is limited to a relatively short stretch of sequence (usually <200 bp) (48). Another inconvenience is the need for several (usually 4) oligonucleotide primers to analyze a genomic region of some dozens of nucleotides. Third, the common protocol requires hybridization of exponentially amplified products to radiolabeled probes to reveal the adducted nucleotides. Finally, the quantitative analysis of adducted bases requires the use of autoradiography/densitometry. We modified the LMPCR protocol, using a total of seven primers (four of them 6-FAM-labeled) and the analysis of fluorescence-amplified PCR products performed with the help of GeneScanTM software (ABI Prism Genetic Analysis instrument, PerkinElmer Life Sciences). With these modifications we analyzed in vitro and in vivo GSCH2CH2Cl-induced adducts on both strands of the complete DNA binding domain (exons 58, total of 1,088 bp) of the human p53 gene. We term this technique FLMPCR (Fig. 1). The complete procedure can be completed in 2 days and renders reproducible results on the localization and (relative) quantification of nucleotide adducts.
Comparison of in Vivo and in Vitro GSCH2CH2Cl-induced Damage Distribution in the p53 GeneWe used concentrations of GSCH2CH2Cl similar to those used to generate the mutants (see above), to have patterns of alkylation that could be compared with the mutant spectra. From our previous (12, 15) and current experience, the level of DNA damage is such that the average exon is damaged no more than once, avoiding a contributing fraction of short fragments resulting from multiple cleavages. Subsequent work (see below) supports the lack of multiple cleavage, in that we did not observe any consistent trends of greater apparent damage at one end of exons (Fig. 2, e.g. panel C). Another technical consideration is variability in the efficiency of the FLMPCR steps among samples. In this work and the subsequent t
DNA damage distribution patterns obtained by FLMPCR from in vivo and in vitro induced GSCH2CH2Cl-treated samples on both strands along exons 58 of p53 are compared in Fig. 2. Although similar, the adduction patterns of in vivo and in vitro treatments showed several reproducible differences, not only regarding the relative patterns of peak heights, but also with the absence of individual bands or even groups of peaks. In vivo treatment with GSCH2CH2Cl produced an apparently higher level of adduction than the in vitro treatment (Fig. 2). Examples of sites preferentially modified following in vivo treatment include nucleotides located in the transcribed strand at codons 127 and 175178 (exon 5), 188193 (exon 6), 243, 248, and 250256 (exon 7), and 278, 297, and 300 (exon 8). In the non-transcribed strand, nucleotides at codons 147, 158, 160, 171, 174, and 180 (exon 5), 188, 199, and 206 (exon 6), and 265 and 287 (exon 8) were also more affected by the in vivo treatment and are considered probably more relevant to mutagenesis. Only nucleotides at codons 129, 137142 (transcribed strand of exon 5), 223 (non-transcribed strand of exon 6), 272 (transcribed strand of exon 8), and 299302 (non-transcribed strand of exon 8) presented higher adduction signals after in vitro treatment than in vivo. These consistently observed differences suggest that chromatin architecture within the cellular environment influences the accessibility of the mutagen to the target DNA. All four bases sustained some DNA damage (adenine could be a N3-alkyl adduct, susceptible to depurination (Ref. 40); limited cleavage at bases other than alkylated A and G has also been reported (Refs. 49 and 50)). As indicated above, the damage distribution was non-random. Some regions, in particular purine-rich tracts and especially those where three or more guanines occur together, showed stronger adduct signals. When each strand of individual exons were examined, it was clear that GSCH2CH2Cl adducts accumulate asymmetrically along the DNA sequence. The non-transcribed strands of exons 6 and 8 and the transcribed strands of exons 5 and 7 were the most affected by the exposure to GSCH2CH2Cl. The asymmetric distribution of guanines along these exons could explain that observation. Although almost all guanine nucleotides present in exons 58 displayed some level of adduction, only a few of them resulted in mutations according to the yeast functional assay (see Table III). Furthermore, no correlation was observed between the extent of initial modification and the resulting mutations. Rather, mutational hotspots were found in nucleotides with somewhat weak damage from GSCH2CH2Cl.
The GSCH2CH2Cl-induced Adduct t
The issue of the molecular basis of mutation spectra (and the relationship to tumor initiation and development) has been considered and reviewed by others (4, 5, 2224, 28, 29, 37, 51). At least five major steps on the pathway from a carcinogen to a mutation may be involved. (i) One step is binding selectivity (of the activated carcinogen) (52). This could be the result of any of at least three factors. One is the influence of the nucleotide sequence on the chemical reactivity at a site, e.g. the enhanced nucleophilicity of a G within a run of G nucleotides (or, more properly, a G with a 5' G neighbor (Refs. 41 and 5356)). Another factor is the influence of the sequence on DNA structure, e.g. favoring intercalation or groove binding (57). Still another factor is the structural influence of bound proteins, which may be a factor in the differences in modification patterns observed in vitro and in vivo (Fig. 2) (58). (ii) A second step in some cases is non-enzymatic conversion of a DNA adduct to another type of adduct, e.g. depurination of an N7-guanyl adduct or conversion to an imidazole ring-opened form (59). However, at the present time there is no evidence that these or other potential transformations (e.g. oxidation of 8-oxo-7,8-dihydroguanine (Ref. 60)) show sequence specificity, although the probability is likely. The product could be either less or more mutagenic (or blocking) than the initial adduct. (iii) The next step considered is DNA repair. Evidence already exists for a role of repair rates in determining the mutation spectra in tumors resulting from exposure of mammalian cells to sunlight (and corresponding to mutations in p53 in human tumors) (27, 28, 61). Some details of how the individual DNA repair systems can be influenced by individual chemical adducts in certain sequences have been published (62, 63), but few generalizations can be made at this time about these systems. (iv) Polymerization is the next step (if the DNA is not repaired). Several possibilities can be considered, depending upon the chemical adduct, the particular polymerase(s) involved, and the DNA sequence. The position at which an adduct is found can influence the blockage of a polymerase, the qualitative and quantitative nature of misincorporation opposite the adduct, and the tendency for the polymerase(s) to extend beyond the site of insertion (opposite the adduct) (64, 65). (v) The last step is the effect on biological activity beyond these steps, i.e. phenotype. The position of the mutation can have a dramatic effect on the phenotype; obviously, some changes will not influence the biological properties of a protein and some will be critical. Further, the phenotypic effects may differ between in vitro assay systems used in the laboratory and in tumors, because different functions may be involved. Because the yeast functional assay identifies mutant p53-expressing clones by a red/white colony color selection based on loss of transactivation ability of an integrated Ade2 gene, it can be assumed that, excepting the clones with double mutations (e.g. clones SPO2, ENU5, ENU13, and GSCH2CH2Cl 49; see Tables II and III), all the mutations obtained with this system affect the transactivating activity of p53 protein. In previous work we compared sites of S-(2-chloroethyl)GSH alkylation of the lacZ gene in bacteriophage DNA with the mutation spectrum (in the bacterium S. typhimurium) (12). That study demonstrated a dominant pattern of G to A transitions, which had previously been suggested in an Escherichia coli assay with ethylene dibromide (9). The results (G to A transitions) argue against a strong influence of depurination, where the apurinic sites that would be expected from depurination of S-[2-(N7-guanyl)ethyl]GSH adducts should produce G to T transversions (66). The predominance of G to A transitions has now been reported for studies with ethylene dibromide in several systems, including Drosophila and Chinese hamster ovary cells (911). In addition to S-[2-(N7-guanyl)ethyl]GSH, some other G adducts have been shown to be produced from S-(2-chloroethyl)GSH, including the N2 and O6-guanyl adducts (12). All three of the known G adducts (N7, N2, O6) have been shown to block replicative polymerases (67), although the G to A transitions have not been specifically attributed to a single G modification yet. Patterns of DNA modification observed in vitro and in vivo were not identical (Fig. 2), presumably because of either the coiling of the DNA or the presence of DNA-bound proteins. Cloutier et al. (58) have also reported differences in alkylation patterns inside and outside of cultured cells. Further, the initial pattern of DNA damage is not identical to that seen after several hours, because of the influence of variable half-lives of DNA repair at individual G sites (Figs. 4, 5, 6, 7). The sites at which mutations were clustered generally corresponded to sites with slower repair. However, some slowly repaired (and phenotypically detectable) sites did not produce mutants with the phenotype under consideration in the assay. In the earlier work with the mutagen S-(2-chloroethyl)GSH in the bacteriophage lacZ complementation system, some overlap was noted between the patterns of piperidine-generated cleavage at Gs and mutations but was not complete (12). As in the current study, the overlap of the patterns had many exceptions, i.e. regions of mutation with low alkylation and regions of extensive alkylation but no mutations. These studies were done in an excision repair-deficient strain of bacteria (to enhance mutations); thus, adduct repair could be discounted as a reason for the lack of correlation between initial adduct sites and mutations. We hypothesized, at that time, that the sequence context of a mutation influences the ability of polymerases to extend beyond the adducted site to a "fix" a mutation in the gene. A similar explanation may contribute to the finding in this study of stably adducted sites that are not adducted and some weakly adducted sites that are mutagenic, although this hypothesis has not been addressed yet. The method used here is focused on N7-guanyl adducts and is a variant of previously described LMPCR methods (4648). The use of a fluorescence end label (FLMPCR) was superior to previous approaches with 32P-end-labeled DNA in that it permitted in vivo studies and direct comparisons of alkylation patterns inside the cells with isolated DNA (Fig. 2). The use of fluorescence methodology was also important in that it allowed the use of a commercial sequencer instead of incubations with radioactive oligonucleotides, conventional gel electrophoresis, and autoradiography. A similar method has been employed in which a near-IR fluorescence label was used (48); our 6-FAM-labeled system was utilized with a conventional sequencer in a facility core sequencing laboratory. The majority of the previous literature on the subject of sequence-related effects of DNA modification comes from work with UV damage (especially cyclobutane dimers and pyrimidinezhy(64)pyrimidone photoproducts) and benzo[a]pyrene (dihydrodiol epoxide)-derived adducts (2631, 61, 6871). Comparisons of the patterns have been made with the sites of p53 mutations in tumors of people who have had exposure to these physical and chemical agents. (In the case of UV light, skin cancers are readily attributed to the agent; with lung cancers, exposure of smokers to benzo[a]pyrene is real but can be complicated by the plethora of other carcinogens.) In the case of UV light damage, the evidence suggests a strong role for DNA repair in influencing which adducts will go on to produce p53 mutations in tumors, as predicted from studies in human fibroblast cells (27). With benzo[a]pyrene-derived adducts, the situation is probably more complex, in part because of questions about the chemical origin of the lung tumors. With regard to benzo[a]pyrene (and its diol epoxides), the mutation spectra appear to be dominated by the sequence selectivity of DNA adduct formation, which is linked to the presence of 5-methyl-C at CpG sites (24, 31, 72). Rates of DNA repair have not been implicated as much as sites of adduct formation in the work by Pfeifers group (70, 71). However, transcription-coupled repair is an issue in terms of which strand is repaired. In other work by Wei and Conney (22, 7375), the mutation spectrum (with other genes) resulting from benzo[a]pyrene diol epoxide treatment of mammalian cells was dependent upon the dose of the carcinogen, and saturation of DNA repair has been proposed as a likely mechanism underlying this phenomenon. Work by Tangs group on human rasK gene codons 12 and 14 has also been interpreted in terms of effects of C methylation on binding, and rates of repair also appear to be a factor (25) (although these rasK codon 12 and 14 adducts have not been extensively evaluated in terms of producing mutations).
Epidemiological evidence for the tumorigenicity of ethylene dibromide in humans is very limited, and its classification as a carcinogen (and restricted industrial and agricultural use (Ref. 8)) is based on animal studies (6). No information is available about ethylene dibromide-generated mutation spectra of p53 or other genes in tumors. The p53 mutant spectra reported here might not predict sites of human p53 mutations in tumors, for several reasons. The influence of C methylation is not reported in the yeast system. Strand bias was not apparent in our work (Figs. 4, 5, 6, 7) and might occur in mammalian systems. Another point is that the mutant spectra are influenced strongly in this system by differential rates of DNA repair. However, the rates may be a factor in yeast more than in mammalian liver (e.g. a t Do these results mean that the S-[2-(N7-guanyl)ethyl]GSH adducts are the ones involved in producing mutations? In earlier work (12) and here, most of the mutations are at G nucleotides (specifically GC to AT transitions) and most of the mutations correspond to sites of apparent guanyl N7-alkylation and slower repair (Figs. 4, 5, 6, 7). However, minor amounts of two other G adducts can be formed from S-(2-chloroethyl)GSH, the N2 and O6 adducts (12). We have not specifically ruled out the possibility that these adducts could contribute to the FLMPCR results on sequence specificity of alkylation, although neither would be expected to on the basis of general chemical principles (39, 40). (However, hot piperidine treatment can lead to some other reactions other than cleavage at N7-guanyl and N3-adenyl alkyl adducts (Refs. 49 and 50), at slower rates.) We have not evaluated the relative efficiency of these reactions relative to N7-guanyl adducts but could not rule out small contributions. A fraction of the mutations produced by S-(2-chloroethyl)GSH in a bacterial system may be attributed to the O6-guanyl adduct, in that this adduct was shown to be a substrate for O6-alkylguanine-DNA alkyltransferase, and expression of the transferase attenuated the mutations produced by S-(2-chloroethyl)GSH (16).2 Mutagenicity of O6-guanyl adducts (which presumably would not be detected in the piperidine cleavage assay) would contribute to a lack of overlap between DNA adduct and mutation spectra (Figs. 4, 5, 6, 7), as would any residual spontaneous mutations (Table II).
The results are consistent with the mutagenicity of N7-guanyl adducts in producing G to A mutations. How could this process occur? Depurination (which would produce G to T transversions) and imidazole ring-opening (not detected at a limit of 1% (Ref. 12)) can be ruled out. The literature generally argues against the mutagenicity of N7-alkylguanine adducts (77). However, closer inspection of the literature and the cited review (77) indicate that most of the evidence, which is rather indirect, was developed with N7-methylguanine. Little direct analysis of misincorporation at N7-alkylguanine adducts has been done and none with the more recently discovered "translesional" polymerases. There is precedent for mutation at an N7-aflatoxin B1 adduct (76). In our earlier studies, we calculated that the frequency of G to A transitions (caused by GSCH2CH2Cl) in bacterial systems could not be attributed to adducts other than the N7-alkylguanine adduct, S-[2-N7-guanyl)ethyl]glutathione (12, 15). S-[2-N7-Guanyl)ethyl]glutathione was misincorporating with E. coli polymerase II (but not polymerase I), with dTTP being inserted (67) (other polymerases have not been examined).2 We have discussed possible mechanisms for misincorporation previously (12, 15, 78), primarily in the context of tautomerism and unusual ionization properties. The concept of base tautomerism as a basis for mispairing is not new (7980) and has theoretical (81) and experimental (8284) backing. The lowered pKa of the N1 atom of guanine resulting from N7-alkylation is well established both with the free base (82) and in double-stranded oligonucleotides (78), lowering the pKa to physiological pH (pH 78). A zwitterionic residue might indeed mispair, particularly if the polymerase is strongly blocked by the N7 substitution, as has been shown (67). These possibilities are being considered in experimental studies.2 In summary, we have characterized a system used as a model for the GSH-dependent mutagenicity of the carcinogen ethylene dibromide. The pattern of DNA alkylation observed inside cells differs from that with isolated DNA, even in the same coding sequence. The mutant spectrum did not correlate well with the pattern of initial DNA alkylation, but the mutations were generally at sites most resistant to repair. These results provide an example of the role of sequence-selective repair in determining the mutation spectrum in this system. The molecular basis of sequence-selective profiles of DNA repair is poorly understood. Our earlier research in this area suggested that nucleotide excision repair is important in preventing the mutagenicity of S-(2-chloroethyl)GSH (12), although more details of the process are unknown (or the potential involvement of other types of repair). Other evidence suggests at least partial roles for sequence selectivity of DNA repair in some other systems, e.g. UV light and PAH damage.
* This work was supported in part by United States Public Health Service Grants R01 ES10546 and P30 ES00267 (to F. P. G.). The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
¶ Present address: Dept. of Neurology, Vanderbilt University School of Medicine, Nashville, TN 37232. || To whom correspondence should be addressed: Dept. of Biochemistry and Center in Molecular Toxicology, Vanderbilt University School of Medicine, 638 Robinson Research Bldg., 23rd and Pierce Aves., Nashville, TN 37232-0146. Tel.: 615-322-2261; Fax: 615-322-3141; E-mail: f.guengerich{at}vanderbilt.edu.
1 The abbreviations used are: PAH, polycyclic aromatic hydrocarbon; ENU, N-ethyl,N-nitrosourea; PBS, phosphate-buffered saline (15 mM potassium phosphate (pH 7.4) containing 150 mM NaCl); LMPCR, ligation-mediated PCR; FLMPCR, fluorescence ligation-mediated PCR; 6-FAM, 6-carboxyfluorescein; ROX-350, 5-carboxy-X-rhodamine.
2 In principle, site-specific mutagenesis methods could be used to resolve the issue of relative mutagenicities of the three G adducts. Precedent now exists with aflatoxin B1 for introducing an unstable guanine N7-alkyl adduct into a cellular system (76). However, a major obstacle is the need to synthesize the adduct by chemical treatment of a G-containing oligonucleotide, followed by separation of all possible products. An issue arises in that the choice of an oligonucleotide target should reflect a site where mutations are known to occur. As shown in Table III, even among pentameric base sequences in mutation sites, only one of the sites identified here (codon 259) has a single G (and that mutation was only found once). Thus, synthesis of an appropriate N7-guanyl adduct for site-specific mutagenesis remains problematic.
We thank N. F. Cariello (University of North Carolina, Chapel Hill, NC) for the statistical program we used to compare patterns and J. A. Krauser (Vanderbilt University, Nashville, TN) for assistance in using this software, K. Bhat of the Vanderbilt sequencing facility for technical contributions, and particularly J. L. Cmarik (Frederick Cancer Center, Frederick, MD) for comments on the manuscript.
This article has been cited by other articles:
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Advertisement | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||