Accurate prediction of electrostatic effects on catalytic activity is an essential component of protein design. Site-directed mutagenesis of charged groups in subtilisin of Bacillus amyloliquefaciens has provided experimental measurements of electrostatic interactions which may be used to test such theoretical methods. The pKa of the histidine of the active site has been perturbed by +0.08 to -1.0 units by modifying one or two residues. Electrostatic effects in proteins can be modelled by the algorithm of Warwicker and Watson, which uses classical electrostatics and considers both the charge position and the shape of the molecule. Here we report that the algorithm can model several pKa shifts in subtilisin to fair accuracy.