How Y-Family DNA polymerase IV is more accurate than Dpo4 at dCTP insertion opposite an N2-dG adduct of benzo[a]pyrene

To bypass DNA damage, cells have Y-Family DNA polymerases (DNAPs). One Y-Family-class includes DNAP κ and DNAP IV, which accurately insert dCTP opposite N(2)-dG adducts, including from the carcinogen benzo[a]pyrene (BP). Another class includes DNAP η and DNAP V, which insert accurately opposite UV-damage, but inaccurately opposite BP-N(2)-dG. To investigate structural differences between Y-Family-classes, regions are swapped between DNAP IV (a κ/IV-class-member) and Dpo4 (a η/V-class-member); the kinetic consequences are evaluated via primer-extension studies with a BP-N(2)-dG-containing template. Four key structural elements are revealed. (1) Y-Family DNAPs have discreet non-covalent contacts between their little finger-domain (LF-Domain) and their catalytic core-domain (CC-Domain), which we call "non-covalent bridges" (NCBs). Arg37 and Arg38 in DNAP IV's CC-Domain near the active site form a non-covalent bridge (AS-NCB) by interacting with Glu251 and Asp252, respectively, in DNAP IV's LF-Domain. Without these interactions dATP/dGTP/dTTP misinsertions increase. DNAP IV's AS-NCB suppresses misinsertions better than Dpo4's equivalent AS-NCB. (2) DNAP IV also suppresses dATP/dGTP/dTTP misinsertions via a second non-covalent bridge, which is ∼8Å from the active site (Distal-NCB). Dpo4 has no Distal-NCB, rendering it inferior at dATP/dGTP/dTTP suppression. (3) dCTP insertion is facilitated by the larger minor groove opening near the active site in DNAP IV versus Dpo4, which is sensible given that Watson/Crick-like [dCTP:BP-N(2)-dG] pairing requires the BP-moiety to be in the minor groove. (4) Compared to Dpo4, DNAP IV has a smaller major groove opening, which suppresses dGTP misinsertion, implying BP-N(2)-dG bulk in the major groove during Hoogsteen syn-adduct-dG:dGTP pairing. In summary, DNAP IV has a large minor groove opening to enhance dCTP insertion, a plugged major groove opening to suppress dGTP misinsertion, and two non-covalent bridges (near and distal to the active site) to suppress dATP/dGTP/dTTP misinsertions; collectively these four structural features enhance DNAP IV's dNTP insertion fidelity opposite a BP-N(2)-dG adduct compared to Dpo4.

A method to accurately quantitate intensities of (32)P-DNA bands when multiple bands appear in a single lane of a gel is used to study dNTP insertion opposite a benzo[a]pyrene-dG adduct by Sulfolobus DNA polymerases Dpo4 and Dbh

Quantitating relative (32)P-band intensity in gels is desired, e.g., to study primer-extension kinetics of DNA polymerases (DNAPs). Following imaging, multiple (32)P-bands are often present in lanes. Though individual bands appear by eye to be simple and well-resolved, scanning reveals they are actually skewed-Gaussian in shape and neighboring bands are overlapping, which complicates quantitation, because slower migrating bands often have considerable contributions from the trailing edges of faster migrating bands. A method is described to accurately quantitate adjacent (32)P-bands, which relies on having a standard: a simple skewed-Gaussian curve from an analogous pure, single-component band (e.g., primer alone). This single-component scan/curve is superimposed on its corresponding band in an experimentally determined scan/curve containing multiple bands (e.g., generated in a primer-extension reaction); intensity exceeding the single-component scan/curve is attributed to other components (e.g., insertion products). Relative areas/intensities are determined via pixel analysis, from which relative molarity of components is computed. Common software is used. Commonly used alternative methods (e.g., drawing boxes around bands) are shown to be less accurate. Our method was used to study kinetics of dNTP primer-extension opposite a benzo[a]pyrene-N(2)-dG-adduct with four DNAPs, including Sulfolobus solfataricus Dpo4 and Sulfolobus acidocaldarius Dbh. Vmax/Km is similar for correct dCTP insertion with Dpo4 and Dbh. Compared to Dpo4, Dbh misinsertion is slower for dATP (∼20-fold), dGTP (∼110-fold) and dTTP (∼6-fold), due to decreases in Vmax. These findings provide support that Dbh is in the same Y-Family DNAP class as eukaryotic DNAP κ and bacterial DNAP IV, which accurately bypass N(2)-dG adducts, as well as establish the scan-method described herein as an accurate method to quantitate relative intensity of overlapping bands in a single lane, whether generated from (32)P-signals or by other means (e.g., staining).

Structural model of the Y-Family DNA polymerase V/RecA mutasome

