Development of a chemometric correlation technique to estimate acid-base descriptors for cationic acids in non-aqueous media

We propose a procedure for estimating acid-base constants in organic solvents or mixture of solvents from the corresponding pK_a values in aqueous medium and from certain properties of the organic solvents that characterize them. To accomplish this, we developed and validated a chemometric correlation for the calculation of the acid-base constants of different cationic acids in a broad variety of non-aqueous solvents. The parameters chosen for building the model were as follows: the acid-base constant of the compound in aqueous medium and those corresponding to the polarity-polarizability, basicity and acidity scales of the solvent. The results of the fitting were significant (p < 0.01), with a root mean error in cross-validation of 18%, with no overfitting. The prediction of the acid-base constants for an external set of compounds had a mean absolute prediction error value of less than 0.8 pK_a units.     

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.     

Benchmarking and validating algorithms that estimate p<Emphasis Type="Italic">K</Emphasis><Subscript>a</Subscript> values of drugs based on their molecular structures

The REGDIA regression diagnostics algorithm in S-Plus is introduced in order to examine the accuracy of pK _a predictions made with four updated programs: PALLAS, MARVIN, ACD/pKa and SPARC. This report reviews the current status of computational tools for predicting the pK _a values of organic drug-like compounds. Outlier predicted pK _a values correspond to molecules that are poorly characterized by the pK _a prediction program concerned. The statistical detection of outliers can fail because of masking and swamping effects. The Williams graph was selected to give the most reliable detection of outliers. Six statistical characteristics (F _exp, R ², , MEP, AIC, and s(e) in pK _a units) of the results obtained when four selected pK _a prediction algorithms were applied to three datasets were examined. The highest values of F _exp, R ², , the lowest values of MEP and s(e), and the most negative AIC were found using the ACD/pK _a algorithm for pK _a prediction, so this algorithm achieves the best predictive power and the most accurate results. The proposed accuracy test performed by the REGDIA program can also be applied to test the accuracy of other predicted values, such as log P, log D, aqueous solubility or certain physicochemical properties of drug molecules.     

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.     

Basicity of some aliphatic amines in mixed H2O?CH3CN solvents

The values of pK_a^ms (K_a^ms represents ionization constant of conjugate acid of amine base in mixed water–acetonitrile solvent) for all amines, except for charged amine bases, show a mild decrease (ca. 0.1–0.4 pK units) with the increase in CH₃CN content from 2 to ∼60% v/v. However, the pK_a^ms values at 70% v/v CH₃CN become nearly equal or slightly larger (by ≤0.7 pK units) than the corresponding pK_a^ms at 2% v/v CH₃CN for all neutral and charged amines. The values of pK_a^ms for phenol increase from 10.17 to 13.38 with the increase in the content of CH₃CN from 2 to 70% v/v in mixed aqueous solvent. Taft reaction constants, ρ*, obtained from the plots of pK_a^ms against ∑σ* for primary and secondary amines decrease by ca. 0.8 ρ* units with the increase in the CH₃CN content from 2 to 70% v/v. The values of pK_a^ms show an empirical linear relationship with the corresponding values of pK_a^w (where pK_a^w represents the pK_a obtained in aqueous solvent containing 2% v/v CH₃CN), which allows the estimation of a pK_a in mixed H₂O CH₃CN solvents from that in water. © 2000 John Wiley & Sons, Inc. Int J Chem Kinet 32: 146–152, 2000     

SAMPL7 blind challenge: quantum–mechanical prediction of partition coefficients and acid dissociation constants for small drug-like molecules

