Synthesis of a duplex oligonucleotide containing a nitrogen mustard interstrand DNA-DNA cross-link.

Many cancer chemotherapeutic agents react with DNA and give adducts that block DNA replication, which is thought to result in cytotoxicity, especially in rapidly proliferating cells such as cancer cells. One class of these agents is bifunctionally reactive (e.g., the nitrogen mustards) and forms DNA-DNA cross-links. It is unknown whether inter- or intrastrand cross-links are more effective at blocking DNA replication. To evaluate this, a DNA shuttle vector is being constructed with an interstrand cross-link at a unique site. In the first step of this project, a duplex oligonucleotide containing an interstrand cross-link is isolated by denaturing polyacrylamide gel electrophoresis from the reaction of nitrogen mustard with two partially complementary oligodeoxynucleotides. The purified oligonucleotide product is characterized and shown to be cross-linked in a 5'-GAC-3' 3'-CTG-5' sequence by a nitrogen mustard moiety that is bound at the N(7)-position of the guanines in the opposing strands; the glycosylic bonds of these guanine adducts are stabilized in their corresponding imidazole ring-opened form. Nitrogen mustard is shown to react with a variety of oligonucleotides and, based upon these results, its preferred targets for interstrand cross-linking are 5'-GXC-3' sequences, where X can be any of the four deoxyribonucleotide bases.

Construction of a human shuttle vector containing a single nitrogen mustard interstrand, DNA-DNA cross-link at a unique plasmid location.

DNA cross-linking reagents are frequently unusually cytotoxic, and many, including the nitrogen mustards, are potent chemotherapeutic agents, presumably because DNA cross-links effectively block DNA replication. Most of these reagents form both inter- and intrastrand DNA cross-links, but it is unknown which is more effective at blocking replication and why. To evaluate the role of interstrand cross-links, a human shuttle vector was constructed that contains a single, nitrogen mustard interstrand cross-link at a unique site. In previous work (J.O. Ojwang, D. A. Grueneberg, and E. L. Loechler, Cancer Res., 49: 6529-6537, 1989) a duplex oligonucleotide was synthesized that had an interstrand cross-link derived from a nitrogen mustard moiety bound at the N(7)- position of the guanines in the opposing strands of a 5'-GAC-3' 3'-CTG-5' sequence. Herein, a procedure is described to incorporate this oligonucleotide into an SV40-based human shuttle vector, which was designed for these experiments. The purified cross-linked vector was characterized and shown: (a) to have a chemical (i.e., a nitrogen mustard) modification at the anticipated genome location; (b) to have a modification that covalently joins the two duplex strands of the vector together; and (c) to contain a single interstrand cross-link per genome. The methodologies described to construct this vector are expected to be generally applicable and, thus, site-specific incorporation of an interstrand cross-link derived from any appropriate chemical should be possible. These procedures complement existing methodologies that permit the incorporation of monoadducts and intrastrand cross-links into vectors in a site-specific manner.

Extrachromosomal probes for mutagenesis by carcinogens: studies on the mutagenic activity of O6-methylguanine built into a unique site in a viral genome.

This work examines the mutagenic activity of O6-methylguanine (O6MeGua), a DNA adduct formed by certain carcinogenic alkylating agents. A tetranucleotide, 5'-HOTpm6GpCpA-3', was synthesized and ligated into a four-base gap in the unique Pst I site of the duplex genome of the E. coli virus, M13mp8. The double-stranded ligation product was converted to single-stranded form and used to transform E. coli to produce progeny phage. The mutation frequency of O6MeGua was defined as the percentage of progeny phage with mutations in their Pst I site, and this value was determined to be 0.4%. To determine the impact of DNA repair on mutagenesis, cellular levels of O6MeGua-DNA methyltransferase (an O6MeGua-repair protein) were depleted by treatment of host cells for virus replication with N-methyl-N'-nitro-N-nitrosoguanidine (MNNG) prior to viral DNA uptake. In these host cells, the mutation frequency due to O6MeGua increased markedly with increasing MNNG dose (the highest mutation frequency observed was 20%). DNA sequence analysis of mutant genomes revealed that in both MNNG treated and untreated cells, O6MeGua induced exclusively G to A transitions.

