Acidity of adenine and adenine derivatives and biological implications. A computational and experimental gas-phase study

The gas-phase acidities of adenine, 9-ethyladenine, and 3-methyladenine have been investigated for the first time, using computational and experimental methods to provide an understanding of the intrinsic reactivity of adenine. Adenine is found to have two acidic sites, with the N9 site being 19 kcal mol(-1) more acidic than the N10 site; the bracketed acidities are 333 +/- 2 and 352 +/- 4 kcal mol(-1), respectively. Because measurement of the less acidic site can be problematic, we benchmarked the adenine N10 measurement by bracketing the acidity of 9-ethyladenine, which has the N9 site blocked and allows for exclusive measurement of the N10 site. The acidity of 9-ethyladenine brackets to 352 +/- 4 kcal mol(-1), comparable to that of the N10 site of the parent adenine. Calculations and experiments with 3-methyladenine, a harmful mutagenic nucleobase, uncovered the surprising result that the most commonly written tautomer of 3-methyladenine is not the most stable in the gas phase. We have found that the most stable tautomer is the "N10 tautomer" 10, as opposed to the imine tautomer 3. The bracketed acidity of 10 is 347 +/- 4 kcal mol(-1). Since 10 is not a viable species in DNA, 3 is a likely tautomer; calculations indicate that this form has an extremely high acidity (320-323 kcal mol(-1)). The biological implications of these results, particularly with respect to enzymes that cleave alkylated bases from DNA, are discussed.     

The acidity of uracil and uracil analogs in the gas phase: four surprisingly acidic sites and biological implications

The gas phase acidities of a series of uracil derivatives (1-methyluracil, 3-methyluracil, 6-methyluracil, 5,6-dimethyluracil, and 1,3-dimethyluracil) have been bracketed to provide an understanding of the intrinsic reactivity of uracil. The experiments indicate that in the gas phase, uracil has four sites more acidic than water. Among the uracil analogs, the N1-H sites have deltaH(acid) values of 331-333 kcal mol(-1); the acidity of the N3 sites fall between 347-352 kcal mol(-1). The vinylic C6 in 1-methyluracil and 3-methyluracil brackets to 363 kcal mol(-1), and 369 kcal mol(-1) in 1,3-dimethyluracil; the C5 of 1,3-dimethyluracil brackets to 384 kcal mol(-1). Calculations conducted at B3LYP/6-31+G* are in agreement with the experimental values. The bracketing of several of these sites involved utilization of an FTMS protocol to measure the less acidic site in a molecule that has more than one acidic site, establishing the generality of this method. In molecules with multiple acidic sites, only the two most acidic sites were bracketable, which is attributable to a kinetic effect. The measured acidities are in direct contrast to in solution, where the two most acidic sites of uracil (N1 and N3) are indifferentiable. The vinylic C6 site is also particularly acidic, compared to acrolein and pyridine. The biological implications of these results, particularly with respect to enzymes for which uracil is a substrate, are discussed.     

Selectivity of purine alkylation by a quinone methide. Kinetic or thermodynamic control?

The alkylation reaction of 9-methyladenine and 9-methylguanine (as prototype substrates of deoxy-adenosine and -guanosine), by the parent o-quinone methide (o-QM), has been investigated in the gas phase and in aqueous solution, using density functional theory at the B3LYP/6-311+G(d,p) level. The effect of the medium on the reactivity, and on the stability of the resulting adducts, has been investigated by using the C-PCM solvation model to assess which adduct arises from the kinetically favorable path, or from an equilibrating process. The calculations indicate that the most nucleophilic site of the methyl-substituted nucleobases in the gas phase is the guanine oxygen atom (O(6)) (DeltaG()(gas) = 5.6 kcal mol(-)(1)), followed by the adenine N1 (DeltaG)(gas) = 10.3 kcal mol(-)(1)), while other centers exhibit a substantially lower nucleophilicity. The bulk effect of water as a solvent is the dramatic reduction of the nucleophilicity of both 9-methyladenine N1 (DeltaG)(solv) = 14.5 kcal mol(-)(1)) and 9-methylguanine O(6) (DeltaG)(solv) = 17.0 kcal mol(-)(1)). As a result there is a reversal of the nucleophilicity order of the purine bases. While O(6) and N7 nucleophilic centers of 9-methylguanine compete almost on the same footing, the reactivity gap between N1 and N7 of 9-methyladenine in solution is highly reduced. Regarding product stability, calculations predict that only two of the adducts of o-QM with 9-methyladenine, those at NH(2) and N1 positions, are lower in energy than reactants, both in the gas phase and in water. However, the adduct at N1 can easily dissociate in water. The adducts arising from the covalent modification of 9-methylguanine are largely more stable than reactants in the gas phase, but their stability is markedly reduced in water. In particular, the oxygen alkylation adduct becomes slightly unstable in water (DeltaG(solv) = +1.4 kcal mol(-)(1)), and the N7 alkylation product remains only moderately more stable than free reactants (DeltaG(solv) = -2.8 kcal mol(-)(1)). Our data show that site alkylations at the adenine N1 and the guanine O(6) and N7 in water are the result of kinetically controlled processes and that the selective modification of the exo-amino groups of guanine N2 and adenine N6 are generated by thermodynamic equilibrations. The ability of o-QM to form several metastable adducts with purine nucleobases (at guanine N7 and O(2), and adenine N1) in water suggests that the above adducts may act as o-QM carriers.     

