Skip to main content
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian
  • Log in
  • My Cart

Main menu

  • Home
  • Articles
    • Current
    • Latest Articles
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • Archive
  • Front Matter
  • News
    • For the Press
    • Highlights from Latest Articles
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Purpose and Scope
    • Editorial and Journal Policies
    • Submission Procedures
    • For Reviewers
    • Author FAQ
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home

Advanced Search

  • Home
  • Articles
    • Current
    • Latest Articles
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • Archive
  • Front Matter
  • News
    • For the Press
    • Highlights from Latest Articles
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Purpose and Scope
    • Editorial and Journal Policies
    • Submission Procedures
    • For Reviewers
    • Author FAQ

New Research In

Physical Sciences

Featured Portals

  • Physics
  • Chemistry
  • Sustainability Science

Articles by Topic

  • Applied Mathematics
  • Applied Physical Sciences
  • Astronomy
  • Computer Sciences
  • Earth, Atmospheric, and Planetary Sciences
  • Engineering
  • Environmental Sciences
  • Mathematics
  • Statistics

Social Sciences

Featured Portals

  • Anthropology
  • Sustainability Science

Articles by Topic

  • Economic Sciences
  • Environmental Sciences
  • Political Sciences
  • Psychological and Cognitive Sciences
  • Social Sciences

Biological Sciences

Featured Portals

  • Sustainability Science

Articles by Topic

  • Agricultural Sciences
  • Anthropology
  • Applied Biological Sciences
  • Biochemistry
  • Biophysics and Computational Biology
  • Cell Biology
  • Developmental Biology
  • Ecology
  • Environmental Sciences
  • Evolution
  • Genetics
  • Immunology and Inflammation
  • Medical Sciences
  • Microbiology
  • Neuroscience
  • Pharmacology
  • Physiology
  • Plant Biology
  • Population Biology
  • Psychological and Cognitive Sciences
  • Sustainability Science
  • Systems Biology

Probing vasoocclusion phenomena in sickle cell anemia via mesoscopic simulations

Huan Lei and George E. Karniadakis
PNAS July 9, 2013 110 (28) 11326-11330; https://doi.org/10.1073/pnas.1221297110
Huan Lei
Division of Applied Mathematics, Brown University, Providence, RI 02912
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
George E. Karniadakis
Division of Applied Mathematics, Brown University, Providence, RI 02912
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: George_Karniadakis@brown.edu
  1. Edited by William R. Schowalter, Princeton University, Princeton, NJ, and approved May 23, 2013 (received for review December 6, 2012)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Abstract

Vasoocclusion crisis is a key hallmark of sickle cell anemia. Although early studies suggest that this crisis is caused by blockage of a single elongated cell, recent experiments have revealed that vasoocclusion is a complex process triggered by adhesive interactions among different cell groups in multiple stages. However, the quantification of the biophysical characteristics of sickle cell anemia remains an open issue. Based on dissipative particle dynamics, we develop a multiscale model for the sickle red blood cells (SS-RBCs), accounting for diversity in both shapes and cell rigidities, to investigate the precise mechanism of vasoocclusion. First, we investigate the adhesive dynamics of a single SS-RBC in shear flow and static conditions, and find that the different cell groups (SS2: young-deformable SS-RBCs, ISCs: rigid-irreversible SS-RBCs) exhibit heterogeneous adhesive behavior due to the diverse cell morphologies and membrane rigidities. We quantify the observed adhesion behavior (in static conditions) in terms of a balance of free energies due to cell adhesion and deformation, and propose a power law that relates the free-energy increase as a function of the contact area. We further simulate postcapillary flow of SS-RBC suspensions with different cell fractions. The more adhesive SS2 cells interact with the vascular endothelium and trap ISC cells, resulting in vasoocclusion in vessels less than Graphic depending on the hematocrit. Under inflammation, adherent leukocytes may also trap ISC cells, resulting in vasoocclusion in even larger vessels.

  • adhesion dynamics
  • free energy analysis
  • disease cell
  • blood flow

Sickle cell anemia (SCA) is a genetic disease originating from the abnormal sickle hemoglobin molecules (HbS). In hypoxic conditions, the intracellular HbS solution transitions into a polymerized state, resulting in a series of alterations in the cell membrane functions. According to the seminal study by Kaul et al. (1), suspensions of sickle red blood cells (SS-RBCs) contain heterogeneous cell density groups, which can be roughly divided into four fractions according to the intracellular mean corpuscular hemoglobin concentration (MCHC). Fractions I (SS1) and II (SS2) with moderate MCHC are mainly composed of reticulocytes and discocytes, respectively, with MCHC similar to healthy cells. On the other hand, fractions III (SS3) and IV (SS4) with high MCHC are mainly composed of rigid discocytes and irreversible sickle cells (ISC). Associated with the heterogeneous cell groups is the abnormal rheology of the SS-RBC suspensions. Compared with healthy blood, the SS-RBC suspensions exhibit higher flow resistance in shear flow (2, 3), as confirmed in microfluidic experiments (4) and in isolated vascular systems (5).

