Improving pKa prediction accuracy for PROTACs