We study the biomechanical interactions between the lipid bilayer and the cytoskeleton in a red blood cell (RBC) by developing a general framework for mesoscopic simulations. We treated the lipid bilayer and the cytoskeleton as two distinct components and developed a unique whole-cell model of the RBC, using dissipative particle dynamics (DPD). The model is validated by comparing the predicted results with measurements from four different and independent experiments. First, we simulated the micropipette aspiration and quantified the cytoskeletal deformation. Second, we studied the membrane fluctuations of healthy RBCs and RBCs parasitized to different intraerythrocytic stages by the malaria-inducing parasite Plasmodium falciparum . Third, we subjected the RBC to shear flow and investigated the dependence of its tank-treading frequency on shear rate. Finally, we simulated the bilayer–cytoskeletal detachment in channel flow to quantify the strength of such interactions when the corresponding bonds break. Taken together, these experiments and corresponding systematic DPD simulations probe the governing constitutive response of the cytoskeleton, elastic stiffness, viscous friction, and strength of bilayer–cytoskeletal interactions as well as membrane viscosities. Hence, the DPD simulations and comparisons with available independent experiments serve as validation of the unique two-component model and lead to useful insights into the biomechanical interactions between the lipid bilayer and the cytoskeleton of the RBC. Furthermore, they provide a basis for further studies to probe cell mechanistic processes in health and disease in a manner that guides the design and interpretation of experiments and to develop simulations of phenomena that cannot be studied systematically by experiments alone.

The cell uptake rate of nanoparticles (NPs) coated with mixed hydrophilic/hydrophobic ligands is known to be strongly influenced by the ligand pattern on the nanoparticle surface. To help reveal the physical mechanism behind this intriguing phenomenon, here we perform dissipative particle dynamics simulations to analyze the evolution of free energy as the ligand-coated NPs pierce through a lipid bilayer. Four characteristic ligand patterns are considered: striated NPs with alternating hydrophilic and hydrophobic groups compared to NPs with randomly mixed ligands at the same hydrophilic to hydrophobic ratio, as well as NPs coated with homogeneous hydrophilic or hydrophobic ligands. The free energy analysis indicates that among the four ligand patterns under study, the striated NP encounters the lowest energy barrier during translocation across the membrane. Further analysis reveals that the translocation of the striated NP is facilitated by the constraint of its rotational degree of freedom by the anisotropic ligand pattern, which prevented the free energy of the system from sinking to a deeper valley as the NP passes through the hydrophobic core of the bilayer. Finally, the critical forces required for almost instant penetration of these patterned NPs across the bilayer are calculated and shown to be consistent with the free energy analysis. These findings provide useful guidelines for the molecular design of patterned NPs for controllable cell penetrability.

We review recent advances in multiscale modeling of the mechanics of healthy and diseased red blood cells (RBCs), and blood flow in the microcirculation. We cover the traditional continuum-based methods but also particle-based methods used to model both the RBCs and the blood plasma. We highlight examples of successful simulations of blood flow including malaria and sickle cell anemia.

This study is partially motivated by the validation of a new two-component multi-scale cell model we developed recently that treats the lipid bilayer and the cytoskeleton as two distinct components...

The motion of a suspension of red blood cells (RBCs) flowing in a Y-shaped bifurcating microfluidic channel is investigated using a validated low-dimensional RBC model based on dissipative particle dynamics. Specifically, the RBC is represented as a closed torus-like ring of ten colloidal particles, which leads to efficient simulations of blood flow in microcirculation over a wide range of hematocrits. Adaptive no-slip wall boundary conditions were implemented to model hydrodynamic flow within a specific wall structure of diverging three-dimensional microfluidic channels, paying attention to controlling density fluctuations. Plasma skimming and the all-or-nothing phenomenon of RBCs in a bifurcating microfluidic channel have been investigated in our simulations for healthy and diseased blood, including the size of a cell-free layer on the daughter branches. The feed hematocrit level in the parent channel has considerable influence on blood-plasma separation. Compared to the blood-plasma separation efficiencies of healthy RBCs, malaria-infected stiff RBCs (iRBCs) have a tendency to travel into the low flow-rate daughter branch because of their different initial distribution in the parent channel. Our simulation results are consistent with previously published experimental results and theoretical predictions.