Early studies of SCA suggested that the major cause of vasoocclusion is the sickling process of the dense SS4 cells during their circulation in capillaries. However, subsequent investigations (5⇓⇓–8) have shown that vasoocclusion mainly occurs in postcapillaries. Moreover, it has been found that there is no direct correlation between the percentage of the dense SS-RBCs and the disease severity (9). Recent studies have revealed that vasoocclusion is a complex process triggered by the interactions between multiple density groups (10, 11), where each group contributes differently to the occlusion crisis. In addition, in an inflammation-activated state, leukocytes adherent to the vessel wall may further trap the SS-RBCs, causing occlusion in venular flow (12).

Despite these findings, understanding the precise mechanism of vasoocclusion and quantifying cell interactions, including interactions with the vascular endothelium, are still open research questions. Several numerical investigations (13, 14) have focused on the effect of membrane viscoelasticity and morphology of SS-RBCs on blood flow resistance. However, the heterogeneous SS-RBC adhesive properties, the physiological contribution to vasoocclusion for individual SS-RBC groups, and the hemodynamic effect of the leukocytes have not been fully understood yet. Numerical investigation of these issues may enhance our understanding of vasoocclusion and potentially reveal new paradigms for therapeutic treatments by targeting specific adhesive conditions. In this paper, we present a systematic simulation study of the biophysical characteristics of vasoocclusion in SCA.

Results and Discussion

The entire simulation framework is based on dissipative particle dynamics (DPD), including models for RBCs, leukocytes, plasma, and the wall boundaries (Figs. S1–S3, SI Text). Unlike the healthy SS-RBCs, the SS-RBC membrane expresses multiple types of adhesive receptors (15), which interact with ligands on the endothelial cells through different pathways. We use a stochastic coarse-grained model to represent the interaction between the SS-RBC and the endothelial cells, where the diverse types of adhesive proteins are modeled by effective coarse-grained particles expressed on the cell membrane and vascular wall (Fig. S4 and SI Text).

First, we examine the adhesive dynamics of single SS-RBCs of different density groups; our simulation results validate the aforementioned density-class dependence of the adhesive properties. Subsequently, we explore the adhesion mechanism by computing the free-energy increase as a response of adhesive interaction and quantify the magnitude of the force between different cell groups and the vascular wall. We then use the same sets of parameters to investigate the hemodynamics of SS-RBC suspensions mixed with different cell groups. By analyzing the distinct contribution of individual cell groups, we identify the specific physiological conditions that trigger vasoocclusion.

Shear Flow System.

Experiments in ref. 16 have used an in vitro shear flow system to show that the light density cells exhibit larger adhesion, whereas the densest ISCs correspond to the least adhesive group. However, the following fundamental question remains unresolved: Are the different adhesive responses a consequence of the different adhesive proteins on the cell membranes, or rather of the biomechanical properties of the different cell groups? Kaul et al. (6) explored this issue by examining the alteration of SS-RBC adhesivity after a dehydration/rehydration treatment on individual cell groups and found that the cell adhesivity of the deformable SS2 and the dense SS4 cells can be reversed after controlled treatment. They hypothesized that both light- and high-density cells have similar “adhesion potential,” whereas the heterogeneous cell adhesivity is mainly due to the differences in cell deformability and shape peculiarities among the multiple cell groups. To investigate the validity of this hypothesis, we simulated the cell adhesive dynamics of SS-RBCs with different cell rigidity and morphologic characteristics.

We consider three distinct SS-RBCs under shear flow, as shown in Fig. 1. Following ref. 6, we assume that the three cells share similar “adhesive potential,” and hence we set identical adhesive parameters in all three simulations. Cell I represents a SS2 deformable discocyte of cell rigidity similar to the healthy RBC. We set the shear modulus Graphic and bending rigidity Graphic, i.e., similar to healthy cell rigidity (17⇓–19). Cell III represents an ISC of the SS4 type. Experimental measurements (20, 21) of cell rigidity report scattered values, which further depend on the deoxygenation–oxygenation cycle. For example, an ektacytometer study (21) showed that the cell rigidity of oxygenated ISC is ∼12 times larger than a healthy cell, whereas a micropipette study (20) showed that the cell rigidity of the SS4 group is about 2–3 times larger than the healthy cell; however, this value increases to 20 times the healthy value after one deoxygenation–reoxygenation cycle. In the present work, we set the shear modulus Graphic, because an ISC undergoes at least one deoxygenation–reoxygenation cycle. The modification to cell bending rigidity by the deoxygenation–reoxygenation process is unknown; here we set Graphic and also conduct sensitivity studies. Cell II represents a rigid discocyte cell (22) of SS3 type with medium MCHC value; we set Graphic and Graphic for the comparative study.

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Sickle cells in shear flow: (A) Successive snapshots of SS-RBCs in shear flow, Movie S1. Labels I, II, and III correspond to a deformable SS2 cell, rigid SS3 cell, and ISC, respectively. Arrow indicates the flow direction. (B and C) Instantaneous contact area and velocity for SS-RBC in shear flow conditions.