The acidity and proton affinity of the damaged base 1,N6-ethenoadenine in the gas phase versus in solution: intrinsic reactivity and biological implications

1,N(6)-ethenoadenine (epsilonA) is a highly mutagenic lesion that is excised from human DNA by the enzyme alkyladenine DNA glycosylase (AAG). In an effort to understand the intrinsic properties of 1,N(6)-ethenoadenine, we examined its gas phase acidity and proton affinity using quantum mechanical calculations and mass spectrometric experimental methods. We measure two acidities for epsilonA: a more acidic site (DeltaH(acid) = 332 kcal mol(-1); DeltaG(acid) = 325 kcal mol(-1)) and a less acidic site (DeltaH(acid) = 367 kcal mol(-1); DeltaG(acid) = 360 kcal mol(-1)). We also find that the proton affinity of the most basic site of 1,N(6)-ethenoadenine is 232-233 kcal mol(-1) (GB = 224 kcal mol(-1)). These measurements, when compared to calculations, establish that, under our experimental conditions, we have only the canonical tautomer of 1,N(6)-ethenoadenine present. We also compare the gas phase acidic properties of epsilonA with that of the normal bases adenine and guanine and find that epsilonA is the most acidic. This supports the theory that AAG and other related enzymes may cleave damaged bases as anions. Furthermore, comparison of the gas phase and aqueous acidities indicates that the nonpolar environment of the enzyme enhances the acidity differences of epsilonA versus adenine and guanine.     

Benzocyclobutene: The Impact of Fusing a Strained Ring onto Benzene

The gas-phase acidities of the two aromatic sites in benzocyclobutene were measured in a Fourier transform mass spectrometer using a kinetic technique (i.e., the DePuy method). Fusion of a cyclobutane ring onto benzene is found to have a slight acidifying effect at the alpha-position (3.2 +/- 1.7 kcal mol(-)(1)) and little, if any, influence on the beta-site (0.8 +/- 1.9 kcal mol(-)(1)). Energetic data (DeltaH degrees (acid) = 386.2 +/- 3.0 kcal mol(-)(1), EA = 0.84 +/- 0.11 eV, and C-H BDE = 92 +/- 4 kcal mol(-)(1)) for the benzylic position were obtained via the bracketing technique and application of a thermodynamic cycle. Differences in the reactivities of the three conjugate bases also were explored. Ab initio and density functional theory calculations were carried out to provide geometries, energies, and insights into the carbanions' electronic structures.     

Equilibrium acidities and homolytic bond dissociation energies of acidic C-H bonds in alpha-arylacetophenones and related compounds

The equilibrium acidities (pK(AH)s) and the oxidation potentials of the congugate anions [E(ox)(A(-))s] were determined in dimethyl sulfoxide (DMSO) for eight ketones of the structure GCOCH(3) and 20 of the structure RCOCH(2)G, (where R = alkyl, phenyl and G = alkyl, aryl). The homolytic bond dissociation energies (BDEs) for the acidic C-H bonds of the ketones were estimated using the equation BDE(AH) = 1.37pK(AH) + 23.1E(ox)(A(-)) + 73.3. While the equilibrium acidities of GCOCH(3) were found to be dependent on the remote substituent G, the BDE values for the C-H bonds remained essentially invariant (93.5 +/- 0.5 kcal/mol). A linear correlation between pK(AH) values and [E(ox)(A(-))s] was found for the ketones. For RCOCH(2)G ketones, both pK(AH) and BDE values for the adjacent C-H bonds are sensitive to the nature of the substituent G. However, the steric bulk of the aryl group tends to exert a leveling effect on BDEs. The BDE of alpha-9-anthracenylacetophenone is higher than that of alpha-2-anthracenylacetophenone by 3 kcal/mol, reflecting significant steric inhibition of resonance in the 9-substituted system. A range of 80.7-84.4 kcal/mol is observed for RCOCH(2)G ketones. The results are discussed in terms of solvation, steric, and resonance effects. Ab initio density functional theory (DFT) calculations are employed to illustrate the effect of steric interactions on radical and anion geometries. The DFT results parallel the trends in the experimental BDEs of alpha-arylacetophenones.     