The physicochemical properties of a drug molecule determine the therapeutic effectiveness of the drug. Thus, the development of fast and accurate theoretical approaches for the prediction of such properties is inevitable. The participation to the SAMPL7 challenge is based on the estimation of logP coefficients and pK_a values of small drug-like sulfonamide derivatives. Thereby, quantum mechanical calculations were carried out in order to calculate the free energy of solvation and the transfer energy of 22 drug-like compounds in different environments (water and n-octanol) by employing the SMD solvation model. For logP calculations, we studied eleven different methodologies to calculate the transfer free energies, the lowest RMSE value was obtained for the M06L/def2-TZVP//M06L/def2-SVP level of theory. On the other hand, we employed an isodesmic reaction scheme within the macro pK_a framework; this was based on selecting reference molecules similar to the SAMPL7 challenge molecules. Consequently, highly well correlated pK_a values were obtained with the M062X/6–311+G(2df,2p)//M052X/6–31+G(d,p) level of theory.

     

Linear Free‐Energy Relationships between a Single Gas‐Phase Ab Initio Equilibrium Bond Length and Experimental pKa Values in Aqueous Solution

Remarkably simple yet effective linear free energy relationships were discovered between a single ab initio computed bond length in the gas phase and experimental pK_a values in aqueous solution. The formation of these relationships is driven by chemical features such as functional groups, meta/para substitution and tautomerism. The high structural content of the ab initio bond length makes a given data set essentially divide itself into high correlation subsets (HCSs). Surprisingly, all molecules in a given high correlation subset share the same conformation in the gas phase. Here we show that accurate pK_a values can be predicted from such HCSs. This is achieved within an accuracy of 0.2 pK_a units for 5 drug molecules.     

Determination of aqueous acid‐dissociation constants of aspartic acid using PCM/DFT method

Determination of acid‐dissociation constants, pK_a, of aspartic acid in aqueous solution, using density functional theory calculations combined with the conductor‐like polarizable continuum model (CPCM) and with integral‐equation‐formalism polarizable continuum model (IEFPCM) based on the UAKS and UAHF radii, was carried out. The computed pK_a values derived from the CPCM and IEFPCM with UAKS cavity model of bare structures of the B3LYP/6‐31+G(d,p)‐optimized tetrahydrated structures of aspartic acid species are mostly close to the experimental pK_a values. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

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     

Potentiometric and spectrophotometric pKa determination of water-insoluble compounds: validation study in a new cosolvent system

In this paper the validation of pK_a determination in MDM-water mixtures is presented. The MDM-water mixture is a new multicomponent cosolvent mixture (consisting of equal volumes of methanol, dioxane and acetonitrile, as organic solvents) that dissolves a wide range of poorly water-soluble compounds. The cosolvent dissociation constants (p_sK_a) of 50 chemically diverse compounds (acids, bases and ampholytes) were measured in 15-56 wt% MDM-water mixtures by potentiometric or spectrophotometric titration and the aqueous pK_a values obtained by extrapolation. Three different extrapolation procedures were compared in order to choose the best extrapolation in MDM-water mixture using a sub-set of 30 water-soluble compounds. The extrapolated results are in good agreement with pK_a values measured in aqueous medium. No significant difference was found among these extrapolation procedures thus the widely used Yasuda-Shedlovsky plot was proposed for MDM cosolvent also. Further we also present that the single point estimation based on measurement in 20%/v MDM-mixture using a general calibration equation may be suitable for rapid pK_a determination in the early phase of drug research.     

Determination of thermodynamic pKa values of pharmaceuticals from five different groups using capillary electrophoresis

A determination of the thermodynamic acid dissociation constants (pK_a) of 22 frequently used pharmaceuticals using capillary electrophoresis in aqueous media is presented in this work. The investigated pharmaceuticals belong to different pharmacological groups: macrolides, fluoroquinolones, sulfonamides, β‐lactams, tetracyclines, and other miscellaneous pharmaceuticals. The electrophoretic mobilities of the investigated analytes were monitored in a pH range from 2.00 to 10.82. The data were fitted with an appropriate mathematical model using a nonlinear regression analysis to obtain pK_a values. Experimentally obtained data were well described by the mathematical model chosen for each analyte that was confirmed by r² values higher than 0.99 for most of the investigated analytes. Extrapolations to zero ionic strength were used to determine the thermodynamic pK_a values. Experimentally obtained acid dissociation constants were interpreted using structural formulae of investigated analytes and the moieties corresponding to specific pK_a were identified.     

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.     