Fig. 1A shows successive snapshots of the three cells. For the same shear rate Graphic and same membrane adhesive characteristics, the three cells exhibit substantially different adhesive dynamics (Movie S1). Cell I shows firm adhesion to the lower plate after the bond formation stage; the corresponding contact area (see SI Text for definition) is about Graphic, as shown in Fig. 1B. Cell II shows transient adhesion to the lower plate initially, but then undergoes a periodic flip movement along the flow direction, characterized by the peak values of the instantaneous cell velocity at 0.30, 0.53 and Graphic s. This indicates a weaker cell adhesivity than cell I. Accordingly, the contact area achieves minimum values at those times. The cell eventually detaches from the plate after two to three flips. Differently from cells I and II, cell III does not show firm/transient adhesion to the plate; instead, it directly detaches from the lower plate and moves freely without adhesive bonds established thereafter. Given the similar adhesive potential, the present results validate the hypothesis of heterogeneous cell adhesive dynamics under shear flow conditions. To further quantify the effect of the cell rigidity and morphology on adhesive properties, next we investigate cell adhesion in static conditions by using free-energy analysis.

Static Conditions.

Fig. 2 shows the instantaneous contact area between the cell and the plate for the three SS-RBCs. As in the shear flow conditions, the contact area shows an inverse relationship with the cell rigidity. Whereas the contact area for all of the three cells increases to Graphic within the initial stage, the contact areas between the plate and cells I and II further increase to Graphic and Graphic, respectively, at the final stage. These results can be understood by analyzing the change of free energy Graphic due to the adhesive interaction, i.e.,Embedded Imagewhere Graphic represents the free-energy change due to the cell deformation, and Graphic represents the free-energy change due to the adhesive bond formation between the cell and adhesive ligands. In the present system Graphic can be approximated by Graphic, where Graphic represents the change of contact area and β is the energy coefficient determined by the adhesive affinity. Fig. 2B shows Graphic for different adhesive affinities (therefore yielding different contact areas). The numerical results are fitted by Graphic, where α depends on the cell rigidity of individual cells. The equilibrium state is determined by the balance between the two energy terms, e.g., Graphic, leading to different contact areas. For SS2 cells (corresponding to small α), the adhesive interaction plays a dominant role in the initial stage, resulting in cell deformation and further increase of the contact area. In contrast, a rigid SS4 cell exhibits “solid”-like properties with a larger energy barrier for cell deformation. After the initial stage, the adhesive interaction driven by the bond formation cannot overcome the free-energy increase induced by the cell deformation and therefore prohibits further increase of contact area.

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Adhesion in static conditions: (A) Instantaneous contact area between the SS-RBC and the plate coated with adhesive ligands. Black curve represents the contact area of a discocyte with cell rigidity similar to the ISC (cell III). (B) Increase of cell free energy as a function of the contact area. (C) Adhesive force between the cell and the plate as a function of the membrane rigidities for cell morphologies. Error bar represents the adhesive force computed from four independent simulations. (Inset) Sketch of the simulation setup, where a uniform lift force is applied on the upper part of an ISC.

Compared with the shear flow system, we note that the equilibrium contact area obtained in the static condition is smaller. This discrepancy is mainly due to the greater cell deformation under the shear flow condition, i.e., part of the deformation free energy Graphic is balanced by the hydrodynamic force exerted on the cell membrane under the shear flow condition. The extended cell membrane facilitates the bond formation and results in a larger contact area. Moreover, compared with the ISC, a discocyte with similar cell rigidity yields a larger contact area in equilibrium state, as shown in Fig. 2A. This result demonstrates that the peculiar elongated and curved characteristics of ISC may further prevent the membrane receptors from interacting with the ligands.

To quantify the effects of shape and rigidity on cell adhesion, we compute the adhesive force between the cell and plate for different cases. Starting from the steady state obtained from the static condition, we apply uniformly a lift force to the upper part of the cell membrane (Fig. 2C, Inset). The adhesive force is determined as the lift force that drives the cell detaching from the plate in a quasistatic way. The deformable SS2 cells exhibit the largest adhesive force of Graphic, whereas this value decreases to 34 pN in the case of rigid cells with Graphic and Graphic. Compared with the discocyte, the ISC with similar cell rigidity yields an adhesive force of 17 pN. For smaller cell rigidity (corresponding to the rehydration treatment of ISC in ref. 6) with Graphic and Graphic, the adhesive force is approximately 24 pN, compared with 43 pN for a discocyte with similar cell rigidity. In a recent experimental study (23), force spectroscopy was used to quantify the cytoadherence of RBC invaded by a Plasmodium falciparum parasite. It was reported that the rupture force between the invaded RBC and the CHO cells is within the range Graphic. This result is similar to the adhesive force magnitude for the RBC–endothelium interaction reported in the present work.

SS-RBC Suspension in Postcapillary Flow.

Single ISCs may occasionally result in occlusion at precapillary junctions (1, 5, 24); however, most of the occlusion sites are in postcapillaries and exhibit a specific cell pattern, where the adherent deformable cells form a sieve-like configuration and selectively trap the ISCs. We model blood circulation by SS-RBC suspensions first in a tube of diameter Graphic, as shown in Fig. 3; small green particles represent the adhesive ligands coated on the tube wall. To quantify the distinct role of different cell groups, we infuse different cell groups into the tube and apply a pressure gradient Graphic.

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

