Ab Initio Study of Guanine Damage by Hydroxyl Radical

High-level ab initio mapping of the multiple H-abstraction pathways of the OH + glycine reaction

We perform a systematic search in the transition-state (TS) and product-channel complex (MIN) regions of the multi-channel OH + glycine → H₂O + H₂N-CH-COOH (CH)/HN-CH₂-COOH (NH)/H₂N-CH₂-COO (COOH) reactions. Geometry optimizations reveal {7, 3, 3} CH-TS, {2, 2, 2} CH-MIN, {17, 10, 5} NH-TS, {35, 19, 19} NH-MIN, and {6, 5, 5} COOH-TS conformers at the {MP2/3-21G, MP2/aug-cc-pVDZ, CCSD(T)-F12b/aug-cc-pVDZ} levels of theory as well as 2 additional CH-TSs based on chemical intuition. The benchmark relative energies of the TS, MIN, and product conformers are obtained by considering basis set effects up to aug-cc-pVQZ using the explicitly-correlated CCSD(T)-F12b method, post-(T) correlation up to CCSDT(Q), core correlation, scalar relativistic effects, spin-orbit coupling, and zero-point energy corrections. All the CH [ΔE_e(ΔH₀) = -38.54(-38.61) kcal mol^-1], NH [ΔE_e(ΔH₀) = -16.72(-17.98) kcal mol^-1], and COOH [ΔE_e = -4.98 kcal mol^-1] reactions are exothermic and proceed via shallow, usually negative, classical(adiabatic) barriers of -0.37(-0.95), -1.91(-2.48), and 1.02(-0.57) kcal mol^-1, respectively. In the entrance channel MRCI/aug-cc-pVTZ computations reveal several complexes with reactive(non-reactive) arrangements and binding energies of 1.0, 1.6, 3.3, (5.2 and 5.9) kcal mol^-1, stabilized by CH⋯OH, NH⋯OH, COOH⋯OH, (OH⋯OC and OH⋯N) hydrogen bonds, respectively.

On the relevance of hydroxyl radical to purine DNA damage

Hydroxyl radical (HO^•) is the most reactive toward DNA among the reactive oxygen species (ROS) generated in aerobic organisms by cellular metabolisms. HO^• is generated also by exogenous sources such as ionizing radiations. In this review we focus on the purine DNA damage by HO^• radicals. In particular, emphasis is given on mechanistic aspects for the various lesion formation and their interconnections. Although the majority of the purine DNA lesions like 8-oxo-purine (8-oxo-Pu) are generated by various ROS (including HO^•), the formation of 5',8-cyclopurine (cPu) lesions in vitro and in vivo relies exclusively on the HO^• attack. Methodologies generally utilized for the purine lesions quantification in biological samples are reported and critically discussed. Recent results on cPu and 8-oxo-Pu lesions quantification in various types of biological specimens associated with the cellular repair efficiency as well as with distinct pathologies are presented, providing some insights on their biological significance.

A theoretical study on the role of stability of cytosine and its tautomers in DNA (deoxyribonucleic acid), and investigation of interactions of Na+, K+, Mg2+, Ca2+, Zn2+ metal ions and OH radical with cytosine tautomers

In the present study, 21 cytosine tautomers were investigated so that some tautomers were reported for the first time in the gas phase and aqueous solution. C₃ tautomer was the most stable tautomer in gas phase but C₁ was the most stable structure in aqueous solution. The potential energy surface of all trajectories was determined for 21 tautomers and 22 transition states. Also, interactions of cytosine tautomers with Na⁺, K⁺, Mg²⁺, Ca²⁺ and Zn²⁺ metal ions were studied in gas phase and aqueous solution. Three types of interactions among metal ions and (N1 and O10), (N3 and O10) and (N3 and N9) of cytosine tautomers were investigated. The study of interaction energies of all complexes showed the stability of complexes in which interactions among Mg²⁺ and Zn²⁺ with tautomers were stronger than interactions among Ca²⁺, Na⁺ and K⁺ with tautomers, respectively. Some interactions of metal ions with cytosine tautomers made the most stable tautomers. So, the stability of rare tutomeric forms had a significant effect on stabilization of anomalous DNA (deoxyribonucleic acid) double helix and spontaneous mutations. Also, one of the most important causes of mutations in DNA (deoxyribonucleic acid) was the reaction of OH radical with nucleotide bases. So, interactions of OH radical with cytosine and its tautomers were investigated in gas phase and aqueous solution.Communicated by Ramaswamy H. Sarma.

DFT Study of the Hydroxyl Radical Addition to 2'-Deoxyguanosine and the Guanine Base in Four Double-Stranded B-Form Dimers

