Modeling of the properties of biochemical components is gaining increasing interest due to its potential for further application within the area of biochemical process development. Generally protein solution properties such as protein solubility are expressed through component activity coefficients which are studied here. The original UNIQUAC model is chosen for the representation of protein activity coefficients and, to the best of our knowledge, this is the first time it has been directly applied to protein solutions. Ten different protein-salt-water systems with four different proteins, serum albumin, alpha-chymotrypsin, beta-lactoglobulin and ovalbumin, are investigated. A root-mean-squared deviation of 0.54% is obtained for the model by comparing calculated protein activity coefficients and protein activity coefficients deduced from osmotic measurements through virial expansion. Model predictions are used to analyze the effect of salt concentrations, pH, salt types, and temperature on protein activity coefficients and also on protein solubility and demonstrate consistency with results from other references.