Improved Treatment of Ligands and Coupling Effects in Empirical Calculation and Rationalization of p K a Values

The new empirical rules for protein pKa predictions implemented in the PROPKA3.0 software package (Olsson et al. J. Chem. Theory Comput.2010, 7, 525-537) have been extended to the prediction of pKa shifts of active site residues and ionizable ligand groups in protein-ligand complexes. We present new algorithms that allow pKa shifts due to inductive (i.e., covalently coupled) intraligand interactions, as well as noncovalently coupled interligand interactions in multiligand complexes, to be included in the prediction. The number of different ligand chemical groups that are automatically recognized has been increased to 18, and the general implementation has been changed so that new functional groups can be added…