The theoretical descriptions based on the mean-field PB theory are limited in scope to systems in the, so-called, weak coupling limit 1 (e. g., weak electrostatic correlations), and are also unable to address the general case of polyelectrolytes or polyampholytes for which the conformation of a molecule is strongly coupled to the state of protonation of its acid/base groups . For such complicated systems, we are forced to rely on computer simulations, which must account for the protonation/deprotonation reactions of polymer's monomers .