Benchmarking and validating algorithms that estimate pK a values of drugs based on their molecular structures

Abstract

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(2), R(P)(2), 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(2), R(P)(2), 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.