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.     

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.     

Proton affinity of β-oxalylaminoalanine (BOAA): Incorporation of direct entropy correction into the single-reference kinetic method

A new version of the single-reference-extended kinetic method is presented in which direct entropy correction is incorporated. Results of calibration experiments with the monodentate base pyridine and the bidentate base ethylenediamine are presented for which the method provides proton affinities in excellent agreement with published values and reasonable predictions for the protonation entropies. The method is then used to determine the proton affinity and protonation entropy of the non-protein amino acid beta-oxalylaminoalanine (BOAA). The PA of BOAA is found to be 933.1 +/- 7.8 kJ/mol and a prediction for the protonation entropy of -39 J mol(-1) K(-1) is also obtained, indicating a significant degree of intramolecular hydrogen bonding in the protonated form. These results are supported by hybrid density functional theory calculations at the B3LYP/6-311++G**//B3LYP/6-31+G* level. They indicate that the preferred site of protonation is the alpha-nitrogen atom (PA = 935.0 kJ/mol) and that protonated BOAA has a strong hydrogen bond between the hydrogen on the alpha-amino group and one of the carbonyl oxygen atoms on the side chain.     

Complexes of DNA bases and Watson-Crick base pairs with small neutral gold clusters

The nature of the DNA-gold interaction determines and differentiates the affinity of the nucleobases (adenine, thymine, guanine, and cytosine) to gold. Our preliminary computational study [Kryachko, E. S.; Remacle, F. Nano Lett. 2005, 5, 735] demonstrates that two major bonding factors govern this interaction: the anchoring, either of the Au-N or Au-O type, and the nonconventional N-H...Au hydrogen bonding. In this paper, we offer insight into the nature of nucleobase-gold interactions and provide a detailed characterization of their different facets, i.e., geometrical, energetic, and spectroscopic aspects; the gold cluster size and gold coordination effects; proton affinity; and deprotonation energy. We then investigate how the Watson-Crick DNA pairing patterns are modulated by the nucleobase-gold interaction. We do so in terms of the proton affinities and deprotonation energies of those proton acceptors and proton donors which are involved in the interbase hydrogen bondings. A variety of properties of the most stable Watson-Crick [A x T]-Au3 and [G x C]-Au3 hybridized complexes are described and compared with the isolated Watson-Crick A x T and G x C ones. It is shown that enlarging the gold cluster size to Au6 results in a rather short gold-gold bond in the Watson-Crick interbase region of the [G x C]-Au6 complex that bridges the G x C pair and thus leads to a significant strengthening of G x C pairing.     

A reevaluation of computed proton affinities for the common α-amino acids

The proton affinities of the 20 common amino acids have been computed at the G3MP2 level using structures derived from broad conformational searches at a variety of levels including G3MP2. In some cases, the conformational surveys identified more stable species than had been used in previous studies of proton affinities, though the differences in energy are sometimes rather small. The present values are likely the most reliable measure of amino acid proton affinities in the gas phase. An analysis of differences between these values and those obtained experimentally via the kinetic method indicates that the extraction of proton affinities from kinetic method data can potentially lead to large errors linked to the estimation of relative protonation entropies.     

Proton affinity of uracil. A computational study of protonation sites

Relative stabilities of uracil tautomers and cations formed by gas-phase protonation were studied computationally with the B3LYP, MP2, QCISD, and QCISD(T) methods and with basis sets expanding from 6-31G(d,p) to 6-311+G(3df,2p). In accordance with a previous density functional theory study, the dioxo tautomer 1a was the most stable uracil isomer in the gas phase. Gibbs free energy calculations using effective QCISD(T)/6-311+G(3df,2p) energies suggested >99.9% of 1a at equilibrium at 523 K. The most stable ion isomer corresponded to N-1 protonated 2,4-dihydroxypyrimidine, which however is not formed by direct protonation of 1a. The topical proton affinities in 1a followed the order O-8 > O-7 > C-5 > N-3 > N-1. The thermodynamic proton affinity of 1a was calculated as 858 kJ mol⁻¹ at 298 K. A revision is suggested for the current estimate included in the ion thermochemistry database.     

Evaluation of the protonation thermochemistry obtained by the extended kinetic method

