Patient-specific modeling of individual sickle cell behavior under transient hypoxia

Abstract

Sickle cell disease (SCD) is a highly complex genetic blood disorder in which red blood cells (RBC) exhibit heterogeneous morphology changes and decreased deformability. We employ a kinetic model for cell morphological sickling that invokes parameters derived from patient-specific data. This model is used to investigate the dynamics of individual sickle cells in a capillary-like microenvironment in order to address various mechanisms associated with SCD. We show that all RBCs, both hypoxia-unaffected and hypoxia-affected ones, regularly pass through microgates under oxygenated state. However, the hypoxia-affected cells undergo sickling which significantly alters cell dynamics. In particular, the dense and rigid sickle RBCs are obstructed thereby clogging blood flow while the less dense and deformable ones are capable of circumnavigating dead (trapped) cells ahead of them by choosing a serpentine path. Informed by recent experiments involving microfluidics that provide in vitro quantitative information on cell dynamics under transient hypoxia conditions, we have performed detailed computational simulations of alterations to cell behavior in response to morphological changes and membrane stiffening. Our model reveals that SCD exhibits substantial heterogeneity even within a particular density-fractionated subpopulation. These findings provide unique insights into how individual sickle cells move through capillaries under transient hypoxic conditions, and offer novel possibilities for designing effective therapeutic interventions for SCD.