Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
J Biomech. Author manuscript; available in PMC 2017 March 27.
Published in final edited form as:
PMCID: PMC5368081

Biomechanics and biorheology of red blood cells in sickle cell anemia


Sickle cell anemia (SCA) is an inherited blood disorder that causes painful crises due to vaso-occlusion of small blood vessels. The primary cause of the clinical phenotype of SCA is the intracellular polymerization of sickle hemoglobin resulting in sickling of red blood cells (RBCs) in deoxygenated conditions. In this review, we discuss the biomechanical and biorheological characteristics of sickle RBCs and sickle blood as well as their implications toward a better understanding of the pathophysiology and pathogenesis of SCA. Additionally, we highlight the adhesive heterogeneity of RBCs in SCA and their specific contribution to vaso-occlusive crisis.

Keywords: Sickle cell disease, Erythrocyte, Biomechanics, Biorheology, Adhesion

1. Introduction

Sickle cell anemia (SCA), the first identified “molecular disease” (Pauling et al., 1949; Strasser, 1999), is one of the most common genetic inherited hematological disorders, which can cause several types of chronic complications such as vaso-occlusive crisis (VOC), hemolytic anemia, and sequestration crisis (Bunn, 1997; Barabino et al., 2010). The pathogenesis of vaso-occlusion involves several processes across multiple time and length scales, from O(101s) to O(103s) for the kinetics of HbS polymerization to the hemodynamics of sickle blood flow, and from O(109m) to O(105m) for the size of the protein to the dimensions of the microcirculatory vessels. During the past few decades, different aspects of this disease have been successfully investigated (Barabino et al., 2010; Steinberg, 1999; Frenette and Atweh, 2007; Yazdani et al., 2016). At the molecular scale, the HbS polymerization process has been characterized by a double nucleation mechanism. At the cellular scale, sickle RBCs are characterized by remarkable heterogeneity in density, morphology, and rigidity. The affected RBCs become more rigid and “sticky” compared to normal RBCs, causing frequent vaso-occlusive episodes and depriving tissues and organs of oxygen. At the microvascular scale, early studies postulated that the HbS polymerization resulted in the entrapment of sickle RBCs in capillaries (Fig. 1) (Kaul et al., 1986), and subsequent studies further revealed the multi-interactional and multi-stage nature of the VOC (Kaul et al., 1989, 1994, 2009; Kaul and Fabry, 2004).

Fig. 1
Vaso-occlusive crisis in SCA

Currently, hydroxyurea (HU) is the only approved medication in widespread use for the treatment of SCA (Ware, 2010). The treatment of SCA patients with HU has the following beneficial effects: (i) increased production of fetal hemoglobin (HbF) and therefore increased delay time of the RBC sickling process (Bridges et al., 1996; Atweh and Schechter, 2001), (ii) reduction of white blood cell (WBC) count and expression pattern of cellular adhesion molecules (Charache et al., 1996), and (iii) reduction in the frequencies of blood transfusion (Ware et al., 1999). These beneficial effects ameliorate the severity of SCA. However, clinical studies report that HU is ineffective for many patients for unclear reasons (Manwani and Frenette, 2013). Moreover, the aforementioned studies indicate that the clinical expression of SCA is heterogeneous, making it hard to predict the risk of VOC, resulting in a serious challenge for disease management. Here, we review experimental studies and predictive simulations related to biomechanical and biorheological properties as well as heterogeneity-related issues associated with SCA.

2. Biomechanical and biorheological properties of sickle RBCs

Quantification of the biomechanical and biorheological characteristics of RBCs can improve our understanding of the etiology of a number of human diseases. In SCA, partial deoxygenation of sickle RBCs in post-capillary venules causes HBS polymerization followed by possible RBC sickling. Repeated RBC sickling can result in the development of defects in the RBC membrane, reduced RBC deformability, increased time of RBC adherence to venules, and eventually in vaso-occlusion.

2.1. Sickle cell biomechanics

Over the past few decades researchers investigated the biomechanics of sickle RBCs as indicators of the severity of the disease. The available experimental methods can measure the biomechanical properties of a large number of sickle RBCs at the same time (Chien et al., 1970; Messmann et al., 1990; Connes et al., 2014), or isolated sickle RBCs (Byun et al., 2012; Maciaszek and Lykotrafitis, 2011). For example, early studies using filtration (Chien et al., 1970) or ektacytometry (Messmann et al., 1990) directly examined the biomechanical properties of the sickle RBC membrane and determined that sickle RBCs are less deformable than normal RBCs. In a recent study, decreased RBC deformability and aggregation, measured using ektacytometry and laser backscatter of Percoll-separated sickle RBCs, have been shown to correlate with hemolysis (Connes et al., 2014). However, these techniques measure properties averaged over all RBCs in a blood sample, without regard to the cell heterogeneity within sickle blood sample. Single-cell experimental methods include micropipette aspiration, optical tweezers, flickering analysis, atomic force microscopy (AFM), diffraction phase microscopy, and recently, microfluidics and ultrasounds. The optical tweezers and micropipette aspiration techniques subject the RBC directly to mechanical deformation and yield shear modulus of sickle RBCs in the range of 8–20 μN m−1 (Fig. 2) (Byun et al., 2012). AFM measurements have found that the Young’s modulus of SCA RBCs are stiffer than normal RBCs with a widely distributed Young’s modulus ranged from 3 kPa to 50 kPa depending on the hypoxic conditions and probably on the clinical severity of the disease (Maciaszek and Lykotrafitis, 2011). The stiffening of sickle RBC membrane may indicate the effect of the polymerization of HbS as well as the possible remodeling of cytoskeleton associated with SCA.

