Theory and simulation of short-range models of globular protein solutions

We report theoretical and simulation studies of phase coexistence in model globular protein solutions, based on short-range, central, pair potential representations of the interaction among macro-particles. After reviewing our previous investigations of hard-core Yukawa and generalized Lennard-Jones potentials, we report more recent results obtained within a DLVO-like description of lysozyme solutions in water and added salt. We show that a one-parameter fit of this model, based on static light scattering and self-interaction chromatography data in the dilute protein regime, yields demixing and crystallization curves in good agreement with experimental protein-rich-protein-poor and solubility envelopes. The dependence of cloud and solubility point temperatures of the model on the ionic strength is also investigated. Our findings highlight the minimal assumptions on the properties of the microscopic interaction sufficient for a satisfactory reproduction of the phase diagram topology of globular protein solutions.