An evaluation of the results obtained by the extended kinetic method for a series of representative bases is presented here. Analysis of the original experimental data is conducted using the orthogonal distance regression (ODR) statistical treatment. A comparison with the proton affinities and protonation entropies obtained from variable temperature equilibrium constant measurements demonstrate deviations, which may be ascribed to random and systematic errors. Considerable random errors are associated with the extended kinetic method if the number of reference bases and the range of effective temperatures are too low. It is also confirmed that large systematic errors on proton affinities and protonation entropies are obtained when large protonation entropy is associated with the considered system. It is, however, encouraging to note that the gas phase basicities obtained by the extended kinetic method are generally comparable to that obtained by other methods within a few kJ/mol.     

Theoretical investigation of the proton affinities of benzazoles in the gas phase and in solution

The ground‐state equilibrium geometries of benzothiazole, benzoxazole, and benzimidazole were optimized at the density functional theory (DFT)/6‐31G** level of theory. Proton affinities on each of the possible sites in the studied series of compounds have been calculated at the DFT/6‐31G**/6‐311++G** level. The results indicate clearly that N‐site protonation is strongly favored over X‐site protonation (X = NH, O, S) for the series studied. Correlation of the computed proton affinities to the energy (E_HOMO) of the highest occupied MO in the gas phase and in solution has been explored and discussed. A comprehensive investigation of the effect of solvent on the process of protonation of the studied compounds has been performed. Different dielectric continuum models (i.e., Onsager, PCM, and IPCM) have been tested; their performance and range of applicability are reported and discussed. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Mononucleotide gas-phase proton affinities as determined by the kinetic method

The goal of this work is to determine the proton affinities of (deoxy)nucleoside 5'- and 3'-monophosphates (mononucleotides) using the kinetic method with fast atom bombardment mass spectrometry. The proton affinities of the (deoxy)nucleoside 5'- and 3'-monophosphates yielded the following trend: (deoxy)adenosine monophosphates > (deoxy)guanosine monophosphates > (deoxy)cytidine monophosphates > deoxythymidine/uridine monophosphates. In all cases the proton affinity decreases or remains the same with the addition of the phosphate group from those values reported for nucleosides. The proton affinity is dependent on the location of the phosphate backbone (5'-vs. 3'-phosphates): the 3'-monophosphates have lower proton affinities than the 5'-monophosphates except for the thymidine/uridine monophosphates where the trend is reversed. Molecular modeling was utilized to determine if multiple protonation sites and intramolecular hydrogen bond formation would influence the proton affinity measurements. Semiempirical calculations of the proton affinities at various locations on each mononucleotide were performed and compared to the experimental results. The possible influence of intramolecular hydrogen bonding between the nucleobases and the phosphate group on the measured and calculated proton affinities is discussed.     

Proton affinities of ketones, vicinal diketones and α-keto esters: a computational study