The precautionary principle in environmental science.

Environmental scientists play a key role in society's responses to environmental problems, and many of the studies they perform are intended ultimately to affect policy. The precautionary principle, proposed as a new guideline in environmental decision making, has four central components: taking preventive action in the face of uncertainty; shifting the burden of proof to the proponents of an activity; exploring a wide range of alternatives to possibly harmful actions; and increasing public participation in decision making. In this paper we examine the implications of the precautionary principle for environmental scientists, whose work often involves studying highly complex, poorly understood systems, while at the same time facing conflicting pressures from those who seek to balance economic growth and environmental protection. In this complicated and contested terrain, it is useful to examine the methodologies of science and to consider ways that, without compromising integrity and objectivity, research can be more or less helpful to those who would act with precaution. We argue that a shift to more precautionary policies creates opportunities and challenges for scientists to think differently about the ways they conduct studies and communicate results. There is a complicated feedback relation between the discoveries of science and the setting of policy. While maintaining their objectivity and focus on understanding the world, environmental scientists should be aware of the policy uses of their work and of their social responsibility to do science that protects human health and the environment. The precautionary principle highlights this tight, challenging linkage between science and policy.

DNA guanine-guanine crosslinking sequence specificity of isophosphoramide mustard, the alkylating metabolite of the clinical antitumor agent ifosfamide.

PURPOSE: The purpose of this investigation was to determine the base sequence specificity of isophosphoramide mustard (IPM), the alkylating metabolite of ifosfamide, by crosslinking of designed DNA oligomers in comparison with the clinical alkylating agents mechlorethamine (ME) (nitrogen mustard) and phosphoramide mustard (PM), the alkylating metabolite of cyclophosphamide. METHODS: IPM, as well as PM and ME were each reacted with three dodecameric duplexes, which were designed to detect interstrand crosslinking between guanines in 5'-GC-3' (I), 5'-GNC-3' (II) or 5'-GNNC-3' (III) sequences (N = A or T). RESULTS: All three agents preferentially react with 5'-GNC-3' target sequences. The 5'-GNNC-3' target sequence is less reactive by a factor of approximately 2.5- to 10-fold, while 5'-GC-3' is of even lower reactivity. CONCLUSION: These results indicate that all three agents show approximately equal preference for reaction with a 5'-GNC-3' target sequence in spite of the fact that IPM yields a 7-atom crosslink, while the other two agents yield 5-atom crosslinks.

Mapping the binding site of aflatoxin B1 in DNA: molecular modeling of the binding sites for the N(7)-guanine adduct of aflatoxin B1 in different DNA sequences.

Aflatoxin B1 (AFB1), a potent mutagen and carcinogen, forms an adduct exclusively at the N(7) position of guanine, but the structure of this adduct in double stranded DNA is not known. Molecular modeling (using the program, PSFRODO) in conjunction with molecular mechanical calculation (using the program, AMBER) are used to assess the binding modes available to this AFB1 adduct. Two modes appear reasonable; in one the AFB1 moiety is intercalated between the base pair containing the adducted guanine and the adjacent base pair on the 5'-side in reference to the adducted guanine, while in the second it is bound externally in the major groove of DNA. Rotational flexibility appears feasible in the latter providing four, potential binding sites. Molecular modeling reveals that the binding sites around the reactive guanine in different sequences are not uniformly compatible for interaction with AFB1. As the sequence is changed, one particular external binding site would be expected to give a pattern of reactivities that is reasonably consistent with the observed sequence specificity of binding that AFB1 shows in its reaction with DNA (Benasutti, M., Ejadi, S., Whitlow, M. D. and Loechler, E. L. (1988) Biochemistry 27, 472-481). The AFB1 moiety is face-stacked in the major groove with its long axis approximately perpendicular to the helix axis. Favorable interactions are formed between exocyclic amino groups that project into the major groove on cytosines and adenines surrounding the reactive guanine, and oxygens in AFB1; unfavorable interactions involve van der Waals contacts between the methyl group on thymine and the AFB1 moiety. "Some of the sequence specificity of binding data can be rationalized more readily if it is assumed that 5'-GG-3' sequences adopt an A-DNA structure." Based upon molecular modeling/potential energy minimization calculation, it is difficult to predict how reactivity would change in different DNA sequences in the case of the intercalative binding mode; however, several arguments suggest that intercalation might not be favored. From these considerations a model of the structure for the transition state in reaction of AFB1 with DNA is proposed involving one particular external binding site.

