Base Substitution Mutagenesis Research Articles

Class-specific recognition and monitoring of environmental steroid estrogens in real water systems utilizing aptamer base substitution mutagenesis approach

The same class of environmental steroid estrogens (SEs) with the highest estrogenic activity share the same chemical core structure and are often found together in the environment, posing significant risks to organismal health and environmental safety due to toxicity accumulation. In this study, a novel method for constructing the group-targeting aptasensor was developed to comprehensively analyze SEs. Through artificial intervention base substitution mutagenesis of adjacent bases T13 and C20 of the aptamer-binding domain recognizing 17β-estradiol, combined with docking calculations, the group-targeting SEs-aptamer for class-specific recognition SEs, such as estrone, estradiol, estriol, and ethinylestradiol were obtained. The binding constant of the SEs-aptamer to the SEs was 108 M−1. The established group-targeting SEs aptasensor exhibited high sensitivity within a concentration range from 0.1 to 10 nM and demonstrated strong interference resistance, as well as high stability and wide pH water applicability. It was further applied to analyze real water samples and monitor changes in SEs concentrations during the removal process by Chlorella pyrenoidosa. These successful applications have demonstrated the excellent ability of this aptasensor to monitor SE in the environment. The method offered a new approach and idea for recognizing and detecting the same class of environmental pollutants in complex systems.

Spontaneous mutagenesis in human cells is controlled by REV1-Polymerase ζ and PRIMPOL

Translesion DNA synthesis (TLS) facilitates replication over damaged or difficult-to-replicate templates by employing specialized DNA polymerases. We investigate the effect on spontaneous mutagenesis of three main TLS control mechanisms: REV1 and PCNA ubiquitylation that recruit TLS polymerases and PRIMPOL that creates post-replicative gaps. Using whole-genome sequencing of cultured human RPE-1 cell clones, we find that REV1 and Polymerase ζ are wholly responsible for one component of base substitution mutagenesis that resembles homologous recombination deficiency, whereas the remaining component that approximates oxidative mutagenesis is reduced in PRIMPOL-/- cells. Small deletions in short repeats appear in REV1-/-PCNAK164R/K164R double mutants, revealing an alternative TLS mechanism. Also, 500-5,000bp deletions appear in REV1-/- and REV3L-/- mutants, and chromosomal instability is detectable in REV1-/-PRIMPOL-/- cells. Our results indicate that TLS protects the genome from deletions and large rearrangements at the expense of being responsible for the majority of spontaneous base substitutions.

The DNA relaxation-dependent OFF-to-ON biasing of the type 1 fimbrial genetic switch requires the Fis nucleoid-associated protein.

The structural genes expressing type 1 fimbriae in Escherichia coli alternate between expressed (phase ON) and non-expressed (phase OFF) states due to inversion of the 314 bp fimS genetic switch. The FimB tyrosine integrase inverts fimS by site-specific recombination, alternately connecting and disconnecting the fim operon, encoding the fimbrial subunit protein and its associated secretion and adhesin factors, to and from its transcriptional promoter within fimS. Site-specific recombination by the FimB recombinase becomes biased towards phase ON as DNA supercoiling is relaxed, a condition that occurs when bacteria approach the stationary phase of the growth cycle. This effect can be mimicked in exponential phase cultures by inhibiting the negative DNA supercoiling activity of DNA gyrase. We report that this bias towards phase ON depends on the presence of the Fis nucleoid-associated protein. We mapped the Fis binding to a site within the invertible fimS switch by DNase I footprinting. Disruption of this binding site by base substitution mutagenesis abolishes both Fis binding and the ability of the mutated switch to sustain its phase ON bias when DNA is relaxed, even in bacteria that produce the Fis protein. In addition, the Fis binding site overlaps one of the sites used by the Lrp protein, a known directionality determinant of fimS inversion that also contributes to phase ON bias. The Fis–Lrp relationship at fimS is reminiscent of that between Fis and Xis when promoting DNA relaxation-dependent excision of bacteriophage λ from the E. coli chromosome. However, unlike the co-binding mechanism used by Fis and Xis at λ attR, the Fis–Lrp relationship at fimS involves competitive binding. We discuss these findings in the context of the link between fimS inversion biasing and the physiological state of the bacterium.

BRCA1 deficiency specific base substitution mutagenesis is dependent on translesion synthesis and regulated by 53BP1

