Relative pKa of some anilinium derivatives in methanol,acetonitrile, and tetrahydrofurane solvents

In this study, the relative pK_a values of nine anilinium derivatives in methanol (MeOH), acetonitrile (AN), and tetrahydrofurane (THF) solutions were successfully calculated with mean absolute deviations of 0.63, 0.68, and 0.75 pK_a units, respectively. To this aim, their gas‐phase basicities were computed using the CBS‐QB3 composite method. Also, conductor‐like polarizable continuum model (CPCM) with UAHF, UAKS and UA0 cavities and SM8 solvation models at HF/6‐31+G(d) level of theory were applied for the calculation of the solvation Gibbs free energies. The obtained results indicate that there is reliable correlation between the experimental and computed pK_a values in the studied solutions. Therefore, to extend the pK_a database for anilines, correlation equations were used to predict the pK_a values in the investigated solvents.     

Theoretical calculation of the pKa values of some drugs in aqueous solution

In this work, calculations of pK_a values have been performed on benzoic acid and its para‐substituted derivatives and some drugs by using Gaussian 98 software package. Gas‐phase energies were calculated with HF/6‐31 G** and B3LYP/6‐31 G** levels of theory. Free energies of solvation have been computed using the polarizable continuum model (PCM), conductor‐like PCM (CPCM), and the integral equation formalism‐PCM at the same levels which have been used for geometry determination in the gas‐phase. The results that show the calculated pK_a values using the B3LYP are better than those using the corresponding HF. In comparison to the other models, the results obtained indicate that the PCM model is a suitable solvation model for calculating pK_a values. For the investigated compounds, a good agreement between the experimental and the calculated pK_a values was also observed. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Accurate pKa Calculation of the Conjugate Acids of Alkanolamines,Alkaloids and Nucleotide Bases by Quantum Chemical Methods

The pK_a of the conjugate acids of alkanolamines, neurotransmitters, alkaloid drugs and nucleotide bases are calculated with density functional methods (B3LYP, M08‐HX and M11‐L) and ab initio methods (SCS‐MP2, G3). Implicit solvent effects are included with a conductor‐like polarizable continuum model (CPCM) and universal solvation models (SMD, SM8). G3, SCS‐MP2 and M11‐L methods coupled with SMD and SM8 solvation models perform well for alkanolamines with mean unsigned errors below 0.20 pK_a units, in all cases. Extending this method to the pK_a calculation of 35 nitrogen‐containing compounds spanning 12 pK_a units showed an excellent correlation between experimental and computational pK_a values of these 35 amines with the computationally low‐cost SM8/M11‐L density functional approach.     

Density functional theoretical study of pKa of RCOOH (RH,CH3, and C2H5) using the combination of the extended clusters‐continuum model

The accurate pK_a determinations for three carboxylic acids have been investigated using the combination of the extended clusters‐continuum model at B3LYP/6‐31+g(d,p) and B3LYP/6‐311++g(d,p) levels. To take into account of the effect of the water combined with carboxylic acids in different positions, eleven molecular clusters were considered. Among these clusters, the one involving the carboxylic acid wrapped up with water molecules and saturated with hydrogen bonds (four hydrogen bonds around ? COOH) leads to the best B3LYP pK_a results compared to the experimental data. For those clusters saturated with hydrogen bonds, when n = 3 (the number of water molecules), the average absolute errors between the calculated pK_a results and experimental data of these three carboxylic acids were 0.19 (0.23) and 0.12 (0.22) pK_a at B3LYP/6‐31+g(d,p)//PCM (IEFPCM) and B3LYP/6‐311++g(d,p)//PCM (IEFPCM) levels, respectively; when n = 4, they are 0.53 (1.23) and 1.09 (1.03) pK_a, respectively. On the basis of the above results, the molecular cluster saturated with four hydrogen bonds formed by three waters and one carboxylic acid molecule was the chief existence in the carboxylic acid solution. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Experimentation with different thermodynamic cycles used for pKa calculations on 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 Barone and Cossi's implementation of the polarizable conductor model (CPCM) continuum solvation methods to calculate pK_a values for six carboxylic acids. Four different thermodynamic cycles were considered in this work. An experimental value of ?264.61 kcal/mol for the free energy of solvation of H⁺, ΔG_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. Thermodynamic cycles that include an explicit water in the cycle are not accurate when the free energy of solvation of a water molecule is used, but appear to become accurate when the experimental free energy of vaporization of water is used. This apparent improvement is an artifact of the standard state used in the calculation. Geometry relaxation in solution does not improve the results when using these later cycles. 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. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Prediction of n-octanol/water partition coefficients and acidity constants (pKa) in the SAMPL7 blind challenge with the IEFPCM-MST model