To synthesize past DNA damaged by chemicals or radiation, cells have lesion bypass DNA polymerases (DNAPs), most of which are in the Y-Family. One class of Y-Family DNAPs includes DNAP η in eukaryotes and DNAP V in bacteria, which have low fidelity when replicating undamaged DNA. In Escherchia coli, DNAP V is carefully regulated to insure it is active for lesion bypass only, and one mode of regulation involves interaction of the polymerase subunit (UmuC) and two regulatory subunits (UmuD') with a RecA-filament bound to ss-DNA. Taking a docking approach, ∼150,000 unique orientations involving UmuC, UmuD' and RecA were evaluated to generate models, one of which was judged best able to rationalize the following published findings. (1) In the UmuD'(2)C/RecA-filament model, R64-UmuC interacts with S117-RecA, which is known to be at the UmuC/RecA interface. (2) At the model's UmuC/RecA interface, UmuC has three basic amino acids (K59/R63/R64) that anchor it to RecA. No other Y-Family DNAP has three basic amino acids clustered in this region, making it a plausible site for UmuC to form its unique interaction with RecA. (3) In the model, residues N32/N33/D34 of UmuC form a second interface with RecA, which is consistent with published findings. (4) Active UmuD' is generated when 24 amino acids in the N-terminal tail of UmuD are proteolyzed, which occurs when UmuD(2)C binds the RecA-filament. When UmuD is included in an UmuD(2)C/RecA-filament model, plausible UmuD/RecA contacts guide the UmuD cleavage site (C24/G25) into the UmuD proteolysis active site (S60/K97). One contact involves E11-UmuD interacting with R243-RecA, where the latter is known to be important for UmuD cleavage. (5) The UmuD(2)C/RecA-filament model rationalizes published findings that at least some UmuD-to-UmuD' cleavage occurs intermolecularly. (6) Active DNAP V is known to be the heterotetramer UmuD'(2)C/RecA, a model of which can be generated by a simple rearrangement of the RecA monomer at the 3'-end of the RecA-filament. The rearranged UmuD'(2)C/RecA model rationalizes published findings about UmuD' residues in proximity to RecA. In summary, docking and molecular simulations are used to develop an UmuD'(2)C/RecA model, whose structure rationalizes much of the known properties of the active form of DNA polymerase V.

Architecture of y-family DNA polymerases relevant to translesion DNA synthesis as revealed in structural and molecular modeling studies

DNA adducts, which block replicative DNA polymerases (DNAPs), are often bypassed by lesion-bypass DNAPs, which are mostly in the Y-Family. Y-Family DNAPs can do non-mutagenic or mutagenic dNTP insertion, and understanding this difference is important, because mutations transform normal into tumorigenic cells. Y-Family DNAP architecture that dictates mechanism, as revealed in structural and modeling studies, is considered. Steps from adduct blockage of replicative DNAPs, to bypass by a lesion-bypass DNAP, to resumption of synthesis by a replicative DNAP are described. Catalytic steps and protein conformational changes are considered. One adduct is analyzed in greater detail: the major benzo[a]pyrene adduct (B[a]P-N²-dG), which is bypassed non-mutagenically (dCTP insertion) by Y-family DNAPs in the IV/κ-class and mutagenically (dATP insertion) by V/η-class Y-Family DNAPs. Important architectural differences between IV/κ-class versus V/η-class DNAPs are discussed, including insights gained by analyzing ~400 sequences each for bacterial DNAPs IV and V, along with sequences from eukaryotic DNAPs kappa, eta and iota. The little finger domains of Y-Family DNAPs do not show sequence conservation; however, their structures are remarkably similar due to the presence of a core of hydrophobic amino acids, whose exact identity is less important than the hydrophobic amino acid spacing.

Amino acid architecture that influences dNTP insertion efficiency in Y-family DNA polymerase V of E. coli

Y-family DNA polymerases (DNAPs) are often required in cells to synthesize past DNA-containing lesions, such as [+ta]-B[a]P-N(2)-dG, which is the major adduct of the potent mutagen/carcinogen benzo[a]pyrene. The current model for the non-mutagenic pathway in Escherichia coli involves DNAP IV inserting deoxycytidine triphosphate opposite [+ta]-B[a]P-N(2)-dG and DNAP V doing the next step(s), extension. We are investigating what structural differences in these related Y-family DNAPs dictate their functional differences. X-ray structures of Y-family DNAPs reveal a number of interesting features in the vicinity of the active site, including (1) the "roof-amino acid" (roof-aa), which is the amino acid that lies above the nucleobase of the deoxynucleotide triphosphate (dNTP) and is expected to play a role in dNTP insertion efficiency, and (2) a cluster of three amino acids, including the roof-aa, which anchors the base of a loop, whose detailed structure dictates several important mechanistic functions. Since no X-ray structures existed for UmuC (the polymerase subunit of DNAP V) or DNAP IV, we previously built molecular models. Herein, we test the accuracy of our UmuC(V) model by investigating how amino acid replacement mutants affect lesion bypass efficiency. A ssM13 vector containing a single [+ta]-B[a]P-N(2)-dG is transformed into E. coli carrying mutations at I38, which is the roof-aa in our UmuC(V) model, and output progeny vector yield is monitored as a measure of the relative efficiency of the non-mutagenic pathway. Findings show that (1) the roof-aa is almost certainly I38, whose beta-carbon branching R-group is key for optimal activity, and (2) I38/A39/V29 form a hydrophobic cluster that anchors an important mechanistic loop, aa29-39. In addition, bypass efficiency is significantly lower both for the I38A mutation of the roof-aa and for the adjacent A39T mutation; however, the I38A/A39T double mutant is almost as active as wild-type UmuC(V), which probably reflects the following. Y-family DNAPs fall into several classes with respect to the [roof-aa/next amino acid]: one class has [isoleucine/alanine] and includes UmuC(V) and DNAP eta (from many species), while the second class has [alanine (or serine)/threonine] and includes DNAP IV, DNAP kappa (from many species), and Dpo4. Thus, the high activity of the I38A/A39T double mutant probably arises because UmuC(V) was converted from the V/eta class to the IV/kappa class with respect to the [roof-aa/next amino acid]. Structural and mechanistic aspects of these two classes of Y-family DNAPs are discussed.

Y-Family DNA polymerases may use two different dNTP shapes for insertion: a hypothesis and its implications

Chemicals and radiation can damage DNA leading to the formation of adducts/lesions, which - if not removed by DNA repair pathways - usually block replicative DNA polymerases (DNAPs). To overcome such potentially lethal blockage, cells have lesion bypass DNAPs, which are often in the Y-Family and include several classes. One class includes human DNAP kappa and E. coli DNAP IV, and they insert dCTP in the non-mutagenic pathway opposite [+ta]-B[a]P-N(2)-dG, which is the major adduct formed by the environmental carcinogen benzo[a]pyrene. Another class includes hDNAP eta and ecDNAP V, and they insert dATP opposite [+ta]-B[a]P-N(2)-dG in the dominant G-->T mutagenic pathway. Herein we develop a hypothesis for why the IV/kappa-class preferentially does cellular dCTP insertion. On the minor groove side of the active site, Y-Family DNAPs have a cleft/hole that can be analyzed based on an analogy to a "chimney." Our models of DNAP IV show a large chimney opening from which the pyrene of [+ta]-B[a]P-N(2)-dG can protrude, which allows canonical adduct-dG:dCTP pairing. In contrast, our models of DNAP V have small chimney openings that forces adduct-dG downward in the active site such that canonical adduct-dG:dCTP pairing is not possible. Based on X-ray structures, sequence alignment and our modeled structures of Y-Family DNAPs, chimney opening size seems primarily controlled by one amino acid ("flue-handle"), which dictates whether nearby amino acids ("flue") plug the chimney or not. Based on this analysis, a correlation is apparent: the flue is closed in V/eta-class DNAPs giving small chimney openings, while the flue is open for the IV/kappa-class giving large chimney openings. Secondarily, a hypothesis is developed for why the V/eta-class might preferentially do cellular dATP insertion opposite [+ta]-B[a]P-N(2)-dG: the small chimney forces adduct-dG lower in the active site, possibly leading to catalysis using a non-canonical dNTP shape that permits syn-adenine:adduct-dG base pairing. In summary, a hypothesize is developed that the pyrene moiety of [+ta]-B[a]P-N(2)-dG protrudes from the large chimney opening of DNAP IV, thus permitting canonical dCTP:adduct-dG pairing, while the small chimney opening of DNAP V forces [+ta]-B[a]P-N(2)-dG lower down in the active site, in which syn-adenine can pair with adduct-dG via a non-canonical dNTP shape.

Environmental and occupational causes of cancer: new evidence 2005-2007

What do we currently know about the occupational and environmental causes of cancer? As of 2007, the International Agency for Research on Cancer (IARC) identified 415 known or suspected carcinogens. Cancer arises through an extremely complicated web of multiple causes, and we will likely never know the full range of agents or combinations of agents. We do know that preventing exposure to individual carcinogens prevents the disease. Declines in cancer rates-such as the drop in male lung cancer cases from the reduction in tobacco smoking or the drop in bladder cancer among cohorts of dye workers from the elimination of exposure to specific aromatic amines-provides evidence that preventing cancer is possible when we act on what we know. Although the overall age-adjusted cancer incidence rates in the United States among both men and women have declined in the last decade, the rates of several types of cancers are on the rise; some of which are linked to environmental and occupational exposures. This report chronicles the most recent epidemiologic evidence linking occupational and environmental exposures with cancer. Peer-reviewed scientific studies published from January 2005 to June 2007 were reviewed, supplementing our state-of-the-evidence report published in September 2005. Despite weaknesses in certain individual studies, we consider the evidence linking the increased risk of several types of cancer with specific exposures somewhat strengthened by recent publications, among them brain cancer from exposure to non-ionizing radiation, particularly radiofrequency fields emitted by mobile telephones; breast cancer from exposure to the pesticide dichlorodiphenyltrichloroethane (DDT) before puberty; leukemia from exposure to 1,3-butadiene; lung cancer from exposure to air pollution; non-Hodgkin's lymphoma (NHL) from exposure to pesticides and solvents; and prostate cancer from exposure to pesticides, polyaromatic hydrocarbons (PAHs), and metal working fluids or mineral oils. In addition to NHL and prostate cancer, early findings from the National Institutes of Health Agricultural Health Study suggest that several additional cancers may be linked to a variety of pesticides. Our report also briefly describes the toxicological evidence related to the carcinogenic effect of specific chemicals and mechanisms that are difficult to study in humans, namely exposures to bis-phenol A and epigenetic, trans-generational effects. To underscore the multi-factorial, multi-stage nature of cancer, we also present a technical description of cancer causation summarizing current knowledge in molecular biology. We argue for a new cancer prevention paradigm, one based on an understanding that cancer is ultimately caused by multiple interacting factors rather than a paradigm based on dubious attributable fractions. This new cancer prevention paradigm demands that we limit exposure to avoidable environmental and occupational carcinogens, in combination with additional important risk factors like diet and lifestyle. The research literature related to environmental and occupational causes of cancer is constantly growing, and future updates will be carried out in light of new biological understanding of the mechanisms and new methods for studying exposures in human populations. The current state of knowledge is sufficient to compel us to act on what we know. We repeat the call of ecologist Sandra Steingraber: "From the right to know and the duty to inquire flows the obligation to act."

Molecular modeling benzo[a]pyrene N2-dG adducts in the two overlapping active sites of the Y-family DNA polymerase Dpo4

The potent, ubiquitous environmental mutagen/carcinogen benzo[a]pyrene (B[a]P) induces a single major adduct [+ta]-B[a]P-N2-dG, whose bypass in most cases results in either no mutation (dCTP insertion) or a G-->T mutation (dATP insertion). Translesion synthesis (TLS) of [+ta]-B[a]P-N2-dG generally requires DNA polymerases (DNAPs) in the Y-family, which exist in cells to bypass DNA damage caused by chemicals and radiation. A molecular dynamics (MD) study is described with dCTP opposite [+ta]-B[a]P-N2-dG in Dpo4, which is the best studied Y-family DNAP from a structural point of view. Two orientations of B[a]P-N2-dG (BPmi5 and BPmi3) are considered, along with two orientations of the dCTP (AS1 and AS2), as outlined next. Based on NMR studies, the pyrene moiety of B[a]P-N2-dG is in the minor groove, when paired with dC, and can point toward either the base on the 5'-side (BPmi5) or the 3'-side (BPmi3). Based on published X-ray structures, Dpo4 appears to have two partially overlapping active sites. The architecture of active site 1 (AS1) is similar to all other families of DNAPs (e.g., the shape of the dNTP). Active site 2 (AS2), however, is non-canonical (e.g., the beta- and gamma-phosphates in AS2 are approximately where the alpha- and beta-phosphates are in AS1). In the Dpo4 models generated herein, using the BPmi3 orientation the pyrene moiety of [+ta]-B[a]P-N2-dG points toward the duplex region of the DNA, and is accommodated without distortions in AS1, but with distortions in AS2. Considering the BPmi5 orientation, the pyrene moiety points toward the ss-region of DNA in Dpo4, and sits in a hole defined by the fingers and little fingers domain ("chimney"); BPmi5 is accommodated in AS2 without significant distortions, but poorly in AS1. In summary, when dCTP is paired with [+ta]-B[a]P-N2-dG in the two overlapping active sites in Dpo4, the pyrene in the BPmi3 orientation is accommodated better in active site 1 (AS1), while the pyrene in the BPmi5 orientation is accommodated better in AS2. Finally, we discuss why Y-family DNAPs might have two catalytic active sites.

Free-energy perturbation methods to study structure and energetics of DNA adducts: results for the major N2-dG adduct of benzo[a]pyrene in two conformations and different sequence contexts

The potent mutagen/carcinogen benzo[a]pyrene (B[a]P) is activated to (+)-anti-B[a]PDE, which induces a variety of mutations (e.g., G --> T, G --> A, etc.) via its major adduct [+ta]-B[a]P-N2-dG. One hypothesis is that adducts (such as [+ta]-B[a]P-N2-dG) induce different mutations via different conformations, probably when replicated by different lesion-bypass DNA polymerases (DNAPs). We showed that Escherichia coli DNAP V was responsible for G --> T mutations with [+ta]-B[a]P-N2-dG in a 5'-TGT sequence (Yin et al., (2004) DNA Repair 3, 323), so we wish to study conformations of this adduct/sequence context by molecular modeling. The development of a CHARMM-based molecular dynamics (MD) simulations protocol with free-energy calculations in the presence of solvent and counterions is described. A representative base-pairing and base-displaced conformation of [+ta]-B[a]P-N2-dG in the 5'-TGT sequence are used: (1) BPmi5, which has the B[a]P moiety in the minor groove pointing toward the base on the 5'-side of the adduct, and (2) Gma5, which has the B[a]P moiety stacked with the surrounding base pairs and the dG moiety displaced into the major groove. The MD output structures are reasonable when compared to known NMR structures. Changes in DNA sequence context dramatically affect the biological consequences (e.g., mutagenesis) of [+ta]-B[a]P-N2-dG. Consequently, we also developed a MD-based free-energy perturbation (FEP) protocol to study DNA sequence changes. FEP involves the gradual "fading-out" of atoms in a starting structure (A) and "fading-in" of atoms in a final structure (B), which allows a realistic assessment of the energetic and structural changes when two structures A and B are closely related. Two DNA sequence changes are described: (1) 5'-TGT --> 5'-TGG, which involves two steps [T:A --> T:C --> G:C], and (2) 5'-TGT --> 5'-TGC, which involves three steps [T:A --> T:2AP --> C:2AP --> C:G], where 2AP (2-aminopurine) is included, because T:2AP and C:2AP retain more-or-less normal pairing orientations between complementary bases. FEP is also used to evaluate the impact that a 5'-TGT to 5'-UGT sequence change might have on mutagenesis with [+ta]-B[a]P-N2-dG. In summary, we developed (1) a CHARMM-based molecular dynamics (MD) simulations protocol with free-energy calculations in the presence of solvent and counterions to study B[a]P-N2-dG adducts in DNA duplexes, and (2) a MD-based free-energy perturbation (FEP) protocol to study DNA sequence context changes around B[a]P-N2-dG adducts.

Genetic effects of oxidative DNA damages: comparative mutagenesis of the imidazole ring-opened formamidopyrimidines (Fapy lesions) and 8-oxo-purines in simian kidney cells

Fapy.dG and 8-oxo-7,8-dihydro-2'-deoxyguanosine (8-oxo-dG) are formed in DNA by hydroxyl radical damage. In order to study replication past these lesions in cells, we constructed a single-stranded shuttle vector containing the lesion in 5'-TGT and 5'-TGA sequence contexts. Replication of the modified vector in simian kidney (COS-7) cells showed that Fapy.dG is mutagenic inducing primarily targeted Fapy.G-->T transversions. In the 5'-TGT sequence mutational frequency of Fapy.dG was approximately 30%, whereas in the 5'-TGA sequence it was approximately 8%. In parallel studies 8-oxo-dG was found to be slightly less mutagenic than Fapy.dG, though it also exhibited a similar context effect: 4-fold G-->T transversions (24% versus 6%) occurred in the 5'-TGT sequence relative to 5'-TGA. To investigate a possible structural basis for the higher G-->T mutations induced by both lesions when their 3' neighbor was T, we carried out a molecular modeling investigation in the active site of DNA polymerase beta, which is known to incorporate both dCTP (no mutation) and dATP (G-->T substitution) opposite 8-oxo-G. In pol beta, the syn-8-oxo-G:dATP pair showed greater stacking with the 3'-T:A base pair in the 5'-TGT sequence compared with the 3'-A:T in the 5'-TGA sequence, whereas stacking for the anti-8-oxo-G:dCTP pair was similar in both 5'-TGT and 5'-TGA sequences. Similarly, syn-Fapy.G:dATP pairing showed greater stacking in the 5'-TGT sequence compared with the 5'-TGA sequence, while stacking for anti-Fapy.G:dCTP pairs was similar in the two sequences. Thus, for both lesions less efficient base stacking between the lesion:dATP pair and the 3'-A:T base pair in the 5'-TGA sequence might cause lower G-->T mutational frequencies in the 5'-TGA sequence compared to 5'-TGT. The corresponding lesions derived from 2'-deoxyadenosine, Fapy.dA and 8-oxo-dA, were not detectably mutagenic in the 5'-TAT sequence, and were only weakly mutagenic (<1%) in the 5'-TAA sequence context, where both lesions induced targeted A-->C transversions. To our knowledge this is the first investigation using extrachromosomal probes containing a Fapy.dG or Fapy.dA site-specifically incorporated, which showed unequivocally that in simian kidney cells Fapy.G-->T substitutions occur at a higher frequency than 8-oxo-G-->T and that Fapy.dA is very weakly mutagenic, as is 8-oxo-dA.

Mirror image stereoisomers of the major benzo[a]pyrene N2-dG adduct are bypassed by different lesion-bypass DNA polymerases in E. coli

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., G-to-T, G-to-A, -1 frameshifts, etc.) via its major adduct [+ta]-B[a]P-N2-dG. We recently showed that the dominant G-to-T mutation depends on DNA polymerase V (DNAP V), but not DNAPs IV or II, when studied in a 5'-TG sequence in E. coli. Herein we investigate what DNAPs are responsible for non-mutagenic bypass with [+ta]-B[a]P-N2-dG, along with its mirror image adduct [-ta]-B[a]P-N2-dG. Each adduct is built into a 5'-TG sequence in a single stranded M13 phage vector, which is then transformed into eight different E. coli strains containing all combinations of proficiency and deficiency in the three lesion-bypass DNAPs II, IV and V. Based on M13 progeny output, non-mutagenic bypass with [-ta]-B[a]P-N2-dG depends on DNAP IV. In contrast, non-mutagenic bypass with [+ta]-B[a]P-N2-dG depends on both DNAPs IV and V, where arguments suggest that DNAP IV is involved in dCTP insertion, while DNAP V is involved in extension of the adduct-G:C base pair. Numerous findings indicate that DNAP II has a slight inhibitory effect on the bypass of [+ta]- and [-ta]-B[a]P-N2-dG in the case of both DNAPs IV and V. In conclusion, for efficient non-mutagenic bypass (dCTP insertion) in E. coli, [+ta]-B[a]P-N2-dG requires DNAPs IV and V, [-ta]-B[a]P-N2-dG requires only DNAP IV, while DNAP II is inhibitory to both, and experiments to investigate these differences should provide insights into the mechanism and purpose of these lesion-bypass DNAPs.

