Mechanisms of Formation of 8-Oxoguanine Due To Reactions of One and Two OH• Radicals and the H2O2 Molecule with Guanine: A Quantum Computational Study

Mutagenicity associated with O6-methylguanine-DNA damage and mechanism of nucleotide flipping by AGT during repair

Methylated guanine damage at O6 position (i.e. O6MG) is dangerous due to its mutagenic and carcinogenic character that often gives rise to G:C-A:T mutation. However, the reason for this mutagenicity is not known precisely and has been a matter of controversy. Further, although it is known that O6-alkylguanine-DNA alkyltransferase (AGT) repairs O6MG paired with cytosine in DNA, the complete mechanism of target recognition and repair is not known completely. All these aspects of DNA damage and repair have been addressed here by employing high level density functional theory in gas phase and aqueous medium. It is found that the actual cause of O6MG mediated mutation may arise due to the fact that DNA polymerases incorporate thymine opposite to O6MG, misreading the resulting O6MG:T complex as an A:T base pair due to their analogous binding energies and structural alignments. It is further revealed that AGT mediated nucleotide flipping occurs in two successive steps. The intercalation of the finger residue Arg128 into the DNA double helix and its interaction with the O6MG:C base pair followed by rotation of the O6MG nucleotide are found to be crucial for the damage recognition and nucleotide flipping.

Quantum theoretical study of cleavage of the glycosidic bond of 2'-deoxyadenosine: base excision-repair mechanism of DNA by MutY

The enzyme adenine DNA glycosylase, also called MutY, is known to catalyze base excision repair by removal of adenine from the abnormal 2'-deoxyadenosine:8-oxo-2'-deoxyguanosine pair in DNA. The active site of the enzyme was considered to consist of a glutamic acid residue along with two water molecules. The relevant reaction mechanism involving different barrier energies was studied theoretically. Molecular geometries of the various molecules and complexes involved in the reaction, e.g., the reactant, intermediate, and product complexes as well as transition states, were optimized employing density functional theory at the B3LYP/6-31G(d,p) level in the gas phase. It was followed by single-point energy calculations at the B3LYP/AUG-cc-pVDZ, BHandHLYP/AUG-cc-pVDZ, and MP2/AUG-cc-pVDZ levels in the gas phase. Single-point energy calculations were also carried out at the B3LYP/AUG-cc-pVDZ and BHandHLYP/AUG-cc-pVDZ levels in aqueous media as well as in the solvents chlorobenzene and dichloroethane. For the solvation calculations, the integral equation formalism of the polarizable continuum model (IEF-PCM) was employed. It is found that glutamic acid along with two water molecules would effectively cleave the glycosidic bond of adenosine by a new two-step reaction mechanism proposed here which is different from the three-step mechanism proposed by other authors earlier regarding the working mechanism of MutY.

Mechanism of scavenging action of N-acetylcysteine for the OH radical: a quantum computational study

N-Acetylcysteine, a precursor of glutathione, is an effective antioxidant present in biological systems. The mechanism of scavenging action of N-acetylcysteine for the OH radical was studied theoretically. For this purpose, reactions of the OH radical at the different sites of N-acetylcysteine were investigated. All the relevant extrema on the potential energy surfaces were located by optimizing the geometries of the reactant and product complexes as well as those of the transition states at the BHandHLYP/AUG-cc-pVDZ level of density functional theory in the gas phase. The solvent effect of aqueous media was treated by performing single point energy calculations at the BHandHLYP/AUG-cc-pVDZ and MP2/AUG-cc-pVDZ levels of theory employing the polarizable continuum model. Correction for basis set superposition error (BSSE) was made by the counterpoise method. Rate constants for all the reaction mechanisms were calculated including the tunneling contributions. Our calculations show that the hydrogen atom of the SH group of N-acetylcysteine would be most efficiently abstracted by the OH group, which is in agreement with experimental observations.

A quantum chemical study of repair of O6-methylguanine to guanine by tyrosine: evaluation of the winged helix-turn-helix model