Within the scope of SAMPL7 challenge for predicting physical properties, the Integral Equation Formalism of the Miertus-Scrocco-Tomasi (IEFPCM/MST) continuum solvation model has been used for the blind prediction of n-octanol/water partition coefficients and acidity constants of a set of 22 and 20 sulfonamide-containing compounds, respectively. The log P and pK_a were computed using the B3LPYP/6-31G(d) parametrized version of the IEFPCM/MST model. The performance of our method for partition coefficients yielded a root-mean square error of 1.03 (log P units), placing this method among the most accurate theoretical approaches in the comparison with both globally (rank 8th) and physical (rank 2nd) methods. On the other hand, the deviation between predicted and experimental pK_a values was 1.32 log units, obtaining the second best-ranked submission. Though this highlights the reliability of the IEFPCM/MST model for predicting the partitioning and the acid dissociation constant of drug-like compounds compound, the results are discussed to identify potential weaknesses and improve the performance of the method.

     

Using polarizable POSSIM force field and fuzzy‐border continuum solvent model to calculate pKa shifts of protein residues

Our Fuzzy‐Border (FB) continuum solvent model has been extended and modified to produce hydration parameters for small molecules using POlarizable Simulations Second‐order Interaction Model (POSSIM) framework with an average error of 0.136 kcal/mol. It was then used to compute pK _a shifts for carboxylic and basic residues of the turkey ovomucoid third domain (OMTKY3) protein. The average unsigned errors in the acid and base pK _a values were 0.37 and 0.4 pH units, respectively, versus 0.58 and 0.7 pH units as calculated with a previous version of polarizable protein force field and Poisson Boltzmann continuum solvent. This POSSIM/FB result is produced with explicit refitting of the hydration parameters to the pK _a values of the carboxylic and basic residues of the OMTKY3 protein; thus, the values of the acidity constants can be viewed as additional fitting target data. In addition to calculating pK _a shifts for the OMTKY3 residues, we have studied aspartic acid residues of Rnase Sa. This was done without any further refitting of the parameters and agreement with the experimental pK _a values is within an average unsigned error of 0.65 pH units. This result included the Asp79 residue that is buried and thus has a high experimental pK _a value of 7.37 units. Thus, the presented model is capable or reproducing pK _a results for residues in an environment that is significantly different from the solvated protein surface used in the fitting. Therefore, the POSSIM force field and the FB continuum solvent parameters have been demonstrated to be sufficiently robust and transferable. © 2016 Wiley Periodicals, Inc.     

Thermodynamic Estimate of pKa Values of the Carboxylic Acids in Aqueous Solution with the Density Functional Theory

