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

Hepatoprotective Effects of Flavonoids against Benzo[a]Pyrene-Induced Oxidative Liver Damage along Its Metabolic Pathways

Benzo[a]pyrene (B[a]P), a highly carcinogenic polycyclic aromatic hydrocarbon primarily formed during incomplete organic matter combustion, undergoes a series of hepatic metabolic reactions once absorbed into the body. B[a]P contributes to liver damage, ranging from molecular DNA damage to the onset and progression of various diseases, including cancer. Specifically, B[a]P induces oxidative stress via reactive oxygen species generation within cells. Consequently, more research has focused on exploring the underlying mechanisms of B[a]P-induced oxidative stress and potential strategies to counter its hepatic toxicity. Flavonoids, natural compounds abundant in plants and renowned for their antioxidant properties, possess the ability to neutralize the adverse effects of free radicals effectively. Although extensive research has investigated the antioxidant effects of flavonoids, limited research has delved into their potential in regulating B[a]P metabolism to alleviate oxidative stress. This review aims to consolidate current knowledge on B[a]P-induced liver oxidative stress and examines the role of flavonoids in mitigating its toxicity.

Carcinogen-induced DNA structural distortion differences in the RAS gene isoforms; the importance of local sequence

BACKGROUND

Local sequence context is known to have an impact on the mutational pattern seen in cancer. The RAS genes and a smoking carcinogen, Benzo[a]pyrene diol epoxide (BPDE), have been utilised to explore these context effects. BPDE is known to form an adduct at the guanines in a number of RAS gene sites, KRAS codons 12, 13 and 14, NRAS codon 12, and HRAS codons 12 and 14.

RESULTS

Molecular modelling techniques, along with multivariate analysis, have been utilised to determine the sequence influenced differences between BPDE-adducted RAS gene sequences as well as the local distortion caused by the adducts.

CONCLUSIONS

We conclude that G:C > T:A mutations at KRAS codon 12 in the tumours of lung cancer patients (who smoke), proposed to be predominantly caused by BPDE, are due to the effect of the interaction methyl group at the C5 position of the thymine base in the KRAS sequence with the BPDE carcinogen investigated causing increased distortion. We further suggest methylated cytosine would have a similar effect, showing the importance of methylation in cancer development.

Base damage, local sequence context and TP53 mutation hotspots: a molecular dynamics study of benzo[a]pyrene induced DNA distortion and mutability

The mutational pattern for the TP53 tumour suppressor gene in lung tumours differs to other cancer types by having a higher frequency of G:C>T:A transversions. The aetiology of this differing mutation pattern is still unknown. Benzo[a]pyrene,diol epoxide (BPDE) is a potent cigarette smoke carcinogen that forms guanine adducts at TP53 CpG mutation hotspot sites including codons 157, 158, 245, 248 and 273. We performed molecular modelling of BPDE-adducted TP53 duplex sequences to determine the degree of local distortion caused by adducts which could influence the ability of nucleotide excision repair. We show that BPDE adducted codon 157 has greater structural distortion than other TP53 G:C>T:A hotspot sites and that sequence context more distal to adjacent bases must influence local distortion. Using TP53 trinucleotide mutation signatures for lung cancer in smokers and non-smokers we further show that codons 157 and 273 have the highest mutation probability in smokers. Combining this information with adduct structural data we predict that G:C>T:A mutations at codon 157 in lung tumours of smokers are predominantly caused by BPDE. Our results provide insight into how different DNA sequence contexts show variability in DNA distortion at mutagen adduct sites that could compromise DNA repair at well characterized cancer related mutation hotspots.

Disturbance in testosterone production in leydig cells by polycyclic aromatic hydevrepocarbons

Polycyclic aromatic hydevrepocarbons (PAHs), which are ubiquitous in the air, are present as volatile and particulate pollutants that result from incomplete combustion. Most PAHs have toxic, mutagenic, and/or carcinogenic properties. Among PAHs, benzo[a]pyrene (B[a]P) and dimethylbenz[a]anthracene (DMBA) are suspected endocrine disruptors. The testis is an important target for PAHs, yet effects on steroidogenesis in Leydig cells are yet to be ascertained. Particularly, disruption of testosterone production by these chemicals can result in serious defects in male reproduction. Exposure to B[a]P reduced serum and intratesticular fluid testosterone levels in rats. Of note, the testosterone level reductions were accompanied by decreased steroidogenic acute regulatory protein (StAR) and 3β-hydevrepoxysteroid dehydevrepogenase isomerase (3β-HSD) expression in Leydig cells. B[a]P exposure can decrease epididymal sperm quality, possibly by disturbing the testosterone level. StAR may be a key steroidogenic protein that is targeted by B[a]P or other PAHs.

