Coarse-grained simulation of the translational and rotational diffusion of globular proteins by dissipative particle dynamics

Abstract

With simplified interactions and degrees of freedom, coarse-grained (CG) simulations have been successfully applied to study the translational and rotational diffusion of proteins in solution. However, in order to reach larger lengths and longer timescales, many CG simulations employ an oversimplified model for proteins or an implicit-solvent model in which the hydrodynamic interactions are ignored, and thus, the real kinetics are more or less unfaithful. In this work, we develop a CG model based on the dissipative particle dynamics (DPD) that can be universally applied to different types of proteins. The proteins are modeled as a group of rigid DPD beads without conformational changes. The fluids (including solvent and ions) are also modeled as DPD beads. The electrostatic interactions between charged species are explicitly considered by including charge distributions on DPD particles. Moreover, a surface friction between the protein and fluid beads is applied to control the slip boundary condition. With this model, we investigate the self-diffusion of a single globular protein in bulk solution. The translational and rotational diffusion coefficients of the protein can be tuned by the surface frictional constant to fit the predictions of the Stokes-Einstein (SE) relation. We find that both translational and rotational diffusion coefficients that meet with the prediction of the SE relation based on experimental results of the hydrodynamic radius are reached at almost the same frictional constant for different types of proteins. Such scaling behavior indicates that the model can be applied to simulate the translational and rotational diffusion together for various types of proteins.