Monte Carlo simulations are used to establish the potential of mean force between two globular proteins in an aqueous electrolyte solution. This potential includes nonelectrostatic contributions arising from dispersion forces first, between the globular proteins, and second, between ions in solution and between each ion and the globular protein. These latter contributions are missing from standard models. The potential of mean force, obtained from simulation, is fitted to an analytic equation. Using our analytic potential of mean force and Barker-Henderson perturbation theory, we obtain phase diagrams for lysozyme solutions that include stable and metastable fluid-fluid and solid-fluid phases when the electrolyte is 0.2 M NaSCN or NaI or NaCl. The nature of the electrolyte has a significant effect on the phase diagram.