Prediction of protonation states in ligand-protein complexes upon ligand binding