Association patterns of platinated purine nucleobases in metal-modified pairs and triples

Blocking of Watson-Crick or Hoogsteen edges in purine nucleobases by a metal entity precludes involvement of these sites in interbase hydrogen bonding, thereby leaving the respective other edge or the sugar edge as potential H bonding sites. In mixed guanine, adenine complexes of trans-a2PtII (a = NH3 or CH3NH2) of composition trans-[(NH3)2Pt(9-EtA-N1)(9-MeGH-N7)](NO3)2 (1a), trans-[(NH3)2Pt(9-EtA-N1)(9-MeGH-N7)](ClO4)2 (1b), and trans,trans-[(CH3NH2)2(9-MeGH-N7)Pt(N1-9-MeA-N7)Pt(9-MeGH-N7)(CH3NH2)2](ClO4)4*2H2O (2) (with 9-EtA = 9-ethyladenine, 9-MeA= 9-methyladenine, 9-MeGH = 9-methylguanine), this aspect is studied. Thus, in 1b pairing of two adenine ligands via Hoogsteen edges and in 2 pairing of two guanine bases via sugar edges is realized. These situations are compared with those found in a series of related complexes.     

Theoretical study on the gas-phase acidity of multiple sites of Cu+-adenine and Cu2+-adenine complexes

The acidities of multiple sites in Cu(+)-adenine and Cu(2+)-adenine complexes have been investigated theoretically. To compare, the acidities of adenine (A) and adenine radical cation (A(*+)) have also been included. The results clearly indicate that the acidities of C-H and N-H groups in Cu(+/2+)-adenine are significantly enhanced relative to the neutral adenine. The acidic order for a given site on adenine and adenine derivatives is as follows: Cu(2+)-adenine > A(*+) > Cu(+)-adenine > A. For Cu(+)-adenine and Cu(2+)-adenine, N3-coordination exhibits N9-H acid, and N1- and N7-coordination exhibits N6-H(a) and N6-H(b) acid, respectively. Additionally, it is found that C2-H group is surprisingly acidic in the coordination complexes. Calculations in aqueous solution reveal that our results can be extrapolated to aqueous solution. Analyses of the electronic properties interpret the highest acidity of Cu(2+)-adenine among the adenine derivatives studied. Also, Electrostatic potential calculations of [A(-H(+))](-) and [A(-H(+))](*) indicate that the removal of H(a) or H(b) from the amino group favors the bidentate coordination, which provides a dative bond from the deprotonated N and the original coordination ligand to copper ion besides the electrostatic interaction between them and thereby stabilizes the [A(-H(+))](-)/[A(-H(+))](*). NBO analysis confirms the electrostatic potential result.     

Equilibrium Acidities and Homolytic Bond Dissociation Enthalpies of the Acidic C-H Bonds in P-(Para-substituted benzyl)triphenylphosphonium Cations and Related Cations