Fig. 2
Biomechanical properties of RBCs in health and in SCA

Advances in mathematical models and computational simulations enable investigation of a broad range of biomechanical problems associated with RBCs in SCA. For example, Dong et al. (1992) developed a mathematical model of RBC flowing in narrow vessels. They showed that the RBCs become stiffer when the amount of intracellular HbS polymer increases. Hemolysis is associated with irreversible structural change of sickle RBCs (Kato et al., 2013). Therefore, the individuals with a higher numbers of irreversibly sickled cells (ISCs) are at greater risk for hemolysis (Serjeant et al., 1969). Fisseha and Katiyar (2012) employed a generalized Voigt-model of nonlinear viscoelastic solids to characterize the viscoelastic properties of sickle RBCs. They found that the ISCs with irreversible alteration in cell membrane structure tend to hemolysis.

2.2. Sickle cell biorheology

Sickle RBCs have increased cell rigidity and decreased cell deformability, causing hemolysis and abnormal hemorheology in SCA (Chien et al., 1970; Usami et al., 1975). The rheological abnormalities are caused primarily by an increase in cytoplasm viscosity due to HbS polymerization upon deoxygenation (DeOxy), as well as biochemical abnormalities in the sickle RBC membrane. The abnormal rheological changes is also related to the cell membrane stiffening due to repeated cycles of polymerization and depolymerization in the circulating RBCs.

Sickle RBCs are heterogeneous in their rheological characteristics. According to Kaul et al. (1983), sickle blood is fractionated into 4 density groups (fractions I–IV). Fractions I (SS1) and II (SS2) are composed primarily of reticulocytes and discocytes, respectively, with mean corpuscular hemoglobin concentration (MCHC) levels similar to healthy RBCs, resulting in a comparable bulk viscosity to that of unseparated healthy blood samples in an oxygenated state (Kaul et al., 1983). Fractions III (SS3) and IV (SS4) are mainly composed of rigid discocytes and ISCs, respectively, with MCHC values considerably higher than those of healthy RBCs, which results in a significant increase in blood viscosity, even under oxygenated state (Kaul et al., 1983).

Recent experimental studies have captured information among all of the aforementioned processes in physiologic regimes, providing insight into the overall dynamics of a vaso-occlusive event (Du et al., 2015; Higgins et al., 2007; Wood et al., 2012; Lu et al., 2016b). In these experiments, whole blood or RBC suspensions from SCA patients flowed through microfluidic channels under constant pressure and low the oxygen concentration, thereby simulating most basic features of a vaso-occlusive event. For example, Higgins and his colleagues found that oxygen tension regulate blood flow for individual SCA patients (Higgins et al., 2007; Wood et al., 2012; Lu et al., 2016b). Recently, Du et al. (2015) developed a high-throughput microfluidics-based model to investigate the sickle cell behavior under transient hypoxia (Fig. 3). Using this microfluidic device to measure blood samples from 25 SCA patients, they quantify the kinetics of cell sickling, unsickling, and individual cell rheology associated with SCA.

Fig. 3
Microfluidic platform for studying sickle cell behavior under transient hypoxic conditions

A change in blood rheological properties is usually linked to hematological diseases; therefore, the viscosity of blood has long been used as an indicator for understanding the implications and treatment routes of these type of diseases. Although various experiments have been performed for blood viscosity measurement in SCA, considerable uncertainty exists with respect to the effect of HU treatment on the viscosity of sickle blood: it has been reported that the viscosity may increase (Fattori et al., 2005), decrease (Ware, 2010) or even remain unchanged (Lemonne et al., 2015) for patients with SCA after treatment with HU.

Computational methods have been used to evaluate biorheological properties, to help explain the complex and multifactorial nature of SCA pathogenesis. For example, Dupin et al. (2008) employed a lattice Boltzmann method to simulate the sickle RBC suspension flowing through an aperture of diameter slightly less than the size of a single RBC. They found that the abnormally shaped and less deformable sickle RBCs get stuck upstream of the constriction. In recent years, a particle-based multiscale RBC (MS-RBC) model has been developed to simulate RBCs in disease like SCA (Lei and Karniadakis, 2013; Li et al., 2014, 2016; Chang et al., 2016). In this MS-RBC model, the cell membrane is represented by a triangulated network in two dimensions, which are connected by Ns edges and Nt triangles (Li et al., 2005; Pivkin and Karniadakis, 2008; Fedosov et al., 2010). The potential energy of the system includes in-plane elastic energy Vs, bending energy Vb, along with energies, Va+v, from surface area and volume conservation. These energies are represented by




where kb is the bending coefficient, lj and lm are the stretched length and contour length of spring j, θj and θ0 are the instantaneous angle and spontaneous angle between two adjacent triangles. Also, Aj is the instantaneous triangle area, and A0 is the initial triangle area. Vtot is the total RBC volume, while V0tot is the initial total RBC volume. The MS-RBC model represents seamlessly the healthy and infected RBCs as well as the plasma flow, therefore, it has been widely used to study the behavior of normal and abnormal RBCs in microfluidics. For example, Lei and Karniadakis (2012b) quantified the hemodynamics of sickle RBC suspensions under various physiological conditions. They showed that the deoxygenated sickle RBCs exhibit “solid” behavior and the viscosity of deoxygenated sickle RBC suspension is nearly shear-independent throughout the entire shear rate regime. Li et al. (2016) examined the effect of the cell morphology on the rheology of sickle RBC suspension under shear. They showed that the abnormal rheological properties of sickle RBCs are correlated with the cell morphology (Fig. 4).

