Theoretical study on protonated and deprotonated 5-substituted uracil derivatives and their complexes with water

Gas-Phase Experimental and Computational Studies of 5-Halouracils: Intrinsic Properties and Biological Implications

The gas-phase acidity and proton affinity (PA) of 5-halouracils (5-fluorouracil, 5-chlorouracil, 5-bromouracil, and 5-iodouracil) have been examined using both theoretical and experimental methods. This work represents a comprehensive study of the thermochemical properties of these nucleobases. Other than 5-fluorouracil acidity, the intrinsic acidity and PA of these halouracils have not been heretofore measured; these new experimental data provide a benchmark for the computational values. Furthermore, we examine these 5-halouracils in the context of the enzyme thymine DNA glycosylase (TDG), which is an enzyme that protects the genome by cleaving these substrates from DNA. Our gas-phase results are compared and contrasted to TDG excision rates to afford insights into the TDG mechanism.

5-Fluorouracil-Complete Insight into Its Neutral and Ionised Forms

5-Fluorouracil (5FU), a common anti-cancer drug, occurs in four tautomeric forms and possesses two potential sites of both protonation and deprotonation. Tautomeric and resonance structures of the ionized forms of 5FU create the systems of connected equilibriums. Since there are contradictory reports on the ionized forms of 5FU in the literature, complex theoretical studies on neutral, protonated and deprotonated forms of 5FU, based on the broad spectrum of DFT methods, are presented. These indicate that the O4 oxygen is more willingly protonated than the O2 oxygen and the N1 nitrogen is more willingly deprotonated than the N3 nitrogen in a gas phase. Such preferences are due to advantageous charge delocalization of the respective ions, which is demonstrated by the NBO and ESP analyses. In an aqueous phase, stability differences between respective protonated and deprotonated forms of 5FU are significantly diminished due to the competition between the mesomeric effect and solvation. The calculated pK_a values of the protonated, neutral and singly deprotonated 5FU indicate that 5FU does not exist in the protonated and double-deprotonated forms in the pH range of 0–14. The neutral form dominates below pH 8 and the N1 deprotonated form dominates above pH 8.

The prediction of intermolecular proton-transfer of guanine-cytosine base pair under the influence of fragments from decomposed MOFs

Metal-organic frameworks (MOFs) can be decomposed into various fragments, including negative/positive charges, Zn⁺ or Cu²⁺ when used as drug delivery materials. To evaluate the safety of MOFs, different mechanisms of intermolecular proton-transfer in guanine-cytosine (GC) base pair under the influence of such fragments were investigated by density functional theory methods. In a vacuum, calculation results show that an excess electron assists proton transfer in the anionic GC radical, and a hole assists proton transfer in the cationic GC radical with small energy barriers. The mechanism for Zn⁺-GC transfer is that the located hole assists proton transfer from G to C. All proton-transfers of Cu²⁺-GC become spontaneous with stable proton-transferred structures, and the driving force is the Cu²⁺ due to its electrostatic and oxidative effects. However, in a micro-water environment, the average energy barrier of all proton-transfer processes increases by 2.8 kcal mol^-1 because of the redistribution of charges. Water molecules play a very important role in buffering, and the influence of fragments on intermolecular proton-transfer processes of GC is reduced.

Theoretical Models for Quantitative Description of the Acid-Base Equilibria of the 5,6-Substituted Uracils

The acidities of 18 5,6-substituted uracils have been numerically estimated as pK_a values in terms of three theoretical models. The first scheme includes the calculation of the gas-phase acidity of uracil with the G3MP2B3 method and taking into account the solvent effect using the polarizable continuum approximation PCM(SMD)-TPSS/aug-cc-pVTZ. The second model is one step and implies calculation of the free Gibbs energies of the hydrate complex of uracil (and its anion) with 5 water molecules by the TPSS/aug-cc-pVTZ method. This model accounts for the solvent effect corresponding to both specific and nonspecific solvation. The third scheme required high time and computational resources and includes the strong features of the two previous schemes. Here, the theoretical estimation of pK_a is performed by the CBS-QB3 composite method. As in the second approach, both specific (as pentahydrate) and nonspecific solvent effects are determined. We have analyzed the advantages and model restrictions of the considered schemes for the pK_a calculations. All models have systematic errors, which have been corrected with the linear empirical regression relations. In the presented model, the absolute mean deviations of the pK_a values of uracils dissociating via the N1-H bonds diminish to 0.25, 0.28, and 0.23 pK_a units (respectively, for I, II, and III models), which corresponds to ∼0.3 kcal/mol on the energy scale. The applicability of our computational schemes to uracils dissociating via N3-H, O-H (orotic acids) and C-H bonds is discussed.

