Abstract
Stomatocyte-discocyte-echinocyte (SDE) transformations in human red blood cells (RBCs) have significant influences on blood dynamics and related disorders. The mechanical properties of the RBC membrane, such as shear modulus and bending elasticity, play crucial roles in determining RBC shapes. Recent biophysical findings reveal that building a comprehensive model capable of describing SDE shape transformations is a challenging problem. Based on dissipative particle dynamics, this study develops a two-component RBC model considering the detachment between the lipid bilayer and cytoskeleton, as well as the cytoskeletal reorganization during echinocyte formation. This model is validated by comparing RBCs' geometric shape and the apparent membrane tension with previous experimental measurements. Results indicate that a complete SDE sequence represented by six typical shapes can be obtained by modulating the model's mechanical and geometric parameters. Furthermore, a phase diagram based on reduced variables is obtained using principal-component analysis, demonstrating the phase transformations among SDE shapes. Our result suggests that the transformation from discocyte to stomatocyte is primarily influenced by dimensionless bending rigidity, whereas, during echinocyte formation, three key variables, i.e., dimensionless bending rigidity, targeting cytoskeleton shrinkage ratio, and connecting pattern, have joint impacts on the formation of spicules or bumps and the development of the cytoskeletal framework. The present two-component RBC model and the associated findings provide a perspective for a deeper understanding of the SDE transformation mechanism. This framework offers new insights into biological science and potential applications in the field of biomedical engineering.