Fig. 4
Shear viscosity of sickle cell suspension with different cell morphology at different shear rate

It is known that the origin of SCA can be traced to a common molecular basis, but individual patients with SCA have a highly variable clinical phenotype. For these reasons, Li et al. (2016) have recently developed a predictive patient-specific model of SCA at the molecular level to quantify the rheological behavior of blood flow in SCA. They determined the shear viscosity of blood from SCA patients with HU treatment and those without HU treatment. Their results demonstrated that treatment with HU can indeed improve the blood flow in SCA.

3. Adhesive properties of sickle RBCs

As noted above, cell sickling is necessary but not sufficient to initiate a VOC. Abnormal adhesive properties of sickle RBCs, including activation of known adhesion receptors and increased interactions with WBCs, platelets, ECs, and extracellular matrix proteins, have drawn intense attention as potential initiating factors in VOC (Hillery et al., 1996; Hebbel et al., 2004). Here we overview significant studies and contributions to this field.

3.1. Adhesive dynamics

Increased adhesive forces between sickle RBCs and ECs have been hypothesized to play a role in the initiation of vaso-occlusion in SCA. In vitro studies have shown that sickle RBCs exhibit heterogeneous cell adhesivity among different density groups (Barabino et al., 2010, 1987). The dense ISCs were found to preferentially adhere to ECs in static assays (Wautier et al., 1985), but other researchers demonstrated that the light-density RBCs adhere more avidly to ECs than ISC-rich dense group under flow conditions mimicking those in postcapillary venules (Barabino et al., 1987). It is likely that the irregular geometry and decreased deformability of dense ISCs prevent the type of contact needed to promote adhesion under shear forces induced by fluid flow. In a recent study, Alapan et al. (2014) presented a microfluidic approach to study the abnormal adhesive properties of individual sickle RBCs in physiological flow conditions. They demonstrated that the most tightly adherent cells are non-deformable sickle RBCs as they exhibit increased adhesion sites compared to the deformable ones.

The adhesion of sickle RBCs is also protein-dependent (Lee et al., 1998; Joneckis et al., 1996). Previous studies have shown that the laminin receptor on sickle RBCs is more available in denser groups (Lee et al., 1998), while thrombospondin supports higher levels of adhesion of light-density sickle RBCs (Joneckis et al., 1996). Epinephrine is a hormone that is increased in the midst of stress. It has been shown that adhesion of sickle RBCs to the endothelium increases in the presence of epinephrine (Telen, 2005; Hines et al., 2003; Zennadi et al., 2004, 2007; Maciaszek et al., 2014). Epinephrine stimulates G-protein coupled receptors which activate the cyclic AMP-protein kinase A (cAMP-PKA) controlled pathway through stimulation of adenylate cyclase (AC) which in turn activates RBC surface adhesion receptors such as BCAM/Lu and ICAM-4 (Hines et al., 2003; Zennadi et al., 2004). Maciaszek et al. (2014) demonstrated that activation of these adhesion receptors is mediated by the scaffolding A-kinase anchoring proteins (AKAPs).

Some studies also probe the effect of HU therapy on the adhesive properties of sickle RBCs (Hillery et al., 2000; Gambero et al., 2007). Therapy with HU has been shown to decrease sickle RBC adhesion and downregulates endothelial adhesion molecules such as thrombospondin, laminin, and fibronectin. For example, Gambero et al. (2007) used a static adhesion assay to analyze the basal adhesion of RBCs, from SCA patients with HU treatment and those without HU treatment, to fibronectin. They found that HU therapy is associated with reduced adhesion molecule gene and protein expression in sickle RBCs with a concomitant reduction in adhesive properties. Odièvre et al. (2008) suggest that HU may modulate the activation of adhesion receptors by directly acting on gene expression and the related signaling cascade. Maciaszek et al. (2014) investigated RBC–EC interactions and their regulation via cAMP using single-molecule AFM experiments. They showed that HU treatment results in a reduced adhesion between sickle RBCs and ECs due to lower receptor expression. Hence, overall, the HU therapy reduce the complications of SCA.

Computational models have been employed to quantify the adhesive properties of sickle RBCs (Lei and Karniadakis, 2012b, 2013). For example, Lei and Karniadakis (2013) examined the adhesive dynamics of sickle RBCs of different density groups. Given the same “adhesive potential”, their results validate the hypothesis that heterogeneous cell adhesive dynamics is mainly due to the different cell rigidities and peculiar cell morphologies (Kaul et al., 1994). A deformable SS2 cell exhibits firm adhesion to the surface with a large contact area, a rigid SS3 cell shows weak adhesivity, while an ISC does not show any adhesion to the surface; instead, it moves freely without adhesive bonds established thereafter. In addition, they computed the adhesive force between the surface and sickle RBCs. They found that the adhesive force exhibits an inverse relationship with the cell rigidity: compared with the discocyte, the ISC exhibits smaller adhesive force given the similar cell rigidity, indicating less adhesivity induced by its peculiar cell morphology.

3.2. Vaso-occlusive crisis

VOC is the hallmark of the biophysical characteristics related to SCA. Typically, it has been viewed as synonymous with the blockage of single sickle/elongated RBC in capillaries (Kaul et al., 1986). Subsequent studies revealed that VOC is a multistep and multicellular paradigm (Kaul et al., 1989, 1994; Kaul and Fabry, 2004). Several models of vaso-occlusion have been proposed. For example, in a two-step model proposed by Kaul and his colleagues, preferential adhesion of deformable sickle RBCs in post-capillary venules that is followed by selective trapping of dense sickled RBCs could result in vaso-occlusion; alternatively, WBC adhesion in inflamed venules could trigger the selective trapping (Kaul and Fabry, 2004). An alternate model is proposed by Frenette (2002), in which sickle cell vaso-occlusion is viewed as a multistep and multicellular process driven by inflammatory stimuli and the vaso-occlusion is mediated by a number of cellular elements in four steps: (i) endothelial activation, (ii) WBC adherence to endothelium, (iii) sickle RBC interactions with adherent WBCs, and (iv) progressive blockage.

