Effects of flexibility in coarse-grained models for bovine serum albumin and immunoglobulin G.

We construct a coarse-grained, structure-based, low-resolution, 6-bead flexible model of bovine serum albumin (BSA, PDB: 4F5S), which is a popular example of a globular protein in biophysical research. The model is obtained via direct Boltzmann inversion using all-atom simulations of a single molecule, and its particular form is selected from a large pool of 6-bead coarse-grained models using two suitable metrics that quantify the agreement in the distribution of collective coordinates between all-atom and coarse-grained Brownian dynamics simulations of solutions in the dilute limit. For immunoglobulin G (IgG), a similar structure-based 12-bead model has been introduced in the literature [Chaudhri et al., J. Phys. Chem. B 116, 8045 (2012)] and is employed here to compare findings for the compact BSA molecule and the more anisotropic IgG molecule. We define several modified coarse-grained models of BSA and IgG, which differ in their internal constraints and thus account for a variation of flexibility. We study denser solutions of the coarse-grained models with purely repulsive molecules (achievable by suitable salt conditions) and address the effect of packing and flexibility on dynamic and static behavior. Translational and rotational self-diffusivity is enhanced for more elastic models. Finally, we discuss a number of effective sphere sizes for the BSA molecule, which can be defined from its static and dynamic properties. Here, it is found that the effective sphere diameters lie between 4.9 and 6.1 nm, corresponding to a relative spread of about ±10% around a mean of 5.5 nm.

Bulk structural information from density functionals for patchy particles.

We investigate bulk structural properties of tetravalent associating particles within the framework of classical density functional theory, building upon Wertheim's thermodynamic perturbation theory. To this end, we calculate density profiles within an effective test-particle geometry and compare to radial distribution functions obtained from computer simulations. We demonstrate that a modified version of the functional proposed by Yu and Wu [J. Chem. Phys. 116, 7094 (2002)] based on fundamental measure theory for hard spheres produces accurate results, although the functional does not satisfy the exactly known low-density limit. In addition, at low temperatures where particles start to form an amorphous tetrahedral network, quantitative differences between simulations and theory emerge due to the absence of geometrical information regarding the patch arrangement in the latter. Indeed, here we find that the theory fits better to simulations of the floating-bond model [E. Zaccarelli et al., J. Chem. Phys. 127, 174501 (2007)], which exhibits a weaker tetrahedral order due to more flexible bonds between particles. We also demonstrate that another common density functional approach by Segura et al. [Mol. Phys. 90, 759 (1997)] fails to capture fundamental structural properties.