Molecular modeling studies of O2-alkylthymines and O4-alkylthymines in DNA: structures that may be pertinent to the incorporation of the corresponding dAlkTTP into DNA by DNA polymerases in vitro

Molecular criteria for mutagenesis by DNA methylation: Some computational elucidations

Alkylating agents and N-nitroso compounds are well-known mutagens and carcinogens which act by alkylating DNA at the nucleobase moieties. Criteria for mutagenicity through DNA alkylation include (a) absence of the Watson-Crick (N1-guanine and N3-thymine) protons, (b) rotation of the alkyl group away from the H-bonding zone, (c) configuration of the alkylated base pair close to the Watson-Crick type. This computational study brings together these three molecular criteria for the first time. Three methylated DNA bases-N7-methylguanine, O6-methylguanine and O4-methylthymine-are studied using computational chemical methods. Watson-Crick proton loss is predicted more feasible for the mutagenic O6-methylguanine and O4-methylthymine than for the non-mutagenic N7-methylguanine in agreement with the observed trend for pKa values. Attainment of a conformer conducive to mutagenesis is more feasible for O6-methylguanine than for O4-methylthymine, though the latter is more mutagenic. These methylated bases yield 9 H-bonded pairs with normal DNA bases. At biological pH, O6-methylguanine and O4-methylthymine would yield stable mutagenic pairs having Watson-Crick type configuration by H-bonded pairing with thymine and guanine respectively, while N7-methylguanine would yield a non-mutagenic pair with cytosine. The three criteria thus well differentiate the non-mutagenic N7-methylguanine from the mutagenic O6-methylguanine and O4-methylthymine in good accord with experimental observations.

DNA Polymerases η and ζ Combine to Bypass O(2)-[4-(3-Pyridyl)-4-oxobutyl]thymine, a DNA Adduct Formed from Tobacco Carcinogens

4-(Methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) and N'-nitrosonornicotine (NNN) are important human carcinogens in tobacco products. They are metabolized to produce a variety 4-(3-pyridyl)-4-oxobutyl (POB) DNA adducts including O(2)-[4-(3-pyridyl)-4-oxobut-1-yl]thymidine (O(2)-POB-dT), the most abundant POB adduct in NNK- and NNN-treated rodents. To evaluate the mutagenic properties of O(2)-POB-dT, we measured the rate of insertion of dNTPs opposite and extension past O(2)-POB-dT and O(2)-Me-dT by purified human DNA polymerases η, κ, ι, and yeast polymerase ζ in vitro. Under conditions of polymerase in excess, polymerase η was most effective at the insertion of dNTPs opposite O(2)-alkyl-dTs. The time courses were biphasic suggesting the formation of inactive DNA-polymerase complexes. The kpol parameter was reduced approximately 100-fold in the presence of the adduct for pol η, κ, and ι. Pol η was the most reactive polymerase for the adducts due to a higher burst amplitude. For all three polymerases, the nucleotide preference was dATP > dTTP ≫ dGTP and dCTP. Yeast pol ζ was most effective in bypassing the adducts; the kcat/Km values were reduced only 3-fold in the presence of the adducts. The identity of the nucleotide opposite the O(2)-alkyl-dT did not significantly affect the ability of pol ζ to bypass the adducts. The data support a model in which pol η inserts ATP or dTTP opposite O(2)-POB-dT, and then, pol ζ extends past the adduct.

Roles of translesion synthesis DNA polymerases in the potent mutagenicity of tobacco-specific nitrosamine-derived O2-alkylthymidines in human cells

The tobacco-specific nitrosamine 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) is a potent human carcinogen. Metabolic activation of NNK generates a number of DNA adducts including O(2)-methylthymidine (O(2)-Me-dT) and O(2)-[4-(3-pyridyl)-4-oxobut-1-yl]thymidine (O(2)-POB-dT). To investigate the biological effects of these O(2)-alkylthymidines in humans, we have replicated plasmids containing a site-specifically incorporated O(2)-Me-dT or O(2)-POB-dT in human embryonic kidney 293T (HEK293T) cells. The bulkier O(2)-POB-dT exhibited high genotoxicity and only 26% translesion synthesis (TLS) occurred, while O(2)-Me-dT was less genotoxic and allowed 55% TLS. However, O(2)-Me-dT was 20% more mutagenic (mutation frequency (MF) 64%) compared to O(2)-POB-dT (MF 53%) in HEK293T cells. The major type of mutations in each case was targeted T → A transversions (56% and 47%, respectively, for O(2)-Me-dT and O(2)-POB-dT). Both lesions induced a much lower frequency of T → G, the dominant mutation in bacteria. siRNA knockdown of the TLS polymerases (pols) indicated that pol η, pol ζ, and Rev1 are involved in the lesion bypass of O(2)-Me-dT and O(2)-POB-dT as the TLS efficiency decreased with knockdown of each pol. In contrast, MF of O(2)-Me-dT was decreased in pol ζ and Rev1 knockdown cells by 24% and 25%, respectively, while for O(2)-POB-dT, it was decreased by 44% in pol ζ knockdown cells, indicating that these TLS pols are critical for mutagenesis. Additional decrease in both TLS efficiency and MF was observed in cells deficient in pol ζ plus other Y-family pols. This study provided important mechanistic details on how these lesions are bypassed in human cells in both error-free and error-prone manner.

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

Base alterations in yeast induced by alkylating agents with differing Swain-Scott substrate constants

The base alterations induced by four alkylating agents, methyl methanesulfonate (MMS), ethyl methanesulfonate (EMS), N-nitroso-N-methylurea (MNU), and N-nitroso-N-ethylurea (ENU), have been determined at the URA3 locus in the yeast Saccharomyces cerevisiae. The mutagen treatment was carried out on yeast cells in the logarithmic phase of growth. The mutants were selected by their resistance to 7.3 mM-5-fluoroorotic acid at pH 3.8. DNA sequence analysis was carried out by the dideoxy chain termination method. The alkylating agents were selected for their widely differing Swain-Scott substrate constants (s values), which are as follows: MMS, s = 0.83; EMS, s = 0.67; MNU, s = 0.42; ENU, s = 0.26. A higher s value is correlated with a higher ratio of 7-alkylguanine to O6-alkylguanine in native DNA in vitro. 125 forward mutations from URA3----ura3 were sequenced with marked differences in the mutational spectra being observed as the s value changed. Five hotspots were recorded for the four alkylating agents. They were all G.C----A.T transition mutations. There was one common hotspot for all of them; there were two additional ones for the two ethylating agents (ENU and EMS) and two different ones for MNU. Four of the five hotspots have the 5'-GG-3' sequence with the 3'-guanine mutated. It was seen that MMS, which has the highest Swain-Scott substrate constant, yielded the widest array of mutational types. As the substrate constants decreased, the types of mutations became more and more restricted to the G.C----A.T transitions and the A.T----T.A transversions. The transitions are consistent with the concept that mutations arise from O6-alkylation of guanine and alkylation of thymine. The transversions are consistent with the notion of N1-alkylation of adenosine or adenylic acid.