Defects in BRCA1, BRCA2 and other genes of the homology-dependent DNA repair (HR) pathway cause an elevated rate of mutagenesis, eliciting specific mutation patterns including COSMIC signature SBS3. Using genome sequencing of knock-out cell lines we show that Y family translesion synthesis (TLS) polymerases contribute to the spontaneous generation of base substitution and short insertion/deletion mutations in BRCA1 deficient cells, and that TLS on DNA adducts is increased in BRCA1 and BRCA2 mutants. The inactivation of 53BP1 in BRCA1 mutant cells markedly reduces TLS-specific mutagenesis, and rescues the deficiency of template switch–mediated gene conversions in the immunoglobulin V locus of BRCA1 mutant chicken DT40 cells. 53BP1 also promotes TLS in human cellular extracts in vitro. Our results show that HR deficiency–specific mutagenesis is largely caused by TLS, and suggest a function for 53BP1 in regulating the choice between TLS and error-free template switching in replicative DNA damage bypass.

From RNA World to SARS-CoV-2: The Edited Story of RNA Viral Evolution.

The current SARS-CoV-2 pandemic underscores the importance of understanding the evolution of RNA genomes. While RNA is subject to the formation of similar lesions as DNA, the evolutionary and physiological impacts RNA lesions have on viral genomes are yet to be characterized. Lesions that may drive the evolution of RNA genomes can induce breaks that are repaired by recombination or can cause base substitution mutagenesis, also known as base editing. Over the past decade or so, base editing mutagenesis of DNA genomes has been subject to many studies, revealing that exposure of ssDNA is subject to hypermutation that is involved in the etiology of cancer. However, base editing of RNA genomes has not been studied to the same extent. Recently hypermutation of single-stranded RNA viral genomes have also been documented though its role in evolution and population dynamics. Here, we will summarize the current knowledge of key mechanisms and causes of RNA genome instability covering areas from the RNA world theory to the SARS-CoV-2 pandemic of today. We will also highlight the key questions that remain as it pertains to RNA genome instability, mutations accumulation, and experimental strategies for addressing these questions.

Two main mutational processes operate in the absence of DNA mismatch repair

The analysis of tumour genome sequences has demonstrated high rates of base substitution mutagenesis upon the inactivation of DNA mismatch repair (MMR), and the resulting somatic mutations in MMR deficient tumours appear to significantly enhance the response to immune therapy. A handful of different algorithmically derived base substitution mutation signatures have been attributed to MMR deficiency in tumour somatic mutation datasets. In contrast, mutation data obtained from whole genome sequences of isogenic wild type and MMR deficient cell lines in this study, as well as from published sources, show a more uniform experimental mutation spectrum of MMR deficiency. In order to resolve this discrepancy, we reanalysed mutation data from MMR deficient tumour whole exome and whole genome sequences. We derived two base substitution signatures using non-negative matrix factorisation, which together adequately describe mutagenesis in all tumour and cell line samples. The two new signatures broadly resemble COSMIC signatures 6 and 20, but perform better than existing COSMIC signatures at identifying MMR deficient tumours in mutation signature deconstruction. We show that the contribution of the two identified signatures, one of which is dominated by C to T mutations at CpG sites, is biased by the different sequence composition of the exome and the whole genome. We further show that the identity of the inactivated MMR gene, the tissue type, the mutational burden or the patient’s age does not influence the mutation spectrum, but that a tendency for a greater contribution by the CpG mutational process is observed in tumours as compared to cultured cells. Our analysis suggest that two separable mutational processes operate in the genomes of MMR deficient cells.

Proteomic Analysis of DNA Synthesis on a Structured DNA Template in Human Cellular Extracts: Interplay Between NHEJ and Replication-Associated Proteins.

