Lactone enols are stable in the gas phase but highly unstable in solution

2-Hydroxyoxol-2-ene (C(5)-1), the enol tautomer of gamma-butyrolactone, was generated in the gas phase as the first representative of the hitherto elusive class of lactone enols and shown by neutralization-reionization mass spectrometry to be remarkably stable as an isolated species. Ab initio calculations by QCISD(T)/6-311+G(3df,2p) provided the enthalpies of formation, proton affinities, and gas-phase basicities for gaseous lactone enols with four- (C(4)-1), five- (C(5)-1), and six-membered rings (C(6)-1). The acid-base properties of C(4)-C(6) lactones and enols and reference carboxylic acid enols CH(2)=C(OH)(2) (3) and CH(2)=C(OH)OCH(3) (4) were also calculated in aqueous solution. The C(4)-C(6) lactone enols show gas-phase proton affinities in the range of 933-944 kJ mol(-)(1) and acidities in the range of 1401-1458 kJ mol(-)(1). In aqueous solution, the lactone enols are 15-20 orders of magnitude more acidic than the corresponding lactones, with enol pK(a) values increasing from 5.6 (C(4)-1) to 14.5 (C(6)-1). Lactone enols are moderately weak bases in water with pK(BH) in the range of 3.9-8.1, whereas the lactones are extremely weak bases of pK(BH) in the range of -10.5 to -17.4. The acid-base properties of lactone enols point to their high reactivity in protic solvents and explain why no lactone enols have been detected thus far in solution studies.     

Ab initio calculations and mass spectrometric determination of the gas-phase proton affinities of 4,4'-disubstituted 2,2'-bipyridines

The gas-phase proton affinities of 4,4'-di(R)-2,2'-bipyridines (R: H, Br, Cl, NO(2), Me) were determined by mass spectrometric measurements and by ab initio calculations at the HF/6-31G and MP2/6-31G levels of theory. The energy barriers for rotation about the central C-C bond were also studied computationally. Two minima were found for both unprotonated and protonated species, the global minima being at the trans planar and cis planar conformations, respectively. Local minima for the unprotonated compounds were at the cis nonplanar conformation and for the protonated compounds at the trans nonplanar. Two different proton affinity values were calculated for each compound by employing different conformations for the protonated species. The computational values were in good agreement with the experimental proton affinities. Substituents affect the proton affinity according to their ability to withdraw or to donate electrons, halogen and nitro-substituted bipyridines having a lower proton affinity and methyl-substituted bipyridine having a higher proton affinity than 2,2'-bipyridine itself.     

First-principles prediction of the pK(a)s of anti-inflammatory oxicams

The gas- and aqueous-phase acidities of a series of oxicams have been computed by combining M05-2X/6-311+G(3df,2p) gas-phase free energies with solvation free energies from the CPCM-UAKS, COSMO-RS, and SMD solvent models. To facilitate accurate gas-phase calculations, a benchmarking study was further carried out to assess the performance of various density functional theory methods against the high-level composite method G3MP2(+). Oxicams are typically diprotic acids, and several tautomers are possible in each protonation state. The direct thermodynamic cycle and the proton exchange scheme have been employed to compute the microscopic pK(a)s on both solution- and gas-phase equilibrium conformers, and these were combined to yield the macroscopic pK(a) values. Using the direct cycle of pK(a) calculation, the CPCM-UAKS model delivered reasonably accurate results with MAD ~ 1, whereas the SMD and COSMO-RS models' performance was less satisfactory with MAD ~ 3. Comparison with experiment also indicates that direct cycle calculations based on solution conformers generally deliver better accuracy. The proton exchange cycle affords further improvement for all solvent models through systematic error cancellation and therefore provides better reliability for the pK(a) prediction of compounds of these types. The latter approach has been applied to predict the pK(a)s of several recently synthesized oxicam derivatives.     

Thermodynamic cycle for the calculation of ab initio pKa values for hydroxamic acids

