Coarse-Grained Molecular Dynamics Simulations for Understanding the Impact of Short-Range Anisotropic Attractions on Structure and Viscosity of Concentrated Monoclonal Antibody Solutions

Understanding protein–protein interactions in concentrated therapeutic monoclonal antibody (mAb) solutions is desirable for improved drug discovery, processing, and administration. Here, we deduce both the net protein charge and the magnitude and geometry of short-ranged, anisotropic attractions of a mAb across multiple concentrations and cosolute conditions by comparing structure factors S(q) obtained from small-angle X-ray scattering experiments with those from molecular dynamics (MD) simulations. The simulations, which utilize coarse-grained 12-bead models exhibiting a uniform van der Waals attraction, uniform electrostatic repulsion, and short-range attractions between specific beads, are versatile enough to fit S(q) of a wide range of protein concentrations and ionic strength with the same charge on each bead and a single anisotropic short-range attraction strength. Cluster size distributions (CSDs) obtained from best fit simulations reveal that the experimental structure is consistent with small reversible oligomers in even low viscosity systems and help quantify the impact of these clusters on viscosity. The ability to systematically use experimental S(q) data together with MD simulations to discriminate between different possible protein–protein interactions, as well as to predict viscosities from protein CSDs, is beneficial for designing mAbs and developing formulation strategies that avoid high viscosities and aggregation at high concentration.