It is established that short inverted repeats trigger base substitution mutagenesis in human cells. However, how the replication machinery deals with structured DNA is unknown. It has been previously reported that in human cell-free extracts, DNA primer extension using a structured single-stranded template is transiently blocked at DNA hairpins. Here, the proteomic analysis of proteins bound to the DNA template is reported and evidence that the DNA-PK complex (DNA-PKcs and the Ku heterodimer) recognizes, and is activated by, structured single-stranded DNA is provided. Hijacking the DNA-PK complex by double-stranded oligonucleotides results in a large removal of the pausing sites and an elevated DNA extension efficiency. Conversely, DNA-PKcs inhibition results in its stabilization on the template, along with other proteins acting downstream in the Non-Homologous End-Joining (NHEJ) pathway, especially the XRCC4-DNA ligase 4 complex and the cofactor PAXX. Retention of NHEJ factors to the DNA in the absence of DNA-PKcs activity correlates with additional halts of primer extension, suggesting that these proteins hinder the progression of the DNA synthesis at these sites. Overall these results raise the possibility that, upon binding to hairpins formed onto ssDNA during fork progression, the DNA-PK complex interferes with replication fork dynamics in vivo.

Pol \u03bc dGTP mismatch insertion opposite T coupled with ligation reveals promutagenic DNA repair intermediate

Incorporation of mismatched nucleotides during DNA replication or repair leads to transition or transversion mutations and is considered as a predominant source of base substitution mutagenesis in cancer cells. Watson-Crick like dG:dT base pairing is considered to be an important source of genome instability. Here we show that DNA polymerase (pol) μ insertion of 7,8-dihydro-8′-oxo-dGTP (8-oxodGTP) or deoxyguanosine triphosphate (dGTP) into a model double-strand break DNA repair substrate with template base T results in efficient ligation by DNA ligase. These results indicate that pol μ-mediated dGTP mismatch insertion opposite template base T coupled with ligation could be a feature of mutation prone nonhomologous end joining during double-strand break repair.

A Polymerase Tautomeric Model for Radiation Induced Bystander Effects: A Model for Untargeted Base Substitution Mutagenesis during Error Prone and SOS Replication of Double Stranded DNA Containing Thymine and Adenine in Rare Tautomeric Forms

A Polymerase Tautomeric Model for Radiation Induced Bystander Effects: A Model for Untargeted Base Substitution Mutagenesis during Error Prone and SOS Replication of Double Stranded DNA Containing Thymine and Adenine in Rare Tautomeric Forms

A Polymerase Tautomeric Model for Radiation Induced Bystander Effects: A Model for Untargeted Base Substitution Mutagenesis during Error Prone and SOS Replication of Double Stranded DNA Containing Thymine and Adenine in Rare Tautomeric Forms

Currently untargeted mutations are studied in the context of bystander effects Up to now the nature and mechanisms of formation of untargeted mutations is poorly understood Untargeted base substitution mutations are base substitution mutations then one or some nucleotides are inserted in DNA molecule on so called undamaged sites of DNA Here the polymerase tautomeric models of ultraviolet mutagenesis and bystander effects are developed The mechanisms of untargeted base substitution mutations formation caused by cis syncyclobutane thymine dimers and bases pairs of adenine thymine or thymine adenine in rare tautomeric forms that are in small neighbor from any cyclobutane dimers are proposed Five rare tautomeric forms may form for thymine and adenine These rare tautomeric forms will be stable if corresponding nucleotides are part of cyclobutane pyrimidine dimers or are in small neighbor of the cyclobutane dimers and during DNA synthesis Structural analysis of bases incorporation shown that adenine in rare tautomeric form A may result in untargeted homologous transversion reg Thymine in rare tautomeric form T may result in untargeted transition reg G C or untargeted homologous transversion reg Molecules of the thymine in rare tautomeric form T may result in transversion A reg G only Thymine in rare tautomeric form T may result in untargeted transversion reg G or untargeted homologous transversion reg Thus molecules of adenine and thymine in definite tautomeric forms are in small neighbor from any cyclobutane pyrimidine dimers may result in untargeted base substitution mutations

Loss of BRCA1 or BRCA2 markedly increases the rate of base substitution mutagenesis and has distinct effects on genomic deletions