Computational modeling and simulations have also been used to tackle the multi-stage and multi-interactional nature of VOC. For example, Lei and Karniadakis (2013) employed a multiscale model of sickle RBC to investigate the biophysical characteristics of the vaso-occlusion in SCA (Fig. 5). They quantified the specific physiological conditions triggering the vaso-occlusion crisis and identified the specific contribution of individual cell groups within the vaso-occlusion process. They found that under physiological conditions similar to microcirculation in post-capillaries, the interplay of deformable SS2 cells and ISCs can potentially trigger full blood occlusion. To probe the WBC recruitment and its effect to blood vaso-occlusion, Lei and Karniadakis employed the multiscale model of sickle RBC and WBC to simulate the cell suspensions (mixed sickle RBCs and WBCs) in a microtube (Fig. 5b). Starting from the initial condition, they found that the WBC migrates toward to the tube wall with firm adhesion, leading to a decrease in blood flow. The adherent WBC further interacts and traps the sickle RBCs, resulting in full occlusion. For venular flow with larger diameter, multiple WBC recruitment may occur in the inflammation regions. They also simulated the sickle RBC suspension with multiple WBCs in a larger microtube. They found that the blood flow undergoes slow down due to the WBC recruitment and the moderate sickle RBC–WBC interaction leads to multiple sickle RBC trapped on the WBCs and the full occlusion (Fig. 5b).

Fig. 5
Vaso-occlusion in post-capillaries

4. Hemoglobin properties

SCA is a molecular disease that affects hemoglobin (Hb), the molecule in RBCs that delivers oxygen to cells throughout the human body. Many different types of Hb exist. The most common ones are normal HbA (α2β2), HbA2 (α2δ2), HbS(α2β2S), HbC(α2β2C), and HbF (α2γ2). Normal RBCs contain HbA composed of 2 α-globin subunits and two β-globin subunits; the pairing of one α subunit and one β subunit produces a Hb dimer (α1β1). Two dimers combine to form a Hb tetramer (α2β2), the functional form of Hb. In SCA, the Glu-6 (β)→ Val mutation induces a hydrophobic interaction two neighboring HbS molecules (Bunn, 1997). HbS is less soluble than other forms of Hb under decreasing oxygen concentrations, leading to polymerization into long polymer fibers that deform the RBC shapes.

4.1. Sickle hemoglobin in SCA

HbS nucleation followed by polymerization and RBC sickling significantly contributes to vaso-occlusion, which is the hallmark of SCA. The polymerization of HbS was modeled as a double nucleation process in Ferrone et al. (1985b) and Vekilov (2007). This nucleation is followed by the growth and alignment of HbS polymer fibers, transforming the cell into the classic sickle shape (Ferrone et al., 1985a,b; Samuel and Briehl, 1990; Vekilov, 2007; Mozzarelli et al., 1987). A different two-step mechanism of HbS polymerization has been suggested by other groups (ten Wolde and Frenkel, 1997; Shiryayev and Gunton, 2004; Lutsko and Nicolis, 2006): (i) formation of dense liquid droplets, and (ii) formation of HbS fiber nuclei within the droplets. Consequently, HbS polymer fibers grow spontaneously and distort the RBCs into sickle shapes.

Once the Hb releases its oxygen, in SCA there is an ensuing molecular chain reaction between HbS molecules to form long polymer chains that are rigid. However, if the RBC is restrained inside a narrow capillary when this reaction begins, its shape is physically restrained from growing outward. In a recent study, Aprelev et al. (2012) employed a microfluidic method to investigate the mechanical interaction between sickled RBCs and capillaries. They constructed a single-cell microfluidic channel to identify the physical forces in RBCs and blood vessels underlying the painful symptoms of SCA. They found that the rigid elongated RBCs do not get stuck in narrow capillaries, and proposed that the timing of polymerization inside sickle RBCs may be important to understanding patients’ susceptibility to symptoms.

Numerical models provide valuable insights into the mechanism of HbS fiber nucleation, polymerization and fiber growth. For example, Turner et al. (2003) proposed a theoretical model to study the thermodynamic stability of HbS of polymer fibers, demonstrating that twist plays an important essential role in stabilizing HbS polymer fibers. Yang et al. (2010) performed dynamic simulations of self-assembled filamentous bundles. They found that chain chirality can be used to control the HbS polymer bundle size. Li et al. (2012b) employed a coarse-grained HbS model to simulate the polymerization of HbS molecules. They demonstrated that the molecular chirality is the main driver for the HbS polymer fiber formation. Li and Lykotrafitis simulated the biomechanical properties of HbS polymer fibers via two different coarse-grained HbS fiber models, demonstrating the significant role of fiber frustration and compression in fiber zipping and unzipping dynamics (Li and Lykotrafitis, 2011; Li et al., 2012a). To gain insight into the nature of the formation of HbS polymer fiber, in a recent study, Lu et al. (2016a) developed a patchy particle HbS model to study the growth dynamics of HbS polymer fibers (Fig. 6a). They demonstrated that the formation process of HbS polymer fiber occurs through monomer addition. In addition, they found that the molecular chirality is the critical determinant of the mechanical and structural properties of HbS polymer fibers.

