A kinetic and thermodynamic study of the glycosidic bond cleavage in deoxyuridine

Preorganized Internal Electric Field Promotes a Double-Displacement Mechanism for the Adenine Excision Reaction by Adenine DNA Glycosylase

Adenine DNA glycosylase (MutY) is a monofunctional glycosylase, removing adenines (A) misinserted opposite 8-oxo-7,8-dihydroguanine (OG), a common product of oxidative damage to DNA. Through multiscale calculations, we decipher a detailed adenine excision mechanism of MutY that is consistent with all available experimental data, involving an initial protonation step and two nucleophilic displacement steps. During the first displacement step, N-glycosidic bond cleavage is accompanied by the attack of the carboxylate group of residue Asp144 at the anomeric carbon (C1'), forming a covalent glycosyl-enzyme intermediate to stabilize the fleeting oxocarbenium ion. After departure of the excised base, water nucleophiles can be recruited to displace Asp144, completing the catalytic cycle with retention of stereochemistry at the C1' position. The two displacement reactions are found to mostly involve the movement of the oxocarbenium ion, occurring with large charge reorganization and thus sensitive to the internal electric field (IEF) exerted by the polar protein environment. Intriguingly, we find that the negatively charged carboxylate group is a good nucleophile for the oxocarbenium ion, yet an unactivated water molecule is not, and that the electric field catalysis strategy is used by the enzyme to enable its unique double-displacement reaction mechanism. A strong IEF, pointing toward 5' direction of the substrate sugar ring, greatly facilitates the second displacement reaction at the expense of elevating the barrier of the first one, thereby allowing both reactions to occur. These findings not only increase our understanding of the strategies used by DNA glycosylases to repair DNA lesions, but also have important implications for how internal/external electric field can be applied to modulate chemical reactions.

Why 6-Iodouridine Cannot Be Used as a Radiosensitizer of DNA Damage? Computational and Experimental Studies

Previous density functional theory (DFT) studies on 6-brominated pyrimidine nucleosides suggest that 6-iodo-2'-deoxyuridine (6IdU) should act as a better radiosensitizer than its 5-iodosubstituted 2'-deoxyuridine analogue. In this work, we show that 6IdU is unstable in an aqueous solution. Indeed, a complete disappearance of the 6IdU signal was observed during its isolation by reversed-phase high-performance liquid chromatography (RP-HPLC). As indicated by the thermodynamic characteristics for the S_N1-type hydrolysis of 6IdU obtained at the CAM-B3LYP/DGDZVP++ level and the polarizable continuum model (PCM) of water, 6-iodouracil (6IU) was already released quantitatively at ambient temperatures. The simulation of the hydrolysis kinetics demonstrated that a thermodynamic equilibrium was reached within seconds for the title compound. To assess the reliability of the calculations carried out, we synthesized 6-iodouridine (6IUrd), which was, unlike 6IdU, sufficiently stable in an aqueous solution at room temperature. The activation barrier for the N-glycosidic bond dissociation in 6IUrd was estimated experimentally using an Arrhenius plot. The stabilities in water calculated for 6IdU, 6IUrd, and 5-iodo-2'-deoxyuridine (5IdU) could be explained by the electronic and steric effects of the 2'-hydroxy group present in the ribose moiety. Our studies highlight the issue of the hydrolytic stability of potentially radiosensitizing nucleotides which, besides having favorable dissociative electron attachment (DEA) characteristics, must be stable in water to have any practical application.

Development of QM/MM (ABEEM polarizable force field) method to simulate the excision reaction mechanism of damaged cytosine

DNA damages are regarded as having harmful effects on cell. The base excision repair mechanism combats these effects by removing damaged bases. The deglycosylation mechanism of excising damaged bases by DNA glycosylase and the state of the leaving base have been controversial. The enzymatic reaction of DNA glycosylase to remove the damaged bases involves not only the formation and breaking of chemical bonds, but also complex polarization effect and charge transfer, which cannot be accurately simulated by the QM/MM method combined with the fixed charge force field. This work has developed the ABEEM fluctuating polarizable force field combining with the QM method, that is (QM/MM[ABEEM]), to accurately simulate the proton transfer, charge transfer and the charge distribution. The piecewise function is used as the valence-state electronegativity in the QM/MM (ABEEM) to realize the accurate fitting of the charge distribution in reaction. And the charge transfer is accurately simulated by the local charge conservation conditions. Four deglycosylation mechanisms including the monofunctional and difunctional mechanisms of four neutral and protonated cytosine derivatives are explored. It is confirmed that the monofunctional mechanism of Asp-activated nucleophile water is a better deglycosylation mechanism and the base is protonated before the reaction occurs. Protonization of the base reduced the activation energy by 10.00-17.00 kcal/mol. Asp provides the necessary charge for the reaction, and DNA glycosylase preferentially cleaves ɛC. This work provides a theoretical basis for the research of excising damaged bases by DNA glycosylase.

DFT Study on the Deglycosylation of Methylated, Oxidized, and Canonical Pyrimidine Nucleosides in Water: Implications for Epigenetic Regulation and DNA Repair

Density functional theory (B3LYP) was used to characterize the kinetics and thermodynamics of the (nonenzymatic) deglycosylation in water for a variety of 2'-deoxycytidine (dC) and 2'-deoxyuridine (dU) nucleoside derivatives that differ in methylation and subsequent oxidation of the C5 substituent. A range of computational models are considered that combine implicit and explicit solvation of the nucleophile and nucleobase. Regardless of the model implemented, our calculations reveal that the glycosidic bond in dC is inherently more stable than that in dU. Furthermore, C5 methylation of either pyrimidine and subsequent oxidation of the methyl group yield overall small changes to the Gibbs reaction energy profiles and thereby preserve lower deglycosylation barriers for the dC compared to those for the dU nucleoside derivatives. However, hydrolytic deglycosylation becomes significantly more energetically favorable when 5-methyl-dC (5m-dC) undergoes two or three rounds of oxidation, with the Gibbs energy barrier decreasing and the reaction becoming more exergonic by up to 40 kJ/mol. In fact, two or three oxidation reactions from 5m-dC result in a deglycosylation barrier similar to that for dU, as well as those for the associated C5-methylated (2'-deoxythymidine) and oxidized (5-hydroxymethyl-dU) derivatives. These predicted trends in the inherent deglycosylation energetics in water directly correlate with the previously reported activity of thymine DNA glycosylase (TDG), which cleaves the glycosidic bond in select dC nucleosides as part of epigenetic regulation and in dU variants as part of DNA repair. Thus, our data suggests that fundamental differences in the intrinsic reactivity of the pyrimidine nucleosides help regulate the function of human enzymes that maintain cellular integrity.

Radiation-induced curcumin release from curcumin–chitosan polymer films

The probability of human exposure to damaging radiation is increased in activities associated with long-term space flight, medical radiation therapies, and responses to nuclear accidents. However, the development of responsive countermeasures to combat radiation damage to biological tissue is lagging behind rates of human exposure. Herein, we report a radiation-responsive drug delivery system that releases doses of curcumin from a chitosan polymer/film in response to low level gamma radiation exposure. As a fibrous chitosan–curcumin polymer, 1 Gy gamma irradiation (¹³⁷Cs) released 5 ± 1% of conjugated curcumin, while 6 Gy exposure releases 98 ± 1% of conjugated curcumin. The same polymer was formed into a film through solvent casting. The films showed similar, albeit attenuated behavior in water (100% released) and isopropyl alcohol (32% released) with statistically significant drug release following 2 Gy irradiation. ATR FT-IR studies confirmed glycosidic bond cleavage in the chitosan–curcumin polymer in response to gamma radiation exposure. Similar behavior was noted upon exposure of the polymer to 20 cGy (1 GeV amu⁻¹, at 20 cGy min⁻¹) high linear energy transfer (LET) ⁵⁶Fe radiation based on FTIR studies. Density Functional Theory calculations indicate homolytic bond scission as the primary mechanism for polymer disintegration upon radiation exposure. Films did not change in thickness during the course of radiation exposure. The successful demonstration of radiation-triggered drug release may lead to new classes of radio-protective platforms for developing countermeasures to biological damage from ionizing radiation.

The probability of human exposure to damaging radiation is increased in activities associated with long-term space flight, medical radiation therapies, and responses to nuclear accidents.

Microsolvated Model for the Kinetics and Thermodynamics of Glycosidic Bond Dissociative Cleavage of Nucleoside D4G

Using the microsolvated model that involves explicit water molecules and implicit solvent in the optimization, two proposed dissociative hydrolysis mechanisms of 2',3'-didehydro-2',3'-dideoxyguanosine (d4G) have been first investigated by means of M06-2X(CPCM, water)/6-31++G(d,p) method. The glycosidic bond dissociation for the generation of the oxacarbenium ion intermediate is the rate-determining step (RDS). The subsequent nucleophilic water attack from different side of the oxacarbenium ion intermediate gives either the α-product [(2S,5S)-5-(hydroxymethyl)-2,5-dihydrofuran-2-ol] or β-product [(2R,5S)-5-(hydroxymethyl)-2,5-dihydrofuran-2-ol] and is thus referred to as α-path (inversion) and β-path (retention). Two to five explicit water molecules (n = 2-5) are considered in the microsolvated model, and n = 3 or 4 is the smallest model capable of minimizing the activation energy for α-path and β-path, respectively. Our theoretical results suggest that α-path (n = 3) is more kinetically favorable with lower free energy barrier (RDS) of 27.7 kcal mol^-1, in contrast to that of 30.7 kcal mol^-1 for the β-path (n = 4). The kinetic preference of the α-path is rationalized by NBO analysis. Whereas thte β-path is more thermodynamically favorable over the α-path, where the formation of β-product and α-product are exergonic and endergonic, respectively, providing theoretical support for the experimental observation that the β-cleavage product was the major one after sufficient reaction time. Comparisons of d4G with analogous cyclo-d4G and dG from kinetic free energy barriers and thermodynamic heterolytic dissociation energies were also carried out. Our kinetic and thermodynamic results manifest that the order of glycosidic bond stability should be d4G < cyclo-d4G < dG, which agrees well with the reported experimental stability order of d4G compounds and analogues and gives further understanding on the influence of 6-cyclopropylamino and unsaturated ribose to the glycosidic bond instability of d4G.

Hydrolytic Glycosidic Bond Cleavage in RNA Nucleosides: Effects of the 2'-Hydroxy Group and Acid-Base Catalysis

Despite the inherent stability of glycosidic linkages in nucleic acids that connect the nucleobases to sugar-phosphate backbones, cleavage of these bonds is often essential for organism survival. The current study uses DFT (B3LYP) to provide a fundamental understanding of the hydrolytic deglycosylation of the natural RNA nucleosides (A, C, G, and U), offers a comparison to DNA hydrolysis, and examines the effects of acid, base, or simultaneous acid-base catalysis on RNA deglycosylation. By initially examining HCOO^-···H₂O mediated deglycosylation, the barriers for RNA hydrolysis were determined to be 30-38 kJ mol^-1 higher than the corresponding DNA barriers, indicating that the 2'-OH group stabilizes the glycosidic bond. Although the presence of HCOO^- as the base (i.e., to activate the water nucleophile) reduces the barrier for uncatalyzed RNA hydrolysis (i.e., unactivated H₂O nucleophile) by ∼15-20 kJ mol^-1, the extreme of base catalysis as modeled using a fully deprotonated water molecule (i.e., OH^- nucleophile) decreases the uncatalyzed barriers by up to 65 kJ mol^-1. Acid catalysis was subsequently examined by selectively protonating the hydrogen-bond acceptor sites of the RNA nucleobases, which results in an up to ∼80 kJ mol^-1 barrier reduction relative to the corresponding uncatalyzed pathway. Interestingly, the nucleobase proton acceptor sites that result in the greatest barrier reductions match sites typically targeted in enzyme-catalyzed reactions. Nevertheless, simultaneous acid and base catalysis is the most beneficial way to enhance the reactivity of the glycosidic bonds in RNA, with the individual effects of each catalytic approach being weakened, additive, or synergistic depending on the strength of the base (i.e., degree of water nucleophile activation), the nucleobase, and the hydrogen-bonding acceptor site on the nucleobase. Together, the current contribution provides a greater understanding of the reactivity of the glycosidic bond in natural RNA nucleosides, and has fundamental implications for the function of RNA-targeting enzymes.

Tautomerization lowers the activation barriers for N-glycosidic bond cleavage of protonated uridine and 2'-deoxyuridine

The gas-phase conformations of protonated uridine, [Urd+H](+), and its 2'-deoxy form, protonated 2'-deoxyuridine, [dUrd+H](+), have been examined in detail previously by infrared multiple photon dissociation action spectroscopy techniques. Both 2,4-dihydroxy tautomers and O4 protonated conformers of [Urd+H](+) and [dUrd+H](+) were found to coexist in the experiments with the 2,4-dihydroxy tautomers dominating the population. In the present study, the kinetic energy dependence of the collision-induced dissociation behavior of [Urd+H](+) and [dUrd+H](+) are examined using a guided ion beam tandem mass spectrometer to probe the mechanisms and energetics for activated dissociation of these protonated nucleosides. The primary dissociation pathways observed involve N-glycosidic bond cleavage leading to competitive elimination of protonated or neutral uracil. The potential energy surfaces (PESs) for these N-glycosidic bond cleavage pathways are mapped out via electronic structure calculations for the mixture of 2,4-dihydroxy tautomers and O4 protonated conformers of [Urd+H](+) and [dUrd+H](+) populated in the experiments. The calculated activation energies (AEs) and heats of reaction (ΔHrxns) for N-glycosidic bond cleavage at both the B3LYP and MP2(full) levels of theory are compared to the measured values. The agreement between experiment and theory indicates that B3LYP provides better estimates of the energetics of the species along the PESs for N-glycosidic bond cleavage than MP2, and that the 2,4-dihydroxy tautomers, which are stabilized by strong hydrogen-bonding interactions, predominantly influence the observed threshold dissociation behavior of [Urd+H](+) and [dUrd+H](+).

Fragmentation of electrospray-produced deprotonated ions of oligodeoxyribonucleotides containing an alkylated or oxidized thymidine

Alkylation and oxidation constitute major routes of DNA damage induced by endogenous and exogenous genotoxic agents. Understanding the biological consequences of DNA lesions often necessitates the availability of oligodeoxyribonucleotide (ODN) substrates harboring these lesions, and sensitive and robust methods for validating the identities of these ODNs. Tandem mass spectrometry is well suited for meeting these latter analytical needs. In the present study, we evaluated how the incorporation of an ethyl group to different positions (i.e., O(2), N3, and O(4)) of thymine and the oxidation of its 5-methyl carbon impact collisionally activated dissociation (CAD) pathways of electrospray-produced deprotonated ions of ODNs harboring these thymine modifications. Unlike an unmodified thymine, which often manifests poor cleavage of the C3'-O3' bond, the incorporation of an alkyl group to the O(2) position and, to a much lesser extent, the O(4) position, but not the N3 position of thymine, led to facile cleavage of the C3'-O3' bond on the 3' side of the modified thymine. Similar efficient chain cleavage was observed when thymine was oxidized to 5-formyluracil or 5-carboxyluracil, but not 5-hydroxymethyluracil. Additionally, with the support of computational modeling, we revealed that proton affinity and acidity of the modified nucleobases govern the fragmentation of ODNs containing the alkylated and oxidized thymidine derivatives, respectively. These results provided important insights into the effects of thymine modifications on ODN fragmentation.

Selecting DFT methods for use in optimizations of enzyme active sites: applications to ONIOM treatments of DNA glycosylases

When using a hybrid methodology to treat an enzymatic reaction, many factors contribute to selecting the method for the high-level region, which can be complicated by the presence of dispersion-driven interactions such as π–π stacking. In addition, the proper treatment of the reaction center often requires a large number of heavy atoms to be included in the high-level region, precluding the use of ab initio methods such as MP2 as well as large basis sets, in the optimization step. In the present work, popular DFT methods were tested to identify an appropriate functional for treating the high-level region in ONIOM optimizations of reactions catalyzed by nonmetalloenzymes. Eight different DFT methods (B3LYP, B97-2, MPW1K, MPWB1K, BB1K, B1B95, M06-2X, and ωB97X-D) in combination with four double-ζ quality Pople basis sets were tested for their ability to optimize noncovalent interactions (hydrogen bonding and π–π) and characterize reactions (proton transfer, SN2 hydrolysis, and unimolecular cleavage). Although the primary focus of this study is accurate structure determination, energetics were also examined at both the optimization level of theory, and with triple-ζ quality basis set and select (M06-2X or ωB97X-D) methods. If dispersion-driven interactions exist within the active site, then MPWB1K/6-31G(d,p) or M06-2X/6-31+G(d,p) are recommended for the optimization step with subsequent triple-ζ quality single-point energies. However, since dispersion-corrected functionals (M06-2X and ωB97X-D) generally require diffuse functions to yield appropriate geometries, the possible size of the high-level region is greatly limited with these methods. In contrast, if the model is large enough to recover steric constraints on π–π interactions, then B3LYP with a small basis set performs comparatively well for the optimization step and is significantly less computationally expensive. Interestingly, the functionals that afford the best geometries often do not yield the best energetics, which emphasizes the importance of structural benchmark studies.

Glycosidic bond cleavage in deoxynucleotides: effects of solvent and the DNA phosphate backbone in the computational model

Density functional theory (B3LYP) was employed to examine the hydrolysis of the canonical 2'-deoxynucleotides in varied environments (gas phase or water) using different computational models for the sugar residue (methyl or phosphate group at C5') and nucleophile (water activated through full or partial proton abstraction). Regardless of the degree of nucleophile activation, our results show that key geometrical parameters along the reaction pathway are notably altered upon direct inclusion of solvent effects in the optimization routine, which leads to significant changes in the reaction energetics and better agreement with experiment. Therefore, despite the wide use of gas-phase calculations in the literature, small model computational work, as well as large-scale enzyme models, that strive to understand nucleotide deglycosylation must adequately describe the environment. Alternatively, although inclusion of the phosphate group at C5' also affects the geometries of important stationary points, the effects cancel to yield unchanged deglycosylation barriers, and therefore smaller computational models can be used to estimate the energy associated with nucleotide deglycosylation, with the 5' phosphate group included if full (geometric) details of the reaction are desired. Hydrogen-bonding interactions with the nucleobase can significantly reduce the barrier to deglycosylation, which supports suggestions that discrete hydrogen-bonding interactions with active-site amino acid residues can play a significant role in enzyme-catalyzed nucleobase excision. Taken together with previous studies, the present work provides vital clues about the components that must be included in future studies of the deglycosylation of isolated noncanonical nucleotides, as well as the corresponding enzyme-catalyzed reactions.

A density functional theory study on the kinetics and thermodynamics of N-glycosidic bond cleavage in 5-substituted 2'-deoxycytidines

B3LYP/6-311+G(2d,p)//B3LYP/6-31+G(d) density functional theory calculations were employed to explore the kinetics and thermodynamics of gas-phase N-glycosidic bond cleavage induced by nucleophilic attack of C1' with a hydroxide ion in 5-substituted 2'-deoxycytidines. The results showed that, among the 5-substituted 2'-deoxycytidine derivatives examined [XdC, where X = H (dC), CH(3) (medC), CH(2)OH (hmdC), CHO (fmdC), COOH (cadC), F (FdC), or Br (BrdC)], fmdC and cadC exhibited the lowest energy barrier and largest exothermicity for N-glycosidic bond cleavage. These results paralleled previously reported nucleobase excision activities of human thymine DNA glycosylase (hTDG) toward duplex DNA substrates harboring a thymine and 5-substituted cytosine derivatives when paired with a guanine. Our study suggests that the inherent chemistry associated with the nucleophilic cleavage of N-glycosidic bond constitutes a major factor contributing to the selectivity of hTDG toward 5-substituted dC derivatives. These findings provided novel insights into the role of TDG in active cytosine demethylation.

Combined effects of π-π stacking and hydrogen bonding on the (N1) acidity of uracil and hydrolysis of 2'-deoxyuridine

M06-2X/6-31+G(d,p) is used to study the simultaneous effects of π-π stacking interactions with phenylalanine (modeled as benzene) and hydrogen bonding with small molecules (HF, H(2)O, and NH(3)) on the N1 acidity of uracil and the hydrolytic deglycosylation of 2'-deoxyuridine (dU) (facilitated by fully (OH(-)) or partially (HCOO(-)···H(2)O) activated water). When phenylalanine is complexed with isolated uracil, the proton affinity of all acceptor sites significantly increases (by up to 28 kJ mol(-1)), while the N1 acidity slightly decreases (by ~6 kJ mol(-1)). When small molecules are hydrogen bound to uracil, addition of the phenylalanine ring can increase or decrease the acidity of uracil depending on the number and nature (acidity) of the molecules bound. Furthermore, a strong correlation between the effects of π-π stacking on the acidity of U and the dU deglycosylation reaction energetics is found, where the hydrolysis barrier can increase or decrease depending on the nature and number of small molecules bound, the nucleophile considered (which dictates the negative charge on U in the transition state), and the polarity of the (bulk) environment. These findings emphasize that the catalytic (or anticatalytic) role of the active-site aromatic amino acid residues is highly dependent on the situation under consideration. In the case of uracil-DNA glycosylase (UNG), which catalyzes the hydrolytic excision of uracil from DNA, the type of discrete hydrogen-bonding interactions with U, the nature of the nucleophile, and the anticipated weak, nonpolar environment in the active site suggest that phenylalanine will be slightly anticatalytic in the chemical step, and therefore experimentally observed contributions to catalysis may entirely result from associated structural changes that occur prior to deglycosylation.

A QM/QM investigation of the hUNG2 reaction surface: the untold tale of a catalytic residue

Human uracil-DNA glycosylase (hUNG2) is a base excision repair enzyme that removes the damaged base uracil from DNA through hydrolytic deglycosylation of the nucleotide. In the present study, the mechanism of hUNG2 is thoroughly investigated using ONIOM(MPWB1K/6-31G(d):PM3) active-site models to generate reaction potential energy surfaces. Active-site models that differ in the hydrogen-bonding arrangement of the nucleophilic water molecule and/or protonation state of His148 are considered. The large barrier calculated using the model with a cationic His148 verifies that this residue is neutral in the early stages of the reaction. The reaction pathways predicted by two models with a neutral His148 are consistent with a wealth of experimental data on the enzyme, including mutational studies, which supports our approach. On the basis of our calculations, we propose a complete mechanism for the chemical step of hUNG2. In the first part of the reaction, His268, Asn204, and a water molecule work together to stabilize the negative charge forming on the uracil moiety. Subsequently, either Asp145 or His148 can act as the general base that activates the water nucleophile depending on the binding orientation of the water molecule in the active site. However, we propose that His148 preferentially acts as the general base. Therefore, in agreement with previous proposals, we assign the primary function of Asp145 to electrostatic stabilization of the positive charge developing on the sugar moiety during the reaction, which is also consistent with a growing theory that the primary function of active-site carboxylate groups present in many glycosylases is transition state stabilization. Most importantly, our work explains, for the first time, the role of His148 in the chemical step and provides additional support for the inclusion of this amino acid in the list of residues (Asp145 and His268) essential to the chemical step of the hUNG2 mechanism.

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.

Effects of nucleophile, oxidative damage, and nucleobase orientation on the glycosidic bond cleavage in deoxyguanosine

Deglycosylation of nucleotides occurs during many essential biological processes, including DNA repair, and is initiated by a variety of nucleophiles. In the present work, density functional theory (B3LYP) was used to investigate the thermodynamics and kinetics of the glycosidic bond cleavage reaction in the model nucleoside forms of guanine and its major oxidation product, 8-oxoguanine. Base excision facilitated by four different nucleophiles (hydroxyl anion (fully activated water), formate-water complex (partially activated water), lysine, and proline) was considered, which spans nucleophiles involved in a collection of spontaneous and enzyme-catalyzed processes. Because some enzymes that catalyze deglycosylation can accommodate more than one orientation of the base with respect to the sugar moiety, the effects of the (anti/syn) base orientation on the barrier height were also considered. We find that the nucleophile has a very large effect on the overall (gas-phase) reaction energetics. Although this effect decreases in different (polar) environments, the nucleophile has the greatest influence on the overall reaction as compared to whether the base is damaged or to the base orientation. Furthermore, the effects are significant in environments that most closely resemble (nonpolar) enzymatic active sites. Our results provide a greater understanding of the relative effects of the nucleophile, damage to the nucleobase, and the nucleobase orientation with respect to the sugar moiety on the deglycosylation pathway, which provide qualitative explanations for relative base excision rates observed in some biological systems.

Modeling the dissociative hydrolysis of the natural DNA nucleosides

Two-dimensional PCM-B3LYP/6-31+G(d) potential energy surfaces for the hydrolysis of the four natural 2'-deoxyribonucleosides (2'-deoxyadenosine, 2'-deoxyguanosine, 2'-deoxycytidine, and thymidine) are characterized using a model that includes both implicit (bulk) solvent effects and (three or four) explicit water molecules in the optimization routine. For the first time, the experimentally predicted dissociative (S(N)1) mechanism is found to be favored over the synchronous (S(N)2) pathway for all nucleosides studied. Due to the success of our model in stabilizing the charge-separated intermediates along the S(N)1 pathway, it is proposed that the new model presented here is the smallest system capable of generating the experimentally predicted oxacarbenium cation intermediate. We therefore stress that dissociative mechanisms should be studied with methodologies that account for the (bulk) environment in the optimization routine, where these effects are often only included as a correction to the energy in the current literature. In addition to accounting for charge stabilization through implicit solvation, nucleophile activation and leaving group stabilization should also be explicitly introduced into the model to further stabilize the system. Our work also emphasizes the importance of studying the Gibbs surface, which in some cases provides a better description of chemically important regions of the reaction surface or changes the calculated trend in the magnitude of dissociative barriers. In addition, it is proposed that the methodology presented in this study can be used to calculate uncatalyzed deglycosylation barriers for a range of DNA nucleosides, which when compared to the corresponding enzyme-catalyzed reactions, will allow the prediction of the rate enhancement (barrier reduction) due to the enzyme.

Glycosidic bond cleavage in deoxynucleotides — A density functional study

Density functional theory was used to study the glycosidic bond cleavage in deoxynucleotides with the main goal to determine the effects of the nucleobase, hydrogen bonding with the nucleobase, and the (bulk) environment on the reaction energetics. Since direct glycosidic bond cleavage is a high-energy process, two nucleophile models were considered (HCOO–···H2O and HO–), which represent different stages of activation of a water nucleophile. The glycosidic bond cleavage barriers were found to decrease, while the reaction exothermicity increases, with an increase in the nucleobase acidity. The gas-phase barriers and reaction energies for bond cleavage in all deoxynucleotides were found to be significantly affected by hydrogen-bonding interactions with the nucleobase (by up to 30 kJ mol–1 depending on the nucleophile). Although the barriers increase and reaction energies become less exothermic in enzymatic and aqueous environments, the effects of the bulk environment are similar in the presence and absence of small molecules bound to the nucleobase. Therefore, the effects of hydrogen bonding with the bases are approximately the same in all environments. Our results suggest that hydrogen bonding with the nucleobase may play an important role in the glycosidic bond cleavage in both pyrimidine and purine nucleotides in a variety of environments.