Disturbance in testosterone production in leydig cells by polycyclic aromatic hydevrepocarbons

Polycyclic aromatic hydrocarbons (PAHs), which are ubiquitous in the air, are present as volatile and particulate pollutants that result from incomplete combustion. Most PAHs have toxic, mutagenic, and/or carcinogenic properties. Among PAHs, benzo[a]pyrene (B[a]P) and dimethylbenz[a]anthracene (DMBA) are suspected endocrine disruptors. The testis is an important target for PAHs, yet effects on steroidogenesis in Leydig cells are yet to be ascertained. Particularly, disruption of testosterone production by these chemicals can result in serious defects in male reproduction. Exposure to B[a]P reduced serum and intratesticular fluid testosterone levels in rats. Of note, the testosterone level reductions were accompanied by decreased steroidogenic acute regulatory protein (StAR) and 3β-hydroxysteroid dehydrogenase isomerase (3β-HSD) expression in Leydig cells. B[a]P exposure can decrease epididymal sperm quality, possibly by disturbing the testosterone level. StAR may be a key steroidogenic protein that is targeted by B[a]P or other PAHs.

Noncovalent interactions between hydroxylated polycyclic aromatic hydrocarbon and DNA: molecular docking and QSAR study

Polycyclic aromatic hydrocarbons (PAHs) can be hydroxylated by CYP450-oxidases (1A1 and 1B1 mainly) and may cause DNA damage and cancer. However, the mechanism of such interactions has not been fully understood. In this study, an integrated molecular docking and QSAR approach was employed to further investigate the binding interactions between hydroxylated PAHs (HO-PAHs) and calf thymus DNA (CT-DNA). Molecular docking, hydrogen-bonding, hydrophobic and π-π interactions were observed to be characteristic interactions between HO-PAHs and DNA. An optimum QSAR model with good robustness and predictability was developed based on the molecular structural parameters calculated by the density function theory and partial least squares. Additionally, the developed QSAR model indicated that the molecular size, polarizability and electrostatic potential of HO-PAHs were related to the binding affinities to DNA.

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

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

Identification and quantitation of benzo[a]pyrene-derived DNA adducts formed at low adduction level in mice lung tissue

The two major metabolic pathways of benzo[a]pyrene (BP) that lead to DNA lesions are monooxygenation that results in diolepoxides (BPDE) and one-electron oxidation that yields a BP radical cation. These pathways result in formation of stable and depurinating DNA adducts, respectively. Most in vivo animal studies with BP, however, have employed dosage/DNA adduct levels several orders of magnitude higher than the DNA damage level expected from environmentally relevant exposures. Presented are results of experiments in which A/J strain mice were intraperitoneally exposed to 50-microg/g doses of BP. It is shown that non-line-narrowed fluorescence and fluorescence line-narrowing spectroscopies possess the selectivity and sensitivity to distinguish between helix-external, base-stacked, and intercalated conformations of DNA-BPDE adducts formed in lung tissue. Concentrations measured by 32P postlabeling 2 and 3 days after intraperitoneal injection were 420-430 and 600-830 amol BPDE-type adducts per microg DNA. The external and base-stacked conformations are attributed mainly to (+)-trans-anti-BPDE-N2dG and the intercalated conformations to (+)-cis-anti adducts. A stable adduct derived from 9-OH-BP-4,5-epoxide was also detected at a concentration about a factor of 10 lower than the above concentrations. The DNA supernatants were analyzed for the presence of depurinating BP-derived adducts by capillary electrophoresis laser-induced fluorescence and high-performance liquid chromatography mass spectrometry.

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

Analysis of DNA and protein adducts of benzo[a]pyrene in human tissues using structure-specific methods

