pKA values are thermodynamic properties, so I don't see why you would require transition states at all. What you need is a thermochemical calculation of the free enthalpy DG of the deprotonation at the desired temperature and then you have DG=-RTlnK. Since that involves the proper inclusion of solvent molecules, that's not easy, either, but transistion states are not required.