Modeling Competitive Reaction Mechanisms of Peroxynitrite Oxidation of Guanine

Tracking Isomerizations of High-Energy Adenine Cation Radicals by UV-Vis Action Spectroscopy and Cyclic Ion Mobility Mass Spectrometry

We report experimental and computational studies of protonated adenine C-8 σ-radicals that are presumed yet elusive reactive intermediates of oxidative damage to nucleic acids. The radicals were generated in the gas phase by the collision-induced dissociation of C-8-Br and C-8-I bonds in protonated 8-bromo- and 8-iodoadenine as well as by 8-bromo- and 8-iodo-9-methyladenine. Protonation by electrospray of 8-bromo- and 8-iodoadenine was shown by cyclic-ion mobility mass spectrometry (c-IMS) to form the N-1-H, N-9-H and N-3-H, N-7-H protomers in 85:15 and 81:19 ratios, respectively, in accordance with the equilibrium populations of these protomers in water-solvated ions that were calculated by density functional theory (DFT). Protonation of 8-halogenated 9-methyladenines yielded single N-1-H protomers, which was consistent with their thermodynamic stability. The radicals produced from the 8-bromo and 8-iodo adenine cations were characterized by UV-vis photodissociation action spectroscopy (UVPD) and c-IMS. UVPD revealed the formation of C-8 σ-radicals along with N-3-H, N-7-H-adenine π-radicals that arose as secondary products by hydrogen atom migrations. The isomers were identified by matching their action spectra against the calculated vibronic absorption spectra. Deuterium isotope effects were found to slow the isomerization and increase the population of C-8 σ-radicals. The adenine cation radicals were separated by c-IMS and identified by their collision cross sections, which were measured relative to the canonical N-9-H adenine cation radical that was cogenerated in situ as an internal standard. Ab initio CCSD(T)/CBS calculations of isomer energies showed that the adenine C-8 σ-radicals were local energy minima with relative energies at 76-79 kJ mol^-1 above that of the canonical adenine cation radical. Rice-Ramsperger-Kassel-Marcus calculations of unimolecular rate constants for hydrogen and deuterium migrations resulting in exergonic isomerizations showed kinetic shifts of 10-17 kJ mol^-1, stabilizing the C-8 σ-radicals. C-8 σ-radicals derived from N-1-protonated 9-methyladenine were also thermodynamically unstable and readily isomerized upon formation.

Reactive Molecular Dynamics Simulation of the Structural Damages of the B-DNA Induced by the Oxidation/Nitration of Guanine

Reactive molecular dynamics simulations (RMD) have been carried out to investigate structural alterations of the dodecamer double-strand B-DNA due to the oxidation/nitration modifications introduced to its guanine bases, including 8-oxoguanine, 8-nitroguanine, and 5-guanidino-4-nitroimidazole, considering two distribution patterns. These modifications may arise in the case of cancer treatment using oxidative/nitrosative reactive nitrogen species as anticancer agents. Results show that these mutations affect structural characteristics of the B-DNA dodecamer in the order 8-nitroguanine > 5-guanidino-4-nitroimidazole ≫ 8-oxoguanine. For instance, the base-pair per turn for these modified B-DNA are changed respectively to 9.79, 10.88 and 10.58 from 10.51 in the native defect-free B-DNA, which is compatible with the experimental value of 10.10. In addition, these mutations allow more water molecules to diffuse into the dodecamer structure and consequently increase the possibility of the penetration of reactive and nonreactive species toward constituting nucleic base-pairs. The largest variation of the B-DNA structure is observed for the mutated B-DNA with 8-nitroguanine modifications applied to its separated CG base-pairs along the dodecamer chain. The structural changes introduced by these nitro-/oxo-modified guanine bases can be considered as a critical step in the damage of the DNA structure and alterations of its function.

Fast and biphasic 8-nitroguanine production from guanine and peroxynitrite

Guanine (Gua), among purines, is a preferred oxidation/nitration target because of its low one-electron redox potential. The reactive oxygen/nitrogen species peroxynitrite (ONOO^-), produced in vivo by the reaction between nitric oxide (^•NO) and superoxide radical (O₂^•‒), is responsible for several oxidative modifications in biomolecules, including nitration, nitrosation, oxidation, and peroxidation. In particular, the nitration of Gua, although detected, as well as its reaction kinetics have been seldom investigated. Thus, we studied the concentration- and temperature-dependent formation of 8-nitroguanine (8-NitroGua) in phosphate buffer (pH 7.40) using stopped-flow spectrophotometry. Traces showed a biexponential behavior, with best-fit rate constants: k_fast = 4.4 s^-1 and k_slow = 0.41 s^-1 (30 °C, 400 μM both Gua and ONOO^-). k_fast increased linearly with the concentration of both reactants whereas k_slow was concentration-independent. Linear regression analysis of k_fast as a function of Gua and ONOO^- concentration yielded values of 2.5-6.3 × 10³ M^-1s^-1 and 1.5-3.5 s^-1 for the second-order (slope) and first-order (ordinate) rate constants, respectively (30 °C). Since ONOO^- is a short-lived species, its decay kinetics was also taken into account for this analysis. The 8-NitroGua product was stable for at least 4 h, so no spontaneous denitration was observed. Stopped-flow assays using antioxidants and free-radical scavengers suggested a mixed direct/indirect reaction mechanism for 8-NitroGua formation. Gua nitration by ONOO^- was also observed in the presence of physiologically relevant CO₂ concentrations. The reaction product identity, its yield (∼4.2%, with 400 μM ONOO^- and 200 μM Gua), and the reaction mechanism were unequivocally determined by HPLC-MS/MS experiments. In conclusion, 8-NitroGua production at physiologic pH reached significant levels in a few hundred milliseconds, suggesting that the process might be kinetically relevant in vivo and can likely cause permanent nitrative damage to DNA bases.

