Structural and biochemical impact of C8-aryl-guanine adducts within the NarI recognition DNA sequence: influence of aryl ring size on targeted and semi-targeted mutagenicity

Structure and function of extreme TLS DNA polymerase TTEDbh from Thermoanaerobacter tengcongensis

Translesion synthesis (TLS) is a kind of DNA repair that maintains the stability of the genome and ensures the normal growth of life in cells under emergencies. Y-family DNA polymerases, as a kind of error-prone DNA polymerase, mainly perform TLS. Previous studies have suggested that the occurrence of tumors is associated with the overexpression of human DNA polymerase of the Y family. And the combination of Y-family DNA polymerase inhibitors is promising for cancer therapy. Here we report the functional and structural characterization of a member of the Y-family DNA polymerases, TTEDbh. We determine TTEDbh is an extreme TLS polymerase that can cross oxidative damage sites, and further identify the amino acids and novel structures that are critical for DNA binding, synthesis, fidelity, and oxidative damage bypass. Moreover, previously unnoticed structural elements with important functions have been discovered and analyzed. These studies provide a more experimental basis for further elucidating the molecular mechanisms of DNA polymerase in the Y family. It could also shed light on the design of drugs to target tumors.

Effect of Guanine Adduct Size, Shape, and Linker Type on the Conformation of Adducted DNA: A DFT and Molecular Dynamics Study

DNA is damaged through various exogenous sources (e.g., automobile exhaust, tobacco smoke, and processed foods), which can yield diverse C8-dG bulky aryl adducts. Adducts are known to induce structural changes to DNA that can lead to various biological outcomes, ranging from cell death to diseases such as cancer. Unfortunately, the relationship between the chemical composition of the damaged product, the adducted DNA structure, and the biological consequences is not well understood, which limits the development of disease detection and prevention strategies. The present study uses density functional theory (DFT) calculations and quintuplicate 1 μs molecular dynamics (MD) simulations to characterize the structure of DNA containing 21 model C8-dG adducts that systematically differ in size (phenyl to pyrenyl), shape (α (2,3), β (3,4) fusion, or ring substitution), and nucleobase-aryl group linkage (N, O, and C-linked). DFT calculations reveal that the inherent structural features of the G nucleobase adducts are impacted by linker type and bulky moiety shape, but not size, with the conformational flexibility reducing with α-ring fusion and linker composition as N > O > C. These structural properties are maintained in nucleoside models, which also reveal an increased propensity for anti-to-syn rotation about the glycosidic bond with N < O < C linker type. Although these diverse chemical features do not influence the global structure of adducted DNA, the adducts differentially impact the conformation local to the adducted site, including the relative populations of structures with the bulky moiety in the major groove (B conformer) and intercalated (stacked) into the helix (S conformer). Specifically, while the smallest phenyl adducts favor the B conformation and the largest pyrenyl-derived adducts stabilize the S conformation, the B/S ratio decreases with an increase in ring size and N > O > C linker composition. The shape and size (length) of the adduct can further finetune the B/S ratio, with β-fused naphthyl or α-fused phenanthryl N-linked adducts and O or C-linked adducts containing ring substitution increasing the prevalence of the S adducted DNA conformation. Overall, this work uncovers the significant effect of bulky moiety size and linker type, as well as the lesser impact of aryl group shape, on adducted DNA structure, which suggests differential replication and repair outcomes, and thereby represents an important step toward rationalizing connections between the structure and biological consequences of diverse DNA adducts.

Effects of a Second Local DNA Damage Event on the Toxicity of the Human Carcinogen 4-Aminobiphenyl: A Molecular Dynamics Study of a Damaged DNA Structure