Vasoocclusion in postcapillaries: instantaneous mean velocity of blood flow in a cylindrical tube of Graphic infused with different SS-RBC suspensions. Red curve represents the resultant velocity infused with SS2 and ISC cell groups. (Inset) Instantaneous snapshots where SS2 cells adhere to vessel wall, consequently trapping the ISCs and resulting in cell blockage. Green curve represents the blood velocity infused with SS2 and ISC cell groups, where adhesive interaction is only applied to the ISC group. The time axis is scaled by 0.5 for better visualization. The Inset represents a snapshot where transient adhesion is established between ISC and the tube wall. Steady flow is recovered as the cell detaches from the tube wall. Blue curve represents the instantaneous velocity of blood flow infused with SS2 and healthy cell groups. Blood flow exhibits a slowdown but not a full occlusion; Movies S2, S3, and S4.

First, we consider suspensions composed of SS2 cells (labeled by blue) and ISC cells (labeled by red) with hematocrit Graphic, similar to experiments (5, 6). Steady flow is achieved by turning off the adhesive interaction during the initial stage, yielding a mean flow velocity of Graphic. Starting from the steady state Graphic, we turn on the adhesive interaction between the SS-RBCs and the ligand particles and compute the instantaneous flow velocity across the tube, as represented by the red curve in Fig. 3. Steady flow is maintained until one of the SS2 cells adheres to the vessel wall, triggering a sharp decrease of blood flow around Graphic. As a positive feedback, the decreased blood flow induces more SS2 cells adherent to the vessel wall, leading to a further decrease of flow rate at Graphic and Graphic. Moreover, these adherent cells decrease the effective vessel lumen near the adherent region, resulting in a secondary trapping of the ISC groups. The final occlusion state is achieved around Graphic with cell patterns similar to the experimental observations (5, 6). We have also performed another three simulations starting from different initial conditions, and in all cases we obtained the same final pattern of full occlusion.

To explore the unique contribution of the deformable cell group, we performed a similar simulation of blood flow where the SS2 cell–ligand interaction is turned off. Starting from the steady flow at Graphic, blood flow is simulated for 6 s, as represented by the green line in Fig. 3. Due to the ISC–ligand interaction, transient adhesive bonds can be formed, resulting in the decrease of blood flow at Graphic. However, firm adhesion cannot be established: blood flow can recover the initial flow rate when the adherent ISCs detach from the tube wall. Blood occlusion is not observed within the simulation time, in agreement with the experimental observation (6) that the ISCs do not result in consistent microvascular blockage when infused alone. This finding reveals the distinct role of the SS2 cell in the vasoocclusion crisis and is consistent with the positive correlation between the SS-RBC deformability index and disease severity reported in clinical observations (9). The ISC cell group, in contrast, contributes differently to the occlusion crisis. We quantify its unique contribution by simulating the blood suspension mixed with the deformable SS2 and healthy cells. Starting from the steady flow at Graphic, we compute the mean velocity with the adhesive interaction between the SS2 and ligands, as represented by the blue line in Fig. 3. The flow velocity decreases due to the cell adhesion at Graphic and Graphic. Although blood flow exhibits a substantial slowdown, full occlusion is not realized due to the high deformability of the healthy RBCs, which enables the cells to squeeze through the sieve-like region formed by adherent SS2 cells. The present result, in turn, reveals the distinct role of the ISC group in the vasoocclusion process. Although the least adhesive, the ISC group, due to its high cell rigidity and peculiar shape, serves as the particular cell group trapped by adherent cells in the postcapillaries. This finding explains the experimental observation that a large number of dense cells accumulate in the occlusion region but are not observed in the peripheral blood circulation (5).

To quantify the effect of the vascular size on occlusion, we have performed similar simulations but in tubes of different diameters. For Graphic, occlusion does not occur with tube diameter larger than Graphic (Fig. S5); the sieve-like pattern formed by SS2 cells cannot fully trap the ISCs due to the larger spatial accommodation near the adhesion region. However, for Graphic, occlusion is observed with Graphic but occlusion is prevented with tube diameter larger than Graphic (see Fig. S6 and SI Text for more details).

Effect of Inflammation-Stimulated Leukocytes.

Recent studies (25) have shown that SCA is often accompanied by an inflammatory endothelial phenotype, resulting in elevated leukocyte recruitment in blood circulation. Moreover, studies by Turhan et al. (12) in transgenic-knockout mice have shown that the inflammation-stimulated (by cytokine TNF-α) adherent leukocytes, upon interactions with SS-RBCs, lead to occlusion in venular flow. We use our computational framework to investigate the vasoocclusion induced by inflammation-stimulated leukocytes (Methods). We first consider blood flow of SS-RBC suspension Graphic with a single leukocyte in a tube of diameter Graphic, as shown in Fig. 4. Due to cell migration, the leukocyte touches the wall at Graphic and firm adhesion is established thereafter. The mean flow velocity drops from Graphic to Graphic during this stage. During the inflammation-stimulated stage, we turn on the cell aggregation interaction between the SS-RBCs and leukocytes (see SI Text for details). Multiple SS-RBCs get trapped by the leukocyte, resulting in full occlusion at Graphic. For venular flow with larger diameter, multiple leukocytes may accumulate in the inflammation region. Similar to the experiment (12), we simulate a SS-RBC suspension Graphic with three leukocytes in a tube of diameter Graphic and compute the mean flow velocity during different time intervals, as shown in Fig. 4. Stage I represents the blood flow during the initial stage of the inflammatory response. Starting from the steady flow at Graphic, the mean flow velocity drops to Graphic due to the attachment of leukoctyes on the vessel wall at Graphic s. Stage II represents the blood flow with moderate RBC–leukocyte interaction; the blood flow exhibits a slowdown due to the adherent leukocytes. Stage III represents the late stage of the inflammatory response, where the RBC–leukoctye interaction is further intensified. Multiple SS-RBCs are trapped on the adherent leukocytes, leading to full occlusion at Graphic (Movies S5, S6, S7, S8, S9, and S10).

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