Homology modeling of four Y-family, lesion-bypass DNA polymerases: the case that E. coli Pol IV and human Pol kappa are orthologs, and E. coli Pol V and human Pol eta are orthologs

Y-family DNA polymerases (DNAPs) are a superfamily of evolutionarily related proteins that exist in cells to bypass DNA damage caused by both radiation and chemicals. Cells have multiple Y-family DNAPs, presumably to conduct translesion synthesis (TLS) on DNA lesions of varying structure and conformation. The potent, ubiquitous environmental mutagen/carcinogen benzo[a]pyrene (B[a]P) induces all classes of mutations with G-->T base substitutions predominating. We recently showed that a G-->T mutagenesis pathway for the major adduct of B[a]P ([+ta]-B[a]P-N2-dG) in Escherichia coli depends on Y-family member DNAP V. Since no X-ray crystal study for DNAP V has been reported, no structure is available to help in understanding the structural basis for dATP insertion associated with G-->T mutations from [+ta]-B[a]P-N2-dG. Herein, we do homology modeling to construct a model for UmuC, which is the polymerase subunit of DNAP V. The sequences of eight Y-family DNAPs were aligned based on the positioning of conserved amino acids and an analysis of conserved predicted secondary structure, as well as insights gained from published X-ray structures of five Y-family members. Starting coordinates for UmuC were generated from the backbone coordinates for the Y-family polymerase Dpo4 for reasons discussed, and were refined using molecular dynamics with CHARMM 27. A survey of the literature revealed that E. coli DNAP V and human DNAP eta show a similar pattern of dNTP insertion opposite a variety of DNA lesions. Furthermore, E. coli DNAP IV and human DNAP kappa show a similar dNTP insertional pattern with these same DNA lesions, although the insertional pattern for DNAP IV/kappa differs from the pattern for DNAPs V/eta. These comparisons prompted us to construct and refine models for E. coli DNAP IV and human DNAPs eta and kappa as well. The dNTP/template binding pocket of all four DNAPs was inspected, focusing on the array of seven amino acids that contact the base of the incoming dNTP, as well as the template base. DNAPs V and eta show similarities in this array, and DNAPs IV and kappa also show similarities, although the arrays are different for the two pairs of DNAPs. Thus, there is a correlation between structural similarities and insertional similarities for the pairs DNAPs V/eta and DNAPs IV/kappa. Although the significance of this correlation remains to be elucidated, these observations point the way for future experimental studies.

