Is FapyG Mutagenic?: Evidence from the DFT Study

Hydroxyl radical-induced C1'-H abstraction reaction of different artificial nucleotides

CONTEXT

Recently, a few antiviral drugs viz Molnupiravir (EIDD-1931), Favipiravir, Ribavirin, Sofosbuvir, Galidesivir, and Remdesivir are shown to be beneficial against COVID-19 disease. These drugs bind to the viral RNA single strand to inhibit the virus genome replication. Similarly, recently, some artificial nucleotides, such as P, J, B, X, Z, V, S, and K were proposed to behave as potent antiviral candidates. However, their activity in the presence of the most reactive hydroxyl (OH) radical is not yet known. Here, the possibility of RNA strand break due to the OH radical-induced C1'-hydrogen (H) abstraction reaction of the above molecules (except Remdesivir) is studied in detail by considering their nucleotide conformation. The results are compared with those of the natural RNA nucleotides (G, C, A, and U). Due to low Gibbs barrier-free energy and high exothermicity, all these nucleotides (except Remdesivir) are prone to OH radical-induced C1'-H abstraction reaction. As Remdesivir contains a C1'-CN bond, the OH radical substitution reactions at the CN and C1' sites would likely liberate the catalytically important CN group, thereby downgrading its activity.

METHOD

Initially, the B3LYP-D3 dispersion-corrected density functional theory method and 6-31 + G* basis set were used to optimize all reactant, transition state, and product complexes in the implicit aqueous medium. Subsequently, the structures of these complexes were further optimized by using the ωB97X-D dispersion-corrected density functional theory method and cc-PVTZ basis set in the aqueous medium. The IEFPCM method was used to model the aqueous medium.

Methodologies for the detection and sequencing of the epigenetic-like oxidative DNA modification, 8-oxo-7,8-dihydroguanine

The human genome is constantly threatened by endogenous and environmental DNA damaging agents that can induce a variety of chemically modified DNA lesions including 8-oxo-7,8-dihydroguanine (OG). Increasing evidence has indicated that OG is not only a biomarker for oxidative DNA damage but also a novel epigenetic-like modification involved in regulation of gene expression in mammalian cells. Here we summarize the recent progress in OG research focusing on the following points: (i) the mechanism of OG production in organisms and its biological consequences in cells, (ii) the accurate identification of OG in low-abundance genomes and complex biological backgrounds, (iii) the development of OG sequencing methods. These studies will be helpful for further understanding of the molecular mechanisms of OG-induced mutagenesis and its potential roles in human development and diseases such as cancer.

FapydG in the Shadow of OXOdG—A Theoretical Study of Clustered DNA Lesions

Genetic information, irrespective of cell type (normal or cancerous), is exposed to a range of harmful factors, which can lead to more than 80 different types of DNA damage. Of these, oxoG and FapyG have been identified as the most abundant in normoxic and hypoxic conditions, respectively. This article considers d[AFapyGAOXOGA]*[TCTCT] (oligo-FapyG) with clustered DNA lesions (CDLs) containing both the above types of damage at the M06-2x/6-31++G** level of theory in the condensed phase. Furthermore, the electronic properties of oligo-FapyG were analysed in both equilibrated and non-equilibrated solvation–solute interaction modes. The vertical/adiabatic ionization potential (VIP, AIP) and electron affinity (VEA, AEA) of the investigated ds-oligo were found as follows in [eV]: 5.87/5.39 and −1.41/−2.09, respectively. The optimization of the four ds-DNA spatial geometries revealed that the transFapydG was energetically privileged. Additionally, CDLs were found to have little influence on the ds-oligo structure. Furthermore, for the FapyGC base-pair isolated from the discussed ds-oligo, the ionization potential and electron affinity values were higher than those assigned to OXOGC. Finally, a comparison of the influence of FapyGC and OXOGC on charge transfer revealed that, in contrast to the OXOGC base-pair, which, as expected, acted as a radical cation/anion sink in the oligo-FapyG structure, FapyGC did not significantly affect charge transfer (electron–hole and excess–electron). The results presented below indicate that 7,8-dihydro-8-oxo-2′-deoxyguanosine plays a significant role in charge transfer through ds-DNA containing CDL and indirectly has an influence on the DNA lesion recognition and repair process. In contrast, the electronic properties obtained for 2,6-diamino-4-hydroxy-5-foramido-2′deoxypyrimidine were found to be too weak to compete with OXOG to influence charge transfer through the discussed ds-DNA containing CDL. Because increases in multi-damage site formation are observed during radio- or chemotherapy, understanding their role in the above processes can be crucial for the efficiency and safety of medical cancer treatment.

Impact of Helicobacter pylori Infection and Its Major Virulence Factor CagA on DNA Damage Repair