Loss-of-function mutations in the BRCA1 and BRCA2 genes increase the risk of cancer. Owing to their function in homologous recombination repair, much research has focused on the unstable genomic phenotype of BRCA1/2 mutant cells manifest mainly as large-scale rearrangements. We used whole-genome sequencing of multiple isogenic chicken DT40 cell clones to precisely determine the consequences of BRCA1/2 loss on all types of genomic mutagenesis. Spontaneous base substitution mutation rates increased sevenfold upon the disruption of either BRCA1 or BRCA2, and the arising mutation spectra showed strong and specific correlation with a mutation signature associated with BRCA1/2 mutant tumours. To model endogenous alkylating damage, we determined the mutation spectrum caused by methyl methanesulfonate (MMS), and showed that MMS also induces more base substitution mutations in BRCA1/2-deficient cells. Spontaneously arising and MMS-induced insertion/deletion mutations and large rearrangements were also more common in BRCA1/2 mutant cells compared with the wild-type control. A difference in the short deletion phenotypes of BRCA1 and BRCA2 suggested distinct roles for the two proteins in the processing of DNA lesions, as BRCA2 mutants contained more short deletions, with a wider size distribution, which frequently showed microhomology near the breakpoints resembling repair by non-homologous end joining. An increased and prolonged gamma-H2AX signal in MMS-treated BRCA1/2 cells suggested an aberrant processing of stalled replication forks as the cause of increased mutagenesis. The high rate of base substitution mutagenesis demonstrated by our experiments is likely to significantly contribute to the oncogenic effect of the inactivation of BRCA1 or BRCA2.

Mechanisms of Base Substitution Mutagenesis in Cancer Genomes

Cancer genome sequence data provide an invaluable resource for inferring the key mechanisms by which mutations arise in cancer cells, favoring their survival, proliferation and invasiveness. Here we examine recent advances in understanding the molecular mechanisms responsible for the predominant type of genetic alteration found in cancer cells, somatic single base substitutions (SBSs). Cytosine methylation, demethylation and deamination, charge transfer reactions in DNA, DNA replication timing, chromatin status and altered DNA proofreading activities are all now known to contribute to the mechanisms leading to base substitution mutagenesis. We review current hypotheses as to the major processes that give rise to SBSs and evaluate their relative relevance in the light of knowledge acquired from cancer genome sequencing projects and the study of base modifications, DNA repair and lesion bypass. Although gene expression data on APOBEC3B enzymes provide support for a role in cancer mutagenesis through U:G mismatch intermediates, the enzyme preference for single-stranded DNA may limit its activity genome-wide. For SBSs at both CG:CG and YC:GR sites, we outline evidence for a prominent role of damage by charge transfer reactions that follow interactions of the DNA with reactive oxygen species (ROS) and other endogenous or exogenous electron-abstracting molecules.

PLEKHG2/FLJ00018, a Rho family-specific guanine nucleotide exchange factor, is tyrosine phosphorylated via the EphB2/cSrc signaling pathway

PLEKHG2/FLJ00018, a Rho family-specific guanine nucleotide exchange factor (RhoGEF), is activated by heterotrimeric GTP-binding protein (G protein) Gβγ subunits, and in turn activates the small G protein Rac and Cdc42, which have been shown to mediate signaling pathways leading to actin cytoskeletal reorganization. In the present study, we show that co-expression of the constitutively active mutant of cSrc, a non-receptor tyrosine kinase, and PLEKHG2 induced the tyrosine phosphorylation of PLEKHG2 in HEK293 cells. Through deletion and base substitution mutagenesis we have identified Tyr489 of PLEKHG2 as the site phosphorylated by cSrc. Furthermore, using a high-throughput src homology 2 (SH2) domain binding assay, the SH2 domain of ABL1 and the PI 3-kinse regulator subunit (PIK3R3) were identified as candidates for the binding partner of tyrosine-phosphorylated PLEKHG2. The interaction between PLEKHG2 and the full-length of PIK3R3, but not ABL1, occurs in a tyrosine-phosphorylation-dependent manner. Furthermore, PLEKHG2 is tyrosine phosphorylated at Tyr489 by ephrinB2 receptor signaling via cSrc. Investigation of the physiological function of tyrosine phosphorylation at Tyr489 in PLEKHG2 remains a subject for future studies.

TaMYB13 is a transcriptional activator of fructosyltransferase genes involved in β‐2,6‐linked fructan synthesis in wheat

