Heuristic pKa calculations with ligands