Microhydration of Deprotonated Nucleobases

Hydration reactions of deprotonated nucleobases (uracil, thymine, 5-fluorouracil,2-thiouracil, cytosine, adenine, and hypoxanthine) produced by electrospray have been experimentally studied in the gas phase at 10 mbar using a pulsed ion-beam high-pressure mass spectrometer. The thermochemical data, ΔH (o) , ΔS (o) , and ΔG (o) , for the monohydrated systems were determined. The hydration enthalpies were found to be similar for all studied systems and varied between 39.4 and 44.8 kJ/mol. A linear correlation was found between water binding energies in the hydrated complexes and the corresponding acidities of the most acidic site of nucleobases. The structural and energetic aspects of the precursors for the hydrated complexes are discussed in conjunction with available literature data. Graphical Abstract ᅟ.

Infrared spectra of protonated uracil, thymine and cytosine

The gas-phase structures of protonated uracil, thymine, and cytosine are probed by using mid-infrared multiple-photon dissociation (IRMPD) spectroscopy performed at the Free Electron Laser facility of the Centre Laser Infrarouge d'Orsay (CLIO), France. Experimental infrared (IR) spectra are recorded for ions that were generated by electrospray ionization, isolated, and then irradiated in a quadrupole ion trap; the results are compared to the calculated infrared absorption spectra of the different low-lying isomers (computed at the B3LYP/6-31++G(d,p) level). For each protonated base, the global energy minimum corresponds to an enolic tautomer, whose infrared absorption spectrum matched very well with the experimental IRMPD spectrum, with the exception of a very weak IRMPD signal observed at about 1800 cm(-1) in the case of the three protonated bases. This signal is likely to be the signature of the second-energy-lying oxo tautomer. We thus conclude that within our experimental conditions, two tautomeric ions are formed which coexist in the quadrupole ion trap.

Investigation of the intermolecular proton transfer in the supersystems adenine-methanol/ethanol/i-propanol: MP2 and DFT levels study

Twelve H-bonded supersystems constructed between the adenine tautomers and methanol, ethanol, and i-propanol were studied at the B3LYP and MP2 levels of theory using 6-311G(d,p) and 6-311++G(d,p) basis functions. The thermodynamic parameters of the complex formations were calculated in order to estimate the exact stability of the supersystems. It was proven that the calculated energy barriers of the alcohol-assisted proton transfers are about 60% lower than those of the intramolecular proton transfers in adenine found earlier (Gu and Leszczynski in J Phys Chem A 103:2744-2750, 1999).

Specificity of human thymine DNA glycosylase depends on N-glycosidic bond stability

Initiating the DNA base excision repair pathway, DNA glycosylases find and hydrolytically excise damaged bases from DNA. While some DNA glycosylases exhibit narrow specificity, others remove multiple forms of damage. Human thymine DNA glycosylase (hTDG) cleaves thymine from mutagenic G.T mispairs, recognizes many additional lesions, and has a strong preference for nucleobases paired with guanine rather than adenine. Yet, hTDG avoids cytosine, despite the million-fold excess of normal G.C pairs over G.T mispairs. The mechanism of this remarkable and essential specificity has remained obscure. Here, we examine the possibility that hTDG specificity depends on the stability of the scissile base-sugar bond by determining the maximal activity (k(max)) against a series of nucleobases with varying leaving-group ability. We find that hTDG removes 5-fluorouracil 78-fold faster than uracil, and 5-chlorouracil, 572-fold faster than thymine, differences that can be attributed predominantly to leaving-group ability. Moreover, hTDG readily excises cytosine analogues with improved leaving ability, including 5-fluorocytosine, 5-bromocytosine, and 5-hydroxycytosine, indicating that cytosine has access to the active site. A plot of log(k(max)) versus leaving-group pK(a) reveals a Brønsted-type linear free energy relationship with a large negative slope of beta(lg) = -1.6 +/- 0.2, consistent with a highly dissociative reaction mechanism. Further, we find that the hydrophobic active site of hTDG contributes to its specificity by enhancing the inherent differences in substrate reactivity. Thus, hTDG specificity depends on N-glycosidic bond stability, and the discrimination against cytosine is due largely to its very poor leaving ability rather than its exclusion from the active site.