Are base substitution and frameshift mutagenesis pathways interrelated? An analysis based upon studies of the frequencies and specificities of mutations induced by the (+)-anti diol epoxide of benzo[a]pyrene.

(+)-anti-B[a]PDE-induced mutagenesis is being investigated, including in a supF gene of the E. coli plasmid pUB3. Based upon various findings a working hypothesis was proposed that the major adduct of (+)-anti-B[a]PDE (formed at N2-Gua) is able to induce different base substitution mutations (e.g., GC-->TA vs. GC-->AT vs. GC-->CG) depending upon its conformation in DNA, which can be influenced by various factors, such as DNA sequence context. Frameshift mutations are also significant and are analyzed herein. In virtually all cases one of three possibilities is observed: (1) some treatments change frameshift and base substitution mutation frequency (MF) in a quantitatively parallel fashion; (2) other treatments, which change frameshift MF, can change base substitution MF in a quantitatively reciprocal fashion; finally, (3) there are treatments that do not change frameshift MF, and also do not change base substitution MF. (Changes can be brought about by SOS induction, differing DNA sequence context, or heating adducted pUB3 prior to transformation. Why different kinds of changes result in (1) vs. (2) vs. (3) is discussed.) Thus, base substitution and frameshift mutagenesis pathways appear to be coupled in some way, which is most easily rationalized if both pathways are interrelated. The simplest mechanism to rationalize this coupling is that a single (+)-anti-B[a]PDE adduct in a single conformation can be bypassed via either a frameshift or a base substitution pathway. The surprising implication is that--although different conformations are likely to be required to induce different base substitution mutations (e.g., GC-->TA vs. GC-->AT; see above)--a single conformation can give rise to either a base substitution or a frameshift mutation. Frameshift and base substitution pathways must eventually diverge, and it is proposed that this is controlled by factors such as DNA sequence context.

An Escherichia coli plasmid-based, mutational system in which supF mutants are selectable: insertion elements dominate the spontaneous spectra.

A new system is described to determine the mutational spectra of mutagens and carcinogens in Escherichia coli; data on a limited number (142) of spontaneous mutants is presented. The mutational assay employs a method to select (rather than screen) for mutations in a supF target gene carried on a plasmid. The E. coli host cells (ES87) are lacI- (am26), and carry the lacZ delta M15 marker for alpha-complementation in beta-galactosidase. When these cells also carry a plasmid, such as pUB3, which contains a wild-type copy of supF and lacZ-alpha, the lactose operon is repressed (off). Furthermore, supF suppression of lacIam26 results in a lactose repressor that has an uninducible, lacIS genotype, which makes the cells unable to grow on lactose minimal plates. In contrast, spontaneous or mutagen-induced supF- mutations in pUB3 prevent suppression of lacIam26 and result in constitutive expression of the lactose operon, which permits growth on lactose minimal plates. The spontaneous mutation frequency in the supF gene is approximately 0.7 and approximately 1.0 x 10(-6) without and with SOS induction, respectively. Spontaneous mutations are dominated by large insertions (67% in SOS-uninduced and 56% in SOS-induced cells), and their frequency of appearance is largely unaffected by SOS induction. These are identified by DNA sequencing to be Insertion Elements; IS1 dominates, but IS4, IS5, gamma-delta and IS10 are also obtained. Large deletions also contribute significantly (19% and 15% for -SOS and +SOS, respectively), where a specific deletion between a 10 base pair direct repeat dominates; the frequency of appearance of these mutations also appears to be unaffected by SOS induction. In contrast, SOS induction increases base pairing mutations (13% and 27% for -SOS and +SOS, respectively). The ES87/pUB3 system has many advantages for determining mutational spectra, including the fact that mutant isolation is fast and simple, and the determination of mutational changes is rapid because of the small size of supF.

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.

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.

