Gas-phase acidity of D-glucose. A density functional theory study

The gas-phase acidity of D-glucopyranose was studied by means of B3LYP calculations combined with 6-31G(d,p) or 6-31+G(d,p) standard basis sets. For each anomer, deprotonation of the various primary and secondary hydroxyl groups was considered. As in solution, the anomeric hydroxyl is found to be the most acidic for both anomers, but only when the 6-31+G(d,p) basis set is used for geometry optimization. Deprotonation of the anomeric hydroxyl induces an important C(1)--O endocyclic bond elongation and subsequently promotes an energetically favored ring-opening process as attested by the very small calculated activation barriers. The results also suggest that interconversion between the various deprotonated alpha- and beta-anomers may easily occur under slightly energetic conditions. B3LYP/6-311+G(2df,2p) calculations led to the an absolute gas-phase acidity of deltaacidGo(298)(alpha-D-glucose) = 1398 kJ mol(-1). This estimate matches well the only experimental value available to date. Finally, this study again confirms that the use of diffuse functions on heavy atoms is necessary to describe anionic systems properly and to achieve good relative and absolute gas-phase acidities.     

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.     

Gas-phase hydrogen/deuterium exchange and conformations of deprotonated flavonoids and gas-phase acidities of flavonoids

Gas-phase hydrogen/deuterium (H/D) exchange was used to probe the conformations, gas-phase acidities, and sites of deprotonation of isomeric flavonoid aglycons and glycosides. The flavonoids in each isomer series were differentiated on the basis of their relative rate constants and total numbers of exchanges. For example, flavonoids that possess neohesperidose-type disaccharides may undergo faster and far more extensive exchange than isomeric rutinoside flavonoids. The structural factors that promote or prevent H/D exchange were identified and correlated with collisionally activated dissociation (CAD) patterns and/or molecular modeling data (both high-level ab initio acidity calculations and conformational analysis with molecular dynamics (MD) simulations), thus providing a framework for the use of H/D exchange reactions in the structural elucidation of new flavonoids.     

Insight into the acidic behavior of oxazolidin-2-one,its thione and selone analogs through computational techniques

Deprotonation thermochemistry of Oxazolidin-2-one (OXA), Oxazolidine-2-thione (OXA-S), and Oxazolidine-2-selone (OXA-Se) has been studied in order to find the most acidic site and relative acidities of these heterocyclics at various sites. The deprotonation enthalpies at MP2/6-311++G**//MP2/6-31+G* and B3LYP/6-31+G* levels, while the free energies for deprotonation process and pKa values at B3LYP/6-31+G* level both in gas and aqueous phase (using PCM continuum model) of the anions of the three heterocyclics have been computed at 298 K. Calculated aqueous phase pKa values of OXA vary by ~6–7 units from the experimental aqueous phase pKa values of OXA and its derivatives. The deprotonation at the nitrogen is favored in OXA over the carbon atoms in contrast to the OXA-S and OXA-Se where in the deprotonation at the carbon attached to the nitrogen is most preferred. Deprotonation at this carbon induces an important C–O bond rupture in OXA-S and OXA-Se promoting an energetically favored ring-opening process. The finding offers a rare case when C–H acidity is able to dominate over the N–H acidity. In order to explain the relative stabilities, relative acidities and deprotonation enthalpies various characteristics of these molecules as well as their anions such as molecular electrostatic potential surface (MEP), frontier molecular orbital (FMO) features, chemical hardness, softness have been governed. The three dimensional MEP maps and HOMO–LUMO orbitals encompassing these molecules yield a reliable relative stability and reactivity (in terms of acidity) map displaying the most probable regions for deprotonation. The differential distribution of the electrostatic potential over the neutral and anionic species of OXA, OXA-S, and OXA-Se molecules is authentically reflected by HOMO–LUMO orbitals and NBO charge distribution analysis. The lone pair occupancies, second order delocalization energies for orbital interactions and the distribution of atomic charges over the entire molecular framework as obtained from natural bond orbital (NBO) analysis are found to faithfully replicate the predictions from the MEP maps and HOMO–LUMO band gaps in respect of explaining the relative stabilities and acidities in most of the cases. Good linear correlations have been obtained between HOMO–LUMO gap and pKa values in the aqueous phase for OXA and OXA-S molecules.     