Effect of leukocytes: instantaneous mean velocity of the blood flow in a tube of Graphic, one leukocyte (A) and Graphic, Graphic, three leukocytes (B). (Inset) Snapshots represent blood cells in free motion, leukocyte adhesion, and blood occlusion states. For Graphic, the Inset plot represents the blood flow velocity of the present study (blue) and the experimental results (red) (12), where measurements are taken on Graphic venules with average diameter Graphic and Graphic before and after inflammation–stimulation; Movies S5, S6, S7, S8, S9, and S10.

Discussion.

Our simulations of individual SS-RBCs in shear flow validated the hypothesis of the importance of cell rigidity and morphology (5, 6, 16), and elucidated the distinct behavior of each individual cell group under adhesive conditions. This behavior was further quantified by computing the adhesive force and free -response for all cell groups. In simulations of postcapillary SS-RBC flow, we identified which combination of cell groups and at what adhesion strengths lead to full or transient occlusion states. Under typical physiological conditions, our sensitivity studies with respect to the tube size and hematocrit show that blood occlusion mainly occurs in postcapillaries with diameter smaller than Graphic, a finding consistent with the experimental results reported in ref. 5. Further free energy analysis on cell adhesion (Fig. S7) and discussion on the important dimensionless groups is included in SI Text. Moreover, under inflammation-stimulated conditions, the present work quantified the influence of the adherent leukocytes, which may further trap the SS-RBC and lead to vasoocclusion in venular flow, especially in larger vessels. In particular, our simulations revealed and quantified the multiinteractional and multistage nature of vasoocclusion, which may be useful in developing effective biomarkers (26) and therapeutic treatments for SCA. In general, our computational framework provides a capability, which, in conjunction with microfluidic experiments, can be used in the preliminary screening of proposed drugs and the interpretation of clinical outcomes.

Methods

Simulation Method.

The DPD method is a Lagrangian-based coarse-grained particle method (27, 28) used to simulate complex fluids (29⇓–31). Each DPD particle can be viewed as a coarse-grained virtual cluster of multiple atomistic particles (32), where the particle motion is governed by pairwise force imposed between the DPD particles (Table S1 and SI Text).

Cell Models.

A healthy RBC (33, 34) is modeled by a triangulated surface (e.g., 500 vertices) of the biconcave RBC shape (35) in equilibrium state. Each vertex on the membrane surface is represented by a DPD particle. Membrane elasticity is imposed by viscoelastic bonds between the DPD particles. Bending rigidity is modeled by the bending resistance between the adjacent triangulated surface. In addition, area and volume constraints are imposed to represent the incompressibility of the lipid bilayer and RBC cytosol. The model has been validated in a number of studies, e.g., for the mechanical properties (33, 34), as well as rheological and hemodynamic properties of blood flow in healthy (36), malaria-infected (37), and SCA conditions (38). The ISC model is constructed by applying surface tension on a healthy RBC membrane (38), representing the distortion of the cell membrane due to the growth of the intracellular sickle hemoglobin polymers. The distorted shape is redefined as the equilibrium state of the ISC with minimum free energy. The leukocyte is modeled by a triangulated surface (e.g., 1,000 vertices) of a sphere of radius Graphic. Similar to the RBC model, proper cell rigidity is imposed by the viscoelastic bonds and bending resistance; total area and volume constraints are imposed for individual cells. More details on the cell models are presented in SI Text.

Cell Interactions.

The adhesive interaction between the SS-RBC/leukocyte and the adhesive ligands on vascular endothelium is modeled by the stochastic bond formation/dissociation according to refs. 39, 40. Cell vertices can interact with the endothelial ligands within interaction distance Graphic. During each time step Graphic, transient bonds can be formed between cell vertices and endothelial ligands with probability Graphic, whereas existing bonds can be ruptured with probability Graphic, where Graphic are, respectively, bond formation and dissociation rate depending on the adhesive affinity between the blood cell and vascular endothelium (Table S2). The interaction between the leukocyte and SS-RBC is modeled by the Morse potential, Graphic, where r is the distance between cell-membrane vertices of adjacent cells, β determines the range of interaction, and Graphic and Graphic define the equilibrium distance and attractive strength, respectively. See SI Text for more details.

Acknowledgments