Mutagenic specificity of alkylated and oxidized DNA bases as determined by site-specific mutagenesis.

This work demonstrates the use of the tools of site-specific mutagenesis to study the mutagenic activity of two DNA adducts, O6-methylguanine and cis-thymine glycol. The former adduct is one of the methylated bases formed by carcinogenic and mutagenic alkylating agents. It was built into the single-stranded genome of bacteriophage M13 and replicated in Escherichia coli (E. coli). The mutation frequency of O6-methylguanine was 0.4% in physiologically normal cells. In cells in which the repair systems for O6-methylguanine were compromised by challenge with an alkylating agent, the mutation frequency rose to approximately 20%. DNA sequencing revealed that O6-methylguanine induced exclusively G----A transitions, which was most consistent with it pairing with thymine during DNA synthesis. The mutagenic effects also were investigated of cis-thymine glycol isomers, which are major, stable products of ionizing radiation and oxidative damage to DNA. By techniques similar to those employed for the study of O6-methylguanine mutagenesis, a single thymine glycol was situated in an M13 phage genome. The genome was replicated in E. coli that were physiologically normal, induced for SOS functions, or deficient in the nth gene product and, in all cases, the mutagenic processing of thymine glycol in vivo yielded mutant progeny phage at a frequency of 0.3-0.4%. All mutations occurred at the site that originally contained thymine glycol, and all were demonstrated by DNA sequencing to have resulted from targeted T----C transitions. These data suggest that thymine glycol pairs with guanine during replication.

Rotation about the C6-O6 bond in O6-methylguanine: the syn and anti conformers can be of similar energies in duplex DNA as estimated by molecular modeling techniques.

O6-Methylguanine (O6MeGua) is generally regarded as the most important premutagenic lesion formed from carcinogenic methylating agents, so its structure and mechanism of action have received considerable attention. Two conformations for the methyl group in O6MeGua are possible: in one the methyl group is syn with respect to the N(1)-position of the purine, points into the helix, and disrupts hydrogen bonding potential; in the second the methyl group is anti with respect to the N1-position of the purine, and points into the major groove. Syn-O6MeGua has been observed when paired with thymine in duplex DNA as determined by NMR, while anti-O6MeGua has been observed when paired with thymine in X-ray diffraction studies. Herein, molecular modeling/computational chemistry is used to evaluate this apparent discrepancy. [N6-Methyladenine (N6MeAde) was also studied, because it is isoelectronic with O6MeGua, and because more information is available about its energetics. Syn-N6MeAde is computed to be favored in small molecules; however, the fraction of the anti-conformer is computed to be approximately 7%, which agrees well with experimentally determined values (4-12%). In contrast, the anti conformation for N6MeAde is calculated to be favored in duplex DNA, which is consistent with what has been observed experimentally using both NMR and X-ray diffraction techniques. The agreement between the calculated and experimentally determined results with N6MeAde suggest that the methods are reasonable.] For O6MeGua, a syn/anti ratio of approximately 10(3) is computed for small molecules. In duplex DNA, syn-O6MeGua is computed to be favored, but the anti-conformer is less than approximately 1 kcal/mol higher in energy. Whether syn- or anti-O6MeGua predominates may depend upon sequence context, as well as environmental factors. The comparison between O6MeGua and N6MeAde suggests a rationale for the puzzling observation that O6MeGua appears to be a cytotoxic lesion in eukaryotic, but not prokaryotic, cells.