Exposure of humans to carcinogenic aromatic amines (AAs) occurs daily. AAs are bioactivated in cells into products that attack DNA, primarily leading to N-linked C8-dG adducts. Previous work on DNA containing a single AA-derived adduct (monoadducted DNA) has shown a structure-function relationship between the damaged DNA conformation and cellular outcomes. However, relatively little is known about the conformation and biological outcomes of DNA containing two bulky adducts (diadducted DNA) in close proximity. To fill this current void in the literature, the present work uses quintuplet 0.5 μs MD simulations to understand the structural impact of DNA exposure to the potent bladder carcinogen 4-aminobiphenyl (ABP), which is found in cigarette smoke and select dyes, and results in the widely studied N-linked ^ABPdG adduct. Specifically, 18 unique DNA duplexes were investigated that contain one or two ^ABPdG adducts in the anti and/or syn glycosidic orientation(s) in all combinations of three G positions in the NarI mutation hotspot for AAs (5'-G₁G₂CG₃CC). Monoadducted DNA displays sequence-dependent conformational heterogeneity, with the G₁ site having the greatest anti preference, and highlights the range of helical structures associated with the syn lesion orientation [i.e., stacked (S), intercalated (I), and wedge (W) conformations]. Diadducted DNA results in interesting lesion separation effects on the conformational heterogeneity, including a greater anti preference for neighboring adducts (G₁G₂) and a greater syn preference for next-nearest neighbor damaged sites (G₂G₃) compared to monoadducted DNA. As a result, an increase in the number of ^ABPdG adducts changes the conformational heterogeneity of ABP-exposed DNA depending on the relative positions of the lesions and thereby could result in increased or decreased toxicity upon human exposure to elevated levels of ABP.

Understanding the ancillary ligand effect on luminescent cyclometalated Ir(III) complex as a reporter for 2-acetylaminofluorene DNA(AAF-dG) adduct

Mutagenic agents such as aromatic amines undergo metabolic activation and produce DNA adducts at C8 position of guanine bases. N-2-acetylaminofluorene (AAF) generates different mutational outcomes when placed at G₁, G₂, and G₃ of a NarI sequence (-G₁G₂CG₃CC/T-). These outcomes are dictated by the conformations adopted by these adducts. Detection of such lesions is of considerable interest owing to their hazardous effects. Here, we report the synthesis of three cyclometalated [Ir(L)₂dppz]⁺ complexes (L = 2-phenylpyridine (ppy) 1; benzo[h]quinoline (bhq) 2; 2-phenylquinoline (pq) 3; dppz = dipyrido[3,2-a:2',3'-c]phenazine) and their interaction with AAF adducted NarI DNA. Remarkably, complexes 1 and 2 displayed dominant ³LC transition characteristic of polar environment despite binding to the adducted sites. On the other hand, complex 3 binds to NarI sequences and behaves as a luminescent reporter for AAF-modified DNA. The results reported here emphasize that molecular light switching phenomenon can be stimulated by switching ancillary ligands and might act as potential probes for covalent-DNA defects.

Screening Internal Donor-Acceptor Biaryl Nucleobase Surrogates for Turn-On Fluorescence Affords an Aniline-Carboxythiophene Probe for Protein Detection by G-Quadruplex DNA

Donor-acceptor biaryls serve as microenvironment fluorescent sensors with highly quenched intramolecular charge transfer (ICT) emission in polar protic solvents that turns on in aprotic media. In DNA, canonical donor-acceptor fluorescent base analogs can be prepared through on-strand Suzuki-Miyaura cross-coupling reactions involving 8-bromo-2'-deoxyguanosine (8-Br-dG) with an acceptor aryboronic acid. Herein, we demonstrate that replacement of 8-Br-dG with N-methyl-4-bromoaniline (4-Br-An) containing an acyclic N-glycol group can be employed in the on-strand Suzuki-Miyaura reaction to afford new donor-acceptor biaryl nucleobase surrogates with a 40-fold increase in emission intensity for fluorescent readout within single-strand oligonucleotides. Screening the best acceptor for turn-on fluorescence upon duplex formation afforded the carboxythiophene derivative [COOTh]An with a 7.4-fold emission intensity increase upon formation of a single-bulged duplex (-1) with the surrogate occupying a pyrimidine-flanked bulge. Insertion of the [COOTh]An surrogate into the lateral TT loops produced by the antiparallel G-quadruplex (GQ) of the thrombin binding aptamer (TBA) afforded a 4.1-fold increase in probe fluorescence that was accompanied by a 20 nm wavelength shift to the blue upon thrombin binding. The modified TBA afforded a limit of detection of 129 nM for thrombin and displayed virtually no emission response to off-target proteins. The fluorescence response of [COOTh]An to thrombin binding highlights the utility of the thienyl-aniline moiety for monitoring DNA-protein interactions.

Impact of DNA Adduct Size, Number, and Relative Position on the Toxicity of Aromatic Amines: A Molecular Dynamics Case Study of ANdG- and APdG-Containing DNA Duplexes

Human exposure to aromatic amines (AAs) can result in carcinogenic DNA adducts. To complement previous work geared toward understanding the mutagenicity of AA-derived adducts, which has almost exclusively studied (monoadducted) DNA containing a single lesion, the present work provides the first in-depth comparison of the structure of monoadducted and diadducted DNA duplexes. Specifically, molecular dynamics (MD) simulations were initially performed on DNA containing the nonmutagenic single-ringed N-(deoxyguanosin-8-yl)-aniline (^ANdG) or the mutagenic four-ringed N-(deoxyguanosin-8-yl)-1-aminopyrene (^APdG) lesion at G₁, G₂, or G₃ in the AA deletion hotspot (5'-G₁G₂CG₃CC) in the anti or syn glycosidic orientation (B/S duplex conformation). Subsequently, diadducted strands were assessed that span each combination of damaged sites (G₁G₂ (nearest neighbors), G₂G₃ (next-nearest neighbors), and G₁G₃ (two intervening nucleotides)) and anti/syn lesion glycosidic orientations. Despite other N-linked C8-dG adducts exhibiting sequence dependence conformational heterogeneity, a single ^ANdG or ^APdG lesion induces helical conformational homogeneity that is exclusively controlled by aryl moiety size. However, the preferred damaged DNA conformation can change upon the addition of a second adduct depending on lesion separation, with neighboring lesions stabilizing a nonmutagenic conformation and next-nearest damaged sites stabilizing a promutagenic conformation regardless of adduct size. As a result, diadducted DNA is found to adopt conformations that are unfavored for the corresponding monoadducted system, pointing to differential replication and repair outcomes for diadducted DNA compared to those for monoadducted DNA. Thus, although the toxicity of monoadducted DNA is most significantly dictated by lesion size, the toxicity can increase or decrease upon a second damaging event depending on lesion size and relative position. Overall, our work adds the number of lesions and their spatial separation to the growing list of factors that determine the structure and biological outcomes of adducted DNA.

Structural insights into the recognition of DNA defects by small molecules

Studies on the binding interaction of small molecules and nucleic acids have been explored for their biological applications. With excellent photophysical/chemical properties, numerous metal complexes have been studied as structural probes for nucleic acids. The recognition of DNA defects is of high importance due to their association with various types of cancers. Small molecules that target DNA defects in a specific and selective manner offer a new avenue for developing novel drugs and diagnostic tools. Transition metal complexes have been studied as probes for abasic sites and DNA/RNA mismatches. By changing the ligand structure or metal center, the probing efficiency of the metal complexes varies towards the defects. In this perspective, we have discussed mainly the structural requirement of metal complexes as probes for abasic sites, mismatches, and covalent DNA adducts, followed by the challenges and future directions.

Phenoxy Radical Reactivity of Nucleic Acids: Practical Implications for Biotinylation

Recent advances in peroxidase-mediated biotin tyramide (BT) signal amplification technology have resulted in high-resolution and subcellular compartment-specific mapping of protein and RNA localization. Horseradish peroxidase (HRP) in the presence of H2 O2 is known to activate phenolic compounds for phenoxy radical reaction with nucleic acids, where biotinylation by BT is a practical example. BT reactivity with RNA and DNA is not understood in detail. We report that BT phenoxy radicals react in a sequence-independent manner with guanosine bases in RNA. In contrast, DNA reactivity with BT cannot be detected by our methods under the same conditions. Remarkably, we show that fluorescein conjugates DNA rapidly and selectively reacts with BT phenoxy radicals, allowing convenient and practical biotinylation of DNA on fluorescein with retention of fluorescence.

Risk assessment of ochratoxin A in food

The European Commission asked EFSA to update their 2006 opinion on ochratoxin A (OTA) in food. OTA is produced by fungi of the genus Aspergillus and Penicillium and found as a contaminant in various foods. OTA causes kidney toxicity in different animal species and kidney tumours in rodents. OTA is genotoxic both in vitro and in vivo; however, the mechanisms of genotoxicity are unclear. Direct and indirect genotoxic and non-genotoxic modes of action might each contribute to tumour formation. Since recent studies have raised uncertainty regarding the mode of action for kidney carcinogenicity, it is inappropriate to establish a health-based guidance value (HBGV) and a margin of exposure (MOE) approach was applied. For the characterisation of non-neoplastic effects, a BMDL 10 of 4.73 μg/kg body weight (bw) per day was calculated from kidney lesions observed in pigs. For characterisation of neoplastic effects, a BMDL 10 of 14.5 μg/kg bw per day was calculated from kidney tumours seen in rats. The estimation of chronic dietary exposure resulted in mean and 95th percentile levels ranging from 0.6 to 17.8 and from 2.4 to 51.7 ng/kg bw per day, respectively. Median OTA exposures in breastfed infants ranged from 1.7 to 2.6 ng/kg bw per day, 95th percentile exposures from 5.6 to 8.5 ng/kg bw per day in average/high breast milk consuming infants, respectively. Comparison of exposures with the BMDL 10 based on the non-neoplastic endpoint resulted in MOEs of more than 200 in most consumer groups, indicating a low health concern with the exception of MOEs for high consumers in the younger age groups, indicating a possible health concern. When compared with the BMDL 10 based on the neoplastic endpoint, MOEs were lower than 10,000 for almost all exposure scenarios, including breastfed infants. This would indicate a possible health concern if genotoxicity is direct. Uncertainty in this assessment is high and risk may be overestimated.

Molecular Dynamics Study of One-Base Deletion Duplexes Containing the Major DNA Adduct Formed by Ochratoxin A: Effects of Sequence Context and Adduct Ionization State on Lesion Site Structure and Mutagenicity

Ochratoxin A (OTA) is a ubiquitous food toxin associated with chronic nephropathy in humans and renal carcinogenicity in rodents. The mutational spectra of cells exposed to OTA reveal that one-base deletions comprise the largest percentage (73%) of the total mutations that occur upon OTA exposure. To contribute toward understanding the prevalence of OTA-induced one-base deletion mutations, the present work uses molecular dynamics (MD) simulations to analyze the conformational preferences of one-base deletion duplexes containing OT-G, the major OTA adduct (addition product) at the C⁸-site of guanine. Specifically, the influence of OT-G in four possible ionization states and three sequence contexts (G¹, G² and G³ in the NarI (5'-G¹G²CG³CC-3'), a prokaryotic mutational hotspot sequence) on the structure of the adducted DNA is investigated. Our data reveal that the damaged helices are stable in two (B-type (B) and stacked (S)) conformations that are structurally similar to those adopted by common N-linked C⁸-guanine lesions. However, the adduct ionization state and sequence context affect the degree of helical distortion and the B/S conformational heterogeneity, which will impact the lesion repair and replication outcomes. This finding correlates with the experimentally reported tissue-specific mutagenicity of OTA exposure. Furthermore, regardless of the adduct conformation, ionization state, or sequence context, more stable lesion-site interactions and lack of disruption of the flanking base pairs in the one-base deletion duplexes compared to the corresponding two-base deletion helices rationalize the greater abundance of OTA induced one-base deletions. Overall, our work provides valuable structural insights that help explain the experimentally observed mutagenicity associated with OTA.

Adduct Fluorescence as a Tool to Decipher Sequence Impact on Frameshift Mutations Mediated by a C-Linked C8-Biphenyl-Guanine Lesion

Aromatic chemicals can undergo metabolic activation to afford electrophilic species that react at the C8-site of 2'-deoxyguanosine (dG) to generate bulky C8-dG adducts as a basis of initiating carcinogenesis. These DNA lesions have served as models to understand the mechanism of frameshift mutagenesis, especially within CG-dinucleotide repeat sequences, such as NarI (5'-GGCXCC-3', where X = C8-dG adduct), however there is still limited capacity to predict the likelihood of mutation arising within particular contexts, and hence chemistry-based strategies are needed for probing relationships between nucleic acid sequence and structure with replication errors. In the NarI sequence, certain C8-dG adducts may trigger in the course of DNA synthesis the formation of a slipped mutagenic intermediate (SMI) that contains a two nucleotide (XC) bulge in the template strand that can form upstream of the polymerase active site. This distortion facilitates polymerization but affords a GC dinucleotide deletion product (-2 frameshift mutation). In the current study, incorporating the fluorescent C-linked 4-fluorobiphenyl-dG (FBP-dG) adduct into two 22-mer templates containing CG-dinucleotide repeats ( NarI: 3'-CXCGGC-5' and CG₃: 3'-CXCGCG-5', X = FBP-dG) and performing primer extension reactions using DNA polymerase I, Klenow fragment exo^- (Kf^-) revealed a dramatic sequence-based difference in polymerase bypass efficiency. Primer extension past FBP-dG within the NarI sequence was strongly blocked, whereas Kf^- extended the primer past FBP-dG within a CG₃ template to afford a full-length product and the GC dinucleotide deletion. To model the nucleotide insertion steps in the fully paired (FP) versus the slipped mutagenic (SM) translesion pathways, adducted template:primer duplexes were constructed and characterized by UV thermal denaturation and fluorescence spectroscopy. The emission intensity of the FBP-dG lesion exhibits sensitivity to SMI formation (turn-on) versus a FP duplex (turn-off), permitting insight into adduct base-pairing within the template:primer duplexes. This fluorescence sensitivity provides a rationale for sequence impact on -2 frameshift mutations mediated by the C-linked FBP-dG lesion.

Manipulating DNA G-Quadruplex Structures by Using Guanosine Analogues

The ability to control the folding topology of DNA G-quadruplexes allows for rational design of quadruplex-based scaffolds for potential use in various therapeutic and technological applications. By exploiting the distinct conformational properties of some base- and sugar-modified guanosine surrogates, conformational transitions can be induced through their judicious incorporation at specific sites in the quadruplex core. Changes may involve tetrad polarity inversions with conservation of the global fold or complete refolding to new topologies. Reliable predictions relating to low-energy conformers formed upon specific chemical perturbations of the system and the rational design of modified sequences suffer from our still limited understanding of the subtle interplay of various favorable and unfavorable interactions within a particular quadruplex scaffold. However, aided by an increasing number of systematic substitution experiments and high-resolution structures of modified quadruplex variants, critical interactions, in addition to glycosidic bond angle propensities, are starting to emerge as important contributors to modification-driven quadruplex refolding.

Loop Length Affects Syn-Anti Conformational Rearrangements in Parallel G-Quadruplexes

A G-quadruplex forming sequence from the MYC promoter region was modified with syn-favoring 8-bromo-2'-deoxyguanosine residues. Depending on the number and position of modifications in the intramolecular parallel G-quadruplex, substitutions with the bromoguanosine analogue at the 5'-tetrad induce conformational rearrangements with concerted all-anti to all-syn transitions for all residues of the modified G-quartet. No unfavorable steric interactions of the C8-substituents in the medium grooves are apparent in the high-resolution structure as determined for a tetrasubstituted MYC quadruplex that exclusively forms the all-syn isomer. In contrast, considerable steric clashes with 5'-phosphate oxygen atoms for those analogues that follow a less flexible 1-nucleotide loop in the native all-anti conformation seem to constitute the major driving force for the tetrad inversion and allow for the rational design of appropriately substituted sequences. Correlations found between the population of species subjected to a tetrad flip and melting temperatures indicate that more effective conformational transitions are compromised by lower thermal stabilities of the modified parallel quadruplexes.

Molecular Mechanism, Dynamics, and Energetics of Protein-Mediated Dinucleotide Flipping in a Mismatched DNA: A Computational Study of the RAD4-DNA Complex

DNA damage alters genetic information and adversely affects gene expression pathways leading to various complex genetic disorders and cancers. DNA repair proteins recognize and rectify DNA damage and mismatches with high fidelity. A critical molecular event that occurs during most protein-mediated DNA repair processes is the extrusion of orphaned bases at the damaged site facilitated by specific repairing enzymes. The molecular-level understanding of the mechanism, dynamics, and energetics of base extrusion is necessary to elucidate the molecular basis of protein-mediated DNA damage repair. The present article investigates the molecular mechanism of dinucleotide extrusion in a mismatched DNA (containing a stretch of three contiguous thymidine-thymidine base pairs) facilitated by Radiation sensitive 4 (RAD4), a key DNA repair protein, on an atom-by-atom basis using molecular dynamics (MD) and umbrella-sampling (US) simulations. Using atomistic models of RAD4-free and RAD4-bound mismatched DNA, the free energy profiles associated with extrusion of mismatched partner bases are determined for both systems. The mismatched bases adopted the most stable intrahelical conformation, and their extrusion was unfavorable in RAD4-free mismatched DNA due to the presence of prohibitively high barriers (>12.0 kcal/mol) along the extrusion pathways. Upon binding of RAD4 to the DNA, the global free energy minimum is shifted to the extrahelical state indicating the key role of RAD4-DNA interactions in catalyzing the dinucleotide base extrusion in the DNA-RAD4 complex. The critical residues of RAD4 contributing to the conformational stability of the mismatched bases are identified, and the energetics of insertion of a β-hairpin of RAD4 into the DNA duplex is examined. The conformational energy landscape-based mechanistic insight into RAD4-mediated base extrusion provided here may serve as a useful baseline to understand the molecular basis of xeroderma pigmentosum C (XPC)-mediated DNA damage repair in humans.

Mutagenicity of Ochratoxin A: Role for a Carbon-Linked C8-Deoxyguanosine Adduct?

Ochratoxin A (OTA) is a fungal toxin that is considered to be a potent kidney carcinogen in rodent models. The toxin produces double strand breaks and has a propensity for deletions, single-base substitutions, and insertions. The toxin reacts covalently with DNA to afford a C8-2'-deoxyguanosine carbon-linked adduct (OT-dG) as the major lesion in animal tissues. Incorporation of model C-linked C8-aryl-dG adducts into the G3 site of the NarI sequence demonstrates a tendency to induce base substitutions and deletion mutations in primer extension assays using model polymerases. The degree of misincorporation induced by the C-linked C8-dG adducts correlates with an ability to adopt the promutagenic syn conformation within the NarI duplex as predicted by molecular dynamics (MD) simulations. MD simulations of the OT-dG adduct within the NarI duplex predict an even greater degree of conformational flexibility, suggesting enhanced in vitro mutagenicity compared to the simpler model C-linked C8-dG adducts. Together these findings support the role of OT-dG in promoting OTA-mediated mutagenicity and carcinogenicity in animal studies.

Molecular Modeling of the Major DNA Adduct Formed from Food Mutagen Ochratoxin A in NarI Two-Base Deletion Duplexes: Impact of Sequence Context and Adduct Ionization on Conformational Preference and Mutagenicity

Exposure to ochratoxin A (OTA), a possible human carcinogen, leads to many different DNA mutations. As a first step toward understanding the structural basis of OTA-induced mutagenicity, the present work uses a robust computational approach and a slipped mutagenic intermediate model previously studied for C⁸-dG aromatic amine adducts to analyze the conformational features of postreplication two-base deletion DNA duplexes containing OT-dG, the major OTA lesion at the C⁸ position of guanine. Specifically, a total of 960 ns of molecular dynamics simulations (excluding trial simulations) were carried out on four OT-dG ionization states in three sequence contexts within oligomers containing the NarI recognition sequence, a known hotspot for deletion mutations induced by related adducts formed from known carcinogens. Our results indicate that the structural properties and relative stability of the competing "major groove" and "stacked" conformations of OTA adducted two-base deletion duplexes depend on both the OTA ionization state and the sequence context, mainly due to conformation-dependent deviations in discrete local (hydrogen-bonding and stacking) interactions at the lesion site, as well as DNA bending. When the structural characteristics of the OT-dG adducted two-base deletion duplexes are compared to those associated with previously studied C⁸-dG adducts, a greater understanding of the effects of the nucleobase-carcinogen linkage, and size of the carcinogenic moiety on the conformational preferences of damaged DNA is obtained. Most importantly, our work predicts key structural features for OT-dG-adducted deletion DNA duplexes, which in turn allow us to develop hypotheses regarding OT-dG replication outcomes. Thus, our computational results are valuable for the design and interpretation of future biochemical studies on the potentially carcinogenic OT-dG lesion.

Molecular Insights into the Translesion Synthesis of Benzyl-Guanine from Molecular Dynamics Simulations: Structural Evidence of Mutagenic and Nonmutagenic Replication

DNA can be damaged by many compounds in our environment, and the resulting damaged DNA is commonly replicated by translesion synthesis (TLS) polymerases. Because the mechanism and efficiency of TLS are affected by the type of DNA damage, obtaining information for a variety of DNA adducts is critical. However, there is no structural information for the insertion of a dNTP opposite an O6-dG adduct, which is a particularly harmful class of DNA lesions. We used molecular dynamics (MD) simulations to investigate structural and energetic parameters that dictate preferred dNTP insertion opposite O6-benzyl-guanine (Bz-dG) by DNA polymerase IV, a prototypical TLS polymerase. Specifically, MD simulations were completed on all possible ternary insertion complexes and ternary -1 base deletion complexes with different Bz-dG conformations. Our data suggests that the purines are unlikely to be inserted opposite anti- or syn-Bz-dG, and dTTP is unlikely to be inserted opposite syn-Bz-dG, because of changes in the active site conformation, including critical hydrogen-bonding interactions and/or reaction-ready parameters compared to natural dG replication. In contrast, a preserved active site conformation suggests that dCTP can be inserted opposite either anti- or syn-Bz-dG and dTTP can be inserted opposite anti-Bz-dG. This is the first structural explanation for the experimentally observed preferential insertion of dCTP and misincorporation of dTTP opposite Bz-dG. Furthermore, we provide atomic level insight into why Bz-dG replication does not lead to deletion mutations, which is in contrast with the replication outcomes of other adducts. These findings provide a basis for understanding the replication of related O6-dG adducts.

In silico prediction of the mutagenicity of nitroaromatic compounds using a novel two-QSAR approach

Certain drugs are nitroaromatic compounds, which are potentially toxic. As such, it is of practical importance to assess and predict their mutagenic potency in the process of drug discovery. A classical quantitative structure-activity relationship (QSAR) model was developed using the linear partial least square (PLS) scheme to understand the underline mutagenic mechanism and a non-classical QSAR model was derived using the machine learning-based hierarchical support vector regression (HSVR) to predict the mutagenicity of nitroaromatic compounds based on a series of mutagenicity data (TA98-S9). It was observed that HSVR performed better than PLS as manifested by the predictions of the samples in the training set, test set, and outlier set as well as various statistical validations. A mock test designated to mimic real challenges also confirmed the better performance of HSVR. Furthermore, HSVR exhibited superiority in predictivity, generalization capabilities, consistent performance, and robustness when compared with various published predictive models. PLS, conversely, revealed some mechanistically interpretable relationships between descriptors and mutagenicity. Thus, this two-QSAR approach using the predictive HSVR and interpretable PLS models in a synergistic fashion can be adopted to facilitate drug discovery and development by designing safer drug candidates with nitroaromatic moiety.

C-Linked 8-aryl guanine nucleobase adducts: biological outcomes and utility as fluorescent probes

Aryl radical species derived from enzymatic transformations of aromatic mutagens preferentially react at the 8-site of the guanine (G) nucleobase to afford carbon-linked 8arylG adducts. The resulting lesions possess altered biophysical and genetic coding properties compared to the precursor G nucleoside in B-form DNA. Unlike other adducts, these lesions also possess useful fluorescent properties, since direct attachment of the 8aryl ring extends the purine π-system to afford G mimics with red-shifted excitation maxima and emission that can be sensitive to the microenvironment of the 8arylG base within nucleic acid structures. In B-form DNA, 8arylG adducts are disruptive to duplex formation because they prefer to adopt the syn-conformation about the bond connecting the nucleobase to the deoxyribose backbone, which perturbs Watson-Crick (WC) H-bonding with the opposing cytosine (C). Thus, in a B-form duplex, the emissive properties of 8arylG adducts can be employed as a tool to provide insight into adduct conformation, which can be related to their biological outcomes. However, since Gs preferentially adopt the syn-conformation in left-handed Z-DNA and antiparallel G-quadruplex (GQ) structures, 8arylG lesions can be inserted into syn-G positions without disrupting H-bonding interactions. In fact, 8arylG lesions can serve as ideal fluorescent probes in an antiparallel GQ because their emission is sensitive to GQ folding. This perspective outlines recent developments in the biological implications of 8arylG formation together with their utility as fluorescent G analogs for use in DNA-based diagnostic systems.

Effect of base sequence context on the conformational heterogeneity of aristolactam-I adducted DNA: structural and energetic insights into sequence-dependent repair and mutagenicity

Aristolochic acids (AAs) are nephrotoxic and potentially carcinogenic plant mutagens that form bulky DNA adducts at the exocyclic amino groups of the purines. The present work utilizes classical molecular dynamics simulations and free energy calculations to investigate the role of lesion site sequence context in dictating the conformational outcomes of DNA containing ALI-N6-dA, the most persistent and mutagenic adduct arising from the AAs. Our calculations reveal that the anti base-displaced intercalated conformer is the lowest energy conformer of damaged DNA in all sequence contexts considered (CXC, CXG, GXC and GXG). However, the experimentally-observed greater mutagenicity of the adduct in the CXG sequence context does not correlate with the relative thermodynamic stability of the adduct in different sequences. Instead, AL-N6-dA adducted DNA is least distorted in the CXG sequence context, which points toward a possible differential repair propensity of the lesion in different sequences. Nevertheless, the structural deviations between adducted DNA with different lesion site sequences are small, and therefore other factors (such as interactions between the adducted DNA and lesion-bypass polymerases during replication) are likely more important for dictating the observed sequence-dependent mutagenicity of ALI-N6-dA.

Synthesis of C8-N-Arylamine-Modified 2'-Deoxyguanosine-5'-Triphosphates and Their Effects on Primer Extension by DNA Polymerases

C8-N-arylamine adducts of 2'-deoxyguanosine (2'-dG) play an important role in the induction of the chemical carcinogenesis caused by aromatic amines. C8-N-acetyl-N-arylamine dG adducts that differ in their substitution pattern in the aniline moiety were converted by cycloSal technology into the corresponding C8-N-acetyl-N-arylamine-2'-deoxyguanosine-5'-triphosphates and C8-NH-arylamine-2'-deoxyguanosine-5'-triphosphates. Their conformation preference has been investigated by NOE spectroscopy and DFT calculations. The substrate properties of the C8-dG adducts were studied in primer-extension assays by using Klenow fragment exo(-) of Escherichia coli DNA polymerase I and human DNA polymerase β. It was shown that the incorporation was independent of the substitution pattern in the aryl moiety and the N-acetyl group. Although the triphosphates were poor substrates for the human polymerases, they were incorporated twice before the termination of the elongation process occurred; this might demonstrate the importance of C8-N-arylamine-2'-deoxyguanosine-5'-triphosphates in chemical carcinogenesis.

Enhancing Bulge Stabilization through Linear Extension of C8-Aryl-Guanine Adducts to Promote Polymerase Blockage or Strand Realignment to Produce a C:C Mismatch

Aryl radicals can react at the C8-site of 2'-deoxyguanosine (dG) to produce DNA adducts with a C8-C linkage (denoted C-linked). Such adducts are structurally distinct from those possessing a flexible amine (N-linked) or ether (O-linked) linkage, which separates the C8-aryl moiety from the guanine nucleobase. In the current study, two model C-linked C8-dG adducts, namely, C8-benzo[b]thienyl-dG ([BTh]G) and C8-(pyren-1-yl)-dG ([Py]G), were incorporated into the NarI (12mer, NarI(12) and 22mer, NarI(22)) hotspot sequence for frameshift mutations in bacteria. For the first time, C-linked C8-dG adducts are shown to stabilize the -2 deletion duplex within the NarI sequence. Primer-elongation assays employing Sulfolobus solfataricus P2 DNA polymerase IV (Dpo4) demonstrates the influence of C8-aryl ring size and shape in promoting Dpo4 blockage or strand realignment to produce a C:C mismatch downstream of the adduct site. Molecular dynamics simulations of the -2 deletion duplex suggest that both anti and syn adduct structures are energetically accessible. These findings provide a rationale for describing the biochemical outcome induced by C-linked C8-dG adducts when processed by Dpo4.

Adenine versus guanine DNA adducts of aristolochic acids: role of the carcinogen-purine linkage in the differential global genomic repair propensity

Computational modeling is employed to provide a plausible structural explanation for the experimentally-observed differential global genome repair (GGR) propensity of the ALII-N(2)-dG and ALII-N(6)-dA DNA adducts of aristolochic acid II. Our modeling studies suggest that an intrinsic twist at the carcinogen-purine linkage of ALII-N(2)-dG induces lesion site structural perturbations and conformational heterogeneity of damaged DNA. These structural characteristics correlate with the relative repair propensities of AA-adducts, where GGR recognition occurs for ALII-N(2)-dG, but is evaded for intrinsically planar ALII-N(6)-dA that minimally distorts DNA and restricts the conformational flexibility of the damaged duplex. The present analysis on the ALII adduct model systems will inspire future experimental studies on these adducts, and thereby may extend the list of structural factors that directly correlate with the propensity for GGR recognition.

Chlorine functionalization of a model phenolic C8-guanine adduct increases conformational rigidity and blocks extension by a Y-family DNA polymerase

Certain phenoxyl radicals can attach covalently to the C8-site of 2'-deoxyguanosine (dG) to afford oxygen-linked C8-dG adducts. Such O-linked adducts can be chemically synthesized through a nucleophilic displacement reaction between a phenolate and a suitably protected 8-Br-dG derivative. This permits the generation of model O-linked C8-dG adducts on scales suitable for insertion into oligonucleotide substrates using solid-phase DNA synthesis. Variation of the C8-aryl moiety provides an opportunity to derive structure-activity relationships on adduct conformation in duplex DNA and replication bypass by DNA polymerases. In the current study, the influence of chlorine C8-dG functionalization on in vitro DNA replication by Klenow fragment exo(-) (Kf(-)) and the Y-family polymerase (Sulfolobus solfataricus P2 DNA polymerase IV (Dpo4)) has been determined. Model O-linked C8-dG adducts derived from the pentachlorophenoxyl radical ([PCP]G) and 2,4,6-trichlorophenoxyl radical ([TCP]G) were inserted into the reiterated G3-position of the NarI sequence (12-mer, NarI(12); and 22-mer, NarI(22)), which is a known hotspot for frameshift mutations mediated by N-linked polycyclic C8-dG adducts in bacterial mutagenesis. Within the NarI(12) duplex, the unsubstituted C8-phenoxy-dG ([PhO]G) adduct adopts a minimally perturbed B-form helix. Chlorination of [PhO]G to afford [PCP]G does not significantly change the adduct conformation within the NarI(12) duplex, as predicted by molecular dynamics simulations. However, when using NarI(22) for DNA synthesis in vitro, the chlorinated [PCP]G and [TCP]G lesions significantly block DNA replication by Kf(-) and Dpo4, whereas [PhO]G is readily bypassed. These findings highlight the impact that chlorine substituents impart to bulky C8-dG lesions.

Positional impact of fluorescently modified G-tetrads within polymorphic human telomeric G-quadruplex structures

Emissive C8-aryl-2'-deoxyguanosines placed within G-tetrads of G-quadruplex structures are useful probes for distinguishing G-quadruplexes from duplex structures using fluorescence spectroscopy. Here, we report the positional impact of C8-furyl-dG ((Fur)dG) on G-quadruplex folding in the human telomere 22-mer oligonucleotide (HTelo22, (d[AG3(T2AG3)3])). The (Fur)dG probe was inserted into four different positions within the three unique G-tetrads of HTelo22, and G-quadruplex folding was monitored by UV-vis thermal denaturation, circular dichroism, and fluorescence spectroscopy. Our studies demonstrate the impacts of C8-aryl-dG adduct formation on G-quadruplex polymorphism in K(+) solution and in the presence of the additives and cosolutes, CH3CN, polyethylene glycol (PEG-600), and N-methyl mesoporphyrin IX (NMM). Our experiments predict that C8-aryl-dG derivatives can serve as useful tools for various in vitro studies aimed at understanding both the implications of DNA adduct formation within G-quadruplex structures and the unique implications imposed by various folding topologies on biological function/recognition.

Dual fluorescent deoxyguanosine mimics for FRET detection of G-quadruplex folding

Replacement of deoxyguanosine (dG) nucleobases within G-tetrads of G-quadruplex folding oligonucleotides with donor (D)/acceptor (A) fluorescent 8aryldG residues provides diagnostic FRET signalling for G-quadruplex detection.