Effects of Hydrogen Bonding on the Acidity of Uracil

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.

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.

Vibrational Energy Flow in the Uracil-H2O Complexes

In the uracil-H₂O complex, the vibrational energy initially stored in the OH(v = 1) stretch efficiently transfers to the first overtone-bending mode under a near-resonant condition. The relaxation of the overtone vibration redistributes its energy to uracil and the two hydrogen bonds in the intermolecular zone, which consists of the OH bond and the bonds between nearby C, N, O, and H atoms of uracil. The uracil NH bond and the hydrogen bond it formed with the H₂O molecule, N-H···O, store the major portion of the energy released by the relaxing bending mode, thus forming a localized hot band in the intermolecular zone. Energy transfer to the bonds beyond the zone is found to be not significant. The excited uracil NH is found to transfer its energy to the bending mode, thus indicating that the hydrogen bond of N-H···O is the principal energy pathway in both directions. In the presence of efficient near-resonant energy transfer pathways, the time evolution of the centers of mass distance shows the phenomenon of beats. One global and two different local minima energy structures are considered. The results of energy transfer do not differ significantly, suggesting that the two hydrogen bonds in all three structures have similar contributions to the energy transfer.

Ground State Destabilization in Uracil DNA Glycosylase: Let's Not Forget "Tautomeric Strain" in Substrates

Enzymes like uracil DNA glycosylase (UDG) can achieve ground state destabilization, by polarizing substrates to mimic rare tautomers. On the basis of computed nucleus independent chemical shifts, NICS(1)_zz, and harmonic oscillator model of electron delocalization (HOMED) analyses, of quantum mechanics (QM) and quantum mechanics/molecular mechanics (QM/MM) models of the UDG active site, uracil is strongly polarized when bound to UDG and resembles a tautomer >12 kcal/mol higher in energy. Natural resonance theory (NRT) analyses identified a dominant O2 imidate resonance form for residue bound 1-methyl-uracil. This "tautomeric strain" raises the energy of uracil, making uracilate a better than expected leaving group. Computed gas-phase S_N2 reactions of free and hydrogen bonded 1-methyl-uracil demonstrate the relationship between the degree of polarization in uracil and the leaving group ability of uracilate.

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.

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 ᅟ.

Structure and Stability of Chemically Modified DNA Bases: Quantum Chemical Calculations on 16 Isomers of Diphosphocytosine

We studied for the first time 16 tautomers/rotamers of diphosphocytosine by four computational methods. Some of these tautomers/rotamers are isoenergetic although they have different structures. High-level electron correlation MP2 and MP4(SDQ) ab initio methods and density functional methods employing a B3LYP and the new M06-2X functional were used to study the structure and relative stability of 16 tautomers/rotamers of diphosphocytosine. The dienol tautomers of diphosphocytosine are shown to be much more stable than the keto-enol and diketo forms. The tautomers/rotamers stability could be ranked as PC3 = PC12 < PC2 = PC11 < PC1 < PC10 < PC8 < PC9 < PC15 < PC16 < PC6 ~ PC7 < PC13 < PC4 ~ PC14 < PC5. This stability order was discussed in the light of stereo and electronic factors. Solvation effect has been modeled in a high dielectric solvent, water using the polarized continuum model (PCM). Consideration of the solvent causes some reordering of the relative stability of diphosphocytosine tautomers: PC3 ~ PC12 ~ PC2 ~ PC11 < PC1 < PC10 < PC8 < PC9 < PC15 ~ PC16 < PC13 < PC6 ~ PC7 ~ PC14 < PC4 ~ PC5.

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.

Effect of the H-bonding on aromaticity of purine tautomers