H.L. acknowledges Dr. Xin Yi (Brown University), Dr. Petia Vlahovska (Brown University), and Dr. Dmitry A. Fedosov (Forschungszentrum Jülich) for useful discussions. This work was sponsored by National Institutes of Health Grant R01HL094270 and the new Collaboratory on Mathematics for Mesocopic Modeling of Materials (CM4) supported by the Department of Energy (DOE). Computations were possible due to a DOE/Innovative and Novel Computational Impact on Theory and Experiment Award.

Footnotes

  • ↵1To whom correspondence should be addressed. E-mail: George_Karniadakis{at}brown.edu.
  • Author contributions: H.L. and G.E.K. designed research, performed research, analyzed data, and wrote the paper.

  • The authors declare no conflict of interest.

  • This article is a PNAS Direct Submission.

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1221297110/-/DCSupplemental.

References

  1. ↵
    1. Kaul DK,
    2. Fabry ME,
    3. Windisch P,
    4. Baez S,
    5. Nagel RL
    (1983) Erythrocytes in sickle cell anemia are heterogeneous in their rheological and hemodynamic characteristics. J Clin Invest 72(1):22–31.
    OpenUrlCrossRefPubMed
  2. ↵
    1. Kaul DK,
    2. Xue H
    (1991) Rate of deoxygenation and rheologic behavior of blood in sickle cell anemia. Blood 77(6):1353–1361.
    OpenUrlAbstract/FREE Full Text
  3. ↵
    1. Kaul DK,
    2. Liu XD
    (1999) Rate of deoxygenation modulates rheologic behavior of sickle red blood cells at a given mean corpuscular hemoglobin concentration. Clin Hemorheol Microcirc 21(2):125–135.
    OpenUrlPubMed
  4. ↵
    1. Higgins JM,
    2. Eddington DT,
    3. Bhatia SN,
    4. Mahadevan L
    (2007) Sickle cell vasoocclusion and rescue in a microfluidic device. Proc Natl Acad Sci USA 104(51):20496–20500.
    OpenUrlAbstract/FREE Full Text
  5. ↵
    1. Kaul DK,
    2. Fabry ME,
    3. Nagel RL
    (1989) Microvascular sites and characteristics of sickle cell adhesion to vascular endothelium in shear flow conditions: pathophysiological implications. Proc Natl Acad Sci USA 86(9):3356–3360.
    OpenUrlAbstract/FREE Full Text
  6. ↵
    1. Kaul DK,
    2. Chen D,
    3. Zhan J
    (1994) Adhesion of sickle cells to vascular endothelium is critically dependent on changes in density and shape of the cells. Blood 83(10):3006–3017.
    OpenUrlAbstract/FREE Full Text
  7. ↵
    1. Kaul DK,
    2. Fabry ME
    (2004) In vivo studies of sickle red blood cells. Microcirculation 11(2):153–165.
    OpenUrlCrossRefPubMed
  8. ↵
    1. Barabino GA,
    2. Platt MO,
    3. Kaul DK
    (2010) Sickle cell biomechanics. Annu Rev Biomed Eng 12:345–367.
    OpenUrlCrossRefPubMed
  9. ↵
    1. Ballas SK,
    2. et al.
    (1988) Rheologic predictors of the severity of the painful sickle cell crisis. Blood 72(4):1216–1223.
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. Ballas SK,
    2. Mohandas N
    (2004) Sickle red cell microrheology and sickle blood rheology. Microcirculation 11(2):209–225.
    OpenUrlCrossRefPubMed
  11. ↵
    1. Chiang EY,
    2. Frenette PS
    (2005) Sickle cell vaso-occlusion. Hematol Oncol Clin North Am 19(5):771–784, v.
    OpenUrlCrossRefPubMed
  12. ↵
    1. Turhan A,
    2. Weiss LA,
    3. Mohandas N,
    4. Coller BS,
    5. Frenette PS
    (2002) Primary role for adherent leukocytes in sickle cell vascular occlusion: A new paradigm. Proc Natl Acad Sci USA 99(5):3047–3051.
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Dong C,
    2. Chadwick RS,
    3. Schechter AN
    (1992) Influence of sickle hemoglobin polymerization and membrane properties on deformability of sickle erythrocytes in the microcirculation. Biophys J 63(3):774–783.
    OpenUrlCrossRefPubMed
  14. ↵
    1. Dupin MM,
    2. Halliday I,
    3. Care CM,
    4. Munn LL
    (2008) Lattice Boltzmann modelling of blood cell dynamics. Int J Comput Fluid Dyn 22(7):481–492.
    OpenUrlCrossRef
  15. ↵
    1. Hebbel RP
    (2008) Adhesion of sickle red cells to endothelium: Myths and future directions. Transfus Clin Biol 15(1-2):14–18.
    OpenUrlCrossRefPubMed
  16. ↵
    1. Barabino GA,
    2. McIntire LV,
    3. Eskin SG,
    4. Sears DA,
    5. Udden M
    (1987) Endothelial cell interactions with sickle cell, sickle trait, mechanically injured, and normal erythrocytes under controlled flow. Blood 70(1):152–157.
    OpenUrlAbstract/FREE Full Text
  17. ↵
    1. Humpert C,
    2. Baumann M
    (2003) Local membrane curvature affects spontaneous membrane fluctuation characteristics. Mol Membr Biol 20(2):155–162.
    OpenUrlCrossRefPubMed
  18. ↵
    1. Mills JP,
    2. Qie L,
    3. Dao M,
    4. Lim CT,
    5. Suresh S
    (2004) Nonlinear elastic and viscoelastic deformation of the human red blood cell with optical tweezers. Mech Chem Biosyst 1(3):169–180.
    OpenUrlPubMed
  19. ↵
    1. Suresh S,
    2. et al.
    (2005) Connections between single-cell biomechanics and human disease states: Gastrointestinal cancer and malaria. Acta Biomater 1(1):15–30.
    OpenUrlCrossRefPubMed
  20. ↵
    1. Itoh T,
    2. Chien S,
    3. Usami S
    (1995) Effects of hemoglobin concentration on deformability of individual sickle cells after deoxygenation. Blood 85(8):2245–2253.
    OpenUrlAbstract/FREE Full Text
  21. ↵
    1. Clark MR,
    2. Mohandas N,
    3. Shohet SB
    (1980) Deformability of oxygenated irreversibly sickled cells. J Clin Invest 65(1):189–196.
    OpenUrlCrossRefPubMed
  22. ↵
    1. Evans EA,
    2. Mohandas N
    (1987) Membrane-associated sickle hemoglobin: A major determinant of sickle erythrocyte rigidity. Blood 70(5):1443–1449.
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Carvalho PA,
    2. Diez-Silva M,
    3. Chen H,
    4. Dao M,
    5. Suresh S
    (2013) Cytoadherence of erythrocytes invaded by Plasmodium falciparum: Quantitative contact-probing of a human malaria receptor. Acta Biomater 9(5):6349–6359, doi:10.1016/j.actbio.2013.01.019.
    OpenUrlCrossRefPubMed
  24. ↵
    1. Lipowsky HH,
    2. Usami S,
    3. Chien S
    (1982) Human SS red cell rheological behavior in the microcirculation of cremaster muscle. Blood Cells 8(1):113–126.
    OpenUrlPubMed
  25. ↵
    1. Kaul DK,
    2. Hebbel RP
    (2000) Hypoxia/reoxygenation causes inflammatory response in transgenic sickle mice but not in normal mice. J Clin Invest 106(3):411–420.
    OpenUrlCrossRefPubMed
  26. ↵
    1. Mohandas N,
    2. Evans E
    (1989) Rheological and adherence properties of sickle cells. Potential contribution to hematologic manifestations of the disease. Ann N Y Acad Sci 565:327–337.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Hoogerbrugge PJ,
    2. Koelman JMVA
    (1992) Simulating microscopic hydrodynamic phenomena with dissipative particle dynamics. Europhys Lett 19(3):155–160.
    OpenUrlCrossRef
  28. ↵
    1. Espanol P,
    2. Warren P
    (1995) Statistical mechanics of dissipative particle dynamics. Europhys Lett 30(4):191–196.
    OpenUrlCrossRef
  29. ↵
    1. Fan X,
    2. Phan-Thien N,
    3. Chen S,
    4. Wu X,
    5. Ng TY
    (2006) Simulating flow of DNA suspension using dissipative particle dynamics. Phys Fluids 18(6):063102.
    OpenUrlCrossRef
  30. ↵
    1. Boek ES,
    2. Coveney PV,
    3. Lekkerkerker HNW,
    4. van der Schoot P
    (1997) Simulating the rheology of dense colloidal suspensions using dissipative particle dynamics. Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics 55(3 Pt B):3124–3133.
    OpenUrlCrossRef
  31. ↵
    1. Groot RD,
    2. Warren PB
    (1997) Dissipative particle dynamics: Bridging the gap between atomistic and mesoscopic simulation. J Chem Phys 107(11):4423–4435.
    OpenUrlCrossRef
  32. ↵
    1. Lei H,
    2. Caswell B,
    3. Karniadakis GE
    (2010) Direct construction of mesoscopic models from microscopic simulations. Phys Rev E Stat Nonlin Soft Matter Phys 81(2 Pt 2):026704.
    OpenUrlCrossRefPubMed
  33. ↵
    1. Pivkin IV,
    2. Karniadakis GE
    (2008) Accurate coarse-grained modeling of red blood cells. Phys Rev Lett 101(11):118105.
    OpenUrlCrossRefPubMed
  34. ↵
    1. Fedosov DA,
    2. Caswell B,
    3. Karniadakis GE
    (2010) A multiscale red blood cell model with accurate mechanics, rheology, and dynamics. Biophys J 98(10):2215–2225.
    OpenUrlCrossRefPubMed
  35. ↵
    1. Evans EA,
    2. Skalak R
    (1980) Mechanics and Thermodynamics of Biomembranes (CRC, Boca Raton, FL).
  36. ↵
    1. Fedosov DA,
    2. Pan W,
    3. Caswell B,
    4. Gompper G,
    5. Karniadakis GE
    (2011) Predicting human blood viscosity in silico. Proc Natl Acad Sci USA 108(29):11772–11777.
    OpenUrlAbstract/FREE Full Text
  37. ↵
    1. Fedosov DA,
    2. Caswell B,
    3. Suresh S,
    4. Karniadakis GE
    (2011) Quantifying the biophysical characteristics of Plasmodium-falciparum-parasitized red blood cells in microcirculation. Proc Natl Acad Sci USA 108(1):35–39.
    OpenUrlAbstract/FREE Full Text
  38. ↵
    1. Lei H,
    2. Karniadakis GE
    (2012) Quantifying the rheological and hemodynamic characteristics of sickle cell anemia. Biophys J 102(2):185–194.
    OpenUrlCrossRefPubMed
  39. ↵
    1. Hammer DA,
    2. Apte SM
    (1992) Simulation of cell rolling and adhesion on surfaces in shear flow: General results and analysis of selectin-mediated neutrophil adhesion. Biophys J 63(1):35–57.
    OpenUrlCrossRefPubMed
  40. ↵
    1. King MR,
    2. Hammer DA
    (2001) Multiparticle adhesive dynamics: Hydrodynamic recruitment of rolling leukocytes. Proc Natl Acad Sci USA 98(26):14919–14924.
    OpenUrlAbstract/FREE Full Text
