| Crystallization of globular proteins has become a very important subject in recent yearn. However there is still no understanding of the particular conditions that lead to the crystallization. Since nucleation of a crystalline droplet is the critical step toward the formation of the solid phase from the supersaturated solution, this is the focus of current studies. In this work we use different approaches to investigate the collective behavior of a system of globular proteins. Especially we focused on the models which have a metastable critical point, because this reflects the properties of solutions of globular proteins.; The first approach is a continuum model of globular proteins. This model was first presented by Talanquer and Oxtoby and is based on the van der Waals theory. The model can have either a stable or a metastable critical point. For the system with the metastable critical point we studied the behavior of the free energy barrier to nucleation; we found that along particular pathways the barrier to nucleation has a minimim around the critical point. As well, the number of molecules in the critical cluster was found to diverge as one approaches the critical point, though most of the molecules are in the fluid tail of the droplet. Our results are an extension of earlier work [17, 7].; The properties of the solvent affect the behavior of the solution. In our second approach, we proposed a model that takes into account the contribution of the solvent free energy to the free energy of the globular proteins. We show that one can map the phase diagram of a repulsive hard core plus attractive square well interacting system to the same system particles in the solvent environment. In particular we show that this leads to phase diagrams with upper critical points, lower critical points and even closed loops with both upper and lower critical points, similar to the one found before [10]. For systems with interaction different from the square well, in the presence of the solvent this mapping procedure can be a first approximation to understand the phase diagram.; The final part of this work is dedicated to the behavior of sickle hemoglobin. While the fluid behavior of the HbS molecules can be approximately explained by the uniform interparticle potential, this model fails to describe the polymerization process and the particular structure of fibers. We develop an anisotropic "patchy" model to describe some features of the HbS polymerization process. To determine the degree of polymerization of the system a "patchy" order parameter was defined. Monte Carlo simulations for the simple two-patch model was performed and reveal the possibility of obtaining chains that can be considered as one dimensional crystals. |