Helicobacter pylori infection induces a plethora of DNA damages. Gastric epithelial cells, in order to maintain genomic integrity, require an integrous DNA damage repair (DDR) machinery, which, however, is reported to be modulated by the infection. CagA is a major H. pylori virulence factor, associated with increased risk for gastric carcinogenesis. Its pathogenic activity is partly regulated by phosphorylation on EPIYA motifs. Our aim was to identify effects of H. pylori infection and CagA on DDR, investigating the transcriptome of AGS cells, infected with wild-type, ΔCagA and EPIYA-phosphorylation-defective strains. Upon RNA-Seq-based transcriptomic analysis, we observed that a notable number of DDR genes were found deregulated during the infection, potentially resulting to base excision repair and mismatch repair compromise and an intricate deregulation of nucleotide excision repair, homologous recombination and non-homologous end-joining. Transcriptome observations were further investigated on the protein expression level, utilizing infections of AGS and GES-1 cells. We observed that CagA contributed to the downregulation of Nth Like DNA Glycosylase 1 (NTHL1), MutY DNA Glycosylase (MUTYH), Flap Structure-Specific Endonuclease 1 (FEN1), RAD51 Recombinase, DNA Polymerase Delta Catalytic Subunit (POLD1), and DNA Ligase 1 (LIG1) and, contrary to transcriptome results, Apurinic/Apyrimidinic Endodeoxyribonuclease 1 (APE1) upregulation. Our study accentuates the role of CagA as a significant contributor of H. pylori infection-mediated DDR modulation, potentially disrupting the balance between DNA damage and repair, thus favoring genomic instability and carcinogenesis.

Role of different tautomers in the base-pairing abilities of some of the vital antiviral drugs used against COVID-19

Base-pair mutations induced by different tautomers of anti-viral drugs are the main reasons for their anti-viral activities.

Recognition of DNA adducts by edited and unedited forms of DNA glycosylase NEIL1

Pre-mRNA encoding human NEIL1 undergoes editing by adenosine deaminase ADAR1 that converts a single adenosine to inosine, and this conversion results in an amino acid change of lysine 242 to arginine. Previous investigations of the catalytic efficiencies of the two forms of the enzyme revealed differential release of thymine glycol (ThyGly) from synthetic oligodeoxynucleotides, with the unedited form, NEIL1 K242 being ≈30-fold more efficient than the edited NEIL1 K242R. In contrast, when these enzymes were reacted with oligodeoxynucleotides containing guanidinohydantoin or spiroiminohydantoin, the edited K242R form was ≈3-fold more efficient than the unedited NEIL1. However, no prior studies have investigated the efficiencies of these two forms of NEIL1 on either high-molecular weight DNA containing multiple oxidatively-induced base damages, or oligodeoxynucleotides containing a bulky alkylated formamidopyrimidine. To understand the extent of changes in substrate recognition, γ-irradiated calf thymus DNA was treated with either edited or unedited NEIL1 and the released DNA base lesions analyzed by gas chromatography-tandem mass spectrometry. Of all the measured DNA lesions, imidazole ring-opened 4,6-diamino-5-formamidopyrimidine (FapyAde) and 2,6-diamino-4-hydroxy-5-formamidopyrimidine (FapyGua) were preferentially released by both NEIL1 enzymes with K242R being ≈1.3 and 1.2-fold more efficient than K242 on excision of FapyAde and FapyGua, respectively. Consistent with the prior literature, large differences (≈7.5 to 12-fold) were measured in the excision of ThyGly from genomic DNA by the unedited versus edited NEIL1. In contrast, the edited NEIL1 was more efficient (≈3 to 5-fold) on release of 5-hydroxycytosine. Excision kinetics on DNA containing a site-specific aflatoxin B₁-FapyGua adduct revealed an ≈1.4-fold higher rate by the unedited NEIL1. Molecular modeling provides insight into these differential substrate specificities. The results of this study and in particular, the comparison of substrate specificities of unedited and edited NEIL1 using biologically and clinically important base lesions, are critical for defining its role in preservation of genomic integrity.

Accurate Base Pair Energies of Artificially Expanded Genetic Information Systems (AEGIS): Clues for Their Mutagenic Characteristics