Gas-Phase Deprotonation of Uracil−Cu2+ and Thiouracil−Cu2+ Complexes

The deprotonation of Cu2+ complexes with uracil, 2-thiouracil, 4-thiouracil, and 2,4-dithiouracil has been investigated by means of B3LYP/ 6-311+G(2df,2p)//6-31G(d) calculations. The most stable [(uracil-H)Cu]+ and [(thiouracil-H)Cu]+ complexes correspond to bidentate structures in which Cu interacts with the deprotonated ring-nitrogen atom and with the oxygen or the sulfur atom of the adjacent carbonyl or thiocarbonyl group. For 2- and 4-thiouracil derivatives, the structures in which the metal cation interacts with the thiocarbonyl group are clearly favored with respect to those in which Cu interacts with the carbonyl group. This is at variance with what was found to be the most stable structure of the corresponding Cu2+ complexes, where association to the carbonyl oxygen was always preferred over the association to the thiocarbonyl group. The [(uracil-H)Cu]+ and [(thiouracil-H)Cu]+ complexes can be viewed as the result of Cu+ attachment to the uracil-H and thiouracil-H radicals formed by the deprotonation of the corresponding uracil+* and thiouracil+* radical cations. As a matter of fact their relative stability is dictated by the intrinsic stability of the corresponding uracil-H and thiouracil-H radical and by the fact that, in general, the N3-deprotonated site is a better electron donor than the N1. In all complexes, the bonding of Cu both to nitrogen and sulfur and to nitrogen and oxygen has a significantly large covalent character.

Effects of Hydrogen Bonding on the Acidity of Adenine, Guanine, and Their 8-Oxo Derivatives

Complexes between ammonia, water, or hydrogen fluoride and adenine, guanine, or their 8-oxo derivatives are investigated using density-functional theory. The binding strengths of the neutral and (N9) anionic complexes are considered for a variety of purine binding sites. The effects of hydrogen-bonding interactions on the (N9) acidity of the purine derivatives are considered as a function of the molecule bound and the binding site. It is found that hydrogen-bonding interactions with one molecule can increase the acidity of purine derivatives by up to 60 kJ mol(-1). The (calculated) simultaneous effects of up to four molecules on the acidity of the purine derivatives are also considered. Our data suggest that the effects of more than one molecule on the acidity of the purines are generally less than the sum of the individual (additive) effects, where the magnitude of the deviation from additivity increases with the number, as well as the acidity, of molecules bound. Nevertheless, the increase in the acidity due to additional hydrogen-bonding interactions is significant, where the effect of two, three, or four hydrogen-bonding interactions can be as large as approximately 95, 115, and 130 kJ mol(-1), respectively. The present study provides a greater fundamental understanding of hydrogen-bonding interactions involving the natural purines, as well as those generated through oxidative DNA damage, which may aid the understanding of important biological processes.

Systematic Study of the Tautomerism of Uracil Induced by Proton Transfer. Exploration of Water Stabilization and Mutagenicity

To systematically investigate all the possible tautomerisms from uracil (U) and its enol form (U) induced by proton transfer, we describe a study of structural tautomer interconversion in the gas phase, in a continuum solvent, and in a microhydrated environment with 1 or 2 explicit water molecules, using density functional theory (DFT) calculations by means of the B3LYP exchange and correlation functions. A total of 62 geometries including 25 transition states were optimized, and the geometrical parameters have been discussed. Some rules of the configuration variation in tautomerization were summarized. The relative stabilities of all the tautomers were established. When a proton transfers from the di-keto form to the keto-enol form, water molecules in different regions show absolutely opposite effects: some assist, whereas others hinder the tautomerization. However, when a proton transfers from the keto-enol form to the di-enol form, water molecules in different regions show similar effects: the Gibbs free energy always increases and the activation energy always decreases. Additionally, some important factors that obviously affect the activation energy and Gibbs free energy were found and discussed in detail. The reasons that water molecules can assist or prevent the proton transfer were given. Furthermore, on the basis of our calculated results, we explain why it is hard to detect the di-enol form of uracil in general experiments.