Protein aggregation is broadly important in diseases and in formulations of biological drugs. Here, we develop a theoretical model for reversible protein-protein aggregation in salt solutions. We treat proteins as hard spheres having square-well-energy binding sites, using Wertheim's thermodynamic perturbation theory. The necessary condition required for such modeling to be realistic is that proteins in solution during the experiment remain in their compact form. Within this limitation our model gives accurate liquid-liquid coexistence curves for lysozyme and γ IIIa-crystallin solutions in respective buffers. It provides good fits to the cloud-point curves of lysozyme in buffer-salt mixtures as a function of the type and concentration of salt. It than predicts full coexistence curves, osmotic compressibilities, and second virial coefficients under such conditions. This treatment may also be relevant to protein crystallization.