A numerical study of a binary Yukawa model in regimes characteristic of globular proteins in solutions

Abstract

The main goal of this paper is to assess the limits of validity, in the regime of low concentration and strong Coulomb coupling (high molecular charges), for a simple perturbative approximation to the radial distribution functions (RDF), based upon a low-density expansion of the potential of mean force and proposed to describe protein-protein interactions in a recent Small-Angle-Scattering (SAS) experimental study. A highly simplified Yukawa (screened Coulomb) model of monomers and dimers of a charged globular protein ($\beta$-lactoglobulin) in solution is considered. We test the accuracy of the RDF approximation, as a necessary complementary part of the previous experimental investigation, by comparison with the fluid structure predicted by approximate integral equations and exact Monte Carlo (MC) simulations. In the MC calculations, an Ewald construction for Yukawa potentials has been used to take into account the long-range part of the interactions in the weakly screened cases. Our results confirm that the perturbative first-order approximation is valid for this system even at strong Coulomb coupling, provided that the screening is not too weak (i.e., for Debye length smaller than monomer radius). A comparison of the MC results with integral equation calculations shows that both the hypernetted-chain (HNC) and the Percus-Yevick (PY) closures have a satisfactory behavior under these regimes, with the HNC being superior throughout. The relevance of our findings for interpreting SAS results is also discussed.

Physical Review E, in press (2005)

References in corpus (5)

• Effective interactions between electric double-layers
• Analytic solutions for Baxter's model of sticky hard sphere fluids within closures different from the Percus_Yevick approximation
• Structure factors for the simplest solvable model of polydisperse colloidal fluids with surface adhesion
• Interaction of proteins in solution from small angle scattering: a perturbative approach
• Small Angle Scattering data analysis for dense polydisperse systems: the FLAC program

Cited by in corpus (5)

• Colloidal systems with competing interactions: from an arrested repulsive cluster phase to a gel
• Effects of patch size and number within a simple model of patchy colloids
• Patchy sticky hard spheres: analytical study and Monte Carlo simulations
• Fluids of spherical molecules with dipolar-like nonuniform adhesion. An analytically solvable anisotropic model
• Effect of Polydispersity and Anisotropy in Colloidal and Protein Solutions: an Integral Equation Approach