A thermodynamic cycle to calculate pK_a values (Minus log of acid dissociation constants) of hydroxamic acids is presented. Hydroxamic acids exist mainly as amide isomers in the aqueous medium. The amide form of hydroxamic acids has two deprotonation sites and may yield either an N-ion or an O-ion upon deprotonation. The thermodynamic cycle proposed includes the gas-phase N–H deprotonation of the hydroxamic acid, the solvent phase transformation of the N-ion to the O-ion and the solvation of the hydroxamic acid molecule and the O-ion in water. The CBS-QB3 method was employed to obtain gas-phase free energy differences between 12 hydroxamic acids and their respective anions. The aqueous solvation Gibbs free energy changes were calculated at the HF/6-31G(d)/CPCM and HF/6-31+G(d)/CPCM levels of theory using HF/6-31+G(d)/CPCM geometries. For the proton, literature values of the gas-phase free energy of formation and the solvation free energy change were used. The free energy change for the transformation of the N-ion to O-ion in the aqueous medium was calculated by employing CBS-QB3/CPCM in the aqueous medium. For this, the hydroxamic acids were divided in two classes according to the substituent at the carbonyl carbon. A common transformation free energy difference for aliphatic substituted hydroxamic acids and a separate common transformation free energy difference for aromatic substituted hydroxamic acids were obtained. The pK_a calculation yielded a root mean square error of 0.32 pK_a units.     

Ab initio procedure for aqueous-phase pKa calculation: the acidity of nitrous acid

We present an ab initio procedure for accurately calculating aqueous-phase pKa values and apply it to study the acidity of nitrous acid (HNO2, or HONO). The aqueous-phase pK(a) of nitrous acid was obtained from calculated gas-phase acidities and solvation free energies via a thermodynamic cycle and the solvation model chemistry of Barone et al. (J. Chem. Phys. 1997, 107, 3210). Solvation free energies were calculated at the HF/6-31G(d) level using the dielectric-polarizable continuum and the integral equation formalism-polarizable continuum solvent models (D-PCM and IEF-PCM, respectively), with the D-PCM model yielding the most accurate pKa values. For HF free energies of solvation, significant improvements in accuracy could be made by moving to the larger 6-311++G(3df,3pd) and aug-cc-pVQZ basis sets. Solvation free energies were also calculated using the density functional theory (DFT) methods B3LYP, TPSS, PBE0, B1B95, VSXC, B98 and O3LYP, with the most accurate methods being TPSS and VSXC, which provided average errors of less than 0.11 pKa units. Solvation free energies calculated with the different DFT methods were relatively insensitive to the basis set used. Our theoretical calculations are compared with experimental results obtained using stopped flow spectrophotometry. The pKa of nitrous acid was measured as 3.16 at 25 degrees C, and the enthalpy and entropy of nitrous acid dissociation were calculated from measurements as 6.7 kJ mol(-1) and -38.4 J mol(-1) K(-1), respectively, between 25 and 45 degrees C. The UV/visible absorption spectra of the nitrite ion and nitrous acid were also examined, and molar extinction coefficients were obtained for each.     

Vibrational study of new Pt(II) and Pd(II) complexes with functionalized nitrogen-containing tertiary phosphine oxides. Ab initio study of ortho-, meta- and para-dimethylphosphinylmethyleneoxyaniline

Vibrational study of new Pt(II) and Pd(II) complexes of functionalized nitrogen-containing tertiary phosphine oxides, namely ortho-, meta- and para-dimethylphosphinylmethyleneoxyaniline (o-, m- and p-dpmoa), (CH3)2P(O)CH2OC6H4NH2, have been presented. Geometry optimization of the ligands was performed at HF/6-31G* and B3LYP/6-31G* levels of the theory. Harmonic frequencies were calculated at HF/6-31G* optimized geometries. Relative gas-phase and solution-phase (H2O and CH3CN) basicities of o-, m- and p-dpmoa ligands have been determined by ab initio calculations at STO-3G level with the Onsager reaction field model. On the basis of the vibrational study, physical and analytical data it was suggested that the ligands in the complexes studied coordinate through the amino group and form square-planar platinum and palladium complexes of the general formula ML2Cl2 (M = Pt, Pd, L = o-, m- and p-dpmoa).     

Determination of the overlapping pK(a) values of chrysin using UV-vis spectroscopy and ab initio methods