Equilibrium acidities (pK(HA)) of six P-(para-substituted benzyl)triphenylphosphonium (p-GC(6)H(4)CH(2)PPh(3)(+)) cations, P-allyltriphenylphosphonium cation, P-cinnamyltriphenylphosphonium cation, and As-(p-cyanobenzyl)triphenylarsonium cation, together with the oxidation potentials [E(ox)(A(-))] of their conjugate anions (ylides) have been measured in dimethyl sulfoxide (DMSO) solution. The acidifying effects of the alpha-triphenylphosphonium groups on the acidic C-H bonds in toluene and propene were found to be ca 25 pK(HA) units (34 kcal/mol). Introduction of an electron-withdrawing group such as 4-NO(2), 4-CN, or 4-Br into the para position of the benzyl ring in p-GC(6)H(4)CH(2)PPh(3)(+) cations resulted in an additional acidity increase, but introduction of the 4-OEt electron-donating group decreases the acidity. The equilibrium acidities of p-GC(6)H(4)CH(2)PPh(3)(+) cations were nicely linearly correlated with the Hammett sigma(-) constants of the substituents (G) with a slope of 4.78 pK(HA) units (R(2) = 0.992) (Figure 1). Reversible oxidation potentials of the P-(para-substituted benzyl)triphenylphosphonium ylides were obtained by fast scan cyclic voltammetry. The homolytic bond dissociation enthalpies (BDEs) of the acidic C-H bonds in these cations, estimated by combining their equilibrium acidities with the oxidation potentials of their corresponding conjugate anions, showed that the alpha-Ph(3)P(+) groups have negligible stabilizing or destabilizing effects on the adjacent radicals. The equilibrium acidity of As-(p-cyanobenzyl)triphenylarsonium cation is 4 pK(HA) units weaker than that of P-(p-cyanobenzyl)triphenylphosphonium cation, but the BDE of the acidic C-H bond in As-(p-cyanobenzyl)triphenylarsonium cation is ca 2 kcal/mol higher than that in P-(p-cyanobenzyl)triphenylphosphonium cation.     

Gas-phase thermochemical properties of pyrimidine nucleobases

The gas-phase acidity and proton affinity of thymine, cytosine, and 1-methyl cytosine have been examined using both theoretical (B3LYP/6-31+G*) and experimental (bracketing, Cooks kinetic) methods. This paper represents a comprehensive examination of multiple acidic sites of thymine and cytosine and of the acidity and proton affinity of thymine, cytosine, and 1-methyl cytosine. Thymine exists as the most stable "canonical" tautomer in the gas phase, with a DeltaH(acid) of 335 +/- 4 kcal mol(-1) (DeltaG(acid) = 328 +/- 4 kcal mol(-1)) for the more acidic N1-H. The acidity of the less acidic N3-H site has not, heretofore, been measured; we bracket a DeltaH(acid) value of 346 +/- 3 kcal mol(-1) (DeltaG(acid) = 339 +/- 3 kcal mol(-1)). The proton affinity (PA = DeltaH) of thymine is measured to be 211 +/- 3 kcal mol(-1) (GB = DeltaG = 203 +/- 3 kcal mol(-1)). Cytosine is known to have several stable tautomers in the gas phase in contrast to in solution, where the canonical tautomer predominates. Using bracketing methods in an FTMS, we measure a DeltaH(acid) for the more acidic site of 342 +/- 3 kcal mol(-1) (DeltaG(acid) = 335 +/- 3 kcal mol(-1)). The DeltaH(acid) of the less acidic site, previously unknown, is 352 +/- 4 kcal mol(-1) (345 +/- 4 kcal mol(-1)). The proton affinity is 228 +/- 3 kcal mol(-1) (GB = 220 +/- 3 kcal mol(-1)). Comparison of these values to calculations indicates that we most likely have a mixture of the canonical tautomer and two enol tautomers and possibly an imine tautomer under our conditions in the gas phase. We also measure the acidity and proton affinity of cytosine using the extended Cooks kinetic method. We form the proton-bound dimers via electrospray of an aqueous solution, which favors cytosine in the canonical form. The acidity of cytosine using this method is DeltaH(acid) = 343 +/- 3 kcal mol(-1), PA = 227 +/- 3 kcal mol(-1). We also examined 1-methyl cytosine, which has fewer accessible tautomers than cytosine. We measure a DeltaH(acid) of 349 +/- 3 kcal mol(-1) (DeltaG(acid) = 342 +/- 3 kcal mol(-1)) and a PA of 230 +/- 3 kcal mol(-1) (GB = 223 +/- 3 kcal mol(-1)). Our ultimate goal is to understand the intrinsic reactivity of nucleobases; gas-phase acidic and basic properties are of interest for chemical reasons and also possibly for biological purposes because biological media can be quite nonpolar.     

Valence anions in complexes of adenine and 9-methyladenine with formic acid: stabilization by intermolecular proton transfer