Acidity of Several Anilinium Derivatives in Pure Tetrahydrofuran

The acidity constant (pK _a) of eleven substituted anilinium ions and the dissociation constants of their perchlorate salts (pK _salt) were determined in pure tetrahydrofuran by potentiometry and conductometry. The pK _a values of the studied aniliniums extend downward the range of previously determined pK _a values. The resolution of acid strength for cationic acids in tetrahydrofuran was compared with those obtained in other amphiprotic and aprotic solvents. It is shown that the resolution in tetrahydrofuran is higher than the ones in water and methanol, similar to those in acetone, dimethyl sulfoxide and isobutylmethylketone, but lower than those in acetonitrile and nitromethane.     

p<Emphasis Type="Italic">K</Emphasis><Subscript>a</Subscript> measurements for the SAMPL6 prediction challenge for a set of kinase inhibitor-like fragments

Determining the net charge and protonation states populated by a small molecule in an environment of interest or the cost of altering those protonation states upon transfer to another environment is a prerequisite for predicting its physicochemical and pharmaceutical properties. The environment of interest can be aqueous, an organic solvent, a protein binding site, or a lipid bilayer. Predicting the protonation state of a small molecule is essential to predicting its interactions with biological macromolecules using computational models. Incorrectly modeling the dominant protonation state, shifts in dominant protonation state, or the population of significant mixtures of protonation states can lead to large modeling errors that degrade the accuracy of physical modeling. Low accuracy hinders the use of physical modeling approaches for molecular design. For small molecules, the acid dissociation constant (pK_a) is the primary quantity needed to determine the ionic states populated by a molecule in an aqueous solution at a given pH. As a part of SAMPL6 community challenge, we organized a blind pK_a prediction component to assess the accuracy with which contemporary pK_a prediction methods can predict this quantity, with the ultimate aim of assessing the expected impact on modeling errors this would induce. While a multitude of approaches for predicting pK_a values currently exist, predicting the pK_as of drug-like molecules can be difficult due to challenging properties such as multiple titratable sites, heterocycles, and tautomerization. For this challenge, we focused on set of 24 small molecules selected to resemble selective kinase inhibitors—an important class of therapeutics replete with titratable moieties. Using a Sirius T3 instrument that performs automated acid–base titrations, we used UV absorbance-based pK_a measurements to construct a high-quality experimental reference dataset of macroscopic pK_as for the evaluation of computational pK_a prediction methodologies that was utilized in the SAMPL6 pK_a challenge. For several compounds in which the microscopic protonation states associated with macroscopic pK_as were ambiguous, we performed follow-up NMR experiments to disambiguate the microstates involved in the transition. This dataset provides a useful standard benchmark dataset for the evaluation of pK_a prediction methodologies on kinase inhibitor-like compounds.     

Direct determination of pKa values of cationic acids conjugated to heterocyclic amine N-oxides in polar aprotic and amphiprotic solvents

The applicability of the direct method of pK_a determination in the case of protonated heterocyclic amine N-oxides in a series of polar non-aqueous (aprotic and amphiprotic) solvents has been tested. The method is based on the pH determination of the non-aqueous solution of complex salt (the semiperchlorate in this case) formed by the N-oxides studied. The direct method not only provides for quick (one data point per each pK_a determined), but also relatively accurate estimates of acidic dissociation constants. It has been experimentally shown on the example of substituted pyridine N-oxides that this method is precise enough in all studied non-aqueous solvents when applied to compounds of not too low basicity (the pK_a being of the order of 5 or higher). To prove this, the pK_a values of protonated monocyclic N-oxides obtained by the direct method have been compared with those resulting from the potentiometric titration curve. The agreement between the results found by using both methods is very good in most cases, the differences being within standard deviations. Based upon this observation it can be inferred that the pK_a values of protonated bicyclic N-oxides in solvents studied determined by using the direct method can be also considered reliable, especially in the case of polar aprotic solvents.     

