The equation of state developed herein is predicated on a hard-sphere reference with perturbations introduced via a potential function to account for electrostatic forces and for attraction between protein particles. During this process, the generalized Lennard-Jones (GLJ) pair potential function is employed. The GLJ pair potential function is employed to represent the protein-protein interaction in two-protein systems. Via the use of the relation between the equation of state and the chemical potential, the phase behavior in the aqueous two-protein system can be estimated. The partition coefficients can be obtained via these processes. The calculated values of the coefficients agree fairly well with the experimental data in the given pH and ionic strength range, with no additional adjustable model parameters.