Photoelectron spectra of adenine-formic acid (AFA(-)) and 9-methyladenine-formic acid (MAFA(-)) anionic complexes have been recorded with 2.540 eV photons. These spectra reveal broad features with maxima at 1.5-1.4 eV that indicate formation of stable valence anions in the gas phase. The neutral and anionic complexes of adenine/9-methyladenine and formic acid were also studied computationally at the B3LYP, second-order M?ller-Plesset, and coupled-cluster levels of theory with the 6-31++G** and aug-cc-pVDZ basis sets. The neutral complexes form cyclic hydrogen bonds, and the most stable dimers are bound by 17.7 and 16.0 kcal/mol for AFA and MAFA, respectively. The theoretical results indicate that the excess electron in both AFA(-) and MAFA(-) occupies a pi* orbital localized on adenine/9-methyladenine, and the adiabatic stability of the most stable anions amounts to 0.67 and 0.54 eV for AFA(-) and MAFA(-), respectively. The attachment of the excess electron to the complexes induces a barrier-free proton transfer (BFPT) from the carboxylic group of formic acid to a N atom of adenine or 9-methyladenine. As a result, the most stable structures of the anionic complexes can be characterized as neutral radicals of hydrogenated adenine (9-methyladenine) solvated by a deprotonated formic acid. The BFPT to the N atoms of adenine may be biologically relevant because some of these sites are not involved in the Watson-Crick pairing scheme and are easily accessible in the cellular environment. We suggest that valence anions of purines might be as important as those of pyrimidines in the process of DNA damage by low-energy electrons.     

Acidity and proton affinity of hypoxanthine in the gas phase versus in solution: intrinsic reactivity and biological implications

Hypoxanthine is a mutagenic purine base that most commonly arises from the oxidative deamination of adenine. Damaged bases such as hypoxanthine are associated with carcinogenesis and cell death. This inevitable damage is counteracted by glycosylase enzymes, which cleave damaged bases from DNA. Alkyladenine DNA glycosylase (AAG) is the enzyme responsible for excising hypoxanthine from DNA in humans. In an effort to understand the intrinsic properties of hypoxanthine, we examined the gas-phase acidity and proton affinity using quantum mechanical calculations and gas-phase mass spectrometric experimental methods. In this work, we establish that the most acidic site of hypoxanthine has a gas-phase acidity of 332 +/- 2 kcal mol-1, which is more acidic than hydrochloric acid. We also bracket a less acidic site of hypoxanthine at 368 +/- 3 kcal mol-1. We measure the proton affinity of the most basic site of hypoxanthine to be 222 +/- 3 kcal mol-1. DFT calculations of these values are consistent with the experimental data. We also use calculations to compare the acidic and basic properties of hypoxanthine with those of the normal bases adenine and guanine. We find that the N9-H of hypoxanthine is more acidic than that of adenine and guanine, pointing to a way that AAG could discriminate damaged bases from normal bases. We hypothesize that AAG may cleave certain damaged nucleobases as anions and that the active site may take advantage of a nonpolar environment to favor deprotonated hypoxanthine as a leaving group versus deprotonated adenine or guanine. We also show that an alternate mechanism involving preprotonation of hypoxanthine is energetically less attractive, because the proton affinity of hypoxanthine is less than that of adenine and guanine. Last, we compare the acidity in the gas phase versus that in solution and find that a nonpolar environment enhances the differences in acidity among hypoxanthine, adenine, and guanine.     

Minor groove site coordination of adenine by platinum group metal ions: effects on basicity,base pairing,and electronic structure