The winged helix-turn-helix model for the repair of O6-MeG to guanine involving the reaction of O6-MeG with a tyrosine residue of the protein O6-alkylguanine-DNA alkyltransferase (AGT) was examined by studying the reaction mechanism and barrier energies. Molecular geometries of the species and complexes involved in the reaction, i.e. the reactant, intermediate and product complexes as well as transition states, were optimized employing density functional theory in gas phase. It was followed by single point energy calculations using density functional theory along with a higher basis set and second order M(phi)ller-Plesset perturbation theory (MP2) along with two different basis sets in gas phase and aqueous media. For the solvation calculations in aqueous media, the integral equation formalism of the polarizable continuum model (IEF-PCM) was employed. Vibrational frequency analysis was performed for each optimized structure and genuineness of transition states was ensured by visualizing the vibrational modes. It is found that tyrosine can repair O6-MeG to guanine by a two-step reaction. The present results have been compared with those obtained considering the helix-turn-helix model where the repair reaction primarily involves cysteine and occurs in a single-step. It is concluded that the repair through tyrosine envisaged in the winged helix-turn-helix model would be less efficient than that through cysteine envisaged in the helix-turn-helix model.

The DNA repair protein NBS1 influences the base excision repair pathway

NBS1 fulfills important functions for the maintenance of genomic stability and cellular survival. Mutations in the NBS1 (Nijmegen Breakage Syndrome 1) gene are responsible for the Nijmegen breakage syndrome (NBS) in humans. The symptoms of this disease and the phenotypes of NBS1-defective cells, especially their enhanced radiosensitivity, can be explained by an impaired DNA double-strand break-induced signaling and a disturbed repair of these DNA lesions. We now provide evidence that NBS1 is also important for cellular survival after oxidative or alkylating stress where it is required for the proper initiation of base excision repair (BER). NBS1 downregulated cells show reduced activation of poly-(adenosine diphosphate-ribose)-polymerase-1 (PARP1) following genotoxic treatment with H(2)O(2) or methyl methanesulfonate, indicating impaired processing of damaged bases by BER as PARP1 activity is stimulated by the single-strand breaks intermediately generated during this repair pathway. Furthermore, extracts of these cells have a decreased capacity for the in vitro repair of a double-stranded oligonucleotide containing either uracil or 8-oxo-7,8-dihydroguanine to trigger BER. Our data presented here highlight for the first time a functional role for NBS1 in DNA maintenance by the BER pathway.

An exploration of mechanisms for the transformation of 8-oxoguanine to guanidinohydantoin and spiroiminodihydantoin by density functional theory

The potential energy surface for formation of 2-amino-5-hydroxy-7,9-dihydropurine-6,8-dione (5-OH-OG), guanidinohydantoin (Gh) and spiroiminodihydantoin (Sp) from 8-oxoguanine (8-oxoG) has been mapped out using B3LYP density functional theory, the aug-cc-pVTZ and 6-31+G(d,p) basis sets and the IEF-polarizable continuum model (PCM) solvation model. Three pathways for formation of 5-OH-OG from 8-oxoG were evaluated: (A) stepwise loss of two electrons and two protons to form the quinonoid intermediate 2-amino-7,9-dihydro-purine-6,8-dione (8-oxoG(ox)) followed by hydration; (B) stepwise loss of two electrons and one proton and net addition of hydroxide, in which the key step is nucleophilic addition to the 8-oxoG radical cation; and (C) stepwise loss of one electron and one proton and addition of hydroxyl radical to the 8-oxoG radical cation. The data suggest that all three pathways are energetically feasible mechanisms for the formation of 5-OH-OG, however, Pathway A may be kinetically favored over Pathway B. Although lower in energy, Pathway C may be of limited biological significance since it depends on the local concentration of hydroxyl radical. Pathways for hydrolysis and decarboxylation of 5-OH-OG to form Gh via either a carboxylic acid or substituted carbamic acid intermediate have been evaluated with the result that cleavage of the N1-C6 bond is clearly favored over that of the C5-C6 bond. Formation of Sp from 5-OH-OG via stepwise proton transfer and acyl migration or ring opening followed by proton transfer and ring closure have also been explored and suggest that deprotonation of the hydroxyl group facilitates a 1,2 acyl shift. Results of the calculations are consistent with experimental studies showing dependence of the Gh/Sp product ratio on pH. Under neutral and basic conditions, the data predict that formation of Sp is kinetically favored over the pathways for formation of Gh. Under acidic conditions, Gh is predicted to be the kinetically favored product.