We review studies which investigate the presence, using structure-specific analytical methods, of DNA or protein adducts of the carcinogen benzo[a]pyrene (BaP) in human tissues. The analytical methods include high performance liquid chromatography with fluorescence detection and gas chromatography-mass spectrometry. Although, for DNA detection these methods are somewhat less sensitive than non-specific techniques such as 32P-postlabeling and immunoassay, they have the distinct advantage of providing reliable structural information. In order to achieve adequate sensitivity, these methods often require the use of fairly large amounts of DNA (>100 microg) or protein (50-100mg). Most studies reviewed here measured tetraols released from DNA or protein by hydrolysis of adducts derived from (7R,8S)-dihydroxy-(9S,10R)-epoxy-7,8,9,10-tetrahydrobenzo[a]pyrene (BPDE), a major ultimate carcinogen of BaP. BPDE-DNA adducts were detected in 39% of 705 samples analyzed. BPDE-protein adducts were found in 59% of 772 samples. There was no single exposure situation that led to an overwhelming presence of detectable adducts. For example, BPDE-DNA adducts were detected in 45% of smokers, 33% of former smokers, 52% of non-smokers, 39% of occupationally exposed individuals, and 34% of environmentally exposed people. Adduct levels were influenced by polymorphisms in carcinogen metabolizing genes such as GSTM1, the presence of which was frequently protective. The relatively high occurrence of non-detectable adducts may result from low levels of BaP exposure and host factors such as genetic polymorphisms. Our analysis demonstrates that the presence of BaP adducts in human tissues cannot be assumed, even in situations where exposure to BaP is relatively high.

Detecting rare mutations associated with cancer risk

For more than a decade, investigators have been searching for a means of determining the risk of individuals developing cancer by detecting rare oncogenic mutations. The accumulation of mutations and the clonal evolvement of tumors provide opportunities for monitoring disease development and intervening prior to the presentation of clinical symptoms, or determining the risk of disease relapse during remission. A number of techniques, mostly polymerase chain reaction (PCR)-based, have been developed that enable the detection of rare oncogenic mutations within the range of 10(-2) to 10(-4) wild-type cells. Only a handful of procedures enable the detection of intragenic single base mutations at one mutant in 10-6 or better. These ultra-sensitive mutation detection techniques have produced some interesting results regarding single base mutation spectra and frequencies in p53, Harvey-ras, N-ras, and other reporter genes and DNA sequences in human tissues. Although there is evidence that some individuals may harbor cells or clones expressing genomic instability, the connection with the processes of carcinogenesis is still tenuous. There remains a need for rigorous epidemiological studies employing these ultra-sensitive mutation detection procedures. Since genomic instability is considered key to tumor development, the relevance of the detection of hypermutable clones in individuals is discussed in the context of cancer risk.

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

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

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

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

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

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

Potent inhibition of werner and bloom helicases by DNA minor groove binding drugs

Maintenance of genomic integrity is vital to all organisms. A number of human genetic disorders, including Werner Syndrome, Bloom Syndrome and Rothmund-Thomson Syndrome, exhibit genomic instability with some phenotypic characteristics of premature aging and cancer predisposition. Presumably the aberrant cellular and clinical phenotypes in these disorders arise from defects in important DNA metabolic pathways such as replication, recombination or repair. These syndromes are all characterized by defects in a member of the RecQ family of DNA helicases. To obtain a better understanding of how these enzymes function in DNA metabolic pathways that directly influence chromosomal integrity, we have examined the effects of non-covalent DNA modifications on the catalytic activities of purified Werner (WRN) and Bloom (BLM) DNA helicases. A panel of DNA-binding ligands displaying unique properties for interacting with double helical DNA was tested for their effects on the unwinding activity of WRN and BLM helicases on a partial duplex DNA substrate. The levels of inhibition by a number of these compounds were distinct from previously reported values for viral, prokaryotic and eukaryotic helicases. The results demonstrate that BLM and WRN proteins exhibit similar sensitivity profiles to these DNA-binding ligands and are most potently inhibited by the structurally related minor groove binders distamycin A and netropsin (K(i) Collapse

Key Words Collapse

MESH Headings

• Adenosine Triphosphatases/antagonists & inhibitors
• Adenosine Triphosphatases/chemistry
• Bloom Syndrome/enzymology
• DNA Helicases/antagonists & inhibitors
• DNA Helicases/chemistry
• Distamycins/pharmacology
• Enzyme Inhibitors/chemistry
• Enzyme Inhibitors/pharmacology
• Exodeoxyribonucleases
• Humans
• Intercalating Agents/chemistry
• Intercalating Agents/pharmacology
• Kinetics
• Ligands
• Netropsin/pharmacology
• RecQ Helicases
• Recombinant Proteins/antagonists & inhibitors
• Topoisomerase I Inhibitors
• Werner Syndrome/enzymology
• Werner Syndrome Helicase

Collapse

Grants Collapse

Affiliation(s)

R M Brosh Laboratory of Molecular Genetics, Box 1, National Institute on Aging, GRC, NIH, 5600 Nathan Shock Drive, Baltimore, MD 21224, USA
J K Karow
E J White
N D Shaw
I D Hickson
V A Bohr

Collapse

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.

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.