Dithioether- or diamine-tethered adenine derivatives react with Pt(II), Pd(II), and Rh(III) ions to give N3-coordinated complexes of the types [MCl(SSN)](+) (M = Pt or Pd), [RhCl(3)(SSN)], or [RhCl(3)(NNN)] (where SSN = 1-(N9-adenine)-3,6-dithia-heptane or 1-(N9-adenine)-4,7-dithia-octane; NNN = ethylenediamine-N,9-ethyladenine). Single-crystal X-ray analysis confirms the nature of the metal-nucleobase interaction and highlights a conserved intermolecular hydrogen-bonding motif for all the complexes, irrespective of the metal-ion geometry. Coordination significantly reduces the basicity of the adeninyl group, as indicated by a pK(a) value of -0.16 for [PtCl(N3-1-(N9-adenine)-3,6-dithia-heptane)]BF(4), compared to a pK(a) value of 4.2 for 9-ethyladenine. The site of proton binding, N1 or N7, could not be unambiguously assigned from the (1)H NMR data, because of the similar effect on the chemical shifts of the H2 and H8 protons. Density functional calculations at the BP-LACVP level suggest N1 as the site of protonation for this type of complex. This is in contrast to the N7-protonation reported for [Pt(dien)(N3-6,6',9-trimethyladenine)](2+), as reported elsewhere (Meiser et al., Chem.-Eur. J. 1997, 3, 388). However, further electronic structure calculations in the gas phase reveal that the preferred site for protonation for N3-bound complexes is conformationally dependent. N3 coordination was also found to reduce the extent of base pairing between adenine and thymine in dimethylsulfoxide for the self-complementary complex [PtCl(L3)](+) (L3 = 1-(N9-adenine)-3,6-dithia-9-(N1-thymine)nonane), compared to that for the uncomplexed ligand.     

Conformation-dependent reaction thermochemistry: study of lactones and lactone enolates in the gas phase

Gas-phase acidities (Delta H degrees (acid)) of lactones with ring sizes from four to seven have been measured on a Fourier transform ion cyclotron resonance mass spectrometer. Electron affinities (EAs) of the corresponding lactone enolate radicals were measured on a continuous-wave ion cyclotron resonance mass spectrometer, and the bond dissociation energies (BDEs) of the alpha C-H bonds were derived. In order of increasing ring size, Delta H degrees (acid) = 368.7 +/- 2., 369.4 +/- 2.2, 367.3 +/- 2.2, and 368.3 +/- 2.2 kcal/mol and BDE = 99.4 +/- 2.3, 94.8 +/- 2.3, 89.2 +/- 2.3, and 92.8 +/- 2.4 kcal/mol for beta-propiolactone, gamma-butyrolactone, delta-valerolactone, and epsilon-caprolactone, respectively. For their corresponding enolate radicals, EA = 44.1 +/- 0.3, 38.8 +/- 0.3, 35.3 +/- 0.3, and 37.9 +/- 0.6 kcal/mol. All of these lactones are considerably more acidic than methyl acetate, consistent with a dipole repulsion model. Both BDEs and EAs show a strong dependence on ring size, whereas Delta H degrees (acid) does not. These findings are discussed, taking into account differential electronic effects and differential strain between the reactant and product species in each reaction.     

Thermochemical studies of N-methylpyrazole and N-methylimidazole

The 351.1 nm photoelectron spectra of the N-methyl-5-pyrazolide anion and the N-methyl-5-imidazolide anion are reported. The photoelectron spectra of both isomers display extended vibrational progressions in the X2A' ground states of the corresponding radicals that are well reproduced by Franck-Condon simulations, based on the results of B3LYP/6-311++G(d,p) calculations. The electron affinities of the N-methyl-5-pyrazolyl radical and the N-methyl-5-imidazolyl radical are 2.054 +/- 0.006 eV and 1.987 +/- 0.008 eV, respectively. Broad vibronic features of the A(2)A' ' states are also observed in the spectra. The gas-phase acidities of N-methylpyrazole and N-methylimidazole are determined from measurements of proton-transfer rate constants using a flowing afterglow-selected ion flow tube instrument. The acidity of N-methylpyrazole is measured to be Delta(acid)G(298) = 376.9 +/- 0.7 kcal mol(-1) and Delta(acid)H(298) = 384.0 +/- 0.7 kcal mol(-1), whereas the acidity of N-methylimidazole is determined to be Delta(acid)G(298) = 380.2 +/- 1.0 kcal mol(-1) and Delta(acid)H(298)= 388.1 +/- 1.0 kcal mol(-1). The gas-phase acidities are combined with the electron affinities in a negative ion thermochemical cycle to determine the C5-H bond dissociation energies, D(0)(C5-H, N-methylpyrazole) = 116.4 +/- 0.7 kcal mol(-1) and D(0)(C5-H, N-methylimidazole) = 119.0 +/- 1.0 kcal mol(-1). The bond strengths reported here are consistent with previously reported bond strengths of pyrazole and imidazole; however, the error bars are significantly reduced.     

Hydroxyl radical reactions with adenine: reactant complexes, transition states, and product complexes

In order to address problems such as aging, cell death, and cancer, it is important to understand the mechanisms behind reactions causing DNA damage. One specific reaction implicated in DNA oxidative damage is hydroxyl free-radical attack on adenine (A) and other nucleic acid bases. The adenine reaction has been studied experimentally, but there are few theoretical results. In the present study, adenine dehydrogenation at various sites, and the potential-energy surfaces for these reactions, are investigated theoretically. Four reactant complexes [A···OH]* have been found, with binding energies relative to A+OH* of 32.8, 11.4, 10.7, and 10.1 kcal mol(-1). These four reactant complexes lead to six transition states, which in turn lie +4.3, -5.4, (-3.7 and +0.8), and (-2.3 and +0.8) kcal mol(-1) below A+OH*, respectively. Thus the lowest lying [A···OH]* complex faces the highest local barrier to formation of the product (A-H)*+H(2)O. Between the transition states and the products lie six product complexes. Adopting the same order as the reactant complexes, the product complexes [(A-H)···H(2)O]* lie at -10.9, -22.4, (-24.2 and -18.7), and (-20.5 and -17.5) kcal mol(-1), respectively, again relative to separated A+OH*. All six A+OH* → (A-H)*+H(2)O pathways are exothermic, by -0.3, -14.7, (-17.4 and -7.8), and (-13.7 and -7.8) kcal mol(-1), respectively. The transition state for dehydrogenation at N(6) lies at the lowest energy (-5.4 kcal mol(-1) relative to A+OH*), and thus reaction is likely to occur at this site. This theoretical prediction dovetails with the observed high reactivity of OH radicals with the NH(2) group of aromatic amines. However, the high barrier (37.1 kcal mol(-1)) for reaction at the C(8) site makes C(8) dehydrogenation unlikely. This last result is consistent with experimental observation of the imidazole ring opening upon OH radical addition to C(8). In addition, TD-DFT computed electronic transitions of the N(6) product around 420 nm confirm that this is the most likely site for hydrogen abstraction by hydroxyl radical.     

Ab initio prediction of the gas- and solution-phase acidities of strong Brønsted acids: the calculation of pKa values less than -10

The intrinsic gas-phase acidities of a series of 21 Br?nsted acids have been predicted with G3(MP2) theory. The G3(MP2) results agree with high level CCSD(T)/CBS acidities for H(2)SO(4), FSO(3)H, CH(3)SO(3)H, and CF(3)SO(3)H to within 1 kcal/mol. The G3(MP2) results are in excellent agreement with experimental gas-phase acidities in the range 342-302 kcal/mol to within <1 kcal/mol for 14 out of 15 acids. Five of the six acids in the range of 302-289 kcal/mol had an average deviation of 5.5 kcal/mol and the strongest acid, (CF(3)SO(2))(3)CH, deviated by 15.0 kcal/mol. These high-level calculations strongly suggest that the experimental acidities in this very acidic part of the scale need to be remeasured. The CCSD(T)/CBS (mixed exponential Gaussian) additive approach for CH(3)CO(2)H, HNO(3), H(2)SO(4), CH(3)SO(3)H, FSO(3)H, and CF(3)SO(3)H gives excellent agreement (+/-1 kcal/mol) with experiment for the DeltaH(f)(0)'s of non-sulfur containing species, and supports the low end of the experimental values for H(2)SO(4) and FSO(3)H. Use of a larger basis set (aug-cc-pV5Z) in the CBS extrapolation improves the agreement with experiment for both H(2)SO(4) and FSO(3)H. The G3(MP2) heats of formation for RSO(3)H molecules tend to be underestimated as compared to the CCSD(T)/CBS approach by 2.5-7.0 kcal/mol. COSMO solvation calculations were used to predict solution free energies and pK(a) values with pK(a)'s up to -17.4. Including the solvation of the proton gives good agreement with experimental pK(a) values in the very acidic regime, whereas it is less reliable for weaker acids. The use of CH(3)CO(2)H and HNO(3) as reference acids in the less acidic and more acidic regions of the scale, respectively, provided improved results to within +/-2 pK(a) units in nearly all cases (+/-3 kcal/mol accuracy).     

Electron attachment to the DNA bases adenine and guanine and dehydrogenation of their anionic derivatives: Density functional study

Density functional theory (DFT) calculations have been used to explore electron attachment to the purines adenine and guanine and their hydrogen atom loss. Calculations show that the dehydrogenation at the N9 site in the adenine and guanine transient anions is the lowest‐cost channel of hydrogen loss, and the N9? H bond scission has Gibbs free energies of dissociation ΔG° of 8.8 kcal mol^?1 for the anionic adenine and 13.9 kcal mol^?1 for the anionic guanine. The relatively high feasibility of low‐energy electron (LEE)‐induced N9? H bond cleavage in the purine nucleobases arises from high electron affinities of their H‐deleted counterparts. Unlike adenine, other N? H bond dissociations are competitive with the N9? H bond fission in the anionic guanine. The replacement of hydrogen in the ring of purine has a significant effect on the N9? H bond fragmentation. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Polar group enhanced gas-phase acidities of carboxylic acids: an investigation of intramolecular electrostatic interaction

We studied the effects of polar groups on the gas-phase acidities of carboxylic acids experimentally and computationally. In this connection, the gas-phase acidities (DeltaH(acid), the enthalpy of deprotonation, and DeltaG(acid), the deprotonation free energy) of borane-complexed methylaminoacetic acid ((CH(3))2N(BH(3))CH(2)CO(2)H) and methylthioacetic acid (CH(3)S(BH(3))CH(2)CO(2)H) were measured using the kinetic method in a flowing afterglow-triple quadrupole mass spectrometer. The values of DeltaH(acid) and DeltaG(acid) of (CH(3))2N(BH(3))CH(2)CO(2)H were determined to be 328.8 +/- 1.9 and 322.1 +/- 1.9 kcal/mol, and those of CH(3)S(BH(3))CH(2)CO(2)H were determined to be 325.8 +/- 1.9 and 319.2 +/- 1.9 kcal/mol, respectively. The theoretical enthalpies of deprotonation of (CH(3))2N(BH(3))CH(2)CO(2)H (329.2 kcal/mol) and CH(3)S(BH(3))CH(2)CO(2)H (325.5 kcal/mol) were calculated at the B3LYP/6-31+G(d) level of theory. The calculated enthalpies of deprotonation of N-oxide-acetic acid (CH(3)NOCH(2)CO(2)H, 329.4 kcal/mol) and S-oxide-acetic acid (CH(3)SOCH(2)CO(2)H, 328.6 kcal/mol) are comparable to the experimental results for borane-complexed methylamino- and methylthioacetic acids. The enthalpy of deprotonation of sulfone-acetic acid (CH(3)SO2CH(2)CO(2)H, 326.1 kcal/mol) is about 2 kcal/mol lower than the S-oxide-acetic acid. The calculated enthalpy of deprotonation of sulfoniumacetic acid, (CH(3))2S+CH(2)CO(2)H, is 243.0 kcal/mol. Compared to the corresponding reference molecules, CH(3)NHCH(2)CO(2)H and CH(3)SCH(2)CO(2)H, the dipolar group and the monopolar group substituted carboxylic acids are stronger acids by 11-14 and 97 kcal/mol, respectively. We correlated the changes of the acidity upon a polar group substitution to the electrostatic free energy within the carboxylate anion. The acidity enhancements in polar group substituted carboxylic acids are the results of the favorable electrostatic interactions between the polar group and the developing charge at the carboxyl group.     

Gas phase SN2 reactions of halide ions with trifluoromethyl halides: front- and back-side attack vs. complex formation

Density functional theory computations and pulsed-ionization high-pressure mass spectrometry experiments have been used to explore the potential energy surfaces for gas-phase S(N)2 reactions between halide ions and trifluoromethyl halides, X(-) + CF(3)Y --> Y(-) + CF(3)X. Structures of neutrals, ion-molecule complexes, and transition states show the possibility of two mechanisms: back- and front-side attack. From pulsed-ionization high-pressure mass spectrometry, enthalpy and entropy changes for the equilibrium clustering reactions for the formation of Cl(-)(BrCF(3)) (-16.5 +/- 0.2 kcal mol(-1) and -24.5 +/- 1 cal mol(-1) K(-1)), Cl(-)(ICF(3)) (-23.6 +/- 0.2 kcal mol(-1)), and Br(-)(BrCF(3)) (-13.9 +/- 0.2 kcal mol(-1) and -22.2 +/- 1 cal mol(-1) K(-1)) have been determined. These are in good to excellent agreement with computations at the B3LYP/6-311+G(3df)//B3LYP/6-311+G(d) level of theory. It is shown that complex formation takes place by a front-side attack complex, while the lowest energy S(N)2 reaction proceeds through a back-side attack transition state. This latter mechanism involves a potential energy profile which closely resembles a condensed phase S(N)2 reaction energy profile. It is also shown that the Cl(-) + CF(3)Br --> Br(-) + CF(3)Cl S(N)2 reaction can be interpreted using Marcus theory, in which case the reaction is described as being initiated by electron transfer. A potential energy surface at the B3LYP/6-311+G(d) level of theory confirms that the F(-) + CF(3)Br --> Br(-) + CF(4) S(N)2 reaction proceeds through a Walden inversion transition state.