The overlapping pK(a) values of 5,7-dihydroxyflavone (chrysin) in EtOH-water solutions were determined by means of a UV-vis spectroscopic method that uses absorbance diagrams, at constant ionic strength (0.050 M) and temperature (25.0+/-0.1 degrees C). It was observed that the pK(a) values increase when the polarity-polarizability and solvation abilities of the reaction medium decrease. In order to calculate the pK(a1) and pK(a2) of chrysin in pure water, various relationships between the determined pK(a) and properties of solvents (relative permittivity, alpha-parameter of Taft and parameter Acity), are proposed. Moreover, with the aim of explaining the first pK(a1) value obtained, the molecular conformations and solute-solvent interactions of the 7(O(-))chrysinate monoanion were also investigated, using ab initio methods. Several ionization reactions and equilibria in water, which possesses a high hydrogen-bond-donor ability, are proposed. These reactions and equilibria constituted the necessary theoretical basis to calculate the first acidity constant of chrysin. The HF/6-31G(d) and HF/6-31+G(d) methods were used for calculations. Tomasi's method was used to analyze the formation of intermolecular hydrogen bonds between the 7(O(-))chrysinate monoanion and water molecules. It was proposed that in alkaline aqueous solutions the monoanion of chrysin is solvated with one water molecule. The agreement between the experimental and theoretical pK(a1) values provides good support for the acid-base reactions proposed in this paper.     

Structural and energetic aspects of the protonation of phenol,catechol, resorcinol,and hydroquinone

The various protonated forms of phenol (1), catechol (2), resorcinol (3), and hydroquinone (4) were explored by ab initio quantum chemical calculations at the MP2/6-31G(d) and B3LYP/6-31G(d) levels. Proton affinities (PA) of 1-4 were calculated by the combined G2(MP2,SVP) method, and their gas-phase basicities were estimated after calculation of the change in entropy on protonation. These theoretical data were compared with the corresponding experimental values determined in a high-pressure mass spectrometer. This comparison confirmed that phenols are essentially carbon bases and that protonation generally occurs in a position para to the hydroxyl group. Resorcinol is the most effective base (PA = 856 kJ mol-1) due to the participation of both oxygen atoms in the stabilization of the protonated form. Since protonation is accompanied by a freezing of the two internal rotations, a significant decrease in entropy is observed. The basicity of catechol (PA = 823 kJ mol-1) is due to the existence of an intramolecular hydrogen bond, which is strengthened upon protonation. The lower basicity of hydroquinone (PA = 808 kJ mol-1) is a consequence of the fact that protonation necessarily occurs in a position ortho to the hydroxyl group. When the previously published data are reconsidered and a corrected protonation entropy is used, a proton affinity value of 820 kJ mol-1 is obtained for phenol.     

Gas-phase basicity of glycine,alanine, proline,serine, lysine,histidine and some of their peptides by the thermokinetic method

The thermokinetic method is applied to a set of six amino acids (glycine, alanine, proline, serine, lysine, histidine) and 30 of their di- and tri-peptides for which experimental proton transfer rate constants were available. The comparison between the presently determined gas-phase basicities, GBs, of the amino acids with values obtained from equilibrium constant determination is generally good (a mean deviation of appoximately 3 kJ mol(-1) is observed). Derived proton affinities values are discussed. The gas-phase basicities of peptides provided by the present study correct several previously estimated values thus offering a more firm basis for structural discussion. Composite reaction efficiency curves indicate the existence, for several peptides, of at least two non-interconverting populations of protonated forms.     

Accurate pK(a) calculations for carboxylic acids using complete basis set and Gaussian-n models combined with CPCM continuum solvation methods

Complete Basis Set and Gaussian-n methods were combined with CPCM continuum solvation methods to calculate pK(a) values for six carboxylic acids. An experimental value of -264.61 kcal/mol for the free energy of solvation of H(+), DeltaG(s)(H(+)), was combined with a value for G(gas)(H(+)) of -6.28 kcal/mol to calculate pK(a) values with Cycle 1. The Complete Basis Set gas-phase methods used to calculate gas-phase free energies are very accurate, with mean unsigned errors of 0.3 kcal/mol and standard deviations of 0.4 kcal/mol. The CPCM solvation calculations used to calculate condensed-phase free energies are slightly less accurate than the gas-phase models, and the best method has a mean unsigned error and standard deviation of 0.4 and 0.5 kcal/mol, respectively. The use of Cycle 1 and the Complete Basis Set models combined with the CPCM solvation methods yielded pK(a) values accurate to less than half a pK(a) unit.     