Reaction mechanism and dynamics for C8-hydroxylation of 9-methylguanine radical cation by water molecules

In contrast to their spontaneous deprotonation in aqueous solution, reactions of guanine and guanosine radical cations with water in the gas phase are exclusively initiated by hydration of the radical cations as reported in recent work (Y. Sun et al., Phys. Chem. Chem. Phys., 2018, 20, 27510). As gas-phase hydration reactions closely mimic the actual scenario for guanine radical cations in double-stranded DNA, exploration of subsequent reactions within their water complexes can provide an insight into the resulting oxidative damage to nucleosides. Herein guided-ion beam mass spectrometry experiment and direct dynamics trajectory simulations were carried out to examine prototype complexes of the 9-methylguanine radical cation with one and two water ligands (i.e., 9MG˙⁺·(H₂O)_1-2) in the gas phase, wherein the complexes were activated by collisional activation in the experiment and by thermal excitation at high temperatures in the simulations. Guided by mass spectroscopic measurements, trajectory results and reaction potential energy surface, three reaction pathways were identified. The first two reaction pathways start with H-atom abstraction from water by the O6 and N7 atoms in 9MG˙⁺ and are referred to as HA_O6 and HA_N7, respectively. The primary products of HA_O6 and HA_N7 reactions, including [9MG + H_O6]⁺/[9MG + H_N7]⁺ and ˙OH, react further to either form [8OH-9MG + H_O6]˙⁺ and [8OH-9MG + H_N7]˙⁺via C8-hydroxylation or form radical cations of 6-enol-guanine (6-enol-G˙⁺) and 7H-guanine (7HG˙⁺) via S_N2-type methanol elimination. The third reaction pathway corresponds to the formation of 8OH-9MG⁺ by H elimination from the complex, referred to as HE. Among these product channels, [8OH-9MG + H_N7]˙⁺ has the most favorable formation probability, especially in the presence of additional water molecules. This product may serve as a preceding structure to the 8-oxo-7,8-dihydroguanine lesion in DNA and has implications for health effects of radiation exposure and radiation therapy.

Deprotonation of Guanine Radical Cation G•+ Mediated by the Protonated Water Cluster

Proton transfer is regarded as a fundamental process in chemical reactions of DNA molecules and continues to be an active research theme due to the connection with charge transport and oxidation damage of DNA. For the guanine radical cation (G^•+) derived from one-electron oxidation, experiments suggest a facile proton transfer within the G^•+:C base pair, and a rapid deprotonation from N1 in free base or single-strand DNA. To address the deprotonation mechanism, we perform a thorough investigation on deprotonation of G^•+ in free G base by combining density functional theory (DFT) and laser flash photolysis spectroscopy. Experimentally, kinetics of deprotonation is monitored at temperatures varying from 280 to 298 K, from which the activation energy of 15.1 ± 1.5 kJ/mol is determined for the first time. Theoretically, four solvation models incorporating explicit waters and the polarized continuum model (PCM), i.e., 3H₂O-PCM, 4H₂O-PCM, 5H₂O-PCM, and 7H₂O-PCM models are used to calculate deprotonation potential energy profile, and the barriers of 5.5, 13.4, 14.4, and 13.7 kJ/mol are obtained, respectively. It is shown that at least four explicit waters are required for properly simulating the deprotonation reaction, where the participation of protonated water cluster plays key roles in facilitating the proton release from G^•+.

Computational Investigation into the Oxidation of Guanine to Form Imidazolone (Iz) and Related Degradation Products

