Erythocyte Membrane Model with Explicit Description of the Lipid Bilayer and the Spectrin Network

The membrane of the red blood cell (RBC) consists of spectrin tetramers connected at actin junctional complexes, forming a two-dimensional (2D) sixfold triangular network anchored to the lipid bilayer. Better understanding of the erythrocyte mechanics in hereditary blood disorders such as spherocytosis, elliptocytosis, and especially, sickle cell disease requires the development of a detailed membrane model. In this study, we introduce a mesoscale implicit-solvent coarse-grained molecular dynamics (CGMD) model of the erythrocyte membrane that explicitly describes the phospholipid bilayer and the cytoskeleton, by extending a previously developed two-component RBC membrane model. We show that the proposed model represents RBC membrane with the appropriate bending stiffness and shear modulus. The timescale and self-consistency of the model are established by comparing our results with experimentally measured viscosity and thermal fluctuations of the RBC membrane. Furthermore, we measure the pressure exerted by the cytoskeleton on the lipid bilayer. We find that defects at the anchoring points of the cytoskeleton to the lipid bilayer (as in spherocytes) cause a reduction in the pressure compared with an intact membrane, whereas defects in the dimer-dimer association of a spectrin filament (as in elliptocytes) cause an even larger decrease in the pressure. We conjecture that this finding may explain why the experimentally measured diffusion coefficients of band-3 proteins are higher in elliptocytes than in spherocytes, and higher than in normal RBCs. Finally, we study the effects that possible attractive forces between the spectrin filaments and the lipid bilayer have on the pressure applied on the lipid bilayer by the filaments. We discover that the attractive forces cause an increase in the pressure as they diminish the effect of membrane protein defects. As this finding contradicts with experimental results, we conclude that the attractive forces are moderate and do not impose a complete attachment of the filaments to the lipid bilayer.