Experimental and theoretical evidence of basic site preference in polyfunctional superbasic amidinazine: N(1),N(1)-dimethyl-N(2)-beta-(2-pyridylethyl)formamidine

The gas-phase basicity (GB) of the flexible polyfunctional N(1),N(1)-dimethyl-N(2)-beta-(2-pyridylethyl)formamidine (1) containing two potential basic sites (the ring N-aza and the chain N-imino) is obtained from proton-transfer equilibrium constant measurements, using Fourier-transform ion-cyclotron resonance mass spectrometry. Comparison of the experimental GB obtained for 1 with those reported for model amidines and azines indicates that the chain N-imino in the amidine group is the favored site of protonation. Semiempirical (AM1) and ab initio calculations (HF, MP2, and DFT), performed for 1 and its protonated forms, confirm this interpretation. These results are in contrast to those found previously for N(1),N(1)-dimethyl-N(2)-azinylformamidines (containing the amidine function directly linked to the azinyl ring), in which the ring N-aza is the most basic site in the gas phase. The separation of the two potential basic sites in 1 by the ethylene chain interrupts the resonance conjugation between the two functions and changes their relative basicities and, thus, the preferable site of protonation. It also increases the chelation effect against the proton and the gas-phase basicity of 1 in such a magnitude that consequently 1 may be classified as a superbase (GB = 241.1 kcal mol(-)(1)). A transition state corresponding to the internal transfer of the proton (ITP) between the ring N-aza and the chain N-imino in 1 is investigated at the DFT(B3LYP)/6-31G level. The energy barrier calculated for the ITP between the two basic sites is small and vanishes when zero-point vibrational terms and thermal corrections are applied to obtain the enthalpy or Gibbs energy of activation for the proton transfer. Additional calculations at the DFT(MPW1K)/6-31G level confirm this behavior. This indicates that the quantum-chemical ITP in 1 has a single-well character. The proton is located on the N-imino site, and the H-bond is formed with the N-aza site.     

Protonation of glycine and alanine: proton affinities, intrinsic basicities and proton transfer path

Proton affinities and intrinsic basicities for nitrogen and oxygen protonation in the gas phase of the amino acids glycine and alanine were calculated using density functional theory (DFT) and ab initio methods at different levels of theory from Hartree-Fock (HF) to G2 approximations. All methods gave good agreement for proton affinities for nitrogen protonation for both amino acids. However, dramatic differences were found between DFT, MP4//MP2, and G2 results on one hand, and MP4//HF results on the other to the calculation of structural and energetic characteristics of oxygen protonation in glycine and alanine. An investigation into the source of these differences revealed that electron correlation effects are chiefly responsible for the differences in calculated oxygen proton affinities between the various methods. It has been found that proton transfer between nitrogen and oxygen protonation sites in both amino acids occurs without a transfer path barrier when correlated methods were used to calculate the path energetics.     

Hydride affinity scale of various substituted arylcarbeniums in acetonitrile