Fructans are soluble fructosyl-oligosaccharides deposited in many cool-season grass species as a carbon reserve; they are synthesised by fructosyltransferases. In wheat and barley fructans can accumulate in mature stems at a very high level and serve as an important carbon source for grain filling. Fructan synthesis in temperate cereals is regulated by sucrose level and developmental signals, and functions as a metabolic adjustment for carbon balance between carbon supply and sink demand. In this study the expression levels of a highly homologous group of Triticum aestivumMYB genes (TaMYB13-1, TaMYB13-2 and TaMYB13-3) were found to be positively correlated with the mRNA levels of sucrose:sucrose 1-fructosyltransferase (1-SST) and sucrose:fructan 6-fructosyltransferase (6-SFT) in wheat stems among recombinant inbred lines with a wide range of fructan concentrations through Affymetrix array expression analysis. This expression correction extended to expression profiles during stem development. TaMYB13 contains an R2R3-type MYB domain. Invitro random DNA-binding site selection followed by base substitution mutagenesis revealed that TaMYB13 bound to a (A/G/T)TT(A/T/C)GGT core sequence, which was present in the promoters of wheat Ta1-SST and Ta6-SFT genes as well as a barley Hv6-SFT gene. Transactivation analysis showed that TaMYB13 was a transcriptional activator and could markedly enhance the expression of 1-SST and 6-SFT promoter-driven reporter genes in wheat. Elimination of TaMYB13-binding sites in Ta6-SFT and Ta1-SST promoters markedly reduced TaMYB13-mediated reporter gene transactivation. These data suggest that TaMYB13 and its orthologues are positive regulators for controlling the expression of major fructosyltransferases involved in the fructan synthetic pathway in temperate cereals.

Regulation of Odor Receptor Genes in Trichoid Sensilla of the Drosophila Antenna

This study concerns the problem of odor receptor gene choice in the fruit fly Drosophila melanogaster. From a family of 60 Odor receptor genes, only one or a small number are selected for expression by each olfactory receptor neuron. Little is known about how an olfactory receptor neuron selects a receptor, or how the nucleotide sequences flanking a receptor gene dictate its expression in a particular neuron. Previous investigation has primarily concerned the maxillary palp, the simpler of the fly's two olfactory organs. Here we focus on genes encoding four antennal receptors that respond to fly odors in an in vivo expression system. To investigate the logic of odor receptor expression, we carry out a genetic analysis of their upstream regulatory sequences. Deletion analysis reveals that relatively short regulatory regions are sufficient to confer expression in the appropriate neurons, with limited if any misexpression. We find evidence for both positive and negative regulation. Multiple repressive functions restrict expression to the antenna, to a region of the antenna, and to neurons. Through deletion and base substitution mutagenesis we identify GCAATTA elements and find evidence that they act in both positive and negative regulation.

ENU-Based Gene-Driven Mutagenesis in the Mouse: A Next-Generation Gene-Targeting System

As a new mouse mutant resource, the RIKEN ENU-based gene-driven mutagenesis system in the mouse has been available to the research community since 2002. By using random base-substitution mutagenesis with ENU, a new reverse genetics infrastructure has been developed as a next-generation gene-targeting system. The construction of a large-scale mutant mouse library and high-throughput mutation discovery systems were the keys making it practically feasible. The RIKEN mutant mouse library consists of ~ 10,000 G1 mice, within which 100-150 mutant strains have been established based on users' requests every year. Use of the system is very simple: users 1) download an application form from our web site and send to us, and 2) design the PCR primers for the target gene. Then, we screen the RIKEN mutant mouse library and report all the detected mutations to the user. From among the allelic series of discovered mutations, users decide which mutant strain(s) to analyze and request the live mutant strain for functional studies of the target gene. Users have been reporting various functional mutations in the RIKEN mutant mouse library: e.g., missense, knockout-type and even functional non-coding mutations. In the near future, next-generation re-sequencing systems should drastically enhance the utility of the ENU-based gene-driven mutagenesis not only for the mouse but also for other species.

RAD51D- and FANCG-dependent base substitution mutagenesis at the ATP1A1 locus in mammalian cells

