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.