Combined with the integral equation formalism polarized continuum model (IEFPCM), the hydride affinities of 96 various acylcarbenium ions in the gas phase and CH(3)CN were estimated by using the B3LYP/6-31+G(d)//B3LYP/6-31+G(d), B3LYP/6-311++G(2df,2p)//B3LYP/6-31+G(d), and BLYP/6-311++G(2df,2p)//B3LYP/6-31+G(d) methods for the first time. The results show that the combination of the BLYP/6-311++G(2df,2p)//B3LYP/6-31+G(d) method and IEFPCM could successfully predict the hydride affinities of arylcarbeniums in MeCN with a precision of about 3 kcal/mol. On the basis of the calculated results from the BLYP method, it can be found that the hydride affinity scale of the 96 arylcarbeniums in MeCN ranges from -130.76 kcal/mol for NO(2)-PhCH(+)-CN to -63.02 kcal/mol for p-(Me)(2)N-PhCH(+)-N(Me)(2), suggesting most of the arylcarbeniums are good hydride acceptors. Examination of the effect of the number of phenyl rings attached to the carbeniums on the hydride affinities shows that the increase of the hydride affinities takes place linearly with increasing number of benzene rings in the arylcarbeniums. Analyzing the effect of the substituents on the hydride affinities of arylcarbeniums indicates that electron-donating groups decrease the hydride affinities and electron-withdrawing groups show the opposite effect. The hydride affinities of arylcarbeniums are linearly dependent on the sum of the Hammett substituent parameters σ(p)(+). Inspection of the correlation of the solution-phase hydride affinities with gas-phase hydride affinities and aqueous-phase pK(R)(+) values reveals a remarkably good correspondence of ΔG(H(-)A)(R(+)) with both the gas-phase relative hydride affinities only if the α substituents X have no large electron-donating or -withdrawing properties and the pK(R)(+) values even though the media are dramatically different. The solution-phase hydride affinities also have a linear relationship with the electrophilicity parameter E, and this dependence can certainly serve as one of the most effective ways to estimate the new E values from ΔG(H(-)A)(R(+)) or vice versa. Combining the hydride affinities and the reduction potentials of the arylcarbeniums, we obtained the bond homolytic dissociation Gibbs free energy changes of the C-H bonds in the corresponding hydride adducts in acetonitrile, ΔG(HD)(RH), and found that the effects of the substituent on ΔG(HD)(RH) are very small. Simple thermodynamic analytic platforms for the three C-H cleavage modes were constructed. It is evident that the present work would be helpful in understanding the nature of the stabilities of the carbeniums and mechanisms of the hydride transfers between carbeniums and other hydride donors.     

水溶液中铵离子pKa值的理论计算

在HF/6-31+G*和B3LYP/6-31+G*水平上, 采用导体极化连续模型(CPCM)及UAKS孔穴计算了11种铵离子在水溶剂中的溶剂化自由能, 与实验值相比较, 平均误差和标准偏差分别为0.17, 12.04和0.96, 10.96 kJ/mol. 结合B3LYP/6-31+G*水平上的11种铵离子气相质子转移反应自由能, 得到了水溶剂中的绝对pKa值, 计算结果与实验数据吻合得很好, 相应的平均误差和标准偏差分别为0.05, 1.50和0.45, 1.40 pKa单位. 可见, 采用CPCM-UAKS模型能够较为精确地计算铵离子型化合物的绝对pKa值.     

Protonation thermochemistry of selected hydroxy- and methoxycarbonyl molecules

The gas-phase basicities of a representative set of hydroxy- and methoxycarbonyl compounds (hydroxyacetone, 1, 3-hydroxybutanone, 2, 3-hydroxy-3-methylbutanone, 3, 1-hydroxy-2-butanone, 4, 4-hydroxy-2-butanone, 5, 5-hydroxy-2-pentanone, 6, methoxyacetone, 7, 3-methoxy-2-butanone, 8, 4-methoxy-2-butanone, 9, and 5-methoxy-2-pentanone, 10) were experimentally determined by the equilibrium method using Fourier transform ion cyclotron resonance and high-pressure mass spectrometry techniques. The latter method allows the measurement of proton transfer equilibrium constants at various temperatures and thus the estimate of both the proton affinities and the protonation entropies of the relevant species. Quantum chemical calculations at the G3 and the B3LYP/6-311+G(3df,2p)//6-31G(d) levels of theory were undertaken in order to find the most stable structures of the neutrals 1-10 and their protonated forms. Conformational and vibrational analyses have been done with the aim of obtaining a theoretical estimate of the protonation entropies.     

Vibrational spectra and ab initio analysis of tert-butyl, trimethylsilyl, and trimethylgermyl derivatives of 3,3-dimethyl cyclopropene V. 3,3-dimethyl-1-(trimethylgermyl)cyclopropene

