A quantum theoretical study of reactions of methyldiazonium ion with DNA base pairs

Role of graphene in scavenging methyl cations: a DFT study

CONTEXT

It is known that methylating agents methylate DNA by transferring a methyl cation (CH3+) to the nucleophilic sites in DNA bases and DNA methylation is implicated in cancer and other pathological conditions. Therefore, it is important to scavenge CH3+ ion in order to protect DNA from methylation. Graphene is considered to be a versatile material for use in a wide variety of fields including sensors, antioxidants, drug delivery and DNA sequencing. In this work, we have theoretically investigated the interaction of CH3+ ions with graphene surface with an aim to understand if pristine graphene can be used as a substrate to adsorb CH3+ cations generated from harmful methylating agents. The computed adsorption energies show that adsorption of one, two and three CH3+ ions on graphene is favourable as the adducts thus formed are found to be substantially stable in both gas phase and aqueous media. The Bader charge transfer analysis and density of states (DOS) calculation also indicate a strong interaction between graphene and CH3+ ions. Thus, our results show that pristine graphene can be used as a substrate to scavenge CH3+ ions.

METHODS

The spin polarised density functional theory (DFT) calculations employing PBE functional, ultrasoft pseudopotentials and plane wave basis set having kinetic energy cut-offs of 40 Ry and 400 Ry, respectively, for wave functions and charge densities were carried out to study the adsorption of CH3+ ion(s) on the pristine graphene surface. The Grimme's DFT-D2 method was used for the estimation of van der Waals interactions. The 'dipole correction' along z-direction was also applied for adsorption study. The Marzari-Vanderbilt smearing and Monkhorst-Pack k-point grid were employed for the Brillouin zone sampling. A 6 × 6 graphene supercell with a vertical cell dimension of 18 Å was considered for the adsorption study. The charge transfer between the CH3+ ion(s) and graphene was estimated using Bader charge analysis. The implicit solvation model (SCCS) was used to estimate the solvent effect of aqueous media. All the calculations were performed using QUANTUM ESPRESSO package.

Free energy profiles for two ubiquitous damaging agents: methylation and hydroxylation of guanine in B-DNA

DNA methylation and hydroxylation are two ubiquitous reactions in DNA damage induction, yet insights are scarce concerning the free energy of activation within B-DNA. We resort to multiscale simulations to investigate the attack of a hydroxyl radical and of the primary diazonium onto a guanine embedded in a solvated dodecamer. Reaction free energy profiles characterize two strongly exergonic processes, yet allow unprecedented quantification of the barrier towards this damage reaction, not higher than 6 kcal mol^-1 and sometimes inexistent, and of the exergonicities. In the case of the [G(C8)-OH]˙ intermediate, we challenge the functional dependence of such simulations: recently-proposed functionals, such as M06-2X and LC-BLYP, agree on a ∼4 kcal mol^-1 barrier, whereas the hybrid GGA B3LYP functional predicts a barrier-less pathway. In the long term, multiscale approaches can help build up a unified panorama of DNA lesion induction. These results stress the importance of DFT/MM-MD simulations involving new functionals towards the sound modelling of biomolecule damage even in the ground state.