Determination of the gas-phase acidities of cysteine-polyalanine peptides using the extended kinetic method

We determined the gas-phase acidities of two cysteine-polyalanine peptides, HSCA3 and HSCA4, using a triple-quadrupole mass spectrometer through application of the extended kinetic method with full entropy analysis. Five halogenated carboxylic acids were used as the reference acids. The negatively charged proton-bound dimers of the deprotonated peptides with the conjugate bases of the reference acids were generated by electrospray ionization. Collision-induced dissociation (CID) experiments were carried out at three collision energies. The enthalpies of deprotonation (Delta(acid)H) of the peptides were derived according to the linear relationship between the logarithms of the CID product ion branching ratios and the differences of the gas-phase acidities. The values were determined to be Delta(acid)H(HSCA3) = 317.3 +/- 2.4 kcal/mol and Delta(acid)H (HSCA4) = 316.2 +/- 3.9 kcal/mol. Large entropy effects (Delta(DeltaS) = 13-16 cal/mol K) were observed for these systems. Combining the enthalpies of deprotonation with the entropy term yielded the apparent gas-phase acidities (Delta(acid)G(app)) of 322.1 +/- 2.4 kcal/mol (HSCA3) and 320.1 +/- 3.9 kcal/mol (HSCA4), in agreement with the results obtained from the CID-bracketing experiments. Compared with that in the isolated cysteine residue, the thiol group in HSCA3,4 has a stronger gas-phase acidity by about 20 kcal/mol. This increased acidity is likely due to the stabilization of the negatively charged thiolate group through internal solvation.     

Origin of the acidity enhancement of formic acid over methanol: resonance versus inductive effects

Density functional theory calculations were employed to study the relative contribution of resonance versus inductive effects toward the 37 kcal/mol enhanced gas-phase acidity (DeltaH degrees (acid)) of formic acid (1) over methanol (2). The gas-phase acidities of formic acid, methanol, vinyl alcohol (5), and their vinylogues (6, 8, and 9) were calculated at the B3LYP/6-31+G level of theory. Additionally, acidities were calculated for the formic acid and vinyl alcohol vinylogues in which the formyl group and the vinyl group, respectively, were perpendicular to the rest of the conjugated system. Comparisons among these calculated acidities suggest that inductive effects are the predominant effects responsible for the enhanced acidity of formic acid over methanol, accounting for between roughly 62% and 65% of the total enhanced acidity; the remaining 38% to 35% of the acidity enhancement appears to be due to resonance effects. Further comparisons suggest that resonance effects are between roughly 58% and 65% of the 26 kcal/mol calculated acidity enhancement of vinyl alcohol over methanol, and the remaining 42% to 35% are due to inductive effects.     

Vinyl hydrogen acidities of two stereoisomers

The gas-phase acidities of the vinyl hydrogens of cis- and trans-2-butene were measured by the silane kinetic method in a Fourier-transform ion cyclotron resonance spectrometer. The acidities of ethene and the secondary vinyl hydrogen of propene were measured by the same method. The method was calibrated using the known acidities of methane and benzene. The vinyl hydrogens of trans-2-butene are more acidic than the vinyl hydrogens of cis-2-butene by 4.5 kcal/mol; the acidities of ethene and the secondary vinyl hydrogen of propene are between those of the two butenes. The acidity of cis-2-butene is 409 +/- 2 kcal/mol, and the acidity of trans-2-butene is 405 +/- 2 kcal/mol. Density functional theory calculations are in good agreement with the experiments. The results are discussed in terms of steric interactions, polarizabilities, dipole-dipole interactions, and charge-dipole interactions.     

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.     

Interpretation of Brønsted acidity by triadic paradigm: a G3 study of mineral acids