View Abstract
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Probing vasoocclusion phenomena in sickle cell anemia via mesoscopic simulations
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
Citation Tools
Quantifying vasoocclusion in sickle cell anemia
Huan Lei, George E. Karniadakis
Proceedings of the National Academy of Sciences Jul 2013, 110 (28) 11326-11330; DOI: 10.1073/pnas.1221297110

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Quantifying vasoocclusion in sickle cell anemia
Huan Lei, George E. Karniadakis
Proceedings of the National Academy of Sciences Jul 2013, 110 (28) 11326-11330; DOI: 10.1073/pnas.1221297110
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley
Proceedings of the National Academy of Sciences: 116 (8)
Current Issue

Submit

Sign up for Article Alerts

Jump to section

  • Article
    • Abstract
    • Results and Discussion
    • Methods
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

News Feature: Cities serve as testbeds for evolutionary change
Urban living can pressure flora and fauna to adapt in intriguing ways. Biologists are starting to take advantage of this convenient laboratory of evolution.
Image credit: Kristin Winchell (Washington University in St. Louis, St. Louis).
Several aspects of the proposal, which aims to expand open access, require serious discussion and, in some cases, a rethink.
Opinion: “Plan S” falls short for society publishers—and for the researchers they serve
Several aspects of the proposal, which aims to expand open access, require serious discussion and, in some cases, a rethink.
Image credit: Dave Cutler (artist).
Featured Profile
PNAS Profile of NAS member and biochemist Hao Wu
 Nonmonogamous strawberry poison frog (Oophaga pumilio).  Image courtesy of Yusan Yang (University of Pittsburgh, Pittsburgh).
