The objective of this project is to develop a validated multiscale modeling methodology for quantifying the biophysical characteristics of sickle cell disease (SCD) -- a hematological disorder that affects tens of thousands of people in US with one in every 500 African-American births resulting in a child with SCD. The pathogenesis of SCD results from (1) irregular red blood cell (RBC) shapes due to hemoglobin polymerization inside the RBCs;(2) stiffening of the RBC membrane;and (3) adhesion of sickle RBCs to the endothelium and the other blood cells. The combination of these phenomena results in vaso-occlusive events or "crises" responsible for the majority of morbidity and mortality associated with SCD but little is certain about the proximal causes or the circumstances in which they occur. The spatio-temporal scales involved in accurately modeling SCD blood flow and vaso-occlusion span at least four orders of magnitude, hence new numerical methods are needed to simulate such multiscale phenomena. We present a general methodology based on 3D dissipative particle dynamics (DPD) to model flow and soft matter seamlessly, i.e., RBCs and other blood cells, blood plasma, cytosol, hemoglobin polymerization, and adhesive dynamics. DPD can be interfaced with molecular dynamics (MD) and with continuum-based description (e.g. Navier-Stokes) based on the triple-decker algorithm we have developed in order to capture molecular details or for computational efficiency in simulating large arteries or networks, respectively. We adopt the same approach here that has proven very effective in our previous work on malaria, namely that models for single RBCs (healthy or sickled), informed and validated from comprehensive single-cell measurements, will be used to predict the collective dynamics and rheology of SCD blood flow. We also present a systematic experimental plan, using microfluidics, nanomechanics and advanced optical techniques, to validate the various stages of the development of our models by targeting individual scales as well as interactions between scales. We will extend the first generation of models to study different modalities of existing and experimental therapeutic interventions for SCD, including simple transfusion, fetal hemoglobin (HbF) induction by hydroxyurea, and RBC hydration. Predictability of multiscale models requires quantifying uncertainty, and, to this end, we will incorporate polynomial chaos methods to model and propagate parametric uncertainties through the multiscale system. We plan to disseminate our models, software tools, and experimental data including the general-purpose triple-decker algorithm, via web-based repositories, existing public open-ware sites, tutorials and through the MSM consortium.

Public Health Relevance

We propose to develop and validate a multiscale modeling methodology for sickle cell disease (SCD) affecting 72,000 people in US. We will model multiscale phenomena across more than four orders of magnitude in spatio-temporal scales. We will develop a new generation of models to study different modalities of therapeutic interventions for SCD, including simple transfusion, HbF induction by hydroxyurea, and RBC hydration.

National Institute of Health (NIH)
National Heart, Lung, and Blood Institute (NHLBI)
Research Project--Cooperative Agreements (U01)
Project #
Application #
Study Section
Special Emphasis Panel (ZEB1-OSR-C (J2))
Program Officer
Qasba, Pankaj
Project Start
Project End
Budget Start
Budget End
Support Year
Fiscal Year
Total Cost
Indirect Cost
Brown University
Biostatistics & Other Math Sci
Schools of Arts and Sciences
United States
Zip Code
Li, Xuejin; Du, E; Lei, Huan et al. (2016) Patient-specific blood rheology in sickle-cell anaemia. Interface Focus 6:20150065
Muenster, Stefan; Beloiartsev, Arkadi; Yu, Binglan et al. (2016) Exposure of Stored Packed Erythrocytes to Nitric Oxide Prevents Transfusion-associated Pulmonary Hypertension. Anesthesiology 125:952-963
Lu, Xinran; Wood, David K; Higgins, John M (2016) Deoxygenation Reduces Sickle Cell Blood Flow at Arterial Oxygen Tension. Biophys J 110:2751-8
Lu, Lu; Li, Xuejin; Vekilov, Peter G et al. (2016) Probing the Twisted Structure of Sickle Hemoglobin Fibers via Particle Simulations. Biophys J 110:2085-93
Hosseini, Poorya; Abidi, Sabia Z; Du, E et al. (2016) Cellular normoxic biophysical markers of hydroxyurea treatment in sickle cell disease. Proc Natl Acad Sci U S A 113:9527-32
Pivkin, Igor V; Peng, Zhangli; Karniadakis, George E et al. (2016) Biomechanics of red blood cells in human spleen and consequences for physiology and disease. Proc Natl Acad Sci U S A 113:7804-9
Guo, Feng; Mao, Zhangming; Chen, Yuchao et al. (2016) Three-dimensional manipulation of single cells using surface acoustic waves. Proc Natl Acad Sci U S A 113:1522-7
Yazdani, Alireza; Li, Xuejin; Em Karniadakis, George (2016) Dynamic and rheological properties of soft biological cell suspensions. Rheol Acta 55:433-449
Li, Peng; Mao, Zhangming; Peng, Zhangli et al. (2015) Acoustic separation of circulating tumor cells. Proc Natl Acad Sci U S A 112:4970-5
Li, Zhen; Bian, Xin; Li, Xiantao et al. (2015) Incorporation of memory effects in coarse-grained modeling via the Mori-Zwanzig formalism. J Chem Phys 143:243128

Showing the most recent 10 out of 24 publications