A characteristic of two-dimensional proteomics gels is a general bimodal distribution of isoelectric (pI) values. Discussion of this feature has focussed on the balance of acidic and basic ionisable residues, and potential relationships between pI distributions and organism classification or protein subcellular location. Electrostatics calculations on a set of protein structures with known subcellular location show that predicted folded state pI are similar to those calculated from sequence alone, but adjusted according to a general stabilising effect from interactions between ionisable groups. Bimodal distributions dominate both pI and the predicted pH of maximal stability. However, there are significant differences between these features. The average pH of maximal stability generally follows organelle pH. Average pI values are well removed from organelle pH in most subcellular environments, consistent with the view that proteins have evolved to carry (on average) net charge in a given subcellular location, and relevant to discussion of solubility in crowded environments. Correlation of the predicted pH of maximum stability with subcellular pH suggests an evolutionary pressure to adjust folded state interactions according to environment. Finally, our analysis of ionisable group contributions to stability suggests that Golgi proteins have the largest such term, although this dataset is small.