How are potent bulky carcinogens able to induce such a diverse array of mutations?

Mutations induced by activated benzo[a]pyrene ((+)-anti-B[a]PDE) in Escherichia coli are being investigated, by using both random and adduct-site-specific mutagenesis approaches. A working hypothesis was proposed that the major adduct of (+)-anti-B[a]PDE (formed at N2-Gua) is able to induce different base-substitution mutations (e.g., GC-->TA vs. GC-->AT) depending upon its conformation in DNA, which can be influenced by various factors, notably DNA sequence context. Frameshift mutations are also common with (+)-anti-B[a]PDE, and other work suggested that the frameshift and base-substitution mutagenesis pathways are coupled. The simplest hypothesis to rationalize this interrelationship is that a single (+)-anti-B[a]PDE adduct in a single conformation can be bypassed via either a frameshift or a base-substitution pathway. This counterintuitive notion can be reconciled if there are two different kinds of conformations on the pathway to mutagenesis: a class I conformation, which is the initial conformation of a DNA adduct in double-stranded DNA before its encounter with a DNA polymerase, and a class II conformation, which is the conformation that forms at a single-strand/double-strand DNA junction during replication by a DNA polymerase. Thus, GC-->TA and GC-->AT mutations may be induced by different class I conformations, whereas base substitution and frameshift mutations may be induced by the same class I conformation but by different class II conformations. The pathway of mutagenesis would be dictated by the relevant class I and II conformations, which in turn would be controlled by various factors, notably DNA sequence context.

Adduct-induced base-shifts: a mechanism by which the adducts of bulky carcinogens might induce mutations.

Most carcinogens have been shown to be mutagens, and DNA adducts are formed when mutagenic/carcinogenic substances react with DNA. It is generally believed these adducts (or their derivatives) induce misreplication events that result in mutations. Many of the more potently mutagenic substances are bulky and three-dimensionally complex, such as the polycyclic aromatic hydrocarbons, aromatic amines, and aflatoxins; little is known about the mechanisms by which they induce mutations. Several theories exist and herein an additional mechanism is proposed by which bulky adducts might induce mutations at GC base pairs. Molecular modeling in conjunction with molecular mechanical calculation is used to assess if the mutagen/carcinogen moiety of the adduct might be able to shift the position of the base moiety of the adduct in such a way that misreplication events might be facilitated. This mechanism is referred to as adduct-induced base-shift, and two classes appeared possible; adduct-induced base-wobble and adduct-induced base-rotation. The latter has been proposed previously. By adduct-induced, base-wobble, the mutagen/carcinogen moiety of the adduct induces a shift in the position of the base moiety of the adduct with respect to the helix axis, which might facilitate mispairing events that are reminisent of non-Watson/Crick pairing that occurs at the wobble base of tRNA during translation. For example, in some guanine adducts, the guanine appears more thymine-like, which might facilitate G.A mispairing and thereby ultimately GC to TA transversion mutations. Adduct-induced base-rotation involves the rotation of the adducted base from the anti to the syn conformation and a variety of mispairing events might result.

The role of adduct site-specific mutagenesis in understanding how carcinogen-DNA adducts cause mutations: perspective, prospects and problems.

