Predicting pK a values from EEM atomic charges

Abstract

The acid dissociation constant p Kais a very important molecular property, and there is a strong interest in the development of reliable and fast methods for p Kaprediction. We have evaluated the p Kaprediction capabilities of QSPR models based on empirical atomic charges calculated by the Electronegativity Equalization Method (EEM). Specifically, we collected 18 EEM parameter sets created for 8 different quantum mechanical (QM) charge calculation schemes. Afterwards, we prepared a training set of 74 substituted phenols. Additionally, for each molecule we generated its dissociated form by removing the phenolic hydrogen. For all the molecules in the training set, we then calculated EEM charges using the 18 parameter sets, and the QM charges using the 8 above mentioned charge calculation schemes. For each type of QM and EEM charges, we created one QSPR model employing charges from the non-dissociated molecules (three descriptor QSPR models), and one QSPR model based on charges from both dissociated and non-dissociated molecules (QSPR models with five descriptors). Afterwards, we calculated the quality criteria and evaluated all the QSPR models obtained. We found that QSPR models employing the EEM charges proved as a good approach for the prediction of p Ka(63% of these models had R2 > 0.9, while the best had R2 = 0.924). As expected, QM QSPR models provided more accurate p Kapredictions than the EEM QSPR models but the differences were not significant. Furthermore, a big advantage of the EEM QSPR models is that their descriptors (i.e., EEM atomic charges) can be calculated markedly faster than the QM charge descriptors. Moreover, we found that the EEM QSPR models are not so strongly influenced by the selection of the charge calculation approach as the QM QSPR models. The robustness of the EEM QSPR models was subsequently confirmed by cross-validation. The applicability of EEM QSPR models for other chemical classes was illustrated by a case study focused on carboxylic acids. In summary, EEM QSPR models constitute a fast and accurate p Kaprediction approach that can be used in virtual screening.Electronic supplementary materialThe online version of this article (doi:10.1186/1758-2946-5-18) contains supplementary material, which is available to authorized users.

Theoretical Study of the Nontraditional Enol‐Based Photoacidity of Firefly Oxyluciferin

“…This type of charge scheme was found to yield good results when used to try to represent variations in pK a . [62,63] Indeed, the partial atomic charge of the oxygen atom decreased upon photoexcitation, which is consistent with the results of Hynes and co-workers. [61] If we compare the FranckCondon charges, it is evident that the atomic charge of the deprotonated oxygen on both benzothiazole and thiazole is similar.…”

A Test of Various Partial Atomic Charge Models for Computations on Diheteroaryl Ketones and Thioketones

“…Partial atomic charges find application in the molecular modeling of several chemically relevant areas, such as explaining structural and reactivity differences between various molecules and their conformers [15][16][17], investigating charge transfers within a single molecule and between several molecules [18][19][20], and predicting pK a variations in series of molecules [21,22]. Since the definition of partial atomic charges is not strict, various models of partial atomic charges can, however, differ significantly in the reliability of their predictions.…”

SAMPL6 challenge results from $$pK_a$$ predictions based on a general Gaussian process model

“…Previous work established partial charges as a useful feature to predict pK a on molecules. These studies used molecule sets with all the same ionizable groups considering the charge on the deprotonated and surrounding atoms [22, 23,36]. In order to apply our model to all identified ionizable groups, we needed a more general approach.…”