Imidazolone (Iz) is one of the many products resulting from oxidative damage to DNA. Three pathways for the formation of Iz and related degradation products have been studied by density functional theory using the ωB97XD functional with the 6-31+G(d,p) basis set and SMD implicit water solvation plus a small number of explicit water molecules positioned to help stabilize charged species and facilitate reaction steps. The first pathway starts with guanine radical and the addition of superoxide at C5. Endoperoxide formation was calculated to have slightly lower barriers than diol formation. The next steps are pyrimidine ring opening and decarboxylation. Ring migration then proceeds via an acyclic intermediate rather than a bicyclic intermediate and is followed by formamide loss to yield Iz. The second pathway starts with 8oxoG and proceeds via C5 superoxide addition and diol formation to a relatively stable intermediate, oxidized guanidinohydantoin (Gh_ox). The barriers for hydroxide ion addition to Gh_ox are much lower than for water addition and should yield more Iz and parabanic acid at higher pH. The third pathway starts with 8-hydroxy guanine radical formed by hydroxyl radical addition to C8 of guanine or water addition to C8 of guanine radical. Superoxide addition at C5 is followed by diol formation, ring opening and decarboxylation similar to pathways 1 and 2, subsequently leading to Iz formation. The calculated pathways are in good agreement with experimental observations.

Computational Study of the pH-Dependent Competition between Carbonate and Thymine Addition to the Guanine Radical

When oligonucleotides are oxidized by carbonate radical, thymine and carbonate can add to guanine radical, yielding either a guanine-thymine cross-link product (G∧T) or 8-oxo-7,8-dehydroguanine (8oxoG) and its further oxidation products such as spiroiminodihydantoin (Sp) and guanidinohydantoin (Gh). The ratio of thymine addition to carbonate addition depends strongly on the pH. Details of the mechanism have been explored by density functional calculations using the ωB97XD/6-31+G(d,p) level of theory with the SMD implicit solvation method, augmented with a few explicit waters. Free energies of intermediates and transition states in aqueous solution have been calculated along the pathways for addition of thymine, CO₃^2-/HCO₃^- and carbonate radical to guanine radical. The pH dependence was examined by using appropriate explicit proton donors/acceptors as computational models for buffers at pH 2.5, 7, and 10. Deprotonation of thymine is required for nucleophilic addition at C8 of guanine radical, and thus is favored at higher pH. The barrier for carbonate radical addition is lower than for bicarbonate or carbonate dianion addition; however, for low concentrations of carbonate radical, the reaction may proceed by addition of bicarbonate/carbonate dianion to guanine radical. Thymine and bicarbonate/carbonate dianion addition are followed by oxidation by O₂, loss of a proton from C8 and decarboxylation of the carbonate adduct. At pH 2.5, guanine radical cation can be formed by oxidization with sulfate radical. Water addition to guanine radical cation is the preferred path for forming 8oxoG at pH 2.5.

Mapping three guanine oxidation products along DNA following exposure to three types of reactive oxygen species

Reactive oxygen and nitrogen species generated during respiration, inflammation, and immune response can damage cellular DNA, contributing to aging, cancer, and neurodegeneration. The ability of oxidized DNA bases to interfere with DNA replication and transcription is strongly influenced by their chemical structures and locations within the genome. In the present work, we examined the influence of local DNA sequence context, DNA secondary structure, and oxidant identity on the efficiency and the chemistry of guanine oxidation in the context of the Kras protooncogene. A novel isotope labeling strategy developed in our laboratory was used to accurately map the formation of 2,2-diamino-4-[(2-deoxy-β-D-erythropentofuranosyl)amino]- 5(2 H)-oxazolone (Z), 8-oxo-7,8-dihydro-2'-deoxyguanosine (OG), and 8-nitroguanine (8-NO₂-G) lesions along DNA duplexes following photooxidation in the presence of riboflavin, treatment with nitrosoperoxycarbonate, and oxidation in the presence of hydroxyl radicals. Riboflavin-mediated photooxidation preferentially induced OG lesions at 5' guanines within GG repeats, while treatment with nitrosoperoxycarbonate targeted 3'-guanines within GG and AG dinucleotides. Little sequence selectivity was observed following hydroxyl radical-mediated oxidation. However, Z and 8-NO₂-G adducts were overproduced at duplex ends, irrespective of oxidant identity. Overall, our results indicate that the patterns of Z, OG, and 8-NO₂-G adduct formation in the genome are distinct and are influenced by oxidant identity and the secondary structure of DNA.

Reaction of hypochlorous acid with imidazole: Formation of 2-chloro- and 2-oxoimidazoles

Reaction of hypochlorous acid (HOCl) with imidazole (Im) taken as a model for the 5-membered ring of guanine, leading to the products 2-chloro- and 2-oxo-imidazoles was investigated at the B3LYP/6-31+G* and B3LYP/AUG-cc-pVDZ levels of density functional theory. For all cases, single point energy calculations were performed at the MP2/AUG-cc-pVDZ level of theory using the geometries optimised at the B3LYP/AUG-cc-pVDZ level. Intrinsic reaction coordinate calculations were performed to ensure genuineness of all the calculated transition states. Effect of aqueous media was investigated by solvating all the species involved in the reactions using the polarizable continuum model. It is found that 2-chloroimidazole (2-ClIm) can be formed following three different reaction schemes while 2-oxoimidazole (2-oxoIm) can be formed following two different reaction schemes. The calculated barrier energies show that formation of 2-oxoIm would be less favored than that of 2-ClIm, which explains the experimental observations on relative yields of 8-chlorodeoxyguanosine and 8-oxodeoxyguanosine.