Four tautomers of purine (1-H, 3-H, 7-H, and 9-H) and their equilibrium H-bonded complexes with F(-) and HF for acidic and basic centers, respectively, were optimized by means of the B3LYP/6-311++G(d,p) level of theory. Purine tautomer stability increases in the following series: 1-H < 3-H < 7-H < 9-H, consistent with increasing aromaticity. Furthermore, the presence of a hydrogen bond with HF does not change this order. For neutral H-bonded complexes, the strongest and the weakest intermolecular interactions occur (-14.12 and -10.49 kcal/mol) for less stable purine tautomers when the proton acceptor is located in the five- and six-membered rings, respectively. For 9-H and 7-H tautomers the order is reversed. The H-bond energy for the imidazole complex with HF amounts to -14.03 kcal/mol; hence, in the latter case, the fusion of imidazole to pyrimidine decreases its basicity. The ionic H-bonds of N(-)···HF type are stronger by ~10 kcal/mol than the neutral N···HF intermolecular interactions. The hydrogen bond N(-)···HF energies in pyrrole and imidazole are -32.28 and -30.03 kcal/mol, respectively, and are substantially stronger than those observed in purine complexes. The aromaticity of each individual ring and of the whole molecule for all tautomers in ionic complexes is very similar to that observed for the anion of purine. This is not the case for neutral complexes and purine as a reference. The N···HF bonds perturb much more the π-electron structure of five-membered rings than that of the six-membered ones. The H-bonding complexes for 7-H and 9-H tautomers are characterized by higher aromaticity and a much lower range of HOMA variability.

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.

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.

A theoretical study of hydration effects on the prototropic tautomerism of selenouracils

The prototropic tautomerism of 2-, 4-selenouracil and 2,4-diselenouracil has been studied using density functional theory (DFT) methods, at the B3LYP/6-311 + G(3df,2p)//B3LYP/6-31G(d,p) level. The relative stability order of selenouracil tautomers does not resemble that of uracil tautomers, but it is similar to that of thiouracils, even though the energy gaps between the different tautomers of selenouracils are smaller than for thiouracils. The tautomerism activation barriers are high enough as to conclude that only the oxo-selenone or the diselenone structures should be found in the gas phase. The specific interaction with one water molecule reduces these barriers by a half, but still the oxo-selenone form is always the most stable tautomer. The addition of a second water molecule has a relatively small effect, as well as bulk effects, evaluated by means of a continuum-polarized model. For isolated 2- and 4-selenouracils, the more favorable tautomerization process corresponds to a hydrogen transfer towards the selenium atom, the activation barriers for transfer towards the oxygen atom being much higher. This situation changes when specific and bulk effects are included, and the latter process becomes the more favorable one. For 2,4-diselenouracil the more favorable tautomerization, in the gas phase, corresponds to the H shift from N1 to the Se atom at C2, while solvation effects favor the transfer from N3 to the Se atom at C4.

Investigation of proton transport tautomerism in clusters of protonated nucleic acid bases (cytosine, uracil, thymine, and adenine) and ammonia by high-pressure mass spectrometry and ab initio calculations

The energetics of the ion-molecule interactions and structures of the clusters formed between protonated nucleic acid bases (cytosine, uracil, thymine, and adenine) and ammonia have been studied by pulsed ionization high-pressure mass spectrometry (HPMS) and ab initio calculations. For protonated cytosine, uracil, thymine, and adenine with ammonia, the measured enthalpies of association with ammonia are -21.7, -27.9, -22.1, and -17.5 kcal mol-1, respectively. Different isomers of the neutral and protonated nucleic acid bases as well as their clusters with ammonia have been investigated at the B3LYP/6-31+G(d,p) level of theory, and the corresponding binding energetics have also been obtained. The potential energy surfaces for proton transfer and interconversion of the clusters of protonated thymine and uracil with ammonia have been constructed. For cytosine, the experimental binding energy is in agreement with the computed binding energy for the most stable isomer, CN01-01, which is derived from the enol form of protonated cytosine, CH01, and ammonia. Although adenine has a proton affinity similar to that of cytosine, the binding energy of protonated adenine to ammonia is much lower than that for protonated cytosine. This is shown to be due to the differing types of hydrogen bonds being formed. Similarly, although uracil and thymine have similar structures and proton affinities, the binding energies between the protonated species and ammonia are different. Strikingly, the addition of a single methyl group, in going from uracil to thymine, results in a significant structural change for the most stable isomers, UN01-01 and TN03-01, respectively. This then leads to the difference in their measured binding energies with ammonia. Because thymine is found only in DNA while uracil is found in RNA, this provides some potential insight into the difference between uracil and thymine, especially their interactions with other molecules.

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