Fig. 6
Multiscale modeling of SCA

4.2. Fetal hemoglobin in SCA

HbF is composed of 2 α-globin subunits and 2 γ-globin subunits. Unlike HbA, HbF actively inhibits the polymerization of HbS, and as a result it reduces the severity of the disease. SCA patients with high HbF levels not only have less severe clinical course, but also show milder clinical complications, as increased production of HbF can reduce the occurrence of sickling-related complications (Akinsheye et al., 2011). Hydroxyurea promotes the production of HbF and can thus be used to treat SCA (Du et al., 2015; Charache et al., 1992; Ferster et al., 2001). In a mixed solution of HbF and HbS, there is a reduction of effective HbS concentration and, more importantly, both HbF and its mixed hybrid tetramer (α2βSγ) cannot enter the HbS polymer, making the RBCs less likely to sickle. In this way a small amount of HbF can reduce the ease with which sickling occurs, and thus ameliorate the disease. In solution, HbF concentration higher than 15% prevents HbS polymerization (Bunn, 1987).

5. Summary and outlook

In this article, we review recent advances in probing and understanding the dynamics of collective processes associated with vaso-occlusion that links together sub-cellular, cellular, and vessel phenomena. We cover the biomechanical, biorheological and adhesive properties of sickle RBCs at cell level and abnormal hemoglobin properties at molecular level. The underlying molecular cause of the disease has been understood for more than half a century; however, progress in developing treatments to prevent painful VOC and the other myriad of associated symptoms has been slow (Parise and Berliner, 2016). Therefore, the need to develop new therapies or improve the existing treatments for SCA remains paramount. For these reasons, it is necessary to have a better understanding of the pathogenesis and pathobiomechanics of SCA and a more accurate evaluation of efficacy and safety of current anti-adhesive and deformability-restoring drugs in subsequent studies. Specifically, it is important to distinguish the effects of the current drug treatments, including HU, statins and adhesion-preventing drugs, from enhanced deformability and reduced adhesion, and examine what adhesive molecular mechanisms are influenced by other potential treatments, through which clinical inventions could potentially be designed and evaluated more effectively. This would potentially facilitate the design and evaluation of new clinical therapeutic approaches and interventions.

Aside from the intensive experimental studies, several computational approaches, including continuum- and particle-based methods, have been developed and applied to investigate a broad range of biomechanical and rheological problems associated with healthy and pathological RBCs at different length and time scales. Traditionally, continuum-based RBC models treat the RBC membrane and intracellular fluids as homogeneous materials, allowing the simulations of large-scale blood flow. However, they cannot describe the structural alteration of the RBC membrane in many hematological disorders such as the uncoupling of cytoskeleton from the lipid bilayer in SCA. In recent years, particle-based RBC models are increasingly popular as a promising tool for modeling of structural, biomechanical and biorheological properties of RBCs in disease. However, they are computationally expensive to scale up to large domains. From the perspective of computational modeling and simulations, the present computational framework can be further extended to investigate the following important questions related to SCA: (i) Development of a hybrid model that integrates the present particle-based RBC models at different levels (Fig. 6) with a continuum description of oxygen transport. Ultimately, computations that encompass all molecular and cellular scales could be used to investigate the entire process from deoxygenation to HbS fiber nucleation and then to vaso-occlusion in SCA. Such simulations would potentially answer questions concerning the links among HbS polymerization, cell sickling, blood flow alteration, and eventually VOC; (ii) Extension of the patient-specific model of SCA by including a macroscopic, continuum description of diffusion of chemicals (e.g., HU, decitabine, erythropoietin) to quantify the therapeutic effects of drug treatments on the microcirculation of the blood flow in SCA. These numerical studies, combined with microfluidic experiments can be used to evaluate and quantify related therapeutic treatment and clinical outcomes.


X.L., M.D., G.K. acknowledge support by the NIH grant U01HL114476. M.D. acknowledges partial support from the Singapore-MIT Alliance for Research and Technology (SMART) Center. G.L. acknowledges support by AHA grant 12SDG12050688, and by NSF grants CMMI-1235025 and PHY-1205910.


Conflict of interest statement

