Block averaging was used to determine the protonation probabilities and standard deviations at each pH value. Fittings to both a single and double Henderson-Hasselbalch (HH) curve were performed via bootstrap; specifically, each protonation probability was normally expanded about its mean and a point was randomly selected from the distribution. The resultant set of points was fit to a single: