The GCN4 leucine zipper is a peptide homodimer that has been the subject of a number of experimental and theoretical investigations into the determinants of affinity and specificity. Here, we utilize this model system to investigate electrostatic effects in protein binding using continuum calculations. A particularly novel feature of the computations made here is that they provide an interaction-by-interaction breakdown of the electrostatic contributions to the free energy of docking that includes changes in the interaction of each functional group with solvent and changes in interactions between all pairs of functional groups on binding. The results show that (1) electrostatic effects disfavor binding by roughly 15 kcal/mol due to desolvation effects that are incompletely compensated in the bound state, (2) while no groups strongly stabilize binding, the groups that are most destabilizing are charged and polar side chains at the interface that have been implicated in determining binding specificity, and (3) attractive intramolecular interactions (e.g., backbone hydrogen bonds) that are enhanced on binding due to reduced solvent screening in the bound state contribute significantly to affinity and are likely to be a general effect in other complexes. A comparison is made between the results obtained in an electrostatic analysis carried out calculationally and simulated results corresponding to idealized data from a scanning mutagenesis experiment. It is shown that scanning experiments provide incomplete information on interactions and, if overinterpreted, tend to overestimate the energetic effect of individual side chains that make attractive interactions. Finally, a comparison is made between the results available from a continuum electrostatic model and from a simpler surface-area dependent solvation model. In this case, although the simpler model neglects certain interactions, on average it performs rather well.