Putative signature of monogamy
A study suggests a putative gene-expression hallmark common to monogamous male vertebrates of some species, namely cichlid fishes, dendrobatid frogs, passeroid songbirds, common voles, and deer mice, and identifies 24 candidate genes potentially associated with monogamy.
Image courtesy of Yusan Yang (University of Pittsburgh, Pittsburgh).
Active lifestyles. Image courtesy of Pixabay/MabelAmber.
Meaningful life tied to healthy aging
Physical and social well-being in old age are linked to self-assessments of life worth, and a spectrum of behavioral, economic, health, and social variables may influence whether aging individuals believe they are leading meaningful lives.
Image courtesy of Pixabay/MabelAmber.

More Articles of This Classification

Biological Sciences

  • Major histocompatibility complex class I diversity limits the repertoire of T cell receptors
  • LINC00116 codes for a mitochondrial peptide linking respiration and lipid metabolism
  • Unnatural verticilide enantiomer inhibits type 2 ryanodine receptor-mediated calcium leak and is antiarrhythmic
Show more

Biophysics and Computational Biology

  • Froghoppers jump from smooth plant surfaces by piercing them with sharp spines
  • Structural hierarchy confers error tolerance in biological materials
  • Structural basis for activity of TRIC counter-ion channels in calcium release
Show more

Physical Sciences

  • Photoexcitation-controlled self-recoverable molecular aggregation for flicker phosphorescence
  • Phosphate graphene as an intrinsically osteoinductive scaffold for stem cell-driven bone regeneration
  • Unnatural verticilide enantiomer inhibits type 2 ryanodine receptor-mediated calcium leak and is antiarrhythmic
Show more

Applied Physical Sciences

  • Observation of chiral surface excitons in a topological insulator Bi2Se3
  • Zintl-phase Eu2ZnSb2: A promising thermoelectric material with ultralow thermal conductivity
  • Amyloid fibril-directed synthesis of silica core–shell nanofilaments, gels, and aerogels
Show more

Related Content

  • No related articles found.
  • Scopus
  • PubMed
  • Google Scholar

Cited by...

  • Simultaneous polymerization and adhesion under hypoxia in sickle cell disease
  • Kinetics of sickle cell biorheology and implications for painful vasoocclusive crisis
  • Scopus (39)
  • Google Scholar

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Latest Articles
  • Archive

PNAS Portals

  • Classics
  • Front Matter
  • Teaching Resources
  • Anthropology
  • Chemistry
  • Physics
  • Sustainability Science

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Press
  • Site Map

Feedback    Privacy/Legal

Copyright © 2019 National Academy of Sciences. Online ISSN 1091-6490