The impact of long range electrostatic interactions on catalysis in the thermolysin-like protease from Bacillus stearothermophilus was studied by analyzing the effects of inserting or removing charges on the protein surface. Various mutations were introduced at six different positions, and double-mutant cycle analysis was used to study the extent to which mutational effects were interdependent. The effects of single point mutations on the k(cat)/K(m) were non-additive, even in cases where the point mutations were located 10 A or more from the active site Zn(2+) and separated from each other by up to 25 A. This shows that catalysis is affected by large electrostatic networks that involve major parts of the enzyme. The interdependence of mutations at positions as much as 25 A apart in space also indicates that other effects, such as active site dynamics, play an important role in determining active site electrostatics. Several mutations yielded a significant increase in the activity, the most active (quadruple) mutant being almost four times as active as the wild type. In some cases the shape of the pH-activity profile was changed significantly. Remarkably, large changes in the pH-optimum were not observed.