I have been recently testing pH_protocol to calculate the pKa of one of my proteins. In the original paper, the authors include the weights function to be used with the protocol:
After replacing hack_elec to fa_elec everything works fine. However, while the results are mostly the same as the ones I get from running on ROSIE, a couple of sidechains get really different Pkas. Any ideas?