The proton affinities of seven different ketones, vicinal diketones, and α-keto esters (acetophenone, 2,2,2-trifluoroacetophenone, 2,3-butanedione, 1-phenyl-1,2-propanedione, methyl pyruvate, ethyl benzoylformate, and ketopantolactone) have been evaluated theoretically using the conventional ab initio HF and several post-HF methods (MP2, MP4, CCSD), density functional methods with the B3LYP hybrid functional, as well as some ab initio model chemistries [CBS-4M, G2(MP2), and G3(MP2)//B3LYP]. The chemical compounds studied are frequently used substrates in the asymmetric hydrogenation over chirally modified platinum catalysts where the protonation properties of the chiral modifier and the substrates are of great interest. In most cases, the proton affinities (PAs) evaluated with the CCSD/6-311+G(d,p)//B3LYP/TZVP and G2(MP2) methods are in good agreement with the existing experimental ones. However, the previously reported PA of 2,3-butanedione seems to be too high by 10-15 kJ mol⁻¹. The B3LYP/TZVP//B3LYP/TZVP and MP2/6-311+G(d,p)//B3LYP/TZVP model chemistries predict proton affinities that are systematically higher and lower than the experimental PAs, respectively. If proton affinities are evaluated as the average of the PAs calculated with these two theoretical methods a very good agreement with the experimental results is obtained. The mean absolute deviation (MAD) from experiment of this combination method for the PAs of 13 test molecules is 4.0 kJ mol⁻¹. For 9 molecules composed only of first-row atoms the MAD is 2.5 kJ mol⁻¹. The B3LYP/TZVP//B3LYP/TZVP and MP2/6-311+G(d,p)//B3LYP/TZVP methods provide significant savings in computational time and disk space compared to the CCSD/6-311+G(d,p)//B3LYP/TZVP and G2(MP2) models. Therefore, it is suggested that if no experimental or highly accurate theoretical data is available (due to computational cost), the proton affinities of similar compounds as investigated in this paper, can be evaluated with the combination method. For the studied molecules, this method gives the following PAs (in kJ mol⁻¹): 788 (2,3-butanedione, exptl 802); 798 (2,2,2-trifluoroacetophenone, exptl 799); 811 (ketopantolactone); 813 (methyl pyruvate); 825 (1-phenyl-1,2-propanedione); 862 (acetophenone, exptl 861); 865 (ethyl benzoylformate).     

Evaluation of MNDO calculated proton affinities

A large set of charged species arising mainly from protonation or deprotonation of hydrocarbons, alcohols, ethers, carboxylic acids, amines, imines, and nitriles has been studied by means of the semiempirical self-consistent-field (SCF ) molecular orbital (MO ) MNDO method. From the calculated heats of formation of such charged species and those of neutral molecules, MNDO -estimated proton affinities have been obtained and the results compared with experimental gas-phase proton affinities. If the small size anions and acetylides, for which the method predicts heats of formation too large, are ruled out, the mean absolute error in calculated proton affinities is ca. 7 kcal/mol for hydrocarbons (22 acid-base pairs) and ca. 8 kcal/mol for oxygen-containing compounds (25 acid-base pairs). For nitrogen-containing molecules it is necessary to discard, in addition, the values corresponding to the protonation of alkylamines and imines in order to achieve a reasonable mean absolute error of 7–8 kcal/mol.     

Structural Effects on the Intrinsic Basicities of alpha,beta-Unsaturated Lactones and Ketones

The proton affinities of 2(5H)-furanone, 1 (836 kJ/mol), 5,6-dihydro-2H-pyran-2-one, 2 (862 kJ/mol), cyclopentenone, 3 (857 kJ/mol), and cyclohexenone, 4 (863 kJ/mol), have been measured by Fourier transform ion cyclotron resonance techniques. A comparison is made with (reexamined) data concerning saturated cyclic and unsaturated aliphatic analogs. Three general observations are made. First, the basicity is found to increase with the size of the ring. Second, unsaturated lactones are more basic than their corresponding aliphatic unsaturated esters. Third, unsaturated and saturated lactones have almost identical gas-phase basicities, while unsaturated and saturated lactones have almost identical gas-phase basicities, while unsaturated cyclic ketones are more basic than their saturated analogs. All these experimental findings have been rationalized by means of ab initio calculations up to the G2(MP2,SVP) level. The basicity trends along the series are the result of two main factors: the different hybridization pattern of the carbonyl carbon as the size of the ring changes and, in the case of lactones, the nonbonding interaction between the proton attached to the carbonyl group and the ether-like oxygen which contributes to the enhanced stability of the protonated form. For unsaturated ketones the C=C double bond participates fully in the change in charge distribution induced by the protonation, while for unsaturated lactones the existence of an oxygen atom within the ring impedes this shift of the electron density.     

Physicochemical and structural properties of bacteriostatic sulfonamides: Theoretical study

A theoretical study at the Hartree–Fock and density functional theory levels is performed on sulfonamide‐type bacteriostatic compounds with the aim to provide an insight into their structure–activity relationship. The basicity of the p‐amino group is analyzed by means of the proton affinities and the protonation energies, showing that molecules presenting bacteriostatic activity are less basic, i.e., they are characterized by larger protonation energies and smaller proton affinities. The acidity of the amide group is analyzed through the deprotonation energy. The results reveal that the more acidic molecules present a larger bacteriostatic activity. This result is also confirmed from a study of bond orders. A bond order analysis of the amide group suggests that the electron attracting group in these molecules is responsible for the increase in acidity. The charge of the SO₂ group is also shown to be affected by the presence of the electron attracting group and consequently related to the acidity of the molecules. A geometric analysis shows that structures in which the amino group is more coplanar with respect to the benzenic ring possess larger bacteriostatic activity. A conformational analysis of these molecules illustrates that active molecules have relatively larger torsion energy barriers. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem 94: 165–172, 2003     

Theoretical calculations of proton affinities of azines. Prediction of the relative basicities and preferred protonation sites

Theoretical calculations at the 3-21G and 3-21 + G ab initio levels and at the MNDO and AM1 semiempirical levels of several six-membered nitrogenated heterocycles and their protonated species have been carried out. The 3–21G calculated proton affinities are systematically too high, in relation to the available experimental data, and it is estimated that inclusion of electron correlation and zero-point corrections is not sufficient to reach the desired agreement; however, additional inclusion of diffuse functions (3-21 + G/3-21G calculations) lowers the calculated proton affinities by 5.4–6.8 kcal/mol, a good agreement being thus obtained, at least for 1–7 . On the other hand, semiempirical methods underestimate the repulsion between each pair of vicinal nitrogens; however, if a correction of ?9 kcal/mol is added to the AM1 results for each pair of neighboring nitrogens containing lone pairs of electrons, the corresponding proton affinities match fairly well the available exoerimental data and corrected 3-21 + G results. As expected, all methods predict that the introduction of additional nitrogens decreases the overall absolute basicity. Futhermore, comparison of the relative basicity of the isomers and of the preferred protonation site for each isomer indicates that nitrogen atoms with (only) one α-nitrogen and without a γ-nitrogen are more basic than any others. In benzazines, MNDO and AM1 suggest that the 2,3-diaza arrangement has a higher intrinsic basicity than the 1,2-diaza arrangement.     

Stereoelectronic effects in negatively and positively (Protonated) charged species. Ab initio studies of the gauche effect in 1,4-dioxa systems

A computational study ab initio of the conformational dependence of proton affinities of 2-methoxyethoxide (MEO), dimethoxyethane (DME), and 1,4-dioxane has been carried out at the MP2/6-31+G level of theory. The results were discussed in comparison with reference systems, from simple alkoxides and ethers to anomeric moieties, in open, cyclic, and bicyclic molecules. The COCCOC species are stronger bases than the COCOC (anomeric) ones and approach regular ethers in their strength. The gauche forms in MEO and DME are altogether stronger bases than the anti forms, and anti (equatorial) protonation is preferred over gauche (axial), unless ditopic protonation is possible, like in aga-DME or cis-tetraoxadecalin. The gauche effect plays a significant role in the formation, relative stability, and reactivity of the charged species.     

Proton affinity of five‐membered heterocyclic amines: Assessment of computational procedures

Ab initio quantum chemical calculations, G3B3, second‐order Møller–Plesset (MP2), and the hybrid density functional method B3LYP were employed to compute the proton affinities of 24 heterocyclic amines. A range of basis sets are employed, starting from double‐ζ polarization quality to triple‐ζ quality basis set with augmented diffuse and polarization function. Experimental values were used to calibrate the performance of various theoretical models. The regioselectivity for the protonation has been unambiguously established by performing B3LYP/6‐31G* calculations on the possible putative sites of attack. For the given series of compounds the performance of B3LYP/6‐31++G** and G3B3 levels of theory have been in excellent agreement with the experimental results with the deviations are of the order comparable with the experimental error. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

The protonation of allene and some heteroallenes, a computational study

Protonation of allene and seven heteroallenes, X = Y = Z, at the terminal and central positions has been studied computationally at the MP2/6-311+G**, B3LYP/6-31+G**, and G3 levels. In all but one case protonation at a terminal position is preferred thermodynamically. The exception is allene, where protonation at C2 giving allyl cation prevails by about 10 kcal/mol over end-protonation, which gives the 2-propenyl cation. In the heteroallenes, protonation at a terminal carbon is strongly favored, activated by electron donation from the other terminal atom. Transition states for identity proton-transfer reactions were found for 10 of the "end-to-end" proton transfers. When the transfer termini are heteroatoms these processes are barrier free. We found no first-order saddle point structures for "center-to-center" proton transfers. An estimate of DeltaH++ for an identity center-to-center proton transfer could be made only for the reaction between the allyl cation and allene; it is approximately 22 kcal/mol higher than DeltaH++ for the end-to-end proton transfer between the 2-propenyl cation and allene. First-order saddle points for the proton transfer from H3S+ to both C1 and C2 of allene were found. The difference in activation enthalpies is 9.9 kcal/mol favoring protonation at C1 in spite of the thermodynamic disadvantage. We infer that protonation of X = Y = Z at central atoms passes through transition states much like primary carbenium (nitrenium, oxenium) cations, poorly conjugated with the attached vinylic or heterovinylic group. Several other processes following upon center protonation were studied and are discussed in the text, special attention being given to comparison of open and cyclic isomers.     

Arrhenius activation parameters for the loss of neutral nucleobases from deprotonated oligonucleotide anions in the gas phase

Arrhenius activation parameters (E(a) and A) for the loss of neutral nucleobase from a series of doubly deprotonated oligodexoynucleotide 10-mers of the type XT(9), T(9)X, and T(5)XT(4), where X = A, C, and G, have been determined using the blackbody infrared radiative dissociation technique. At temperatures of 120 to 190 degrees C, the anions dissociate exclusively by the loss of a neutral nucleobase (XH), followed by cleavage of the sugar 3' C-O bond leading to (a-XH) and w type ions or, in the case of the T(9)X(2-) ions, the loss of H(2)O. The dissociation kinetics and energetics are sensitive to the nature and position of X. Over the temperature range investigated, the kinetics for the loss of AH and GH were similar, but approximately 100 times faster than for the loss of CH. For the loss of AH and GH, the values of E(a) are sensitive to the position of the base. The order of the E(a)s for the loss of XH from the 5' and 3' termini is: C > G > A; while for T(5)XT(4) the order is: C > A > G. The trends in the values of E(a) do not parallel the trend in deprotonation enthalpies or proton affinities of the nucleobases in the gas phase, indicating that the energetic differences do not simply reflect differences in their gas phase acidity or basicity. The pre-exponential factors (A) vary from 10(10) to 10(15) s(-1), depending on the nature and position of X. These results suggest that the reactivity of individual nucleobases is influenced by stabilizing intramolecular interactions.     

Carbon-to-carbon identity proton transfer from allene, ketene, ketenimine, and thioketene to their respective conjugate anions in the gas phase. An ab initio study.

Gas-phase acidities of CH2=C=X (X = CH2, NH, O, and S) and barriers for the identity proton transfers (X=C=CH2 + HC triple bond C-X- right harpoon over left harpoon -X-C triple bond CH + CH2=C=X) as well as geometries and charge distributions of CH2=C=X, HC triple bond C-X- and the transition states of the proton transfer were determined by ab initio methods at the MP2/6-311+G(d,p)//MP2/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory. The acidities were also calculated at the CCSD(T)/6-311+G(2df,p) level. A major objective of this study was to examine how the enhanced unsaturation of CH2=C=X compared to that of CH3CH=X may affect acidities, transition state imbalances, and intrinsic barriers of the identity proton transfer. The results show that the acidities are all higher while the barriers are lower than for the corresponding CH3CH=X series. The transition states are all imbalanced but less so than for the reactions of CH3CH=X.     

Theoretical study of the symmetry of the (OH...O)- hydrogen bonds in vinyl alcohol-vinyl alcoholate systems

The interactions between substituted vinyl alcohols and vinyl alcoholates (X = NH(2), H, F, Cl, CN) are studied at the B3LYP/6-311++G(d,p) level of theory. In a first step, the conformation of the monomers is investigated and the proton affinities (PA(A(-))) of the enolates are calculated. The enols and enolates are held together by strong (OH...O)(-) hydrogen bonds, the hydrogen bond energies ranging from 19.1 to 34.6 kcal mol(-1). The optimized O...O distances are between 2.414 and 2.549 A and the corresponding OH distances from 1.134 and 1.023 A. The other geometry parameters such as C[double bond]C or CO distances also indicate that, in the minimum energy configuration, the hydrogen bonds are characterized by a double well potential. The Mulliken charges on the different atoms of the proton donors and proton acceptors and the frequencies of the nu(OH) stretching vibrations agree with this statement. All the data indicate that the hydrogen bonds are the strongest in the homomolecular complexes. The transition state for hydrogen transfer is located with the transition barrier estimated to be about zero. Upon addition of the zero-point vibration energies to the total potential energy, the barrier vanishes. This is a characteristic feature of low-barrier hydrogen bonds (LBHBs). The hydrogen bond energies are correlated to the difference 1.5 PA(AH) - PA(A(-)). The correlation predicts different energies for homomolecular hydrogen bonds, in agreement with the theoretical calculations. Our results suggest that a PA (or pK(a)) match is not a necessary condition for forming LBHBs in agreement with recent data on the intramolecular hydrogen bond in the enol form of benzoylacetone (J. Am. Chem. Soc. 1998, 120, 12117).