Density functional theory was used to study the thermodynamics and kinetics for the glycosidic bond cleavage in deoxyuridine. Two reaction pathways were characterized for the unimolecular decomposition in vacuo. However, these processes are associated with large reaction barriers and highly endothermic reaction energies, which is in agreement with experiments that suggest a (water) nucleophile is required for the nonenzymatic glycosidic bond cleavage. Two (S(N)1 and S(N)2) reaction pathways were characterized for direct hydrolysis of the glycosidic bond by a single water molecule; however, both pathways also involve very large barriers. Activation of the water nucleophile via partial proton abstraction steadily decreases the barrier and leads to a more exothermic reaction energy as the proton affinity of the molecule interacting with water increases. Indeed, our data suggests that the barrier heights and reaction energies range from that for hydrolysis by water to that for hydrolysis by the hydroxyl anion, which represents the extreme of (full) water activation (deprotonation). Hydrogen bonds between small molecules (hydrogen fluoride, water, or ammonia) and the nucleobase were found to further decrease the barrier and overall reaction energy but not to the extent that the same hydrogen-bonding interactions increase the acidity of the nucleobase. Our results suggest that the nature of the nucleophile plays a more important role in reducing the barrier to glycosidic bond cleavage than the nature of the small molecule bound, and models with more than one hydrogen fluoride molecule interacting with the nucleobase provide further support for this conclusion. Our results lead to a greater fundamental understanding of the effects of the nucleophile, activation of the nucleophile, and interactions with the nucleobase for this important biological reaction.

Environmental Effects on the Enhancement in Natural and Damaged DNA Nucleobase Acidity Because of Discrete Hydrogen-Bonding Interactions

The present study uses density functional theory to carefully consider the effects of the environment on the enhancement in (natural and damaged) DNA nucleobase acidities because of multiple hydrogen-bonding interactions. Although interactions with one small molecule can increase the acidity of the nucleobases by up to 60 kJ mol-1 in the gas phase, the maximum increase in enzymatic-like environments is expected to be approximately 40 kJ mol-1, which reduces to approximately 30 kJ mol-1 in water. Furthermore, the calculated (simultaneous) effects of two, three, or four molecules are increasingly less than the sum of the individual (additive) effects with an increase in the number and acidity of the small molecules bound or the dielectric constant of the solvent. Regardless of these trends, our calculations reveal that additional hydrogen-bonding interactions will have a significant effect on nucleobase acidity in a variety of environments, where the exact magnitude of the effect depends on the properties of the small molecule bound, the nucleobase binding site, and the solvent. The maximum increase in nucleobase acidity because of interactions with up to four small molecules is approximately 80 kJ mol-1 in enzymatic-like environments (or 65 kJ mol-1 in water). These results suggest that hydrogen-bonding interactions likely play an important role in many biological processes by changing the physical and chemical properties of the nucleobases.

Effects of Hydrogen-Bonding and Stacking Interactions with Amino Acids on the Acidity of Uracil

The effects of hydrogen-bonding interactions with amino acids on the (N1) acidity of uracil are evaluated using (B3LYP) density functional theory. Many different binding arrangements of each amino acid to three uracil binding sites are considered. The effects on the uracil acidity are found to significantly depend upon the nature of the amino acid and the binding orientation, but weakly depend on the binding site. Our results reveal that in some instances small models for the amino acids can be used, while for other amino acids larger models are required to properly describe the binding to uracil. The gas-phase acidity of uracil is found to increase by up to approximately 60 kJ mol(-1) due to discrete hydrogen-bonding interactions. Although (MP2) stacking interactions with aromatic amino acids decrease the acidity of uracil, unexpected increases in the acidity are found when any of the aromatic amino acids, or the backbone, hydrogen bond to uracil. Consideration of enzymatic and aqueous environments leads to decreases in the effects of the amino acids on the acidity of uracil. However, we find that the magnitude of the decrease varies with the nature of the molecule bound, as well as the (gas-phase) binding orientations and strengths, and therefore solvation effects should be considered on a case-by-case basis in future work. Nevertheless, the effects of amino acid interactions within enzymatic environments are as much as approximately 35 kJ mol(-1). The present study has general implications for understanding the nature of active site amino acids in enzymes, such as DNA repair enzymes, that catalyze reactions involving anionic nucleobase intermediates.

