The physical properties of chemicals are strongly influenced by their protonation state, affecting, for example, solubility or hydrogen-bonding characteristics. The ability to accurately calculate protonation states in the form of pK as is, therefore, desirable. Calculations of pK a changes in a series of substituted pyridines are presented. Computations were performed using both ab initio and semiempirical approaches, including free energies of solvation via reaction-field models. The selected methods are readily accessible with respect to both software and computational feasibility. Comparison of calculated and experimental pK as shows the experimental trends to be reasonably reproduced by the computations with root-mean-square differences ranging from 1.22 to 4.14 pK a units. Of the theoretical methods applied the best agreement occurred using the second-order Møller–Plesset/6-31G(d)/isodensity surface polarized continuum solvation model, while the more computationally accessible Austin model 1/Solvent model 2 (SM2) approach yielded results similar to the ab initio methods. Analysis of component contributions to the calculated pK as indicates the largest source of error to be associated with the free energies of solvation of the protonated species followed by the gas-phase protonation energies; while the latter may be improved via the use of higher levels of theory, enhancements in the former require improvements in the solvation models. The inclusion of alternate minimum in the computation of pK as is also indicated to contribute to differences between experimental and calculated pK a values.