Mutagenesis studies with four stereoisomeric N2-dG benzo[a]pyrene adducts in the identical 5′-CGC sequence used in NMR studies: G→T mutations dominate in each case

Benzo[a]pyrene (B[a]P) is a polycyclic aromatic hydrocarbon (PAH) and a potent mutagen/carcinogen found ubiquitously in the environment. B[a]P is primarily metabolized to diol epoxides, which react principally at N2-dG in DNA. B[a]P-N2-dG adducts have been shown to induce a variety of mutations, notably G-->T, G-->A, G-->C and -1 frameshifts. Four stereoisomers of B[a]P-N2-dG (designated: [+ta]-;, [+ca]-, [-ta] and [-ca]) were studied by NMR in duplex 11mers in a 5'-CGC sequence context, and each adopted a different adduct conformation (Geacintov, et al. (1997) Chem. Res. Toxicol., 10, 111). Herein these four identical B[a]P-containing 11mers are built into duplex plasmid genomes and mutagenesis studied in Escherichia coli following SOS-induction. In nucleotide excision repair (NER) proficient E.coli, no adduct-derived mutants are detected. In NER deficient E.coli, G-->T mutations dominate for all four stereoisomers [+ta]-, [+ca]-, [-ta] and [-ca]-B[a]P-N(2)-dG, and mutation frequency is similar. Thus, the mutagenic pattern for these four B[a]P-N2-dG stereoisomers is the same, in spite of the fact that they adopt dramatically different conformations in ds-oligonucleotides as determined by NMR. These findings suggest that adduct conformation must be fluid enough in the 5'-CGC sequence that the duplex DNA conformation can interconvert to mutagenic and non-mutagenic conformations during lesion-bypass. A comparison of all published studies with these four B[a]P-N2-dG stereoisomers in E.coli reveals that B[a]P-N2-dG adduct stereochemistry tends to have a lesser impact on mutagenic pattern (e.g. G-->T versus G-->A mutations) than does DNA sequence context, which is discussed.

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.

A role for DNA polymerase V in G --> T mutations from the major benzo[a]pyrene N2-dG adduct when studied in a 5'-TGT sequence in E. coli

Benzo[a]pyrene (B[a]P), a potent mutagen/carcinogen, is metabolically activated to (+)-anti-B[a]PDE, which induces a full spectrum of mutations (e.g. GC --> TA, GC --> AT, etc.) principally via its major adduct [+ta]-B[a]P-N2-dG. Recent findings suggest that different lesion bypass DNA polymerases may be involved in different mutagenic pathways, which is the subject of this report. [+ta]-B[a]P-N2-dG built into a plasmid in a 5'-TGT sequence gives approximately equal numbers of G --> T and G --> A mutations when host E. coli are UV irradiated prior to transformation, so this sequence context was chosen to investigate what DNA polymerases are involved in G --> T versus G --> A mutations. G --> T mutations decline (>10-fold) if E. coli either are not UV-irradiated or are deficient in DNA polymerase V ((delta)umuD/C), demonstrating a role for damage-inducible DNA Pol V in a G --> T pathway. G --> T mutations are not affected by transformation into E. coli deficient in either DNA polymerases II or IV. While the work herein was in progress, Lenne-Samuel et al. [Mol. Microbiol. 38 (2000) 299] built the same adduct into a plasmid in a 5'-GGA sequence, and showed that the frequency of G --> T mutations was similar in UV-irradiated and unirradiated host E. coli cells, suggesting no involvement by damage-inducible, lesion bypass DNA polymerases (i.e., not II, IV or V); furthermore, a role for DNA Pol V was explicitly ruled out. The easiest way to reconcile the findings of Lenne-Samuel et al. with the findings herein is if two G --> T mutagenic pathways exist for [+ta]-B[a]P-N2-dG, where sequence context dictates which pathway is followed. In contrast to the G --> T mutations, herein G --> A mutations from [+ta]-B[a]P-N2-dG in the 5'-TGT sequence context are shown not to be affected by UV-irradiation of host E. coli, and are not dependent on DNA Pol V, or Pol II, Pol IV, or the damage-inducible, but SOS-independent UVM system. Published studies, however, have shown that G --> A mutations are usually enhanced by UV-irradiation of host E. coli prior to the introduction of plasmids either site-specifically modified with [+ta]-B[a]P-N2-dG or randomly adducted with (+)-anti-B[a]PDE; both findings imply the involvement of a lesion-bypass DNA polymerase. These disparate results suggest the existence of two G --> A mutagenic pathways for [+ta]-B[a]P-N2-dG as well, although confirmation of this awaits further study. In conclusion, a comparison between the evidence presented herein and published findings suggests the existence of two distinct mutagenic pathways for both G --> T and G --> A mutations from [+ta]-B[a]P-N2-dG, where in each case one pathway is not damage-inducible and not dependent on a lesion-bypass DNA polymerase, while the second pathway is damage-inducible and dependent on a lesion-bypass DNA polymerase. Furthermore, DNA sequence context appears to dictate which pathway (as defined by the involvement of different DNA polymerases) is followed in each case.

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