Usually, a particular mutagen/carcinogen forms adducts at many sites in DNA, making it impossible to determine which type of adduct causes which mutation and why. Adduct site-specific mutagenesis studies, in which a single adduct is built into a vector, can be used to overcome this problem. The adduct can be situated in double-stranded DNA, single-stranded DNA or in a single-stranded gap, and the benefit and concerns associated with each are addressed. An adduct site-specific study is most useful when it is compared to a mutagenesis study with its corresponding mutagen/carcinogen. Mutations induced by a particular mutagen/carcinogen can be influenced by DNA sequence context, mutagen/carcinogen dose (and other changes in conditions), level of SOS induction, cell type and other factors. Thus, it is important to match the conditions of the adduct study versus the mutagen/carcinogen study as closely as possible. DNA sequence context can profoundly affect the quantitative and qualitative pattern of adduct mutagenesis, which is addressed. In vitro studies with DNA polymerases, frameshift mutagenesis and semi-targeted mutagenesis, whereby a mutation is induced near but not at the site of the adduct, are each discussed. Finally, the relationship between structural studies on adducts and mutagenesis is considered.

Molecular modeling of the conformational complexity of (+)-anti-B[a]PDE-adducted DNA using simulated annealing.

Benzo[a]pyrene (B[a]P), a potent mutagen/carcinogen, reacts with DNA following metabolism to its corresponding (+)-anti-7,8-diol-9,10-epoxide [(+)-anti-B[a]PDE], giving a major adduct (+)-trans-anti-B[a]P-N2-dG. Evidence suggests that this adduct is responsible for most of the different kinds of mutations (e.g. G-->T, G-->A, etc.) induced by (+)-anti-B[a]PDE, raising the question of how can a single adduct cause many different kinds of mutations? One hypothesis is that different mutations are induced depending upon the conformation of this adduct when bypassed during DNA replication. If true, then it becomes imperative to explore different reasonable conformations for this adduct. Herein a simulated annealing protocol is employed to study the conformation of (+)-trans-anti-B[a]P-N2-dG with the B[a]P moiety in the minor groove and pointing toward the base on its 5'-side in a 5'-CGC-3' sequence context in duplex DNA. This conformation and sequence were chosen because there is a structure derived from NMR constraints for comparison. A four step procedure is followed: the adduct is docked in canonical B-DNA, after which the structure is subjected to an initial conjugate gradient minimization, followed by simulated annealing and a final conjugate gradient minimization. The quality and final energy of structures is assessed as a function of changes in six parameters, including the length of the DNA helix, the initial annealing temperature (T0), the annealing time (t), the molecular dynamics time step (tau) and two other parameters. While there is no single set of optimum parameters, reasonable low energy structures were obtained using the values t approximately 40 ps (or longer), T0 approximately 750 K and tau approximately 1.0 fs with a helix length of 7 bp. The structures that emerge all retain the basic features of the input structure, being B-DNA-like with the B[a]P moiety in the minor groove pointing toward the base on the 5'-side. However, within this broad category there are at least six subclasses of structures, of which four have lowest energy members that differ by < approximately 5 kcal/mol. The fact that a variety of distinct but related structures emerge from a single starting structure as this parameter set is varied suggests that the use of a large but manageable number of simulated annealing runs should be considered in the search for a cohort of related structures. This is especially important given that this breadth of potentially relevant structures of approximately the same energy may indeed be relevant to the hypothesis that different mutations arise from a single adduct in different conformations.

Mutagenesis by the (+)-anti-diol epoxide of benzo[a]pyrene: what controls mutagenic specificity?

