Development of a QM/MM(ABEEM) method for the deprotonation of neutral and cation radicals in the G-tetrad and GGX(8-oxo-G) tetrad

The rapid deprotonation of G˙+ in the DNA strand impedes positive charge (hole) transfer, whereas the slow deprotonation rate of G˙+ in the G-tetrad makes it a more suitable carrier for hole conduction. The QM/MM(ABEEM) combined method, which involves the integration of QM and the ABEEM polarizable force field (ABEEM PFF), was developed to investigate the deprotonation of neutral and cation free radicals in the G-tetrad and GGX(8-oxo-G) tetrad (xanthine and 8-oxoguanine dual substituted G-tetrad). By incorporating valence-state electronegativity piecewise functions χ*(r) and implementing charge local conservation conditions, QM/MM(ABEEM) possesses the advantage of accurately simulating charge transfer and polarization effect during deprotonation. The activation energy calculated by the QM method of X˙ is the lowest among other bases in the GGX(8-oxo-G) tetrad, which is supported by the computation of the average electronegativity calculated by ABEEM PFF. By utilizing QM/MM(ABEEM) with a two-way free energy perturbation method, the deprotonation activation energy of X˙ in the GGX(8-oxo-G) tetrad is determined to be 33.0 ± 2.1 kJ mol-1, while that of G˙+ in the G-tetrad is 20.7 ± 0.6 kJ mol-1, consistent with the experimental measurement of 20 ± 1.0 kJ mol-1. These results manifest that X˙ in the GGX(8-oxo-G) tetrad exhibits a slower deprotonation rate than G˙+ in the G-tetrad, suggesting that the GGX(8-oxo-G) tetrad may serve as a more favorable hole transport carrier. Furthermore, the unequal average electronegativities of bases in the GGX(8-oxo-G) tetrad impede the deprotonation rate. This study provides a potential foundation for investigating the microscopic mechanism of DNA electronic devices.

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.

DNA Nucleobase under Ionizing Radiation: Unexpected Proton Transfer by Thymine Cation in Water Nanodroplets

The effect of ionizing radiation on DNA constituents is a widely studied fundamental process using experimental and computational techniques. In particular, radiation effects on nucleobases are usually tackled by mass spectrometry in which the nucleobase is embedded in a water nanodroplet. Here, we present a multiscale theoretical study revealing the effects and the dynamics of water droplets towards neutral and ionized thymine. In particular, by using both hybrid quantum mechanics/molecular mechanics and full ab initio molecular dynamics, we reveal an unexpected proton transfer from thymine cation to a nearby water molecule. This leads to the formation of a neutral radical thymine and a Zundel structure, while the hydrated proton localizes at the interface between the deprotonated thymine and the water droplet. This observation opens entirely novel perspectives concerning the reactivity and further fragmentation of ionized nucleobases.

Reactions of water with radical cations of guanine, 9-methylguanine, 2′-deoxyguanosine and guanosine: keto–enol isomerization, C8-hydroxylation, and effects of N9-substitution

The reactions of D₂O with guanine radical cations in nucleobases and nucleosides were studied in the gas phase using the guided-ion-beam experiment and computational modeling.