The human erythrocyte membrane skeleton may be an ionic gel. II. Numerical analyses of cell shapes and shape transformations.

In the first paper in this series (Stokke et al. Eur Biophys J 1986, 13:203-218) we developed the general theory of the mechanochemical properties and the elastic free energy of the protein gel--lipid bilayer membrane model. Here we report on an extensive numerical analysis of the human erythrocyte shapes and shape transformations predicted by this new cell membrane model. We have calculated the total elastic free energy of deformation of four different cell shape classes: disc-shaped cells, cup-shaped cells, crenated cells, and cells with membrane invaginations. We find that which of these shape classes is favoured depends strongly on the spectrin gel osmotic tension, IIGu, and the surface tensions, IIEu and IIPu, of the extracellular and protoplasmic halves of the membrane lipid bilayer, respectively. For constant ratio IIEu/IIPu greater than O large negative or positive values of IIGu favour respectively the crenated and invaginated cell shape classes. For small absolute values of IIGu, IIEu, and IIPu, biconcave or cup-shaped cells are the stable ones. Our numerical analysis shows that the higher the membrane skeleton compressibility is, the smaller are the values of IIGu needed to induce cell shape transformation. We find that the stable and metastable shapes of discocytes and stomatocytes generally depend both on the shape of the stressfree membrane skeleton and the membrane skeleton compressibility.