Ground- and Excited-State Proton Transfer and Rotamerism in 2-(2-Hydroxyphenyl)-5-phenyl-1,3,4-oxadiazole and Its O/“NH or S”-Substituted Derivatives

The intramolecular proton-transfer process, rotational process, and optical properties of 2-(2-hydroxyphenyl)-5-phenyl-1,3,4-oxadiazole (HOXD) and its O/"NH"- and O/"S"-substituted derivatives, 2-(2-hydroxyphenyl)-5-phenyl-1,3,4-triazole (HOT) and 2-(2-hydroxyphenyl)-5-phenyl-1,3,4-thiadiazole (HOTD), respectively, have been studied. DFT (B3LYP/6-31+G**) single-point energy calculations were performed using HF- and DFT-optimized geometries in the ground state (S0). TD-B3LYP/6-31+G** calculations using CIS-optimized geometries were carried out to investigate the properties of the first singlet excited state (S1) and first triplet excited state (T1). The computational results revealed that a high-energy barrier inhibits the proton transfer from cis-enol (Ec) to keto (K) form in S0, whereas the proton transfer in S1 can take place through a very-low-energy barrier. The rotation between Ec and trans-enol (Et) can occur in S0 through a low-energy barrier, whereas it is prohibited because of the high-energy barrier in S1 for each of the three molecules. The vertical excitation energies were calculated using the TD-B3LYP/6-31+G** method based on the HF- and CIS-optimized geometries. Absorption and fluorescence wavelengths of HOT show a hypsochromic shift (6-15 nm) relative to HOXD, while those of HOTD show a bathochromic shift (21-29 nm). The phosphorescence wavelength of HOTD shows a significant bathochromic shift relative to that of HOXD.

H2O3 as a Reactive Oxygen Species: Formation of 8-Oxoguanine from Its Reaction with Guanine

Reaction of guanine with H2O3 in the absence and presence of a water molecule leading to the formation of 8-oxoguanine (8-oxoG) was investigated. Initial calculations were performed using imidazole (Im) as a model for the five-membered ring of guanine. The reactant, intermediate, and product complexes as well as transition states were obtained in gas phase at the B3LYP/6-31+G* and B3LYP/AUG-cc-pVDZ levels of theory. In all the cases, except for the reactions involving imidazole, single-point energy calculations were performed in gas phase at the MP2/AUG-cc-pVDZ level of theory. Solvation calculations in aqueous media were carried out using the polarizable continuum model (PCM) of the self-consistent reaction field (SCRF) theory. Vibrational frequency analysis and intrinsic reaction coordinate (IRC) calculations were performed to ensure that the transition states connected the reactant and product complexes properly. Zero-point energy (ZPE)-corrected total energies and Gibbs free energies at 298.15 K in gas phase and aqueous media were obtained. When a reaction of H2O3 in place of H2O2 with guanine is considered, the major barrier energy which is encountered at the first step is almost halved showing that H2O3 would be much more reactive than H2O2. Considering the reaction schemes investigated here and the observed fact that H2O3 is dissociated easily under ambient conditions, it appears that H2O3 would serve as an effective reactive oxygen species.

The role of nucleobase carboradical and carbanion on DNA lesions: a theoretical study

DNA base release induced by H and OH radical addition to thymine and their corresponding electron adducts is studied at the DFT B3LYP/6-31+G(d,p) level in gas phase and in solution. H atom transfer after radical formation from C2' on the sugar to the C6 site on the base is shown to be prohibited for the radical species. Their corresponding electron adducts, albeit minor events in cellular systems, show excellent capabilities to proton transfer from C2' on the sugar to the C6 site on the base. The barriers for subsequent N-glycosidic bond dissociation range from 0.1 to 1.6 kcal mol(-1) at the B3LYP level and around 5 kcal mol(-1) using the BB1K functional, implying that these reactions can serve as a source to abasic sites. Analysis of bond dissociation energies show that all the reactions are exothermic, which is consistent with the changes in N-glycosidic bond lengths during the proton-transfer reactions. Bulk solvation plays a reverse influence on proton transfer and the bond rupture reactions. Molecular orbitals, NPA charges, and electron affinities are calculated to shed further light on the properties leading up to the intramolecular reactions.

