Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
Comparative Study
. 2009 Dec;49(12):2801-12.
doi: 10.1021/ci900289x.

Comparison of nine programs predicting pK(a) values of pharmaceutical substances

Affiliations
Comparative Study

Comparison of nine programs predicting pK(a) values of pharmaceutical substances

Chenzhong Liao et al. J Chem Inf Model. 2009 Dec.

Abstract

Knowledge of the possible ionization states of a pharmaceutical substance, embodied in the pK(a) values (logarithm of the acid dissociation constant), is vital for understanding many properties essential to drug development. We compare nine commercially available or free programs for predicting ionization constants. Eight of these programs are based on empirical methods: ACD/pK(a) DB 12.0, ADME Boxes 4.9, ADMET Predictor 3.0, Epik 1.6, Marvin 5.1.4, Pallas pKalc Net 2.0, Pipeline Pilot 5.0, and SPARC 4.2; one program is based on a quantum chemical method: Jaguar 7.5. We compared their performances by applying them to 197 pharmaceutical substances with 261 carefully determined and highly reliable experimental pK(a) values from a literature source. The programs ADME Boxes 4.9, ACD/pK(a) DB 12.0, and SPARC 4.2 ranked as the top three with mean absolute deviations of 0.389, 0.478, and 0.651 and r(2) values of 0.944, 0.908, and 0.894, respectively. ACD/pK(a) DB 12.0 predicted all sites, whereas ADME Boxes 4.9 and SPARC 4.2 failed to predict 5 and 18 sites, respectively. The performance of the quantum chemical-based program Jaguar 7.5 was not as expected, with a mean absolute deviation of 1.283 and an r(2) value of 0.579, indicating the potential for further development of this type of approach to pK(a) prediction.

PubMed Disclaimer

Figures

Figure 1.
Figure 1.
Part 3 of 3. Structure set used for this study. All structures were extracted from ref 41. Calculated acidic and basic atoms are indicated in red and green, respectively. Structures of compounds 23, 74, 140, and 154 in the original reference are wrong and are shown here in their corrected form.
Figure 1.
Figure 1.
Part 3 of 3. Structure set used for this study. All structures were extracted from ref 41. Calculated acidic and basic atoms are indicated in red and green, respectively. Structures of compounds 23, 74, 140, and 154 in the original reference are wrong and are shown here in their corrected form.
Figure 1.
Figure 1.
Part 3 of 3. Structure set used for this study. All structures were extracted from ref 41. Calculated acidic and basic atoms are indicated in red and green, respectively. Structures of compounds 23, 74, 140, and 154 in the original reference are wrong and are shown here in their corrected form.
Figure 2.
Figure 2.
Distribution of (A) molecular weights of the calculated 197 compounds and (B) experimental pKa values of the 261 protonation sites used in the study.
Figure 3.
Figure 3.
Predicted pKa values calculated by nine different programs versus experimental pKa values for 261 protonation sites.
Figure 4.
Figure 4.
Summary of the pKa value predictions of 261 protonation sites. The predictions are binned by prediction accuracy with a resolution of 0.5 log unit.

Similar articles

Cited by

References

    1. Comer J; Tam K In Pharmacokinetic Optimization in Drug Research: Biological, Physicochemical, and Computational Strategies; Testa B, Waterbeemd H, Folkers G, Guy R, Eds.; Wiley-VCH: Weinheim; New York, 2001; pp 275–304.
    1. Comer JEA In Comprehensive Medicinal Chemistry II. Testa B, Waterbeemd H, Eds.; Elsevier: Oxford, UK, 2007; Vol. 5, pp 357–398.
    1. Mitra R; Shyam R; Mitra I; Miteva MA; Alexov E Calculating the protonation states of proteins and small molecules: Implications to ligand-receptor interactions. Curr. Comput-Aided Drug Des. 2008, 4, 169–179.
    1. Prue JE Ionic Equilibria. Pergamon Press: Oxford, New York, 1966.
    1. Poole SK; Patel S; Dehring K; Workman H; Poole CF Determination of acid dissociation constants by capillary electrophoresis. J. Chromatogr., A 2004, 1037, 445–454. - PubMed

Publication types

MeSH terms

Substances