None declared.


  • Akinsheye I, Alsultan A, Solovieff N, Ngo D, Baldwin CT, Sebastiani P, Chui DH, Steinberg MH. Fetal hemoglobin in sickle cell anemia. Blood. 2011;118:19–27. [PubMed]
  • Alapan Y, Little JA, Gurkan UA. Heterogeneous red blood cell adhesion and deformability in sickle cell disease. Sci Rep. 2014;4:7173. [PMC free article] [PubMed]
  • Aprelev A, Stephenson W, Noh H, Meier M, Ferrone FA. The physical foundation of vasoocclusion in sickle cell disease. Biophys J. 2012;103:L38–L40. [PubMed]
  • Atweh GF, Schechter AN. Pharmacologic induction of fetal hemoglobin: raising the therapeutic bar in sickle cell disease. Curr Opin Hematol. 2001;8:123–130. [PubMed]
  • Barabino GA, McIntire LV, Eskin SG, Sears DA, Udden M. Endothelial cell interactions with sickle cell, sickle trait, mechanically injured, and normal erythrocytes under controlled flow. Blood. 1987;70:152–157. [PubMed]
  • Barabino GA, Platt MO, Kaul DK. Sickle cell biomechanics. Annu Rev Biomed Eng. 2010;12:345–367. [PubMed]
  • Bridges KR, Barabino GD, Brugnara C, Cho MR, Christoph GW, Dover G, Ewenstein BM, Golan DE, Guttmann CR, Hofrichter J, Mulkern RV, Zhang B, Eaton WA. A multiparameter analysis of sickle erythrocytes in patients undergoing hydroxyurea therapy. Blood. 1996;88:4701–4710. [PubMed]
  • Bunn HF. Subunit assembly of hemoglobin: an important determinant of hematologic phenotype. Blood. 1987;69:1–6. [PubMed]
  • Bunn HF. Pathogenesis and treatment of sickle cell disease. N Engl J Med. 1997;337:762–769. [PubMed]
  • Byun H, Hillman TR, Higgins JM, Diez-Silva M, Peng Z, Dao M, Dasari RR, Suresh S, Park Y. Optical measurement of biomechanical properties of individual erythrocytes from a sickle cell patient. Acta Biomater. 2012;8:4130–4138. [PMC free article] [PubMed]
  • Chang HY, Li XJ, Li H, Karniadakis GE. MD/DPD multiscale framework for predicting morphology and stresses of red blood cells in health and disease. PLOS Comput Biol. 2016;12:e1005173. [PMC free article] [PubMed]
  • Charache S, Barton FB, Moore RD, Terrin ML, Steinberg MH, Dover GJ, Ballas SK, McMahon RP, Castro O, Orringer EP. Hydroxyurea and sickle cell anemia – clinical utility of a myelosuppressive “switching” agent: the multicenter study of hydroxyurea in sickle cell anemia. Medicine. 1996;75:300–326. [PubMed]
  • Charache S, Dover GJ, Moore RD, Eckert S, Ballas SK, Koshy M, Milner P, Orringer EP, Phillips G, Platt O, Thomas GH. Hydroxyurea: effects on hemoglobin F production in patients with sickle cell anemia. Blood. 1992;79:2555–2565. [PubMed]
  • Chien S, Usami S, Bertles JF. Abnormal rheology of oxygenated blood in sickle cell anemia. J Clin Invest. 1970;49:623–634. [PMC free article] [PubMed]
  • Chien S, Usami S, Taylor HM, Lundberg JL, Gregersen MI. Effects of hematocrit and plasma proteins on human blood rheology at low shear rates. J Appl Phys. 1966;21:81–87. [PubMed]
  • Connes P, Lamarre Y, Waltz X, Ballas SK, Lemonne N, Etienne-Julan M, Hue O, Hardy-Dessources MD, Romana M. Haemolysis and abnormal hae-morheology in sickle cell anaemia. Br J Haematol. 2014;165:564–572. [PubMed]
  • Dong C, Chadwick R, Schechter A. Influence of sickle hemoglobin poly-merization and membrane properties on deformability of sickle erythrocytes in the microcirculation. Biophys J. 1992;63:774–783. [PubMed]
  • Du E, Diez-Silva M, Kato GJ, Dao M, Suresh S. Kinetics of sickle cell biorheology and implications for painful vasoocclusive crisis. Proc Natl Acad Sci USA. 2015;112:1422–1427. [PubMed]
  • Dupin MM, Halliday I, Care CM, Munn LL. Lattice Boltzmann modelling of blood cell dynamics. Int J Comput Fluid D. 2008;22:481–492.
  • Fattori A, de Souza RA, Saad ST, Costa FF. Acute myocardial infarction in sickle cell disease: a possible complication of hydroxyurea treatment. Hematol J. 2005;5:589–590. [PubMed]
  • Fedosov DA, Caswell B, Karniadakis GE. A multiscale red blood cell model with accurate mechanics, rheology, and dynamics. Biophys J. 2010;98:2215–2225. [PubMed]
  • Ferrone FA, Hofrichter J, Eaton WA. Kinetics of sickle hemoglobin polymerization. I. Studies using temperature-jump and laser photolysis techniques. J Mol Biol. 1985a;183:591–610. [PubMed]
  • Ferrone FA, Hofrichter J, Eaton WA. Kinetics of sickle hemoglobin polymerization. II. A double nucleation mechanism. J Mol Biol. 1985b;183:611–631. [PubMed]
  • Ferster A, Tahriri P, Vermylen C, Sturbois G, Corazza F, Fondu P, Devalck C, Dresse MF, Feremans W, Hunninck K, Toppet M, Philippet P, Van Geet C, Sariban E. Five years of experience with hydroxyurea in children and young adults with sickle cell disease. Blood. 2001;97:3628–3632. [PubMed]
  • Fisseha D, Katiyar VK. Analysis of mechanical behavior of red cell membrane in sickle cell disease. Appl Math. 2012;2:40–46.
  • Frenette PS. Sickle cell vaso-occlusion: multistep and multicellular paradigm. Curr Opin Hematol. 2002;9:101–106. [PubMed]
  • Frenette PS, Atweh GF. Sickle cell disease: old discoveries, new concepts, and future promise. J Clin Invest. 2007;117:850–858. [PMC free article] [PubMed]
  • Gambero S, Canalli AA, Traina F, Albuquerque DM, Saad ST, Costa FF, Con-ran N. Therapy with hydroxyurea is associated with reduced adhesion molecule gene and protein expression in sickle red cells with a concomitant reduction in adhesive properties. Eur J Haematol. 2007;78:144–151. [PubMed]
  • Hebbel RP, Osarogiagbon R, Kaul D. The endothelial biology of sickle cell disease: inflammation and a chronic vasculopathy. Microcirculation. 2004;11:129–151. [PubMed]
  • Higgins JM, Eddington DT, Bhatia SN, Mahadevan L. Sickle cell vasoocclusion and rescue in a microfluidic device. Proc Natl Acad Sci USA. 2007;104:20496–20500. [PubMed]
  • Hillery C, Du M, Montgomery R, Scott J. Increased adhesion of erythrocytes to components of the extracellular matrix: isolation and characterization of a red blood cell lipid that binds thrombospondin and laminin. Blood. 1996;87:4879–4886. [PubMed]
  • Hillery CA, Du MC, Wang WC, Scott JP. Hydroxyurea therapy decreases the in vitro adhesion of sickle erythrocytes to thrombospondin and laminin. Br J Haematol. 2000;109:322–327. [PubMed]
  • Hines PC, Zen Q, Burney SN, Shea DA, Ataga KI, Orringer EP, Telen MJ, Parise LV. Novel epinephrine and cyclic AMP-mediated activation of BCAM/Lu-dependent sickle RBC adhesion. Blood. 2003;101:3281–3287. [PubMed]
  • Joneckis C, Shock D, Cunningham M, Orringer E, Parise L. Glycoprotein IV-independent adhesion of sickle red blood cells to immobilized thrombos-pondin under flow conditions. Blood. 1996;87:4862–4870. [PubMed]
  • Kato GJ, Nouraie SM, Gladwin MT. Lactate dehydrogenase and hemolysis in sickle cell disease. Blood. 2013;122:1091–1092. [PubMed]
  • Kaul D, Chen D, Zhan J. Adhesion of sickle cells to vascular endothelium is critically dependent on changes in density and shape of the cells. Blood. 1994;83:3006–3017. [PubMed]
  • Kaul D, Fabry M, Nagel R. Vaso-occlusion by sickle cells: evidence for selective trapping of dense red cells. Blood. 1986;68:1162–1166. [PubMed]
  • Kaul DK, Fabry M, Windisch P, Baez S, Nagel RL. Erythrocytes in sickle cell anemia are heterogeneous in their rheological and hemodynamic characteristics. J Clin Invest. 1983;72:22–31. [PMC free article] [PubMed]
  • Kaul DK, Fabry ME. In vivo studies of sickle red blood cells. Microcirculation. 2004;11:153–165. [PubMed]
  • Kaul DK, Fabry ME, Nagel RL. Microvascular sites and characteristics of sickle cell adhesion to vascular endothelium in shear flow conditions: pathophysiological implications. Proc Natl Acad Sci USA. 1989;86:3356–3360. [PubMed]
  • Kaul DK, Finnegan E, Barabino GA. Sickle red cell-endothelium interactions. Microcirculation. 2009;16:97–111. [PMC free article] [PubMed]
  • Kaul DK, Xue H. Rate of deoxygenation and rheologic behavior of blood in sickle cell anemia. Blood. 1991;77:1353–1361. [PubMed]
  • Lee SP, Cunningham ML, Hines PC, Joneckis CC, Orringer EP, Parise LV. Sickle cell adhesion to laminin: potential role for the α5 chain. Blood. 1998;92:2951–2958. [PubMed]
  • Lei H, Karniadakis GE. Predicting the morphology of sickle red blood cells using coarse-grained models of intracellualr aligned hemoglobin polymers. Soft Matter. 2012a;8:4507–4516. [PMC free article] [PubMed]
  • Lei H, Karniadakis GE. Quantifying the rheological and hemodynamic characteristics of sickle cell anemia. Biophys J. 2012b;102:185–194. [PubMed]
  • Lei H, Karniadakis GE. Probing vasoocclusion phenomena in sickle cell anemia via mesoscopic simulations. Proc Natl Acad Sci USA. 2013;110:11326–11330. [PubMed]
  • Lemonne N, Charlot K, Waltz X, Ballas SK, Lamarre Y, Lee K, Hierso R, Connes C, Etienne-Julan M, Romana M, Connes P. Hydroxyurea treatment does not increase blood viscosity and improves red blood cell rheology in sickle cell anemia. Haematologica. 2015;100:e383–e386. [PubMed]
  • Li H, Ha V, Lykotrafitis G. Modeling sickle hemoglobin fibers as one chain of coarse-grained particles. J Biomech. 2012a;45:1947–1951. [PubMed]
  • Li H, Lykotrafitis G. A coarse-grain molecular dynamics model for sickle hemoglobin fibers. J Mech Behav Biomed Mater. 2011;4:162–173. [PubMed]
  • Li J, Dao M, Lim CT, Suresh S. Spectrin-level modeling of the cytoskeleton and optical tweezers stretching of the erythrocyte. Biophys J. 2005;88:3707–3719. [PubMed]
  • Li XJ, Caswell B, Karniadakis GE. Effect of chain chirality on the self-assembly of sickle hemoglobin. Biophys J. 2012b;103:1130–1140. [PubMed]
  • Li XJ, Peng ZL, Lei H, Dao M, Karniadakis GE. Probing red blood cell mechanics, rheology and dynamics with a two-component multiscale model. Philos Trans R Soc A. 2014;372:20130389. [PMC free article] [PubMed]
  • Li XJ, Du E, Lei H, Tang YH, Dao M, Suresh S, Karniadakis GE. Patient-specific blood rheology in sickle-cell anemia. Interface Focus. 2016;6:20150065. [PMC free article] [PubMed]
  • Lu LJ, Li X, Vekilov PG, Karniadakis GE. Probing the twisted structure of sickle hemoglobin fibers via particle simulations. Biophys J. 2016a;110:2085–2093. [PubMed]
  • Lu X, Wood DK, Higgins JM. Deoxygenation reduces sickle cell blood flow at arterial oxygen tension. Biophys J. 2016b;110:2751–2758. [PubMed]
  • Lutsko JF, Nicolis G. Theoretical evidence for a dense fluid precursor to crystallization. Phys Rev Lett. 2006;96:46102. [PubMed]
  • Maciaszek JL, Andemariam B, Abiraman K, Lykotrafitis G. AKAP-dependent modulation of BCAM/Lu adhesion on normal and sickle cell disease rbcs revealed by force nanoscopy. Biophys J. 2014;106:1258–1267. [PubMed]
  • Maciaszek JL, Lykotrafitis G. Sickle cell trait human erythrocytes are significantly stiffer than normal. J Biomech. 2011;44:657–661. [PubMed]
  • Manwani D, Frenette PS. Vaso-occlusion in sickle cell disease: pathophysiology and novel targeted therapies. Blood. 2013;122:3892–3898. [PubMed]
  • Messmann R, Gannon S, Sarnaik S, Johnson R. Mechanical properties of sickle cell membranes. Blood. 1990;75:1711–1717. [PubMed]
  • Mozzarelli A, Hofrichter J, Eaton WA. Delay time of hemoglobin S polymerization prevents most cells from sickling in vivo. Science. 1987;237:500–506. [PubMed]
  • Odièvre MH, Bony V, Benkerrou M, Lapouméroulie C, Alberti C, Ducrocq R, Jacqz-Aigrain E, Elion J, Cartron JP. Modulation of erythroid adhesion receptor expression by hydroxyurea in children with sickle cell disease. Haematologica. 2008;93:502–510. [PubMed]
  • Parise LV, Berliner N. Sickle cell disease: challenges and progress. Blood. 2016;127:789. [PubMed]
  • Pauling L, Itano HA, Singer SJ, Wells IC. Sickle cell anemia, a molecular disease. Science. 1949;110:543–548. [PubMed]
  • Pivkin IV, Karniadakis GE. Accurate coarse-grained modeling of red blood cells. Phys Rev Lett. 2008;101:118105. [PubMed]
  • Samuel RE, Briehl RW. Nucleation and growth of fibres and gel formation in sickle cell hemoglobin. Nature. 1990;345:833–835. [PubMed]
  • Serjeant GR, Serjeant BE, Milner PF. The irreversibly sickled cell; a determinant of haemolysis in sickle cell anemia. Br J Haematol. 1969;17:527–533. [PubMed]
  • Shiryayev A, Gunton JD. Crystal nucleation for a model of globular proteins. J Chem Phys. 2004;120:8318–8326. [PubMed]
  • Strasser BJ. Molecular medicine – “Sickle cell anemia, a molecular disease” Science. 1999;286:1488–1490. [PubMed]
  • Rees DC, Williams TN, Gladwin MT. Sickle-cell disease. Lancet. 2010;376:2018–2031. [PubMed]
  • Telen MJ. Erythrocyte adhesion receptors: blood group antigens and related molecules. Transfus Med Rev. 2005;19:32–44. [PubMed]
  • ten Wolde PR, Frenkel D. Enhancement of protein crystal nucleation by critical density fluctuations. Science. 1997;277:1975–1978. [PubMed]
  • Turner MS, Briehl RW, Ferrone FA, Josephs R. Twisted protein aggregations and disease: the stability of sickle hemoglobin fibers. Phys Rev Lett. 2003;90:128103. [PubMed]
  • Usami S, Chein S, Scholtz PM, Bertles FJ. Effect of deoxygenation on blood rheology in sickle cell disease. Microvasc Res. 1975;9:324–334. [PubMed]
  • Vekilov PG. Sickle-cell hemoglobin polymerization: is it the primary pathogenic event of sickle-cell anaemia? Br J Haematol. 2007;139:173–184. [PubMed]
  • Ware RE. How I use hydroxyurea to treat young patients with sickle cell anemia. Blood. 2010;115:5300–5311. [PubMed]
  • Ware RE, Zimmerman SA, Schultz WH. Hydroxyurea as an alternative to blood transfusions for the prevention of recurrent stroke in children with sickle cell disease. Blood. 1999;94:3022–3026. [PubMed]
  • Wautier JL, Galacteros F, Wautier MP, Pintigny D, Beuzard Y, Rosa J, Caen JP. Clinical manifestations and erythrocyte adhesion to endothelium in sickle cell syndrome. Am J Hematol. 1985;19:121–130. [PubMed]
  • Wood DK, Soriano A, Mahadevan L, Higgins JM, Bhatia SN. A biophysical indicator of vaso-occlusive risk in sickle cell disease. Sci Transl Med. 2012;4:123ra26. [PMC free article] [PubMed]
  • Yang YS, Meyer RB, Hagan MF. Self-limited self-assembly of chiral filaments. Phys Rev Lett. 2010;104:258102. [PubMed]
  • Yazdani A, Li XJ, Karniadakis GE. Dynamic and rheological properties of soft biological cell suspensions. Rheol Acta. 2016;55:433–449. [PMC free article] [PubMed]
  • Zennadi R, Hines PC, De Castro LM, Cartron JP, Parise LV, Telen MJ. Epinephrine acts through erythroid signaling pathways to activate sickle cell adhesion to endothelium via LW-αvβ3 interactions. Blood. 2004;104:3774–3781. [PubMed]
  • Zennadi R, Moeller BJ, Whalen EJ, Batchvarova M, Xu K, Shan S, Delahunty M, Dewhirst MW, Telen MJ. Epinephrine-induced activation of LW-mediated sickle cell adhesion and vaso-occlusion in vivo. Blood. 2007;110:2708–2717. [PubMed]