Molecular modeling of four stereoisomers of the major B[a]PDE adduct (at N(2)-dG) in five cases where the structure is known from NMR studies: molecular modeling is consistent with NMR results

The potent mutagen/carcinogen benzo[a]pyrene (B[a]P) is metabolically activated to (+)-anti-B[a]PDE, which is known to induce a variety 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-N(2)-dG when bypassed during DNA replication (perhaps by different DNA polymerases). Our previous molecular modeling studies have suggested that conformational complexity might be extensive in that B[a]P-N(2)-dG adducts appeared capable of adopting at least sixteen potential conformational classes in ds-DNA [e.g., Kozack and Loechler (1999) Carcinogenesis 21, 1953], although only eight seemed likely to be relevant to base substitution mutagenesis. Such molecular modeling studies are only likely to be valuable for the interpretation of mutagenesis results if global minimum energy conformations for adducts are found and if the differences in the energies of these different conformations can be computed reasonably accurately. One approach to assessing the reliability of our molecular modeling techniques is considered herein. Using a five-step molecular modeling protocol, which importantly included a molecular dynamics version of simulated annealing, eight conformations are studied in each of five cases. (The five cases are listed below, and were chosen because in each case the preferred solution conformation is known from a NMR study.) Of the eight conformations studied, the one computed to be lowest in energy is the same conformation as the one observed by NMR in four of the five cases: 5'-CGC sequence with [+ta]-, [-ta]-, and [+ca]-B[a]P-N(2)-dG, and 5'-TGC sequence with [+ta]-B[a]P-N(2)-dG. In the fifth case (5'-CGC sequence with [-ca]-B[a]P-N(2)-dG), the known NMR conformation is computed to be second lowest in energy, but it is within approximately 1.7 kcal of the computed lowest energy conformation. These results suggest that molecular modeling is surprisingly accurate in computing lowest energy conformations and that it should be useful in assessing the relative energies of different conformations. This is especially important given that currently molecular modeling is the only means available to study the energetics of minor conformations of DNA adducts.

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.

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.

DNA polymerase II (polB) is involved in a new DNA repair pathway for DNA interstrand cross-links in Escherichia coli

DNA-DNA interstrand cross-links are the cytotoxic lesions for many chemotherapeutic agents. A plasmid with a single nitrogen mustard (HN2) interstrand cross-link (inter-HN2-pTZSV28) was constructed and transformed into Escherichia coli, and its replication efficiency (RE = [number of transformants from inter-HN2-pTZSV28]/[number of transformants from control]) was determined to be approximately 0.6. Previous work showed that RE was high because the cross-link was repaired by a pathway involving nucleotide excision repair (NER) but not recombination. (In fact, recombination was precluded because the cells do not receive lesion-free homologous DNA.) Herein, DNA polymerase II is shown to be in this new pathway, since the replication efficiency (RE) is higher in a polB+ ( approximately 0. 6) than in a DeltapolB (approximately 0.1) strain. Complementation with a polB+-containing plasmid restores RE to wild-type levels, which corroborates this conclusion. In separate experiments, E. coli was treated with HN2, and the relative sensitivity to killing was found to be as follows: wild type < polB < recA < polB recA approximately uvrA. Because cells deficient in either recombination (recA) or DNA polymerase II (polB) are hypersensitive to nitrogen mustard killing, E. coli appears to have two pathways for cross-link repair: an NER/recombination pathway (which is possible when the cross-links are formed in cells where recombination can occur because there are multiple copies of the genome) and an NER/DNA polymerase II pathway. Furthermore, these results show that some cross-links are uniquely repaired by each pathway. This represents one of the first clearly defined pathway in which DNA polymerase II plays a role in E. coli. It remains to be determined why this new pathway prefers DNA polymerase II and why there are two pathways to repair cross-links.

DNA polymerase II (polB) is involved in a new DNA repair pathway for DNA interstrand cross-links in Escherichia coli

DNA-DNA interstrand cross-links are the cytotoxic lesions for many chemotherapeutic agents. A plasmid with a single nitrogen mustard (HN2) interstrand cross-link (inter-HN2-pTZSV28) was constructed and transformed into Escherichia coli, and its replication efficiency (RE = [number of transformants from inter-HN2-pTZSV28]/[number of transformants from control]) was determined to be approximately 0.6. Previous work showed that RE was high because the cross-link was repaired by a pathway involving nucleotide excision repair (NER) but not recombination. (In fact, recombination was precluded because the cells do not receive lesion-free homologous DNA.) Herein, DNA polymerase II is shown to be in this new pathway, since the replication efficiency (RE) is higher in a polB+ ( approximately 0. 6) than in a DeltapolB (approximately 0.1) strain. Complementation with a polB+-containing plasmid restores RE to wild-type levels, which corroborates this conclusion. In separate experiments, E. coli was treated with HN2, and the relative sensitivity to killing was found to be as follows: wild type < polB < recA < polB recA approximately uvrA. Because cells deficient in either recombination (recA) or DNA polymerase II (polB) are hypersensitive to nitrogen mustard killing, E. coli appears to have two pathways for cross-link repair: an NER/recombination pathway (which is possible when the cross-links are formed in cells where recombination can occur because there are multiple copies of the genome) and an NER/DNA polymerase II pathway. Furthermore, these results show that some cross-links are uniquely repaired by each pathway. This represents one of the first clearly defined pathway in which DNA polymerase II plays a role in E. coli. It remains to be determined why this new pathway prefers DNA polymerase II and why there are two pathways to repair cross-links.

The major, N2-dG adduct of (+)-anti-B[a]PDE induces G-->A mutations in a 5'-AGA-3' sequence context

Previously, in a random mutagenesis study, the (+)-anti diol epoxide of benzo[a]pyrene [(+)-anti-B[a]PDE] was shown to induce a complex mutational spectrum in the supF gene of an Escherichia coli plasmid, which included insertions, deletions and base substitution mutations, notably a significant fraction of GC-->TA, GC-->AT and GC-->CG mutations. At some sites, a single type of mutation dominated and to understand individual mutagenic pathways these sites were chosen for study by site-specific means to determine whether the major adduct, [+ta]-B[a]P-N2-dG, was responsible. [+ta]-B[a]P-N2-dG was shown to induce approximately 95% G-->T mutations in a 5'-TGC-3' sequence context and approximately 80% G-->A mutations in a 5'-CGT-3' sequence context. (+)-anti-B[a]PDE induced principally GC-->CG mutations in the G133 sequence context (5'-AGA-3') in studies using both SOS-uninduced or SOS-induced E. coli. Herein, [+ta]-B[a]P-N2-dG is shown to induce principally G-->A mutations (>90%) either without or with SOS induction in a closely related 5'-AGA-3' sequence context (identical over 7 bp). This is the first time that there has been a discrepancy between the mutagenic specificity of (+)-anti-B[a]PDE versus [+ta]-B[a]P-N2-dG. Eight explanations for this discordance are considered. Four are ruled out; e.g. the second most prevalent adduct [+ca]-B[a]P-N2-dG also induces a preponderance of G-->A mutations (>90%), so it also is not responsible for (+)-anti-B[a]PDE-induced G133-->C mutations. The four explanations not ruled out are discussed and include that another minor adduct might be responsible and that the 5'-AGA-3' sequence context differed slightly in the studies with [+ta]-B[a]P-N2-dG versus (+)-anti-B[a]PDE. In spite of the discordance, [+ta]-B[a]P-N2-dG induces G-->A mutations in the context studied herein and this result has proven useful in generating a hypothesis for what conformations of [+ta]-B[a]P-N2-dG are responsible for G-->T versus G-->A mutations.

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

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

Molecular modeling (simulated annealing) was used to study the conformations in dsDNA of [+ta]-B[a]P-N2-dG (R.E. Kozack and E.L.Loechler, accompanying paper), which is the major benzo[a]pyrene (B[a]P) adduct. Sixteen classes of conformations were identified, and are analyzed herein vis-a-vis the two most prominent B[a]P mutations, G-->T and G-->A base substitutions. Eight conformations seem more relevant to frameshift mutagenesis, so they are excluded, leaving eight conformations as follows. Two conformations (BPmi5 and BPmi3) retain Watson-Crick G:C base pairing having the B[a]P moiety of the adduct in the minor groove. Two conformations (BPma5 and BPma3) have the Hoogsteen orientation with B[a]P in the major groove. Four conformations are base displaced and have B[a]P stacked in the helix with the dG moiety of the adduct displaced into either the major groove (Gma5 and Gma3) or the minor groove (Gmi5 and Gmi3). Three of these eight conformations (BPma5, BPma3 and Gma3) are universally high in energy. The two conformations that retain G:C base pairing potential (BPmi5 and BPmi3) are likely to be non-mutagenic. Of the three remaining conformations, Gmi5 can be relatively low in energy, but is distorted. A correlation exists between the calculated energies for the remaining two base displaced conformations and mutagenesis for [+ta]-B[a]P-N2-dG, leading to the hypothesis that Gma5 is responsible for G-->T mutations and Gmi3 is responsible for G-->A mutations. Gma5 and Gmi3 resemble each other, except that dG is in the major and minor grooves, respectively. An incipient rationale for this hypothesis is discussed: DNA polymerase might be triggered to follow a different mutagenic pathway depending upon whether a non-informational lesion has bulk protruding into the major or minor groove. A pathway for interconversion between these eight conformations is also proposed and its implications are discussed; e.g. four steps are required to interconvert between Gma5 and Gmi3.

How stereochemistry affects mutagenesis by N2-deoxyguanosine adducts of 7,8-dihydroxy-9,10-epoxy-7,8,9,10-tetrahydrobenzo[a]pyrene: configuration of the adduct bond is more important than those of the hydroxyl groups

Previous work has shown that the major adduct from the (+)-anti diol epoxide of benzo[a]pyrene (B[a]P), which forms at N2-deoxyguanosine [(+)-trans-anti-B[a]P-N2-dG], is capable of inducing either predominantely G --> T mutations ( approximately 95%) in a 5'-TGC-3 sequence context or predominantly G --> A mutations ( approximately 80%) in a 5'-CGT-3' sequence context. This is likely to be attributable to the major adduct being in a different mutagenic conformation in each case. In the next phase of this work, the questions to be addressed are what conformation is associated with what mutation and why? To help define what aspect of adduct structure is important to mutagenesis, the work herein reports on the mutations induced in a single sequence context by four stereoisomers of B[a]P-N2-dG: (+)-trans-, (+)-cis-, (-)-trans-, and (-)-cis-. The (+)-trans- and (-)-cis-adducts show a remarkably similar mutational pattern with G --> A mutations predominating ( approximately 80%). The (-)-trans- and (+)-cis-adducts also show a similar mutational pattern with a more even mixture of G --> T, G --> A, and G --> C mutations. Each of these adducts has an adduct bond and three hydroxyl groups at four consecutive saturated carbons in the B[a]P moiety of the adduct; the stereochemistry at these four positions differs in each of the adducts. The (+)-trans- and (-)-cis-adducts are a pair sharing the S configuration for the adduct bond, although they are a mirror image vis-a-vis the hydroxyl groups. The (-)-trans- and (+)-cis-adducts share the opposite adduct bond stereochemistry (R) but differ in the stereochemistry of their hydroxyl groups. Thus, there is a correlation suggesting that anti-B[a]P-N2-dG adduct mutagenesis is more dependent on the stereochemistry of the adduct bond than on the stereochemistry of the hydroxyl groups.

The major, N2-dG adduct of (+)-anti-B[a]PDE shows a dramatically different mutagenic specificity (predominantly, G --> A) in a 5'-CGT-3' sequence context

Mutations induced by the (+)-anti diol epoxide of benzo[a]pyrene [(+)-anti-B[a]PDE] were described previously in the supF gene of the Escherichia coli plasmid pUB3 [Rodriguez et al.(1993) Biochemistry, 32, 1759]. (+)-anti-B[a]PDE induced a complex pattern of mutations, including insertions, deletions, frameshifts, as well as base substitution mutations, which for G:C base pairs alone included a significant fraction of G:C --> T:A, A:T and C:G mutations. A variety of results suggest that most of these mutations arise from the major adduct ([+ta]-B[a]P-N2-dG), raising the question how can a single adduct induce different kinds of mutations? Our working hypothesis in this regard is that (1) an adduct can adopt multiple conformations; (2) different conformations cause different mutations; and (3) adduct conformation is controlled by various factors, such as DNA sequence context. To investigate what conformation is associated with what mutation, it is essential to find examples where [+ta]-B[a]P-N2-dG induces principally one kind of mutation as a prelude to the study in that same context of the conformation(s) potentially relevant to mutagenesis. Earlier work indicated that (+)-anti-B[a]PDE gave a preponderance of G --> A mutations in a 5'-CGT-3 sequence context, and herein it is shown that these mutations are likely to be attributable to the major adduct, since in this same sequence context [+ta]-B[a]P-N2-dG studied site specifically also induces principally G --> A mutations ( approximately 82%). Previously, [+ta]-B[a]P-N2-dG was shown to induce principally G --> T mutations (approximately 97%) in a 5'-TGC-3' sequence context. Thus, by simply altering its surrounding sequence context this adduct can give a preponderance of either G --> A or G --> T mutations. This is the most dramatic change in base substitution mutagenic specificity for an adduct described to date and illustrates that the qualitative pattern of mutagenesis by a bulky adduct can be remarkably diverse.

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.

Sequence specific mutagenesis of the major (+)-anti-benzo[a]pyrene diol epoxide-DNA adduct at a mutational hot spot in vitro and in Escherichia coli cells

In the supF gene, most (+)-anti-benzo[a]pyrene diol epoxide ((+)-anti-B[a]PDE) mutagenesis hot spots in Escherichia coli are in 5'-GG sequences [Rodriguez and Loechler (1993) Carcinogenesis 14, 373-383]. A major hot spot was detected at G1 in the sequence 5'-GCG1G2-CCAAAG, whereas G2 yielded very few mutants. In order to investigate the details of such sequence context effects of (+)-anti-B[a]PDE mutagenesis, we have constructed 25-mer oligonucleotides and single-stranded M13 genomes containing the above decamer sequence, in which the trans-N2-dG adduct induced by (+)-anti-B[a]PDE [(+)-trans-anti-B[a]P-N2-dG] at G1 or G2 was introduced. In vitro DNA synthesis on the adducted 25-mers was strongly blocked at each site, although the 3'-->5' exonuclease-deficient Klenow fragment could incorporate a nucleotide opposite the adduct in the presence of Mn2+. For both sites purine nucleotides were preferred. The ratio Vmax/K(m) indicated that the efficiency of incorporation of dGTP opposite these sites was very similar, but dATP incorporation opposite the adduct at G1 was five-fold more efficient than that at G2. For each site, further extension beyond the adducted nucleotide was investigated by annealing four different primers, in which only the nucleotide opposite the adducted deoxyguanosine was altered. Significant extension was only observed when deoxyadenosine was located opposite adducted G1. When the M13 genomes containing the (+)-trans-anti-B[a]P-N2-dG were replicated in E. coli, survival of each adducted genome was less than 1% as compared to the unadducted genome. Upon induction of SOS, viability increased 2-6-fold. DNA sequencing showed no base substitutions in the progeny from SOS-uninduced cells, although small deletions in a quasipalindromic sequence occurred with the adduct being located at either site. However, following SOS induction, up to 40% targeted base substitutions were detected when the adduct was located at G1, while approximately 12% of the progeny were mutants with the adduct at G2. Most base substitutions were targeted G-->T transversions. We conclude that (+)-trans-anti-B[a]P-N2-dG is a highly mutagenic and replication blocking lesion. In addition, the biological consequence of this adduct depends on whether it is located at G1 or G2, suggesting that sequence context plays a major role in the mutagenic processing of this adduct.

Evidence for a recombination-independent pathway for the repair of DNA interstrand cross-links based on a site-specific study with nitrogen mustard

DNA-DNA interstrand cross-links are thought to be important for the cytotoxicity of many chemotherapeutic agents. To study this more definitively, adduct site-specific methods are used to construct a plasmid with a single nitrogen mustard interstrand cross-link (inter-HN2-pTZSV28). Replication efficiency (RE = [colonies from (inter-HN2-pTZSV28)/(control with no cross-link)]) is approximately 0.3 following transformation into Escherichia coli, implying that the cross-link is repaired. The commonly accepted pathway for cross-link repair, which involves both nucleotide excision repair (NER) and recombination, is ruled out since RE is approximately 0.3 in a delta recA strain. Non-RecA-directed recombination such as copy-choice is also unlikely. However, NER is involved since RE was approximately 0.02 in strains deficient in NER. Base excision repair is not important since RE is approximately 0.3 in strains deficient in 3-methyladenine DNA glycosylases I and II, FAPY DNA glycosylase, both known apurinic/apyrimidinic endonucleases, or DNA deoxyribophosphodiesterase. Another hypothetical repair pathway hinging on a 5' --> 3' exonuclease activity is unlikely since RE is approximately 0.3 in cells deficient in either the 5' --> 3' exonuclease activities of DNA polymerase I, exonuclease VII, or RecJ. Thus, aside from NER, it is unclear what else participates in this recombination-independent repair pathway, although a pathway involing NER followed by replicative bypass of the lesion is the current working hypothesis. Psoralen interstrand cross-links appear not to be repairable by this second pathway, which may have implications for the relative cytotoxicity of interstrand cross-links from different agents.

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.

Targeted A --> T and G --> T mutations induced by site-specific deoxyadenosine and deoxyguanosine adducts, respectively, from the (+)-anti-diol epoxide of dibenz[a,j]anthracene in M13mp7L2

The studies described in this report directly examined the mutagenicity in Escherichia coli of both a deoxyadenosine (dAdo) and a deoxyguanosine (dGuo) adduct derived from (+)-anti-dibenz[a,j]-anthracene-3,4-diol 1,2-epoxide [(+)anti-DB[a,j]A-DE] that were site-specifically placed in a single-stranded M13mp7L2 replication vector. An 11-base oligonucleotide (5'-CTC ACG CTT CT-3') containing either a single (+)anti-DB[a,j]A-DE--trans-N2-dGuo or (+)anti-DB[a,j]A-DE--trans-N6dAdo adduct was successfully incorporated into single-stranded M13mp7L2 plasmid via ligation. In vitro studies using E. coli DNA polymerase I (Klenow fragment)indicated that both adducts were effective blocks for polymerase action. E. coli strains JM103 and JM103 uvrA6 were subsequently transformed with control (unadducted) and adduct-containing M13mp7L2 constructs followed by analysis of progeny DNA. In both JM103 and JM103 uvrA6 cells, plaque yields were markedly reduced with adduct containing vectors compared to control vectors. Activation of the inducible bacterial DNA repair system (SOS) by UV light only slightly increased the number of plaques recovered from either bacterial strain transformed with adduct-containing vectors. Targeted mutations were obtained with both adduct-containing vectors in both bacterial strains, whereas no mutations were detected in plaques recovered from control M13mp7L2 vectors. In JM103 cells, (+)anti-DB[a,j]A-DE--N6-dAdo induced exclusively A --> t transversions and (+)anti-DB[a,j]A-DE--N2-dGuo induced exclusively G --> T transversions. In JM103 uvrA6 cells, similar targeted transversion mutations were also obtained except that a few C deletions (i.e., aprroximately 10% of the mutations) were detected immediately 3' to the dAdo adduct. While mutagenesis was SOS dependent in JM103 cells [<0.15% (-SOS) vs approximately 1.3% (+SOS)], it appeared to be SOS independent in JM103 uvrA6 cells (approximately 1-2% in the presence or absence of SOS induction). It is argued that adduct-induced G --> T mutations can be rationalized by either misinformational or noninformational mechanisms. In contrast, A --> T mutations are unlikely to arise via a misinformational pathway, which provides the strongest support to date that bulky DNA adducts can induce mutations via a noninformational pathway.

The major, N2-Gua adduct of the (+)-anti-benzo[a]pyrene diol epoxide is capable of inducing G-->A and G-->C, in addition to G-->T, mutations

Mutations induced by the (+)-anti-diol epoxide of benzo[a]pyrene [(+)-anti-B[a]PDE] were collected in the supF gene of the Escherichia coli plasmid pUB3. pUB3 was reacted with (+)-anti-B[a]-PDE and then either (1) transformed immediately into E. coli or (2) heated at 80 degrees C for 10 min and then cooled prior to transformation--the latter to probe mechanism [Rodriguez & Loechler (1993) Biochemistry 32, 1759]. Qualitatively, heating did not affect the mutagenic pattern, except at the major base substitution hotspot in supF, G115, where principally G-->T mutations were obtained prior to heating, while after heating, G-->A and G-->C mutations became statistically significantly more prevalent. Several studies have suggested that a heat-induced chemical transformation of a (+)-anti-B[a]PDE adduct at G115 (e.g., into an apurinic site) is not likely to explain the change in mutational pattern. The most likely model is that (+)-anti-B[a]P-N2-Gua is initially trapped in a metastable conformation giving principally G-->T mutations, while heating induces a change to a stable conformation(s) resulting in G-->T, A, and C mutations. This suggests that adduct conformational complexity is at the root of adduct mutational complexity. To investigate this model, a plasmid (B[a]P-G115-pRE1) with (+)-anti-B[a]P-N2-Gua in the G115 sequence context is constructed using adduct site-specific techniques.(ABSTRACT TRUNCATED AT 250 WORDS)

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.

The major, N2-Gua adduct of the (+)-anti-benzo[a]pyrene diol epoxide can be unstable in double-stranded DNA

The mechanisms of mutagenesis by the (+)-anti diol epoxide of benzo[a]pyrene [(+)-anti-B[a]PDE] was investigated in supF of the Escherichia coli plasmid pUB3 [Rodriguez & Loechler (1993) Biochemistry 32, 1759]. pUB3 was reacted with (+)-anti-B[a]PDE, then either (1) transformed immediately into E. coli or (2) heated at 80 degrees C for 10 min prior to transformation--the latter to probe mechanism. Qualitatively, heating did not have a statistically significant effect on the mutagenic pattern, except at the major base substitution hot spot, G115, in supF; principally, G115-->T mutations were obtained prior to heating, while after heating, G115-->A and G115-->C mutations became more prevalent. Quantitatively, heating caused an approximately 2-fold decrease in mutation frequency. Heating released a small fraction of adducts (approximately 5%), and the chemistry and implications of this reaction are investigated herein. Although the major adduct of (+)-anti-B[a]PDE (formed at N2-Gua) is generally regarded to be heat stable, it can be quite unstable in double-stranded (but not single-stranded) DNA at low [Mg2+]. Heating releases the corresponding tetraols from (+)-anti-B[a]P-N2-Gua in approximately the same ratio as for simple hydrolysis of (+)-anti-B[a]PDE itself. This and other results suggest that guanine remains in DNA when (+)-anti-B[a]P-N2-Gua adducts are hydrolyzed. [No evidence for any other chemical change in (+)-anti-B[a]PDE adducts was found.] While no general acid/base or nucleophilic catalysis was observed, adduct hydrolysis was specific acid catalyzed down to pH approximately 5.6, where the pH-rate profile showed a break to a slope of approximately 0.0. This break probably indicates the pKa of (+)-anti-B[a]P-N2-Gua protonated at the N2-position, which is higher than expected.(ABSTRACT TRUNCATED AT 250 WORDS)

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.

Construction of Escherichia coli vectors containing deoxyadenosine and deoxyguanosine adducts from (+)-anti-dibenz[a,j]anthracene diol epoxide at a defined site

Dibenz[a,j]anthracene (DB[a,j]A) is a carcinogenic polycylic aromatic hydrocarbon, which is metabolically activated through the formation of bay region diol epoxides. Site-specifically modified M13mp19-based vectors containing a single (+)-anti-dibenz[a,j] anthracene diol epoxide [(+)-anti-DB[a,j[A-DE]-deoxyguanosine (dGuo) or -deoxyadenosine (dAdo) adduct were constructed. Four-base oligonucleotides, 5'-HOTGCA-3' and 5'-HOCATG-3', corresponding to the central four base pairs in the PstI and SphI restriction endonuclease sites, respectively, in the multiple cloning region of M13mp19, were reacted in solution with (+/-)-anti-DB[a,j]A-DE. The resulting adducted oligonucleotides were separated and purified using reverse-phase HPLC. Several different singly adducted oligonucleotides were isolated, consisting of the various cis and trans addition products of the (+) and (-) enantiomers of the diol epoxide bound to dGuo or dAdo in the oligonucleotides. 5'-HOTGCA-3' containing the (+)-anti-DB[a,j]A-trans-N2-dGuo adduct [T(DB[a,j]A-N2)GCA] and 5'-HOCATG-3' containing the (+)-anti-DB[a,j]A-trans-N6-dAdo adduct [C(DB[a,j]A-N6)ATG) were selected for subsequent ligation into M13mp19 vectors that had been constructed with a corresponding four base gap in the minus strand. Both unmodified and adducted oligonucleotides were successfully ligated into the M13mp19 vectors, [yields: unmodified -TGCA-M13mp19 (approximately 32%) and -CATG- M13mp19 (approximately 42%); adducted T(DB[a,j]A-N2)GCA-M13mp19 (approximately 13%) and C(DB[a,j]A-N6)ATG-M13mp19 (approximately 12%)]. The dAdo adduct-containing vector was characterized. The presence of a dAdo-DNA adduct at the recognition site of SphI inhibited restriction by SphI.(ABSTRACT TRUNCATED AT 250 WORDS)

AP sites are not significantly involved in mutagenesis by the (+)-anti diol epoxide of benzo[a]pyrene: the complexity of its mutagenic specificity is likely to arise from adduct conformational polymorphism

In previous work, mutations induced by the (+)-anti diol epoxide of benzo[a]pyrene [(+)-anti-B[a]PDE] were scored in the supF gene of the Escherichia coli plasmid pUB3 [Rodriguez & Loechler (1993) Biochemistry 32, 1759]. pUB3 was reacted with (+)-anti-B[a]PDE and then either (1) transformed immediately into E. coli or (2) heated at 80 degrees C for 10 min prior to transformation. Heating only released a small fraction of adducts (approximately 5%) and did not significantly affect the mutagenic pattern at most sites in supF. However, at the major base substitution hotspot, G115, principally G-->T mutations (87%) were obtained prior to heating, while after heating, G-->T mutations decreased (45%) and G-->A (21%) and G-->C (33%) mutations became more prevalent. One model for this result is that prior to heating a heat-labile adduct at G115 causes one pattern of mutagenesis, but after heating the labile adduct is hydrolyzed to an apurinic site (AP site), which causes a second mutational pattern. To test this, a role for AP sites generated from labile adducts by heating at 80 degrees C for 10 min is investigated. It is shown that when plasmid pUB3 contains 22 (+)-anti-B[a]PDE adducts, 0.6% (or fewer) are converted to AP sites as determined in an assay based upon the action of an AP-endonuclease. In a separate line of investigation not involving (+)-anti-B[a]PDE adducts, mutation frequency (MF) per AP site is estimated. (In these experiments, AP sites were introduced into pUB3 by the classic procedure of heating at 70 degrees C/pH 5.0 to hydrolyze purines.(ABSTRACT TRUNCATED AT 250 WORDS)

Evidence from in vitro replication that O6-methylguanine can adopt multiple conformations

The effect of O6-methylguanine (m6G) on replication, in a partially double-stranded defined 25-base oligonucleotide, has been studied under nonlimiting conditions of unmodified dNTPs and over an extended time period, using the Klenow fragment of Escherichia coli DNA polymerase I. The sequence surrounding m6G has flanking cytosines (C-m6G-C), and the initial steady-state kinetics have been reported. When the primer was annealed so that the first base to be replicated was m6G, replication was virtually complete in approximately 5 min, although the reaction appears biphasic. When annealed with a primer where thymine or cytosine is paired opposite template m6G, about half the molecules were replicated in the first 15 sec, and no significant further replication was seen over a 1-hr period. When m6G was dealkylated by DNA-O6-methylguanine-methyltransferase, replication was rapid with no blockage. These data suggest that there can be two (or more) conformations of m6G. In these studies the term syn refers to conformers interfering with base-pairing, whereas anti refers to those allowing such base-pairing. Previous physical studies by others indicate that syn- and anti-conformers of the methyl group relative to the N1 of guanine are possible. Here molecular modeling/computational studies are described, suggesting that syn- and anti-m6G can be of similar energy in DNA, and, therefore, these two conformers may explain the two types of species observed during in vitro replication. An alternative explanation could be the possibility that the different species may manifest differential interactions of m6G with Klenow fragment. These results may provide a rationale for why m6G lesions in vivo have been reported to be lethal as well as mutagenic.

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.

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)

Mutagenic specificity of the (+)anti-diol epoxide of dibenz[a,j]anthracene in the supF gene of an Escherichia coli plasmid

This study was designed to examine the mutagenic specificity of (+)anti-dibenz[a,j]anthracene 3,4-diol-1,2-epoxide ((+)anti-DB[a,j]A-DE) in SOS-induced repair-proficient Escherichia coli ES87 (delta pro-lac, strA)/F' (pro+, lac1Q, lac1am26, lacZ delta M15). The plasmid pUB3, which contains the mutation target gene, supF, was modified with (+)anti-DB[a,j]A-DE in vitro (two to five adducts/plasmid) and then transformed into bacteria by electroporation. The spontaneous mutation frequency for unmodified pUB3 in uninduced cells was about 2 x 10(-6) and for SOS-induced cells, about 8 x 10(-6). The spontaneous supF- mutations were primarily insertions, deletions, and frameshifts. The mutation frequency for (+)anti-DB[a,j]A-DE-modified pUB3 was about 8 x 10(-6) and about 32 x 10(-6) for uninduced cells and SOS-induced cells, respectively. (+)anti-DB[a,j]A-DE induced primarily point mutations in supF in SOS-induced cells. GC-->AT transitions were the major mutations observed in SOS-induced cells (37%). GC-->TA (21%) and GC-->CG (8.6%) transversion mutations were also observed, whereas mutations at AT base pairs were rare (1.9%). Furthermore, a large number of tandem GC/GC-->AT/AT transition mutations were also observed (about 15% of all mutations in SOS-induced cells). Taken together, single and tandem GC-->AT mutations accounted for slightly over half (about 51%) of the mutations observed in SOS-induced cells. These results demonstrated that (+)anti-DB[a,j]A-DE was mutagenic in repair-proficient E. coli; however, unlike other polycyclic aromatic hydrocarbons that induce primarily transversion mutations, (+)anti-DB[a,j]A-DE caused mostly GC-->AT transitions.

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.

Mutagenesis by (+)-anti-B[a]P-N2-Gua, the major adduct of activated benzo[a]pyrene, when studied in an Escherichia coli plasmid using site-directed methods

The suspected major mutagenic adduct of benzo[a]pyrene, (+)-anti-B[a]P-N2-Gua, is built into the unique PstI recognition site of the Escherichia coli plasmid, pUC19, in order to study its mutagenic potential. The adduct can either be at G437, which is replicated during leading strand DNA synthesis, or at G438, which is replicated during lagging strand DNA synthesis. The DNA strand complementary to the strand containing the (+)-anti-B[a]P-N2-Gua adduct is saturated with UV lesions to minimize its potential to generate progeny. Although all in-frame mutations could have been detected, a G437----T transversion mutation is virtually exclusively obtained at a frequency of approximately 0.04% per adduct following transformation into Uvr+ E. coli when SOS is not induced, and approximately 0.18% when SOS is induced. The mutation frequency of the adduct in a Uvr- background is estimated to be approximately 0.2% when SOS is not induced, and approximately 0.9% when SOS is induced. The absence of G438----T mutations is rationalized. G----T mutations from (+)-anti-B[a]P-N2-Gua are compared to the mutational specificity of the ultimate mutagenic form of activated benzo[a]pyrene.