The 96 pK_a values of 85 carboxylic acids in aqueous solution were calculated with the density functional theory method at the level of B3LYP/6‐31+G(d,p) and the polarizable continuum model (PCM) was used to describe the solvent. In the calculations of pK_a values, the dissociation Gibbs free energies were directly calculated using carboxylic acid dissociation reactions in aqueous solution, i. e., no thermodynamic cycle was employed, which is different from the previous literatures. A highly significant correlation of R²=0.95 with a standard deviation (SD) of 0.36 between the experimental pK_a values and the calculated dissociation Gibbs free energies [ΔG(calc.)] was found. The slope of pK_a vs. (G(calc.)/(20303RT) is only 47.6% of the theoretically expected value, which implies that the ΔG(calc.) value from the theoretical calculation is larger than the actual one for all 85 carboxylic acids studied. Thus, by adding the 0.476 scaling‐factor into the slope, we can derive a reliably procedure that can reproduce the experimental pK_a values of carboxylic acids. The pK_a values furnished by this procedure are in good agreement with the experimental results for carboxylic acids in aqueous solution.     

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.     

Calculation of acidic dissociation constants in water: solvation free energy terms. Their accuracy and impact

Three polarizable continuum models, DPCM, CPCM, and IEFPCM, have been applied to calculate free energy differences for nine neutral compounds and their anions. On the basis of solvation free energies, the pK_a values were obtained for the compounds in question by using three thermodynamic cycles: one, involving the combined experimental and calculated data, as well as two other cycles solely with calculated data. This paper deals with the influence of factors such as the SCRF model applied, choice of a particular thermodynamic cycle, atomic radii used to build a cavity in the solvent (water), optimization of geometry in water, inclusion of electron correlation, and the dimension of the basis set on the solvation free energies and on the calculated pK_a values. Electronic supplementary material The online version of this article (doi: ) contains supplementary material, which is available to authorized users.     

Experimental and computational study of the pKa of coumaric acid derivatives

The dissociation constant (pK_a) of a drug is a key parameter in drug discovery and pharmaceutical formulation. The hydroxy substituent has a significant effect on the acidity of hydroxycinnamic acid. In this work, the acidic constants of coumaric acids are obtained experimentally by spectrophotometry using the chemometric method and calculated theoretically using ab initio quantum mechanical method at the HF/6‐31G* level of theory in combination with the SMD continuum solvation method. Rank annihilation factor analysis (RAFA) is an efficient chemometric technique based on the elimination of the contribution of one of the chemical components from the data matrix. RAFA cannot be performed because the pure spectrum of HA^? is not available. So, two‐rank annihilation factor analysis (TRAFA) is proposed for the determination of the pK_a OF H₂A. A comparison between the pK_a values obtained previously by TRAFA for the molecules o‐coumaric acid (4.13, 9.58), m‐coumaric acid (4.48, 10.35), and p‐coumaric acid (4.65, 9.92) makes it clear that there is good agreement between the results obtained by TRAFA and ab initio quantum mechanical method.     

Avoiding gas-phase calculations in theoretical p<Emphasis Type="Italic">K</Emphasis><Subscript>a</Subscript> predictions

CBS-QB3, two simplified and less computationally demanding versions of CBS-QB3, DFT-B3LYP, and HF quantum chemistry methods have been used in conjunction with the CPCM continuum solvent model to calculate the free energies of proton exchange reactions in water solution following an isodesmic reaction approach. According to our results, the precision of the predicted pK _a values when compared to experiment is equivalent to that of the thermodynamic cycles that combine gas-phase and solution-phase calculations. However, in the aqueous isodesmic reaction schema, the accuracy of the results is less sensitive to the presence of explicit water molecules and to the global charges of the involved species since the free energies of solvation are not required. In addition, this procedure makes easier the prediction of pK _a values for molecules that undergo large conformational changes in solvation process and makes possible the pK _a prediction of unstable species in gas-phase such as some zwitterionic tautomers. The successive pK _a values of few amino acids corresponding to the ionization of the α-carboxylic acid and α-amine groups, which is one of the problematic cases for thermodynamic cycles, were successfully calculated by employing the aqueous isodesmic reaction yielding mean absolute deviations of 0.22 and 0.19 pK _a units for the first and second ionization processes, respectively.     

Prediction of the reduction potential in transition‐metal containing complexes: How expensive? For what accuracy?

Accurate computationally derived reduction potentials are important for catalyst design. In this contribution, relatively inexpensive density functional theory methods are evaluated for computing reduction potentials of a wide variety of organic, inorganic, and organometallic complexes. Astonishingly, SCRF single points on B3LYP optimized geometries with a reasonably small basis set/ECP combination works quite well‐‐B3LYP with the BS1 [modified‐LANL2DZ basis set/ECP (effective core potential) for metals, LANL2DZ(d,p) basis set/LANL2DZ ECP for heavy nonmetals (Si, P, S, Cl, and Br), and 6‐31G(d') for other elements (H, C, N, O, and F)] and implicit PCM solvation models, SMD (solvation model based on density) or IEFPCM (integral equation formalism polarizable continuum model with Bondi atomic radii and α = 1.1 reaction field correction factor). The IEFPCM‐Bondi‐B3LYP/BS1 methodology was found to be one of the least expensive and most accurate protocols, among six different density functionals tested (BP86, PBEPBE, B3LYP, B3P86, PBE0, and M06) with thirteen different basis sets (Pople split‐valence basis sets, correlation consistent basis sets, or Los Alamos National Laboratory ECP/basis sets) and four solvation models (SMD, IEFPCM, IPCM, and CPCM). The MAD (mean absolute deviation) values of SCRF‐B3LYP/BS1 of 49 studied species were 0.263 V for SMD and 0.233 V for IEFPCM‐Bondi; and the linear correlations had respectable R ² values (R ² = 0.94 for SMD and R ² = 0.93 for IEFPCM‐Bondi). These methodologies demonstrate relatively reliable, convenient, and time‐saving functional/basis set/solvation model combinations in computing the reduction potentials of transition metal complexes with moderate accuracy. © 2017 Wiley Periodicals, Inc.     

Improving the Capture of CO2 by Substituted Monoethanolamines: Electronic Effects of Fluorine and Methyl Substituents

The influence of electronic and steric effects on the reaction between CO₂ and monoethanolamine (MEA) absorbents is investigated using computational methods. The pK_a of the alkanolamine, the reaction enthalpy for carbamate formation, and the hydrolytic carbamate stability are important factors for the efficiency of CO₂ capture. The steric and electronic effects of CH₃, CH₂F, CHF₂, CF₃, F, dimethyl, difluoro, and bis(2‐trifluoromethyl) substituents at the α carbon of MEA on this reaction are investigated. Density functional theory (DFT) (B3LYP, M06‐2X, M08‐HX and M11‐L) and ab initio methods [spin component‐scaled second‐order Møller‐Plesset theory (SCS‐MP2), G3], each coupled with solvent models [conductor‐like polarizable continuum model (CPCM) and universal solvation models (SM8 and SMD)], are shown to yield accurately calculated pK_a values of the substituted MEAs. Specifically, G3, SCS‐MP2, and M11‐L methods coupled with the SMD and SM8 solvation models perform well with a mean unsigned error (MUE) of only 0.15, 0.24 and 0.25 pK_a units, respectively. SCS‐MP2 is used to calculate the reaction enthalpy for carbamate formation and the carbamate stability towards hydrolysis. With the introduction of β‐fluoro substituents (especially the CH₂F moiety) the reaction enthalpy for the formation of carbamates can be fine‐tuned to be less exothermic than that using the unsubstituted MEA. This implies a reduced energy requirement for the solvent‐regeneration step in the post‐combustion carbon‐capture method, which is currently the energy‐limiting step in efficient CO₂ capture. β‐Fluoro‐substituted MEAs are also shown to form less stable carbamates than MEA. Thus, β‐fluoro‐substituted MEAs display a great potential for the use in the post‐combustion carbon‐capture process. Finally, a clear correlation is observed between the gas‐phase basicity and the tendency to form carbamates. This allows for the rapid prediction of which species will be formed experimentally, and thus the CO₂‐absorbing capacities of alkanolamines can be estimated.     

Calculation of dissociation constants and chemical hardness of some biologically important molecules: A theoretical study

The gas‐phase geometries of neutral, protonated, and deprotonated forms of some biologically important molecules, alanine (Ala), glycine (Gly), phenylalanine (Phe), and tyrosine (Tyr), were optimized using density functional theory at B3LYP/6‐311++G(d) and the ab initio HF/6‐311++G(d) level of theories. The neutral and different stable ionic states of Ala, Gly, Phe, and Tyr have also been solvated in aqueous medium using polarizable continuum model for the determination of solvation free energies in the aqueous solution. The gas‐phase acidity constants of above four molecules have been also calculated at both levels of theories and found that the values calculated at HF/6‐311++G(d) method are in good agreement with experimental results. A thermodynamic cycle was used to determine the solvation free energies for the proton dissociation process in aqueous solution and the corresponding pK_a values of these molecules. The pK_a values calculated at B3LYP/6‐311++G(d) method are well supported by the experimental data with a mean absolute deviation 0.12 pK_a units. Additionally, the chemical hardness and the ionization potential (IP) for these molecules have been also explored at both the level of theories. The Tyr has less value of chemical hardness and IP at both levels of theories compared with other three molecules, Ala, Gly, and Phe. The calculated values of chemical hardness and IP are decreasing gradually with the substitution of the various functional groups in the side chain of the amino acids. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Evaluation of Acid Dissociation Constants in DMSO and DMF by Quantum-Chemical Methods

We have calculated total electronic energies (E) and Gibbs energies (G) of a large number of acids and their anions in water, dimethylsulfoxide, and dimethylformamide using the hybrid B3LYP functional DFT method in the 6-31++G(d,p) basis set, taking into account the solvent effect by the conductor-like polarizable continuum model method. A linear correlation has been found between the experimental values of acid dissociation constants (pK_a) of different nature and the difference between anion and acid E values, and between pK_a and the difference between anion and acid G values. The obtained correlations allowed us to evaluate the pK_a values of both inorganic and organic acids. Such an evaluation is of special importance for nonaqueous solvents as it is quite problematic to determine these dissociation constants.     

The salt–cocrystal spectrum in salicylic acid–adenine: the influence of crystal structure on proton‐transfer balance

At one extreme of the proton‐transfer spectrum in cocrystals, proton transfer is absent, whilst at the opposite extreme, in salts, the proton‐transfer process is complete. However, for acid–base pairs with a small ΔpK_a (pK_a of base ? pK_a of acid), prediction of the extent of proton transfer is not possible as there is a continuum between the salt and cocrystal ends. In this context, we attempt to illustrate that in these systems, in addition to ΔpK_a, the crystalline environment could change the extent of proton transfer. To this end, two compounds of salicylic acid (SaH) and adenine (Ad) have been prepared. Despite the same small ΔpK_a value (≈1.2), different ionization states are found. Both crystals, namely adeninium salicylate monohydrate, C₅H₆N₅⁺·C₇H₅O₃^?·H₂O, I , and adeninium salicylate–adenine–salicylic acid–water (1/2/1/2), C₅H₆N₅⁺·C₇H₅O₃^?·2C₅H₅N₅·C₇H₆O₃·2H₂O, II , have been characterized by single‐crystal X‐ray diffraction, IR spectroscopy and elemental analysis (C, H and N) techniques. In addition, the intermolecular hydrogen‐bonding interactions of compounds I and II have been investigated and quantified in detail on the basis of Hirshfeld surface analysis and fingerprint plots. Throughout the study, we use crystal engineering, which is based on modifications of the intermolecular interactions, thus offering a more comprehensive screening of the salt–cocrystal continuum in comparison with pure pK_a analysis.     

A new method of accurate pKb determinations for some organic amines

The accurate pK_b determinations for some amines have been investigated using the combination of the extended clusters‐continuum model with the polarizable continuum solvation model. The formation of molecular clusters by means of the amines wrapped up with water molecules leads to the weakness of the mutual function between the polar solvents and the amines, and, hence, the accuracy of pK_b has been enhanced by using a coherent and well‐defined approach without external approximations or experimental data. The calculations are performed at the HF/6‐311++G(d, p) level and agreed well with experimental data because electron correlation effect cancels mutually in the calculated value of ΔG which is not an absolute value, but a relative value. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Evaluation of DFT Methods and Implicit Solvation Models for Anion‐Binding Host‐Guest Systems

Although supramolecular chemistry is traditionally an experimental discipline, computations have emerged as important tools for the understanding of supramolecules. We have explored how well commonly used density functional theory quantum mechanics and polarizable continuum solvation models can calculate binding affinities of host‐guest systems. We report the calculation of binding affinities for eight host–guest complexes and compare our results to experimentally measured binding free energies that span the range from ?2.3 to ?6.1 kcal mol^?1. These systems consist of four hosts (biotin[6]uril, triphenoxymethane, cryptand, and bis‐thiourea) with different halide ions (F^?, Cl^?, Br^?) in various media including organic and aqueous. The mean average deviation (MAD) of calculated from measured ΔG_a is 2.5 kcal mol^?1 when using B3LYP‐D3 with either CPCM or PCM. This MAD value lowers even more by eliminating two outliers: 1.1 kcal mol^?1 for CPCM and 1.2 kcal mol^?1 for PCM. The best DFT and implicit solvation model combination that we have studied is B3LYP?D3 with either CPCM or PCM.     

Investigation of Solution p<Emphasis Type="Italic">K</Emphasis><Subscript>a</Subscript> and Thermodynamic Values of Lamivudine and Pefloxacin Drugs by Ab initio and DFT Methods

In this work we investigate the thermodynamic properties and pK_a value of lamivudine and pefloxacin drugs, in aqueous solutions, by ab initio and density functional theory (DFT) methods at different temperatures. Molecular structures and solute–solvent effects of the anions, cations, and neutral molecules of lamivudine and pefloxacin were studied by the polarizable continuum model (PCM). The calculation was done at the DFT-B3LYP/6-31+G(d) level of theory using Tomasi’s method to analyze the formation of intermolecular hydrogen bonds (IHB) in aqueous solution. The pK_a1 values of lamivudine and pK_a2 values of pefloxacin increase with temperature increase. In contrast, the pK_a1 values of pefloxacin decrease when the temperature increases. Further, the thermodynamic properties of the ionization processes (?H, ?S and ?G) of the drugs in aqueous solution were determined and discussed. The results of this work are in good agreement with the literature data at 298.15 K.