Unimolecular Reactivity of Uracil–Cu2+ Complexes in the Gas Phase

The gas-phase interaction of copper(II) ions with uracil are studied by means of mass spectrometry and B3LYP/6-311+G(2df,2p)//B3LYP/6-31G(d) calculations. Positive-ion electrospray spectra show that the reaction of uracil with copper(II) gives rise to singly charged species, whereby the [Cu(uracil--H)](+) complex is the most intense ion in the spectra at low concentration. Mass spectrometry/mass spectrometry (MS/MS) experiments show that the loss of HNCO and NCO are the dominant fragmentation processes, accompanied by a minor loss of CO. A systematic study of the spectra obtained with different labeled species, namely, 2-(13)C- (m/z 175), 2-(13)C,1,3-(15)N(2)- (m/z 177) and 3-(15)N-uracil (m/z 175), concludes unambiguously that both the loss of HNCO and NCO involve exclusively C2 and N3, whereas only C4 is involved in the loss of CO. Suitable mechanisms for these fragmentation processes are proposed through a theoretical survey of the corresponding potential energy surface. In these mechanisms, pi complexes, which lie high in energy with respect to the global minimum, play a significant role in the loss of NCO; this explains why both products, HNCO and NCO involve the same atoms of the ring.

Hydration and stability of nucleic acid bases and base pairs

Empirical, quantum chemical calculations and molecular dynamics simulations of the role of a solvent on tautomerism of nucleic acid bases and structure and properties of nucleic acid base pairs are summarized. Attention was paid to microhydrated (by one and two water molecules) complexes, for which structures found by scanning of empirical potential surfaces were recalculated at a correlated ab initio level. Additionally, isolated as well as mono- and dihydrated H-bonded, T-shaped and stacked structures of all possible nucleic acid base pairs were studied at the same theoretical levels. We demonstrate the strong influence of a solvent on the tautomeric equilibrium between the tautomers of bases and on the spatial arrangement of the bases in a base pair. The results provide clear evidence that the prevalence of either the stacked or hydrogen-bonded structures of the base pairs in the solvent is not determined only by its bulk properties, but rather by specific hydrophilic interactions of the base pair with a small number of solvent molecules.

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.

Microsolvation Effect, Hydrogen-Bonding Pattern, and Electron Affinity of the Uracil−Water Complexes U−(H2O)n (n = 1, 2, 3)

To achieve a systematic understanding of the influence of microsolvation on the electron accepting behaviors of nucleobases, the reliable theoretical method (B3LYP/DZP++) has been applied to a comprehensive conformational investigation on the uracil-water complexes U-(H(2)O)(n) (n = 1, 2, 3) in both neutral and anionic forms. For the neutral complexes, the conformers of hydration on the O2 of uracil are energetically favored. However, hydration on the O4 atom of uracil is more stable for the radical anions. The electron structure analysis for the H-bonding patterns reveal that the CH...OH(2) type H-bond exists only for di- and trihydrated uracil complexes in which a water dimer or trimer is involved. The electron density structure analysis and the atoms-in-molecules (AIM) analysis for U-(H(2)O)(n) suggest a threshold value of the bond critical point (BCP) density to justify the CH...OH(2) type H-bond; that is, CH...OH(2) could be considered to be a H-bond only when its BCP density value is equal to or larger than 0.010 au. The positive adiabatic electron affinity (AEA) and vertical detachment energy (VDE) values for the uracil-water complexes suggest that these hydrated uracil anions are stable. Moreover, the average AEA and VDE of U-(H(2)O)(n) increase as the number of the hydration waters increases.

