Modeling sickle hemoglobin fibers as one chain of coarse-grained particles

Sickle cell disease (SCD) is caused by a single point mutation in the beta-chain hemoglobin gene, resulting in the presence of abnormal hemoglobin S (HbS) in the patients' red blood cells (RBCs). In the deoxygenated state, the defective hemoglobin tetramers polymerize forming stiff fibers which distort the cell and contribute to changes in its biomechanical properties. Because the HbS fibers are essential in the formation of the sickle RBC, their material properties draw significant research interests. Here, a solvent-free coarse-grain molecular dynamics (CGMD) model is introduced to simulate single HbS fibers as a chain of particles. First, we show that the proposed model is able to efficiently simulate the mechanical behavior of single HbS fibers. Then, the zippering process between two HbS fibers is studied and the effect of depletion forces is investigated. Simulation results illustrate that depletion forces play a role comparable to direct fiber-fiber interaction via Van der Waals forces. This proposed model can greatly facilitate studies on HbS polymerization, fiber bundle and gel formation as well as interaction between HbS fiber bundles and the RBC membrane.