Density functional theory (DFT) calculations of reactions between 2'-deoxyguanosine (dR-Gua) and hydroxyl radical (HO^•) with water molecules (H₂O)_n, n = 0, 1, and 2, were carried out. The HO^• addition to three carbon sites, C(4), C(5), and C(8), and the subsequent ring cleavage of the three HO adducts were investigated. The addition to C(5) is of the smallest activation energy according to the largest lobe of the dR-Gua highest occupied molecular orbital (HOMO) at C(5). However, its adduct has small stability, and the C(8) adduct has the largest one. The C(8) adduct and the ring-opened amide have similar stability, which would lead to the apparent small yield of the former. Calculations were also performed on HO^• additions to the C(4) and C(8) sites of the guanine moiety of four dimer sequence models of B-form DNA with nucleotide moieties (a) 5'-GA-3', (b) 5'-GG-3', (c) 5'-GT-3', and (d) 5'-GC-3'. For instance, the (a) 5'-GA-3' model has a molecular formula C₃₉H₅₀N₁₅Na₂O₂₁P₂. The HO^• attack to C(4) is ruled out owing to the reinforced deformation of the parallel stacking of base pairs. The clear selectivity that the (b) 5'-GG-3' sequence is most reactive was found with the inclusion of the water dimer.

Targeted DNA oxidation and trajectory of radical DNA using DFT based QM/MM dynamics

Molecular insight into electronic rearrangements and structural trajectories arising from oxidative damages to DNA backbone is of crucial importance in understanding the effect of ionizing radiation, developing DNA biosensors and designing effective DNA cleaving molecules. Employing a Density Functional Theory based multi-scale Quantum-Mechanical-Molecular-Mechanical (QM/MM) simulation and a suitable partitioning of the Hamiltonian on solvated nucleotide, and single-, and double-stranded DNA, we mimic hydrogen transfer reactions from the backbone by OH radicals and report structural trajectories arising from on-the-fly electronic charge- and spin-density redistribution in these three different structural topologies of DNA. Trajectories reveal that H4′ abstraction can disrupt the deoxyribose moiety through the formation of C4′=O4′ ketone and a π-bond with base at C1′-N9 in a nucleotide versus only partial ketone formation in single- and double-stranded DNA, where the orientation of the base is topologically restrained. However, H5′ abstraction can lead DNA cleavage at 5′ end through the formation of C5′=O5′ ketone and breakage of P-O5′ bond. Results demonstrate that structural damages from oxidative reactions are restrained by base stacking and base-pair hydrogen bonding. The methodology can be suitably used to study targeted DNA and RNA damages from radicals and radiomimetic drugs to design DNA cleaving molecules for chemotherapy.

Catalytic Effect of Aqueous Solution in Water-Assisted Proton-Transfer Mechanism of 8-Hydroxy Guanine Radical

Water-assisted proton-transfer process is a key step in guanine damage reaction by hydroxyl radical in aqueous solution. In this article, we quantitatively determine the solvent effect in water-assisted proton-transfer mechanism of 8-hydroxy guanine radical using combined quantum mechanics and molecular mechanism with an explicit solvation model. Atomic-level reaction pathway was mapped, which shows a synchronized two-proton-transfer mechanism between the assistant water molecule and 8-hydroxy guanine radical. The transition-state dipole moment is the largest along the reaction pathway, which electrostatically stabilizes the proton-transfer transition-state complex. The free-energy reaction barrier for this water-assisted proton-transfer reaction was calculated at 19.2 kcal/mol with the density functional theory/M08-SO/cc-pVTZ+/molecular mechanics level of theory. The solvent effect not only has a big impact on geometries, but also dramatically changes the energetics along the reaction pathway. Among the solvent effect contributions to the transition state, the solvent energy contribution is -28.5 kcal/mol and the polarization effect contribution is 19.9 kcal/mol. In total, the solvent effect contributes -8.6 kcal/mol to the free-energy barrier height, which means that the presence of aqueous solution has a catalytic effect on the reaction mechanism and enhances the proton-transfer reactivity in aqueous solution.

Multi-level Quantum Mechanics and Molecular Mechanics Study of Ring Opening Process of Guanine Damage by Hydroxyl Radical in Aqueous Solution

Combining multi-level quantum mechanics theories and molecular mechanics with an explicit water model, we investigated the ring opening process of guanine damage by hydroxyl radical in aqueous solution. The detailed, atomic-level ring-opening mechanism along the reaction pathway was revealed in aqueous solution at the CCSD(T)/MM levels of theory. The potentials of mean force in aqueous solution were calculated at both the DFT/MM and CCSD(T)/MM levels of the theory. Our study found that the aqueous solution has a significant effect on this reaction in solution. In particular, by comparing the geometries of the stationary points between in gas phase and in aqueous solution, we found that the aqueous solution has a tremendous impact on the torsion angles much more than on the bond lengths and bending angles. Our calculated free-energy barrier height 31.6 kcal/mol at the CCSD(T)/MM level of theory agrees well with the one obtained based on gas-phase reaction profile and free energies of solvation. In addition, the reaction path in gas phase was also mapped using multi-level quantum mechanics theories, which shows a reaction barrier at 19.2 kcal/mol at the CCSD(T) level of theory, agreeing very well with a recent ab initio calculation result at 20.8 kcal/mol.