Mutagenesis by (+)-anti-benzo[a]pyrene-7,8-dihydrodiol-9,10-epoxide [(+)-anti-B[a]PDE], an important mutagenic/carcinogenic metabolite of benzo[a]pyrene (B[a]P), is being studied in order to understand the factors that influence mutagenesis both quantitatively and qualitatively. A new mutational system, which permits the selection of supF- mutations in an Escherichia coli plasmid, pUB3, was used. The work described herein is an extension of previous work, which involved plasmid adduction and then immediate transformation (Rodriguez & Loechler, 1993), and began with the observation that mutation frequency (MF) decreased approximately 2-fold when the (+)-anti-B[a]PDE-adducted plasmid pUB3 is either (1) frozen and then thawed prior to transformation or (2) heated at 80 degrees C for 10 min prior to transformation. Several results suggest that this decrease is not due to the loss of labile adducts. To begin to understand this phenomenon, the mutagenic spectra are compared for (+)-anti-B[a]PDE in supF for the unheated (187 mutants), the freeze/thawed (134 mutants), and the heated (254 mutants) samples. In general, freeze/thawing and heating cause a decrease in all classes of mutations. Considering substitution mutations at G.C base pairs, which predominate, the mutagenic specificity for the combined data sets is GC-->TA (57%), GC-->AT (23%), and GC-->CG (20%). This raises the question, how does (+)-anti-B[a]PDE generate this complex mutagenic specificity, which contrasts with the situation for, e.g., simple methylating agents? One factor is that mutagenic specificity at a particular guanine residue can be influenced by the base on its immediate 5'-side, most notably where mutations are virtually exclusively restricted to GC-->TA in 5'-TG-3' sequence contexts. One unexpected finding may provide additional insight. G115 in supF, which is the major hot spot for base-pairing mutagenesis, is the only site where the qualitative pattern of mutagenesis is significantly affected by heating the (+)-anti-B[a]PDE-adducted plasmid prior to transformation. Without heating, G115-->T mutations predominate, but following heating there is a statistically significant increase in the fraction of G115-->A and G115-->C mutations. The most likely model to explain this and other results is (1) a particular DNA adduct can adopt multiple conformations, (2) the conformation adopted by an adduct can be influenced by various factors, including DNA sequence context, as well as heating and freeze/thawing, and (3) each of these conformations can cause a different pattern of mutation.(ABSTRACT TRUNCATED AT 400 WORDS)

Mutational specificity of the (+)-anti-diol epoxide of benzo[a]pyrene in a supF gene of an Escherichia coli plasmid: DNA sequence context influences hotspots, mutagenic specificity and the extent of SOS enhancement of mutagenesis.

Mutagenesis by the suspected major mutagenic metabolite of activated benzo[a]pyrene, which is (+)-anti-BPDE, was analyzed with a new system, which permits the selection of supF- mutations in an Escherichia coli plasmid using lactose minimal plates. (+)-anti-BPDE enhances base pairing mutations--principally at G:C base pairs, frameshift mutations and large deletions. Frameshift mutagenesis principally involves deletions and insertions of a single G:C base pair in runs of G:C base pairs. Base pairing mutations are significantly enhanced by SOS induction, especially GC-->TA mutations. Nearest neighbor analysis was performed assuming that a guanine (underlined) is being mutated, and (+)-anti-BPDE base pairing mutagenesis is enhanced by SOS induction in 5'-(A/T)G-3' sequences approximately 4-fold more than in 5'-(G/C)G-3' sequences, and in 5'-G(C/G)-3' sequences approximately 4-fold more than in 5'-G(A/T)-3' sequences; this is discussed. The influence of sequence context on quantitative aspects of (+)-anti-BPDE mutagenesis is considered, and hotspots are found at most, but not all, 5'-GG-3' sequences. The influence of sequence context on qualitative aspects of (+)-anti-BPDE mutagenesis (i.e. mutagenic specificity) is also considered. For example, the sequences, 5'-AG-3', 5'-CG-3' and 5'-GG-3', all have examples of G-->T, G-->A and G-->C mutations, while in the sequence, 5'-TG-3', only G-->T mutations have been detected. (The latter finding correlates with a recent site-specific study on the major adduct of (+)-anti-BPDE formed at N2-Gua in a 5'-TG-3' context, where G-->T mutations predominated [Carcinogenesis (1992) 13, 1415-1425].) These results suggest that sequence context plays a role in defining the kind of mutation (i.e. GC-->TA versus GC-->AT versus GC-->CG) induced by (+)-anti-BPDE, where the base on the 5' side of the guanine undergoing mutation seems to be influential. The most likely model for this is that sequence context influences adduct conformation, which controls mutagenic specificity.