Formation of 8-nitroguanine and 8-oxoguanine due to reactions of peroxynitrite with guanine

Reactions of peroxynitrite with guanine were investigated using density functional theory (B3LYP) employing 6-31G** and AUG-cc-pVDZ basis sets. Single point energy calculations were performed at the MP2/AUG-cc-pVDZ level. Genuineness of the calculated transition states (TS) was tested by visually examining the vibrational modes corresponding to the imaginary vibrational frequencies and applying the criterion that the TS properly connected the reactant and product complexes (PC). Genuineness of all the calculated TS was further ensured by intrinsic reaction coordinate (IRC) calculations. Effects of aqueous media were investigated by solvating all the species involved in the reactions using the polarizable continuum model (PCM). The calculations reveal that the most stable nitro-product complex involving the anion of 8-nitroguanine and a water molecule i.e. 8NO(2)G(-) + H(2)O can be formed according to one reaction mechanism while there are two possible reaction mechanisms for the formation of the oxo-product complex involving 8-oxoguanine and anion of the NO(2) group i.e. 8OG + NO(2)(-). The calculated relative stabilities of the PC, barrier energies of the reactions and the corresponding enthalpy changes suggest that formation of the complex 8OG + NO(2)(-) would be somewhat preferred over that of the complex 8NO(2)G(-) + H(2)O. The possible biological implications of this result are discussed.

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.

Successive Attachment of Electrons to Protonated Guanine: (G+H)• Radicals and (G+H)- Anions

The structures, energetics, and vibrational frequencies of nine hydrogenated 9H-keto-guanine radicals (G+H)(*) and closed-shell anions (G+H)(-) are predicted using the carefully calibrated (Chem. Rev. 2002, 102, 231) B3LYP density functional method in conjunction with a DZP++ basis set. These radical and anionic species come from consecutive electron attachment to the corresponding protonated (G+H)(+) cations in low pH environments. The (G+H)(+) cations are studied using the same level of theory. The proton affinity (PA) of guanine computed in this research (228.1 kcal/mol) is within 0.7 kcal/mol of the latest experiment value. The radicals range over 41 kcal/mol in relative energy, with radical r1, in which H is attached at the C8 site of guanine, having the lowest energy. The lowest energy anion is a2, derived by hydride ion attachment at the C2 site of guanine. No stable N2-site hydride should exist in the gas phase. Structure a9 was predicted to be dissociative in this research. The theoretical adiabatic electron affinities (AEA), vertical electron affinities, and vertical detachment energies were computed, with AEAs ranging from 0.07 to 3.12 eV for the nine radicals.

A theoretical study of some new analogues of the anti-cancer drug camptothecin

The enzyme topoisomerase I (topo I), which is essential for cell replication, transiently causes a DNA single strand break and makes a complex with it. The anti-cancer agent camptothecin (CPT) binds to the topo I-DNA complex and stabilizes it, preventing resealing of the broken DNA strand and cell growth. Considering the structural factors of CPT that are believed to be involved in stabilizing the topo I-DNA complex via hydrogen bonding and stacking interactions, designs of two new analogues of CPT (topo I inhibitors) have been suggested. The molecular geometries of CPT, two of its analogues and certain other related molecules included in the study were fully optimized in both gas phase and aqueous media at the B3LYP/6-311++G(d,p) level of density functional theory. Solvation effects of aqueous media were treated using the polarizable continuum model (PCM). Net CHelpG charges and surface molecular electrostatic potentials (MEP) near the atomic sites of the molecules were studied. Structural analogy and surface MEP values suggests that the two new CPT analogues studied here would be potent topoisomerase I inhibitors.