Recently, several artificial nucleobases, such as B, S, J, V, X, K, P, and Z, have been proposed to help in the expansion of the genetic information system and diagnosis of diseases. Among these bases, P and Z were identified to form stable DNA and to participate in the replication. However, the stabilities of P:Z and other artificial base pairs are not fully understood. The abilities of these unnatural nucleobases in mispairing with themselves and with natural bases are also not known. Here, the ωB97X-D dispersion-corrected density functional theoretical and complete basis set (CBS-QB3) methods are used to obtain accurate structural and energetic data related to base pair interactions involving these unnatural nucleobases. The roles of protonation and deprotonation of certain artificial bases in inducing mutations are also studied. It is found that each artificial purine has a complementary artificial pyrimidine, the base pair interactions between which are similar to those of the natural Watson-Crick base pairs. Hence, these base pairs will function naturally and would not impart mutagenicity. Among these base pairs, the J:V complex is found to be the most stable and promising artificial base pair. Remarkably, the noncomplementary artificial nucleobases are found to form stable mispairs, which may generate mutagenic products in DNA. Similarly, the misinsertions of natural bases opposite artificial bases are also found to be mutagenic. The mechanisms of these mutations are explained in detail. These results are in agreement with earlier biochemical studies. It is thus expected that this study would aid in the advancement of the synthetic biology to design more robust artificial nucleotides.

Analogues of P and Z as Efficient Artificially Expanded Genetic Information System

To artificially expand the genetic information system and to realize artificial life, it is necessary to discover new functional DNA bases that can form stable duplex DNA and participate in error-free replication. It is recently proposed that the 2-amino-imidazo[1,2- a]-1,3,5-triazin-4(8 H)one (P) and 6-amino-5-nitro-2(1 H)-pyridone (Z) would form a base pair complex, which is more stable than that of the normal G-C base pair and would produce an unperturbed duplex DNA. Here, by using quantum chemical calculations in aqueous medium, it is shown that the P and Z molecules can be modified with the help of electron-withdrawing and -donating substituents mainly found in B-DNA to generate new bases that can produce even more stable base pairs. Among the various bases studied, P3, P4, Z3, and Z5 are found to produce base pairs, which are about 2-15 kcal/mol more stable than the P-Z base pair. It is further shown that these base pairs can be stacked onto the G-C and A-T base pairs to produce stable dimers. The consecutive stacking of these base pairs is found to yield even more stable dimers. The influence of charge penetration effects and backbone atoms in stabilizing these dimers are also discussed. It is thus proposed that the P3, P4, Z3, and Z5 would form promiscuous artificial genetic information system and can be used for different biological applications. However, the evaluations of the dynamical effects of these bases in DNA-containing several nucleotides and the efficacy of DNA polymerases to replicate these bases would provide more insights.

Error-prone replication bypass of the imidazole ring-opened formamidopyrimidine deoxyguanosine adduct

Addition of hydroxyl radicals to the C8 position of 2'-deoxyguanosine generates an 8-hydroxyguanyl radical that can be converted into either 8-oxo-7,8-dihydro-2'-deoxyguanosine or N-(2-deoxy-d-pentofuranosyl)-N-(2,6-diamino-4-hydroxy-5-formamidopyrimidine) (Fapy-dG). The Fapy-dG adduct can adopt different conformations and in particular, can exist in an unnatural α anomeric configuration in addition to canonical β configuration. Previous studies reported that in 5'-TGN-3' sequences, Fapy-dG predominantly induced G → T transversions in both mammalian cells and Escherichia coli, suggesting that mutations could be formed either via insertion of a dA opposite the 5' dT due to primer/template misalignment or as result of direct miscoding. To address this question, single-stranded vectors containing a site-specific Fapy-dG adduct were generated to vary the identity of the 5' nucleotide. Following vector replication in primate cells (COS7), complex mutation spectra were observed that included ∼3-5% G → T transversions and ∼14-21% G → A transitions. There was no correlation apparent between the identity of the 5' nucleotide and spectra of mutations. When conditions for vector preparation were modified to favor the β anomer, frequencies of both G → T and G → A substitutions were significantly reduced. Mutation frequencies in wild-type E. coli and a mutant deficient in damage-inducible DNA polymerases were significantly lower than detected in COS7 and spectra were dominated by deletions. Thus, mutagenic bypass of Fapy-dG can proceed via mechanisms that are different from the previously proposed primer/template misalignment or direct misinsertions of dA or dT opposite to the β anomer of Fapy-dG. Environ. Mol. Mutagen. 58:182-189, 2017. © 2017 Wiley Periodicals, Inc.

Conformational stabilities of iminoallantoin and its base pairs in DNA: implications for mutagenicity

Under acidic conditions, insertion of G opposite Ia may lead to G to C mutations in DNA.

Normal and reverse base pairing of Iz and Oz lesions in DNA: structural implications for mutagenesis

During replication, incorporation of G opposite Oz lesion is mainly responsible for G to C mutations in DNA.

The R- and S-diastereoisomeric effects on the guanidinohydantoin-induced mutations in DNA

Although, Gh (Gh1 or Gh2) in DNA would induce mainly G to C mutations, other mutations cannot be ignored.