Non-Aqueous Solvation of n-Octanol and Ethanol: Spectroscopic and Computational Studies

Raman spectroscopy was used to examine the interactions of the free O-H bonds in n-octanol and ethanol with the organic solvents carbon tetrachloride (CCl(4)), cyclohexane, and benzene. These spectra reveal that the solvents CCl(4) and cyclohexane have a small effect on the free O-H peak of alcohols, whereas benzene as a solvent significantly red-shifts the free O-H band. Calculated spectra were generated via MP2/6-31G* calculations and the B3LYP/6-31+G**//MP2/6-31G*-derived Boltzmann populations of each ethanol complex and are consistent with the experimental results. Additional spectra were calculated using Boltzmann populations derived from single-point energies at the polarizable continuum model (PCM) level with the B3LYP/6-31+G** level of theory to take overall solvent effects into account, and these simulated spectra are also largely consistent with the experimental results. Analysis of the computational results reveals a lengthening of the O-H bond from the O-H interaction with the delocalized electronic structure of benzene as well as a bimodal distribution of the free O-H peak of the alcohol/benzene mixtures due to two distinctly different types of alcohol/benzene complexes.

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.

Valence- and Dipole-Bound Anions of the Thymine−Water Complex: Ab Initio Characterization of the Potential Energy Surfaces

The potential energy surfaces of the neutral and anionic thymine-water complexes are investigated using high-level ab initio calculations. Both dipole-bound (DB) and valence-bound (VB) anionic forms are considered. Four minima and three first-order stationary points are located, and binding energies are computed. All minima, for both anions, are found to be vertically and adiabatically stable. The binding energies are much higher for valence-bound than for dipole-bound anions. Adiabatic electron affinities are in the 66-287 meV range for VB anions and the 4-60 meV range for DB anions, and vertical detachment energies are in the 698-977 meV and 10-70 meV range for VB and DB anions, respectively. For cases where literature data are available, the computed values are in good agreement with previous experimental and theoretical studies. It is observed that electron attachment modifies the shape of the potential energy surfaces of the systems, especially for the valence-bound anions. Moreover, for both anions the size of the energy barrier between the two lowest energy minima is strongly reduced, rendering the coexistence of different structures more probable.

The Hydrogen Bonding Properties of Cytosine: A Computational Study of Cytosine Complexed with Hydrogen Fluoride, Water, and Ammonia

Density functional theory is used to study the hydrogen bonding pattern in cytosine, which does not contain alternating proton donor and acceptor sites and therefore is unique compared with the other pyrimidines. Complexes between various small molecules (HF, H(2)O, and NH(3)) and four main binding sites in (neutral and (N1) anionic) cytosine are considered. Two complexes (O2(N1) and N3(N4)) involve neighboring cytosine proton acceptor and donor sites, which leads to cooperative interactions and bidendate hydrogen bonds. The third (less stable) complex (N4) involves a single cytosine donor. The final (O2-N3) complex involves two cytosine proton acceptors, which leads to an anticooperative hydrogen bonding pattern for H(2)O and NH(3). On the neutral surface, the anticooperative O2-N3 complex is less stable than those involving bidentate hydrogen bonds, and the H(2)O complex cannot be characterized when diffuse functions are included in the (6-31G(d,p)) basis set. On the contrary, the anionic O2-N3 structure is the most stable complex, while the HF and H(2)O N3(N4) complexes cannot be characterized with diffuse functions. B3LYP and MP2 potential energy surface scans are used to consider the relationship between the water N3(N4) and O2-N3 complexes. These calculations reveal that diffuse functions reduce the conversion barrier between the two complexes on both the neutral and anionic surfaces, where the reduction leads to a (O2-N3) energy plateau on the neutral surface and complete (N3(N4)) complex destabilization on the anionic surface. From these complexes, the effects of hydrogen bonds on the (N1) acidity of cytosine are determined, and it is found that the trends in the effects of hydrogen bonds on the (N1) acidity are similar for all pyrimidines.

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.