Elaborate processes act at the DNA replication fork to minimize the generation of chromatid discontinuity when lesions are encountered. To prevent collapse of stalled replication forks, mutagenic translesion synthesis (TLS) polymerases are recruited temporarily to bypass DNA lesions. When a replication-associated (one-ended) double-strand break occurs, homologous recombination repair (HRR) can restore chromatid continuity in what has traditionally been regarded as an “error-free” process. Our previous mutagenesis studies show an important role for HRR in preventing deletions and rearrangements that would otherwise result from error-prone nonhomologous end joining (NHEJ) after fork breakage. An analogous, but distinct, role in minimizing mutations is attributed to the proteins defective in the cancer predisposition disease Fanconi anemia (FA). Cells from FA patients and model systems show an increased proportion of gene-disrupting deletions at the hprt locus as well as decreased mutation rates in the hprt assay, suggesting a role for the FANC proteins in promoting TLS, HRR, and possibly also NHEJ. It remains unclear whether HRR, like the FANC pathway, impacts the rate of base substitution mutagenesis. Therefore, we measured, in isogenic rad51d and fancg CHO mutants, mutation rates at the Na +/K +-ATPase α-subunit ( ATP1A1) locus using ouabain resistance, which specifically detects base substitution mutations. Surprisingly, we found that the spontaneous mutation rate was reduced ∼2.5-fold in rad51d knockout cells, an even greater extent than observed in fancg cells, when compared with parental and isogenic gene-complemented control lines. A ∼2-fold reduction in induced mutations in rad51d cells was seen after treatment with the DNA alkylating agent ethylnitrosurea while a lesser reduction occurred in fancg cells. Should the model ATP1A1 locus be representative of the genome, we conclude that at least 50% of base substitution mutations in this mammalian system arise through error-prone polymerase(s) acting during HRR-mediated restart of broken replication forks.

TaNAC69 from the NAC superfamily of transcription factors is up-regulated by abiotic stresses in wheat and recognises two consensus DNA-binding sequences.

NAC proteins are one of the largest families of plant transcription factors and have recently been implicated in diverse physiological processes. To elucidate their role in gene regulation, we determined the DNA-binding specificity of a drought- and cold-inducible NAC protein, TaNAC69 from wheat, and analysed its homologues from other species. Two consensus DNA-binding sequences (spanning 23-24 bp) of TaNAC69 were identified through binding site selection and both consisted of two half sites. Comprehensive data on the DNA-binding specificity of TaNAC69 were generated through extensive base substitution mutagenesis. TaNAC69 and its homologue in Arabidopsis, NAP, sharing 75% sequence identity in the NAC domain, exhibited similar DNA-binding specificity. TaNAC69 was able to homodimerise through its NAC domain. The NAC domain consists of five conserved subdomains. Subdomain mutation showed that a loss or reduction in TaNAC69 dimerisation capacity was accompanied with abolition or decrease in its DNA-binding activity. These data suggest that all subdomains are necessary to maintain a functional NAC domain structure required for interaction with DNA and dimerisation.

Identification of mutator genes and mutational pathways in Escherichia coli using a multicopy cloning approach.

We searched for genes that create mutator phenotypes when put on to a multicopy plasmid in Escherichia coli. In many cases, this will result in overexpression of the gene in question. We constructed a random shotgun library with E. coli genomic fragments between 3 and 5 kbp in length on a multicopy plasmid vector that was transformed into E. coli to screen for frameshift mutators. We identified a total of 115 independent genomic fragments that covered 17 regions on the E. coli chromosome. Further studies identified 12 genes not previously known as causing mutator phenotypes when overproduced. A striking finding is that overproduction of the multidrug resistance transcription regulator, EmrR, results in a large increase in frameshift and base substitution mutagenesis. This suggests a link between multidrug resistance and mutagenesis. Other identified genes include those encoding DNA helicases (UvrD, RecG, RecQ), truncated forms of the DNA mismatch repair protein (MutS) and a primosomal component (DnaT), a negative modulator of initiation of replication/GATC-binding protein (SeqA), a stationary phase regulator AppY, a transcriptional regulator PaaX and three putative open reading frames, ycgW, yfjY and yjiD, encoding hypothetical proteins. In addition, we found three genes encoding proteins that were previously known to cause mutator effects under overexpression conditions: error-prone polymerase IV (DinB), DNA methylase (Dam) and sigma S factor (RpoS). This genomic strategy offers an approach to identify novel mutator effects resulting from the multicopy cloning (MCC) of specific genes and therefore complementing the conventional gene inactivation approach to finding mutators.

Molecular modeling of the major benzo[ a]pyrene N 2-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–N 2-dG when bypassed during DNA replication (probably by different DNA polymerases). Previous molecular modeling studies suggested that B[ a]P–N 2-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–N 2-dG adducts of differing stereoisomers and DNA sequence contexts, as well as in five cases where the conformation is known by NMR. Of the ∼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–N 2-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).