Deprotonation enthalpies and the gas-phase acidities of 24 inorganic acids are calculated by using composite G3 and G2 methodologies. The computed values are in very good accordance with available measured data. It is found that the experimental DeltaH(acid) values of the FSO(3)H and CF(3)SO(3)H are too high by some 6 and 7 kcal mol(-1), respectively. Furthermore, a new DeltaH(acid) value for HClO(4) of 300 kcal mol(-1) is recommended and suggested as a threshold of superacidicity in the gas phase. The calculated deprotonation enthalpies are interpreted by employing the trichotomy paradigm. Taking into account that the deprotonation enthalpy is a measure of acidity, it can be safely stated that the pronounced acidities of mineral acids are to a very large extent determined by Koopmans' term with very few exceptions, one of them being H(2)S. To put it in another way, acidities are predominantly a consequence of the ability of the conjugate bases to accommodate the excess electron charge, since Koopmans' term in trichotomy analysis is related to conjugate base anion. The final state is decisive in particular for superacids like ClSO(3)H, CF(3)SO(3)H, HClO(4), HBF(4), HPF(6), HAlCl(4), and HAlBr(4). However, in the latter two molecules the bond dissociation energy of the halogen-H bond substantially contributes to their high acidity too. Therefore, acidity of these two most powerful superacids studied here is determined by cooperative influence of both initial and final state effects. It should be emphasized that acidity of hydrogen halides HCl and HBr is a result of concerted action of all three terms included in triadic analysis. A byproduct of the triadic analysis are the first adiabatic ionization energies of the anionic conjugate bases. They are in fair to good agreement with the experimental data, which are unfortunately sparse. A fairly good qualitative correlation is found between the gas-phase deprotonation enthalpies of six mineral O-H acids and available Hammett-Taft sigma(p)- constants of the corresponding substituent groups.     

7-Ethyl-2,3,5,6,8-pentahydroxy-1,4-naphthoquinone (echinochrome A): a DFT study of the antioxidant mechanism 2.* The structure of monosodium salts of echinochrome A and their reactions with the hydroperoxyl radical

The relative gas-phase acidities of all OH groups of 7-ethyl-2,3,5,6,8-pentahydroxy-1,4-naphthoquinone (echinochrome A, 1) were evaluated by the B3LYP/6-311G(d) and B3LYP/6-311G(d,p) methods. Calculations predict that gB-OH groups at the C(2) and C(6) atoms are the most acidic in molecule 1and their acidity is higher than that of o-nitrophenol. Conformational analysis of undissociated monosodium salts of 1and their radicals was performed. It was shown that gas-phase quenching reactions of the hydroperoxyl radical by mono-sodium salts of 1are exothermic.     

Regioselectivity of pyridine deprotonation in the gas phase

The regioselective deprotonation of pyridine in the gas phase has been investigated by using chemical reactivity studies. The mixture of regioisomers, trapped as carboxylates, formed in an equilibrium mixture is determined to result from 70-80% deprotonation in the 4-position, and 20-30% deprotonation at the 3-position. The ion formed by deprotonation in the 2-position is not measurably deprotonated at equilibrium because the ion is destabilized by lone-pair repulsion. From the composition of the mixture, the gas-phase acidities (DeltaH degrees acid) at the 4-, 3-, and 2-positions are determined to be 389.9 +/- 2.0, 391.2-391.5, and >391.5 kcal/mol, respectively. The relative acidities of the 4- and 3-positions are explained by using Hammett-Taft parameters, derived by using the measured gas-phase acidities of pyridine carboxylic acids. The values of sigmaF and sigmaR are -0.18 and 0.74, respectively, showing the infused nitrogen in pyridine to have a strong pi electron-withdrawing effect, but with little sigma-inductive effect.     

Gas-phase protonation and deprotonation of acrylonitrile derivatives N[triple chemical bond]C--CH==CH--X (X=CH(3), NH(2), PH(2), SiH(3))

A combined experimental and theoretical study on the gas-phase basicity and acidity of a series of cyanovinyl derivatives is presented. The gas-phase basicities and acidities of (N[triple chemical bond]C--CH==CH--X, X=CH(3), NH(2)) were obtained by means of Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometry techniques. The corresponding calculated values were obtained at the G3B3 level of theory. The effects of exchanging CH(3) for SiH(3), and NH(2) for PH(2), were analyzed at the same level of theory. For the neutral molecules, the Z isomer is always the dominant species under standard gas-phase conditions at 298 K. The loss of the proton from the substituent X was found systematically to be much more favorable than deprotonation of the HC==CH linking group. The corresponding isomeric E ion is much more stable than the Z ion, so that only the former should be found in the gas phase. The most significant structural changes upon deprotonation occur for the methyl and amino derivatives because, in both cases, deprotonation of X leads to a significant charge delocalization in the corresponding anion. Protonation takes place systematically at the cyano group, whereby the isomeric E ion is again more stable than the Z ion. Push-pull effects explain the preference of aminoacrylonitrile to be protonated at the cyano group, which also explains the high basicity of this derivative relative to other members of the analyzed series that present rather similar gas-phase basicities, GB approximately 780 kJ mol(-1), indicating that the different nature of the substituents has only a weak effect on the intrinsic basicity of the cyano group. The cyanovinyl derivatives have a significantly stronger gas-phase acidity than that of the corresponding vinyl compounds CH(2)==CH--X. This acidity-strengthening effect of the cyano group is attributed to the greater stabilization of the anion with respect to the corresponding neutral compound.     

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

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.     

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.     