Determination of acidity constants of pyridines,imidazoles, and oximes by capillary electrophoresis

The acidity constant in the form of pK_a is one of the most important physicochemical quantities. There are prediction tools available for calculating the pK_a, but they only deliver precise calculated values for a relatively small set of chemicals. For complex structures with multiple functional groups in particular, the error in the predicted pK_a is high due to the application domain of the corresponding models. Thus, we aim to enlarge the dataset of experimentally determined pK_a values using capillary electrophoresis. We, therefore, selected various pyridines, imidazoles, and oximes to determine the pK_a values using the internal standard approach and the classical method. Especially oximes were not investigated in the past, and predictions for them include larger errors. Thus, our experimentally determined values could contribute to an improved understanding of various functional groups impacting the pK_a values and serve as additional datasets to develop improved pK_a prediction tools.     

Absolute and relative p<Emphasis Type="Italic">K</Emphasis><Subscript>a</Subscript> predictions via a DFT approach applied to the SAMPL6 blind challenge

In this work, quantum mechanical methods were used to predict the microscopic and macroscopic pK_a values for a set of 24 molecules as a part of the SAMPL6 blind challenge. The SMD solvation model was employed with M06-2X and different basis sets to evaluate three pK_a calculation schemes (direct, vertical, and adiabatic). The adiabatic scheme is the most accurate approach (RMSE?=?1.40 pK_a units) and has high correlation (R²?=?0.93), with respect to experiment. This approach can be improved by applying a linear correction to yield an RMSE of 0.73 pK_a units. Additionally, we consider including explicit solvent representation and multiple lower-energy conformations to improve the predictions for outliers. Adding three water molecules explicitly can reduce the error by 2–4 pK_a units, with respect to experiment, whereas including multiple local minima conformations does not necessarily improve the pK_a prediction.     

Polynuclear Branched Tetrazole Systems. 2. New 2-(5-Tetrazolyl)ethylpodands and Their NH-Acidity

Nine new polynuclear 2-(5-tetrazolyl)ethyl podands have been obtained by the azidation of the corresponding nitriles. Using Bjerrum distribution functions, the values of pK _a ¹, pK _a ², pK _a ³, and pK _a ⁴ have been determined by a potentiometric method for 14 polynuclear tetrazoles in aqueous and aqueous methanolic solution. The found values lie in the range from 3.5 to 7.5 pH units. The overall rules and the sequence of the ionization of the spatially separated tetrazole fragments in these podand systems are discussed.     

Physicochemical Properties of a New Multicomponent Cosolvent System for the pKa Determination of Poorly Soluble Pharmaceutical Compounds

A mixture of cosolvents is described that significantly improves the solubility of most pharmaceutical compounds. The mixture consists of equal volumes of MeOH, 1,4‐dioxane, and MeCN, thereby containing polar and nonpolar solvents, and is referred to as MDM (from MeOH, dioxane, and MeCN). MDM is mixed with H₂O until the required composition is reached. The utility of this system is that it enables analytical measurements to be performed on a wide range of compounds where measurements would be impaired in aqueous solution. We present the physicochemical characteristics of MDM/H₂O mixtures (density, dielectric constant, p_sK_w) and the principles of pK_a measurement in this solvent/H₂O mixture. We also present pK_a values in H₂O of several drug compounds determined from values measured in MDM/H₂O mixtures.     

Calculations of pKa values of carboxylic acids in aqueous solution using density functional theory

With the specific aim of calculating the acidity equilibrium constant (K_a) of carboxylic acids in aqueous solution we investigated the solute-solvent interactions of these acids and their corresponding anions. The pK_a (−lg K_a) values have been calculated using density functional theory (DFT). The polarized continuum model (PCM) is used to describe the solvent. Using these methods, we successfully predicted the pK_as of 66 carboxylic acids in aqueous with the average error of 0.5 in pK_a units. Two different thermodynamic cycles have been studied. The theoretical values are in better agreement with the experimental results for those acids with moderate strength of acidity with the pK_a value higher than 3.