3,3-dimethyl-1-(trimethylgermyl)cyclopropene (I) was synthesised using a standard procedure. The IR and Raman spectra of I in the liquid phase were measured. The molecular geometry of I was optimised completely at the HF/6-31G* level. The HF/6-31G*//HF/6-31G* force field was calculated and scaled using the set of scale factors transferred from those determined previously for scaling the theoretical force fields of 3,3-dimethylbutene-1 and 1-methyl-, 1,2-dimethyl-, and 3,3-dimethylcyclopropene. The assignments of the observed vibrational bands were performed using the theoretical frequencies calculated from the scaled HF/6-31G*//HF/6-31G* force field and the ab initio values of the IR intensities, Raman cross-sections and depolarisation ratios. The theoretical spectra are given. The completely optimised structural parameters of I and its vibrational frequencies are compared with corresponding data of related molecules.     

Ab initio calculations of proton affinities of glycine, proline, cysteine and phenylalanine: comparison with the experimental values obtained using an electrospray ionisation ion trap mass spectrometer

The proton affinities of four characteristic amino acids, glycine, proline, cysteine and phenylalanine, have been calculated using three different types of wave function, HF, DFT and MP2. These wave functions were combined with several basis sets using the valence double- or triple-xi with polarisation functions, and included or not a diffuse d function on heavy atoms. Calculations were carried out using Gaussian 98 on a LINUX system (2 GHz, 2 GB of RAM). The calculated results have been compared with the experimental values obtained using Cooks' kinetic method, in particular on an electrospray ionisation ion trap mass spectrometer. An excellent agreement was found between the experimental values and the theoretical results obtained using the B3P86/6-31+G*//B3LYP/6-31G* level.     

Protonation studies of modified adenine and adenine nucleotides by theoretical calculations and (15)N NMR

The acid/base character of nucleobases affects phenomena such as self-association, interaction with metal ions, molecular recognition by proteins, and nucleic acid base-pairing. Therefore, the investigation of proton-transfer equilibria of natural and synthetic nucleos(t)ides is of great importance to obtain a deeper understanding of these phenomena. For this purpose, a set of ATP prototypes was investigated using (15)N NMR spectroscopy, and the corresponding adenine bases were investigated by theoretical calculations. (15)N NMR measurements provided not only acidity constants but also information on the protonation site(s) on the adenine ring and regarding the ratio of the singly protonated species in equilibrium. Substituents of different nature and position on the adenine ring did not change the preferred protonation site, which remained N1. However, for 2-thioether-ATP derivatives a mixed population of N1 and N7 singly protonated species was observed. Reduction of basicity of 0.4-1 pK(a) units relative to ATP was also observed for all evaluated ATP derivatives, except for 2-Cl-ATP, for which K(a) was ca. 10,000-fold lower. To explain the substitution-dependent variations in the experimental pK(a) values of the ATP analogues, gas-phase proton affinities (PA), Delta Delta G(hyd), and pK(a) values of the corresponding adenine bases were calculated using quantum mechanical methods. The computed PA and Delta Delta G(hyd) values successfully explained the experimental pK(a) values. A computational procedure for the prediction of accurate pK(a) values was developed using density functional theory and polarizable continuum model calculations. In this procedure, we developed a set of parameters for the polarizable continuum model that was fitted to reproduce experimental pK(a) values of nitrogen heterocycles. This method is proposed for the prediction of pK(a) values and protonation site(s) of purine analogues that have not been synthesized or analyzed.     

Theoretical evaluation of pK(a) in phosphoranes: implications for phosphate ester hydrolysis

Knowledge of the pK(a) of phosphoranes is important for the interpretation of phosphate ester hydrolysis. Calculated pK(a)'s of the model phosphorane, ethylene phosphorane, are reported. The method of calculation is based on the use of dimethyl phosphate as a reference state for evaluating relative pK(a) values, and on the optimization of the oxygen and acidic hydrogen van der Waals radii to give reasonable pK(1)(a), pK(2)(a), and pK(3)(a) for phosphoric acid in solution. Density functional theory is employed to calculate the gas-phase protonation energies, and continuum dielectric methods are used to determine the solvation corrections. The calculated pK(1)(a) and p(2)(a) for the model phosphorane are 7.9 and 14.3, respectively. These values are within the range of proposed experimental values, 6.5-11.0 for pK(1)(a), and 11.3-15.0 for pK(2)(a). The mechanistic implications of the calculated pK(a)'s are discussed.