Theoretical Study of Gas-Phase Acidities of Selected Angiotensin-Converting Enzyme Inhibitors

The structure and gas-phase acidities of six members of large angiotensin-converting enzyme (ACE) inhibitor family [cilazaprilat, silanediol, fosinoprilat, AcSDKP, angiotensin_I (terminal part), and RXP_407] have been studied using the ONIOM Becke3LYP/6-311+G(d,p):HF/3-21G^∗ method. The investigated ACE inhibitors are weak acids with calculated acidity of about 1270–1650 kJ mol⁻¹. Of acids studied the highest gas-phase acidity (1273 kJ mol⁻¹) possesses experimental ACE inhibitor RXP_407. This drug, according to the computed pK_a value (3.2), is also in water solution the most acidic compound of the ACE inhibitors investigated.     

Towards highly powerful neutral organic superacids—a DFT study of some polycyano derivatives of planar hydrocarbons

Density functional theory (DFT) calculations at B3LYP/6-311+G(2d,p)//B3LYP/6-31G(d) level have been carried out on indene, cyclopentaphenanthrene and 1H-phenalene and their heptacyano and nonacyano derivatives, respectively, in order to examine their acidities in the gas-phase and DMSO. It is found that polycyano derivatives represent powerful organic superacids, the most acidic being nonacyano-1H-phenalene. The origin of the highly pronounced acidity is identified as a strong anionic resonance in the resulting conjugate base. Comparison of the calculated ΔH_acid value for 1H-phenalene with the experimental NIST value shows that the latter is too large by 8-11 kcal mol⁻¹. A possible reason for this error is briefly discussed.     

A comparison of the gas phase acidities of phospholipid headgroups: Experimental and computational studies

Proton-bound dimers consisting of two glycerophospholipids with different headgroups were prepared using negative ion electrospray ionization and dissociated in a triple quadrupole mass spectrometer. Analysis of the tandem mass spectra of the dimers using the kinetic method provides, for the first time, an order of acidity for the phospholipid classes in the gas phase of PE < PA < PG < PS < PI. Hybrid density functional calculations on model phospholipids were used to predict the absolute deprotonation enthalpies of the phospholipid classes from isodesmic proton transfer reactions with phosphoric acid. The computational data largely support the experimental acidity trend, with the exception of the relative acidity ranking of the two most acidic phospholipid species. Possible causes of the discrepancy between experiment and theory are discussed and the experimental trend is recommended. The sequence of gas phase acidities for the phospholipid headgroups is found to (1) have little correlation with the relative ionization efficiencies of the phospholipid classes observed in the negative ion electrospray process, and (2) correlate well with fragmentation trends observed upon collisional activation of phospholipid [M - H](-) anions.     

The acidity of chloro-substituted benzenes: a comparison of gas phase, ab initio, and kinetic data

The deprotonation energies of benzene, chlorobenzene, all di-, tri-, tetrachlorobenzenes, and pentachlorobenzene have been determined in the gas phase using a Fourier transform ion cyclotron resonance mass spectrometer. The values measured differ only slightly, though significantly, from the corresponding data for oligofluorobenzenes. The heavier halogen acidifies orthopositions slightly less and meta-positions slightly more than fluorine does. Moreover, the contributions of three or more chloro substituents are not perfectly additive. In fact the accumulation attenuates the contributions somewhat. Quantum chemical calculations at the MP2/6-311+G* level reproduce the gas-phase acidities fairly well, but reveal special effects when extended to experimentally not observable benzenides carrying the halogens at anion-remote positions. Competition experiments have been performed to assess the relative reactivity of nine oligochlorobenzenes towards sec-butyllithium in tetrahydrofuran at -100 degrees C. An almost exact linear correlation between logarithmic rates and gas-phase acidities has been found.