Patient-specific blood rheology in sickle-cell anaemia

Biomechanics of phagocytosis of red blood cells by macrophages in the human spleen

The clearance of senescent and altered red blood cells (RBCs) in the red pulp of the human spleen involves sequential processes of prefiltration, filtration, and postfiltration. While prior work has elucidated the mechanisms underlying the first two processes, biomechanical processes driving the postfiltration phagocytosis of RBCs retained at interendothelial slits (IES) are still poorly understood. We present here a unique computational model of macrophages to study the role of cell biomechanics in modulating the kinetics of phagocytosis of aged and diseased RBCs retained in the spleen. After validating the macrophage model using in vitro phagocytosis experiments, we employ it to probe the mechanisms underlying the kinetics of phagocytosis of mechanically altered RBCs, such as heated RBCs and abnormal RBCs in hereditary spherocytosis (HS) and sickle cell disease (SCD). Our simulations show pronounced deformation of the flexible and healthy RBCs in contrast to minimal shape changes in altered RBCs. Simulations also show that less deformable RBCs are engulfed faster and at lower adhesive strength than flexible RBCs, consistent with our experimental measurements. This efficient sensing and engulfment by macrophages of stiff RBCs retained at IES are expected to temper splenic congestion, a common pathogenic process in malaria, HS, and SCD. Altogether, our combined computational and in vitro experimental studies suggest that mechanical alterations of retained RBCs may suffice to enhance their phagocytosis, thereby adapting the kinetics of their elimination to the kinetics of their mechanical retention, an equilibrium essential for adequately cleaning the splenic filter to preserve its function.

Gene-environmental influence of space and microgravity on red blood cells with sickle cell disease

A fundamental question in human biology and for hematological disease is how do complex gene-environment interactions lead to individual disease outcome? This is no less the case for sickle cell disease (SCD), a monogenic disorder of Mendelian inheritance, both clinical course, severity, and treatment response, is variable amongst affected individuals. New insight and discovery often lie between the intersection of seemingly disparate disciplines. Recently, opportunities for space medicine have flourished and have offered a new paradigm for study. Two recent Nature papers have shown that hemolysis and oxidative stress play key mechanistic roles in erythrocyte pathogenesis during spaceflight. This paper reviews existing genetic and environmental modifiers of the sickle cell disease phenotype. It reviews evidence for erythrocyte pathology in microgravity environments and demonstrates why this may be relevant for the unique gene-environment interaction of the SCD phenotype. It also introduces the hematology and scientific community to methodological tools for evaluation in space and microgravity research. The increasing understanding of space biology may yield insight into gene-environment influences and new treatment paradigms in SCD and other hematological disease phenotypes.

A Compact Differential Dynamic Microscopy-based Device (cDDM): An Approach Tool for Early Detection of Hypercoagulable State in Transfusion-Dependent-β-Thalassemia Patients

β-Thalassemia especially transfusion-dependent thalassemia (TDT) associates with a hypercoagulable state, which is the main cause of thromboembolic events (TEE). Plasma viscosity and rheological parameters could be essential markers for determining hypercoagulable state in β-thalassemia patients. The traditional methods for measuring viscosity are often limited by large sample volumes and are impractical for routine clinical monitoring. The compact differential dynamic microscopy-based device (cDDM), an optical microscopy for quantitative rheological assessment, was developed and applied for prognosis of the hypercoagulable state in β-TDT with and without splenectomy. The device was performed plasma viscosity measurement using low plasma volume (8 μL) and revealed a value as modulus of complex viscosity |η(ω)| in 7 min. We also parallelly demonstrated the correlation of the viscosity and related-coagulable parameters: complete blood count, prothrombin time (PT), activated partial thromboplastin time (APTT), protein C (PC), protein S (PS), CD62P and CD63 expression, and platelet aggregation test. The thalassemia plasma exhibited a higher value of |η(ω)| than healthy plasma, which can represent a different viscoelastic property among the groups. Even all related-coagulable parameters indicated hypercoagulable state in both nonsplenectomies and splenectomies β-TDT patients when compared to control, only high platelet numbers significantly correlated to high plasma viscosity in the splenectomy group. However, the other coagulable parameters have shown a trend of positive relationship with high plasma viscosity in all β-1thalassemia TDT patients. The relative results suggested that our device would be an approach tool for early detection of hypercoagulable state in transfusion-dependent-β-thalassemia patients, which can help to prevent TEE and the critical consequent-complications.

In silico modeling of patient-specific blood rheology in type 2 diabetes mellitus

Increased blood viscosity in type 2 diabetes mellitus (T2DM) is a risk factor for the development of insulin resistance and diabetes-related vascular complications; however, individuals with T2DM exhibit heterogeneous hemorheological properties, including cell deformation and aggregation. Using a multiscale red blood cell (RBC) model with key parameters derived from patient-specific data, we present a computational study of the rheological properties of blood from individual patients with T2DM. Specifically, one key model parameter, which determines the shear stiffness of the RBC membrane (μ) is informed by the high-shear-rate blood viscosity of patients with T2DM. At the same time, the other, which contributes to the strength of the RBC aggregation interaction (D0), is derived from the low-shear-rate blood viscosity of patients with T2DM. The T2DM RBC suspensions are simulated at different shear rates, and the predicted blood viscosity is compared with clinical laboratory-measured data. The results show that the blood viscosity obtained from clinical laboratories and computational simulations are in agreement at both low and high shear rates. These quantitative simulation results demonstrate that the patient-specific model has truly learned the rheological behavior of T2DM blood by unifying the mechanical and aggregation factors of the RBCs, which provides an effective way to extract quantitative predictions of the rheological properties of the blood of individual patients with T2DM.

Microfluidic study of retention and elimination of abnormal red blood cells by human spleen with implications for sickle cell disease

The spleen clears altered red blood cells (RBCs) from circulation, contributing to the balance between RBC formation (erythropoiesis) and removal. The splenic RBC retention and elimination occur predominantly in open circulation where RBCs flow through macrophages and inter-endothelial slits (IESs). The mechanisms underlying and interconnecting these processes significantly impact clinical outcomes. In sickle cell disease (SCD), blockage of intrasplenic sickled RBCs is observed in infants splenectomized due to acute splenic sequestration crisis (ASSC). This life-threatening RBC pooling and organ swelling event is plausibly triggered or enhanced by intra-tissular hypoxia. We present an oxygen-mediated spleen-on-a-chip platform for in vitro investigations of the homeostatic balance in the spleen. To demonstrate and validate the benefits of this general microfluidic platform, we focus on SCD and study the effects of hypoxia on splenic RBC retention and elimination. We observe that RBC retention by IESs and RBC-macrophage adhesion are faster in blood samples from SCD patients than those from healthy subjects. This difference is markedly exacerbated under hypoxia. Moreover, the sickled RBCs under hypoxia show distinctly different phagocytosis processes from those non-sickled RBCs under hypoxia or normoxia. We find that reoxygenation significantly alleviates RBC retention at IESs, and leads to rapid unsickling and fragmentation of the ingested sickled RBCs inside macrophages. These results provide unique mechanistic insights into how the spleen maintains its homeostatic balance between splenic RBC retention and elimination, and shed light on how disruptions in this balance could lead to anemia, splenomegaly, and ASSC in SCD and possible clinical manifestations in other hematologic diseases.

Cardiovascular consequences of sickle cell disease

Sickle cell disease (SCD) is an inherited blood disorder caused by a single point mutation within the beta globin gene. As a result of this mutation, hemoglobin polymerizes under low oxygen conditions causing red blood cells to deform, become more adhesive, and increase in rigidity, which affects blood flow dynamics. This process leads to enhanced red blood cell interactions with the endothelium and contributes to vaso-occlusion formation. Although traditionally defined as a red blood cell disorder, individuals with SCD are affected by numerous clinical consequences including stroke, painful crisis episodes, bone infarctions, and several organ-specific complications. Elevated cardiac output, endothelium activation along with the sickling process, and the vaso-occlusion events pose strains on the cardiovascular system. We will present a review of the cardiovascular consequences of sickle cell disease and show connections with the vasculopathy related to SCD. We will also highlight biophysical properties and engineering tools that have been used to characterize the disease. Finally, we will discuss therapies for SCD and potential implications on SCD cardiomyopathy.

Microfluidic Microcirculation Mimetic as a Tool for the Study of Rheological Characteristics of Red Blood Cells in Patients with Sickle Cell Anemia

Sickle cell disorder (SCD) is a multisystem disease with heterogeneous phenotypes. Although all patients have the mutated hemoglobin (Hb) in the SS phenotype, the severity and frequency of complications are variable. When exposed to low oxygen tension, the Hb molecule becomes dense and forms tactoids, which lead to the peculiar sickled shapes of the affected red blood cells, giving the disorder its name. This sickle cell morphology is responsible for the profound and widespread pathologies associated with this disorder, such as vaso-occlusive crisis (VOC). How much of the clinical manifestation is due to sickled erythrocytes and what is due to the relative contributions of other elements in the blood, especially in the microcapillary circulation, is usually not visualized and quantified for each patient during clinical management. Here, we used a microfluidic microcirculation mimetic (MMM), which has 187 capillary-like constrictions, to impose deformations on erythrocytes of 25 SCD patients, visualizing and characterizing the morpho-rheological properties of the cells in normoxic, hypoxic (using sodium meta-bisulfite) and treatment conditions (using hydroxyurea). The MMM enabled a patient-specific quantification of shape descriptors (circularity and roundness) and transit time through the capillary constrictions, which are readouts for morpho-rheological properties implicated in VOC. Transit times varied significantly (p < 0.001) between patients. Our results demonstrate the feasibility of microfluidics-based monitoring of individual patients for personalized care in the context of SCD complications such as VOC, even in resource-constrained settings.

Rheological Impact of GBT1118 Cessation in a Sickle Mouse Model

In sickle cell disease (SCD), higher whole blood viscosity is a risk factor for vaso-occlusive crisis, avascular necrosis, and proliferative retinopathy. Blood viscosity is strongly impacted by hemoglobin (Hb) levels and red blood cell (RBC) deformability. Voxelotor is a hemoglobin S (HbS) polymerization inhibitor with anti-sickling properties that increases the Hb affinity for oxygen, thereby reducing HbS polymerization. In clinical trials, voxelotor increased Hb by an average of 1g/dl, creating concern that this rise in Hb could increase viscosity, particularly when the drug was cleared. To investigate this potential rebound hyperviscosity effect, we treated SCD mice with GBT1118, a voxelotor analog, and stopped the treatment to determine the effect on blood viscosity and RBC deformability under a range of oxygen concentrations. GBT1118 treatment increased Hb, improved RBC deformability by increasing the elongation index under normoxic (EImax) and hypoxic conditions (EImin), and decreased the point of sickling (PoS) without increasing blood viscosity. The anti-sickling effects and improvement of RBC deformability balanced the effect of increased Hb such that there was no increase in blood viscosity. Forty-eight hours after ceasing GBT1118, Hb declined from the rise induced by treatment, viscosity did not increase, and EImin remained elevated compared to control animals. Hb and PoS were not different from control animals, suggesting a return to native oxygen affinity and clearance of the drug. RBC deformability did not return to baseline, suggesting some residual rheological improvement. These data suggest that concerns regarding viscosity rise above pre-treatment levels upon sudden cessation of voxelotor are not warranted.

A portable impedance microflow cytometer for measuring cellular response to hypoxia

This article presents the development and testing of a low-cost (<$60), portable, electrical impedance-based microflow cytometer for single-cell analysis under a controlled oxygen microenvironment. The system is based on an AD5933 impedance analyzer chip, a microfluidic chip, and an Arduino microcontroller operated by a custom Android application. A representative case study on human red blood cells (RBCs) affected by sickle cell disease is conducted to demonstrate the capability of the cytometry system. Impedance values of sickle blood samples exhibit remarkable deviations from the common reference line obtained from two normal blood samples. Such deviation is quantified by a conformity score, which allows for the measurement of intrapatient and interpatient variations of sickle cell disease. A low conformity score under oxygenated conditions or drastically different conformity scores between oxygenated and deoxygenated conditions can be used to differentiate a sickle blood sample from normal. Furthermore, an equivalent circuit model of a suspended biological cell is used to interpret the electrical impedance of single flowing RBCs. In response to hypoxia treatment, all samples, regardless of disease state, exhibit significant changes in at least one single-cell electrical property, that is, cytoplasmic resistance and membrane capacitance. The overall response to hypoxia is less in normal cells than those affected by sickle cell disease, where the change in membrane capacitance varies from -23% to seven times as compared with -17% in normal cells. The results reported in this article suggest that the developed method of testing demonstrates the potential application for a low-cost screening technique for sickle cell disease and other diseases in the field and low-resource settings. The developed system and methodology can be extended to analyze cellular response to hypoxia in other cell types.

A comprehensive review of hydroxyurea for β-haemoglobinopathies: the role revisited during COVID-19 pandemic

Background

Hydroxyurea is one of the earliest drugs that showed promise in the management of haemoglobinopathies that include β-thalassaemia and sickle cell disease. Despite this, many aspects of hydroxyurea are either unknown or understudied; specifically, its usefulness in β-thalassaemia major and haemoglobin E β-thalassaemia is unclear. However, during COVID-19 pandemic, it has become a valuable adjunct to transfusion therapy in patients with β-haemoglobinopathies. In this review, we aim to explore the available in vitro and in vivo mechanistic data and the clinical utility of hydroxyurea in β-haemoglobinopathies with a special emphasis on its usefulness during the COVID-19 pandemic.

Main body

Hydroxyurea is an S-phase-specific drug that reversibly inhibits ribonucleoside diphosphate reductase enzyme which catalyses an essential step in the DNA biosynthesis. In human erythroid cells, it induces the expression of γ-globin, a fetal globin gene that is suppressed after birth. Through several molecular pathways described in this review, hydroxyurea exerts many favourable effects on the haemoglobin content, red blood cell indices, ineffective erythropoiesis, and blood rheology in patients with β-haemoglobinopathies. Currently, it is recommended for sickle cell disease and non-transfusion dependent β-thalassaemia. A number of clinical trials are ongoing to evaluate its usefulness in transfusion dependent β-thalassaemia. During the COVID-19 pandemic, it was widely used as an adjunct to transfusion therapy due to limitations in the availability of blood and logistical disturbances. Thus, it has become clear that hydroxyurea could play a remarkable role in reducing transfusion requirements of patients with haemoglobinopathies, especially when donor blood is a limited resource.

Conclusion

Hydroxyurea is a well-tolerated oral drug which has been in use for many decades. Through its actions of reversible inhibition of ribonucleoside diphosphate reductase enzyme and fetal haemoglobin induction, it exerts many favourable effects on patients with β-haemoglobinopathies. It is currently approved for the treatment of sickle cell disease and non-transfusion dependent β-thalassaemia. Also, there are various observations to suggest that hydroxyurea is an important adjunct in the treatment of transfusion dependent β-thalassaemia which should be confirmed by randomised clinical trials.

Standardized microfluidic assessment of red blood cell-mediated microcapillary occlusion: Association with clinical phenotype and hydroxyurea responsiveness in sickle cell disease

OBJECTIVES

We present a standardized in vitro microfluidic assay and Occlusion Index (OI) for the assessment of red blood cell (RBC)-mediated microcapillary occlusion and its clinical associations in sickle cell disease (SCD).

METHODS

Red blood cell mediated microcapillary occlusion represented by OI and its clinical associations were assessed for seven subjects with hemoglobin-SC disease (HbSC), 18 subjects with homozygous SCD (HbSS), and five control individuals (HbAA).

RESULTS

We identified two sub-populations with HbSS based on the OI distribution. HbSS subjects with relatively higher OIs had significantly lower hemoglobin levels, lower fetal hemoglobin (HbF) levels, and lower mean corpuscular volume (MCV), but significantly higher serum lactate dehydrogenase levels and absolute reticulocyte counts, compared to subjects with HbSS and lower OIs. HbSS subjects who had relatively higher OIs were more likely to have had a concomitant diagnosis of intrapulmonary shunting (IPS). Further, lower OI associated with hydroxyurea (HU) responsiveness in subjects with HbSS, as evidenced by significantly elevated HbF levels and MCV.

CONCLUSIONS

We demonstrated that RBC-mediated microcapillary occlusion and OI associated with subject clinical phenotype and HU responsiveness in SCD. The presented standardized microfluidic assay may be useful for evaluating clinical phenotype and assessing therapeutic outcomes in SCD, including emerging targeted and curative treatments that aim to improve RBC deformability and microcirculatory health.

A microfluidic computational fluid dynamics model for cellular interaction studies of sickle cell disease vaso-occlusions

Individuals with sickle cell disease are plagued with vaso-occlusions, chronic blockages within the vasculature. Several factors including stiffer sickle red blood cells and increased cell aggregation contribute to vaso-occlusion formation; however much remains to be understood. We present a computational fluid dynamics blood flow simulation within a microfluidic platform using the Carreau model and Murray's law. Vaso-occlusions form preferentially near bifurcations within 60 s in the sickle cell disease simulation. Velocity profiles and shear rates align with clinical and experimental reports. We assert that results from this study can be utilized to inform experimental investigations and microfluidic system design decisions.

Multi-GPU Immersed Boundary Method Hemodynamics Simulations

Large-scale simulations of blood flow that resolve the 3D deformation of each comprising cell are increasingly popular owing to algorithmic developments in conjunction with advances in compute capability. Among different approaches for modeling cell-resolved hemodynamics, fluid structure interaction (FSI) algorithms based on the immersed boundary method are frequently employed for coupling separate solvers for the background fluid and the cells within one framework. GPUs can accelerate these simulations; however, both current pre-exascale and future exascale CPU-GPU heterogeneous systems face communication challenges critical to performance and scalability. We describe, to our knowledge, the largest distributed GPU-accelerated FSI simulations of high hematocrit cell-resolved flows with over 17 million red blood cells. We compare scaling on a fat node system with six GPUs per node and on a system with a single GPU per node. Through comparison between the CPU- and GPU-based implementations, we identify the costs of data movement in multiscale multi-grid FSI simulations on heterogeneous systems and show it to be the greatest performance bottleneck on the GPU.

Factors Influencing Motivation and Engagement in Mobile Health Among Patients With Sickle Cell Disease in Low-Prevalence, High-Income Countries: Qualitative Exploration of Patient Requirements

Background

Sickle cell disease (SCD) is a hematological genetic disease affecting over 25 million people worldwide. The main clinical manifestations of SCD, hemolytic anemia and vaso-occlusion, lead to chronic pain and organ damages. With recent advances in childhood care, high-income countries have seen SCD drift from a disease of early childhood mortality to a neglected chronic disease of adulthood. In particular, coordinated, preventive, and comprehensive care for adults with SCD is largely underresourced. Consequently, patients are left to self-manage. Mobile health (mHealth) apps for chronic disease self-management are now flooding app stores. However, evidence remains unclear about their effectiveness, and the literature indicates low user engagement and poor adoption rates. Finally, few apps have been developed for people with SCD and none encompasses their numerous and complex self-care management needs.

Objective

This study aimed to identify factors that may influence the long-term engagement and user adoption of mHealth among the particularly isolated community of adult patients with SCD living in low-prevalence, high-income countries.

Methods

Semistructured interviews were conducted. Interviews were audiotaped, transcribed verbatim, and analyzed using thematic analysis. Analysis was informed by the Braun and Clarke framework and mapped to the COM-B model (capability, opportunity, motivation, and behavior). Results were classified into high-level functional requirements (FRs) and nonfunctional requirements (NFRs) to guide the development of future mHealth interventions.

Results

Overall, 6 males and 4 females were interviewed (aged between 21 and 55 years). Thirty FRs and 31 NFRs were extracted from the analysis. Most participants (8/10) were concerned about increasing their physical capabilities being able to stop pain symptoms quickly. Regarding the psychological capability aspects, all interviewees desired to receive trustworthy feedback on their self-care management practices. About their physical opportunities, most (7/10) expressed a strong desire to receive alerts when they would reach their own physiological limitations (ie, during physical activity). Concerning social opportunity, most (9/10) reported wanting to learn about the self-care practices of other patients. Relating to motivational aspects, many interviewees (6/10) stressed their need to learn how to avoid the symptoms and live as normal a life as possible. Finally, NFRs included inconspicuousness and customizability of user experience, automatic data collection, data shareability, and data privacy.

Conclusions

Our findings suggest that motivation and engagement with mHealth technologies among the studied population could be increased by providing features that clearly benefit them. Self-management support and self-care decision aid are patients’ major demands. As the complexity of SCD self-management requires a high cognitive load, pervasive health technologies such as wearable sensors, implantable devices, or inconspicuous conversational user interfaces should be explored to ease it. Some of the required technologies already exist but must be integrated, bundled, adapted, or improved to meet the specific needs of people with SCD.

Blood rheology biomarkers in sickle cell disease

Sickle cell disease (SCD) is the most common inherited blood disorder, affecting approximately 100,000 patients in the U.S. and millions more worldwide. Patients with SCD experience a wide range of clinical complications, including frequent pain crises, stroke, and early mortality, all originating from a single-point mutation in the β-globin subunit. The RBC changes resulting from the sickle mutation lead to a host of rheological abnormalities that diminish microvascular blood flow, and produce severe anemia due to RBC hemolysis, and ischemia from vaso-occlusion initiated by sticky, rigid sickle RBCs. While the pathophysiology and mechanisms of SCD have been investigated for many years, therapies to treat the disease are limited. In addition to RBC transfusion, there are only two US Food and Drug Administration (FDA)-approved drugs to ameliorate SCD complications: hydroxyurea (HU) and L-glutamine (Endari™). The only curative therapy currently available is allogeneic hematopoietic stem cell transplantation (HSCT), which is generally reserved for individuals with a matched related donor, comprising only 10–15% of the total SCD population. Potentially curative advanced gene therapy approaches for SCD are under investigation in ongoing clinical trials. The ultimate goal of any curative treatment should be to repair the hemorheological abnormalities caused by SCD, and thus normalize blood flow and prevent clinical complications. Our mini-review highlights a set of key hemorheological biomarkers (and the current and emerging technologies used to measure them) that may be used to guide the development of novel curative and palliative therapies for SCD, and functionally assess outcomes.

Impact statement

Severe impairment of blood rheology is the hallmark of SCD pathophysiology, and one of the key factors predisposing SCD patients to pain crises, organ damage, and early mortality. As novel therapies emerge to treat or cure SCD, it is crucial that these treatments are functionally evaluated for their effect on blood rheology. This review describes a comprehensive panel of rheological biomarkers, their clinical uses, and the technologies used to obtain them. The described technologies can produce highly sensitive measurements of the ability of current treatments to improve blood rheology of SCD patients. The goal of curative therapies should be to achieve blood rheology biomarkers measurements in the range of sickle cell trait individuals (HbAS). The use of the panel of rheological biomarkers proposed in this review could significantly accelerate the development, optimization, and clinical translation of novel therapies for SCD.

Deformation of a Red Blood Cell in a Narrow Rectangular Microchannel

The deformability of a red blood cell (RBC) is one of the most important biological parameters affecting blood flow, both in large arteries and in the microcirculation, and hence it can be used to quantify the cell state. Despite numerous studies on the mechanical properties of RBCs, including cell rigidity, much is still unknown about the relationship between deformability and the configuration of flowing cells, especially in a confined rectangular channel. Recent computer simulation techniques have successfully been used to investigate the detailed behavior of RBCs in a channel, but the dynamics of a translating RBC in a narrow rectangular microchannel have not yet been fully understood. In this study, we numerically investigated the behavior of RBCs flowing at different velocities in a narrow rectangular microchannel that mimicked a microfluidic device. The problem is characterized by the capillary number C a , which is the ratio between the fluid viscous force and the membrane elastic force. We found that confined RBCs in a narrow rectangular microchannel maintained a nearly unchanged biconcave shape at low C a , then assumed an asymmetrical slipper shape at moderate C a , and finally attained a symmetrical parachute shape at high C a . Once a RBC deformed into one of these shapes, it was maintained as the final stable configurations. Since the slipper shape was only found at moderate C a , measuring configurations of flowing cells will be helpful to quantify the cell state.

Faster Sickling Kinetics and Sickle Cell Shape Evolution during Repeated Deoxygenation and Oxygenation Cycles

Kinetics of cell sickling and morphological change have been recognized as important parameters that are correlated closely with altered blood rheology and vasoocclusion in microcirculation. A microfluidic transient hypoxia assay was developed to create repeated hypoxia-normoxia cycles for real time observation of repetitive sickling and unsickling of freely suspended red blood cells (RBCs) from sickle cell disease patients. Cell sickling behavior and kinetics were found to be influenced by its previous sickling-unsickling processes accumulatively, where those sickled RBCs that had a history of sickling in a previous hypoxia cycle would sickle again in subsequent hypoxia/sickling cycles and the collective sickling kinetics became progressively faster (with reduced delay time and higher sickled fraction versus deoxygenation time). Individual sickled RBCs would sickle into drastically different shapes randomly in subsequent hypoxia/sickling cycles, however, the collective shape distribution retained similar characteristics. These observations indicate a gradual worsening trend in sickling kinetics over repeated hypoxia cycles, as well as a relatively stable collective shape characteristics within a limited number of hypoxia-normoxia cycles.

Quantifying Shear-Induced Deformation and Detachment of Individual Adherent Sickle Red Blood Cells

Vaso-occlusive crisis, a common painful complication of sickle cell disease, is a complex process triggered by intercellular adhesive interactions among blood cells and the endothelium in all human organs (e.g., the oxygen-rich lung as well as hypoxic systems such as liver and kidneys). We present a combined experimental-computational study to quantify the adhesive characteristics of sickle mature erythrocytes (SMEs) and irreversibly sickled cells (ISCs) under flow conditions mimicking those in postcapillary venules. We employed an in vitro microfluidic cell adherence assay, which is coated uniformly with fibronectin. We investigated the adhesion dynamics of SMEs and ISCs in pulsatile flow under well-controlled hypoxic conditions, inferring the cell adhesion strength by increasing the flow rate (or wall shear stress (WSS)) until the onset of cell detachment. In parallel, we performed simulations of individual SMEs and ISCs under shear. We introduced two metrics to quantify the adhesion process, the cell aspect ratio (AR) as a function of WSS and its rate of change (the dynamic deformability index). We found that the AR of SMEs decreases significantly with the increase of WSS, consistent between the experiments and simulations. In contrast, the AR of ISCs remains constant in time and independent of the flow rate. The critical WSS value for detaching a single SME in oxygenated state is in the range of 3.9-5.5 Pa depending on the number of adhesion sites; the critical WSS value for ISCs is lower than that of SMEs. Our simulations show that the critical WSS value for SMEs in deoxygenated state is above 6.2 Pa (multiple adhesion sites), which is greater than their oxygenated counterparts. We investigated the effect of cell shear modulus on the detachment process; we found that for the same cell adhesion spring constant, the higher shear modulus leads to an earlier cell detachment from the functionalized surface. These findings may aid in the understanding of individual roles of sickle cell types in sickle cell disease vaso-occlusion.

Shape oscillations of single blood drops: applications to human blood and sickle cell disease

Sickle cell disease (SCD) is an inherited blood disorder associated with severe anemia, vessel occlusion, poor oxygen transport and organ failure. The presence of stiff and often sickle-shaped red blood cells is the hallmark of SCD and is believed to contribute to impaired blood rheology and organ damage. Most existing measurement techniques of blood and red blood cell physical properties require sample contact and/or large sample volume, which is problematic for pediatric patients. Acoustic levitation allows rheological measurements in a single drop of blood, simultaneously eliminating the need for both contact containment and manipulation of samples. The technique shows that the shape oscillation of blood drops is able to assess blood viscosity in normal and SCD blood and demonstrates an abnormally increased viscosity in SCD when compared with normal controls. Furthermore, the technique is sensitive enough to detect viscosity changes induced by hydroxyurea treatment, and their dependence on the total fetal hemoglobin content of the sample. Thus this technique may hold promise as a monitoring tool for assessing changes in blood rheology in sickle cell and other hematological diseases.

GBT440 improves red blood cell deformability and reduces viscosity of sickle cell blood under deoxygenated conditions

BACKGROUND

In sickle cell disease (SCD), polymerization of hemoglobin S (HbS) leads to the formation of rigid, non-deformable sickled RBCs. Loss of RBC deformability, sickling and irreversible membrane damage causes abnormal blood rheology, and increases viscosity which contributes to vasoocclusion and other SCD pathophysiology. GBT440 (generic name voxelotor) is a novel anti-polymerization and anti-sickling agent currently undergoing clinical evaluation for the treatment of SCD.

OBJECTIVE

The purpose of this study was to determine the effects of GBT440 on deformability of sickle RBCs (SS RBCs) and the hyperviscosity of sickle cell blood (SS blood).

METHODS

The mechanical and rheological properties of GBT440-treated SS RBCs were measured using micropipette and filtration techniques. The viscosity of sickle blood was measured using a Wells-Brookfield cone/plate viscometer.

RESULTS

GBT440 restored movement of deoxygenated SS RBCs through a gel filtration column and reduced the pressure required to pass SS RBCs through a polycarbonate filter. Moreover, GBT440 decreased the membrane shear elastic modulus of SS RBCs assessed via micropipette aspiration and reduced the hyperviscosity of SS blood under deoxygenated conditions.

CONCLUSIONS

GBT440 maintains SS RBC deformability and improves SS blood viscosity by inhibiting HbS polymerization under deoxygenated conditions. These results further support development of GBT440 as a disease-modifying agent in SCD patients.

Synergistic Integration of Laboratory and Numerical Approaches in Studies of the Biomechanics of Diseased Red Blood Cells

In red blood cell (RBC) disorders, such as sickle cell disease, hereditary spherocytosis, and diabetes, alterations to the size and shape of RBCs due to either mutations of RBC proteins or changes to the extracellular environment, lead to compromised cell deformability, impaired cell stability, and increased propensity to aggregate. Numerous laboratory approaches have been implemented to elucidate the pathogenesis of RBC disorders. Concurrently, computational RBC models have been developed to simulate the dynamics of RBCs under physiological and pathological conditions. In this work, we review recent laboratory and computational studies of disordered RBCs. Distinguished from previous reviews, we emphasize how experimental techniques and computational modeling can be synergically integrated to improve the understanding of the pathophysiology of hematological disorders.

Modeling erythrocyte electrodeformation in response to amplitude modulated electric waveforms

We present a comprehensive theoretical-experimental framework for quantitative, high-throughput study of cell biomechanics. An improved electrodeformation method has been developed by combing dielectrophoresis and amplitude shift keying, a form of amplitude modulation. This method offers a potential to fully control the magnitude and rate of deformation in cell membranes. In healthy human red blood cells, nonlinear viscoelasticity of cell membranes is obtained through variable amplitude load testing. A mathematical model to predict cellular deformations is validated using the experimental results of healthy human red blood cells subjected to various types of loading. These results demonstrate new capabilities of the electrodeformation technique and the validated mathematical model to explore the effects of different loading configurations on the cellular mechanical behavior. This gives it more advantages over existing methods and can be further developed to study the effects of strain rate and loading waveform on the mechanical properties of biological cells in health and disease.

Shear dependent red blood cell adhesion in microscale flow

Non-adherence and deformability are the key intrinsic biomechanical features of the red blood cell (RBC), which allow it to tightly squeeze and pass through even the narrowest of microcirculatory networks. Blockage of microcirculatory flow, also known as vaso-occlusion, is a consequence of abnormal cellular adhesion to the vascular endothelium. In sickle cell disease (SCD), an inherited anaemia, even though RBCs have been shown to be heterogeneous in adhesiveness and deformability, this has not been studied in the context of physiologically relevant dynamic shear gradients at the microscale. We developed a microfluidic system that simulates physiologically relevant shear gradients of microcirculatory blood flow at a constant single volumetric flow rate. Using this system, shear dependent adhesion of RBCs from 28 subjects with SCD and from 11 healthy subjects was investigated using vascular endothelial protein functionalized microchannels. We defined a new term, RBC Shear Gradient Microfluidic Adhesion (SiGMA) index to assess shear dependent RBC adhesion in a subject-specific manner. We have shown for the first time that shear dependent adhesion of RBCs is heterogeneous in a microfluidic flow model, which correlates clinically with inflammatory markers and iron overload in subjects with SCD. This study reveals the complex dynamic interactions between RBC-mediated microcirculatory occlusion and clinical outcomes in SCD. These interactions may also be relevant to other microcirculatory disorders and microvascular diseases.

Electrical impedance microflow cytometry with oxygen control for detection of sickle cells

Polymerization of intracellular sickle hemoglobin induced by low oxygen tension has been recognized as a primary determinant of the pathophysiologic manifestations in sickle cell disease. Existing flow cytometry techniques for detection of sickle cells are typically based on fluorescence markers or cellular morphological analysis. Using microfluidics and electrical impedance spectroscopy, we develop a new, label-free flow cytometry for non-invasive measurement of single cells under controlled oxygen level. We demonstrate the capability of this new technique by determining the electrical impedance differential of normal red blood cells obtained from a healthy donor and sickle cells obtained from three sickle cell patients, under normoxic and hypoxic conditions and at three different electrical frequencies, 156 kHz, 500 kHz and 3 MHz. Under normoxia, normal cells and sickle cells can be separated completely using electrical impedance at 156 kHz and 500 kHz but not at 3 MHz. Sickle cells, intra-patient and inter-patient show significantly different electrical impedance between normoxia and hypoxia at all three frequencies. This study shows a proof of concept that electrical impedance signal can be used as an indicator of the disease state of a red blood cell as well as the cell sickling events in sickle cell disease. Electrical impedance-based microflow cytometry with oxygen control is a new method that can be potentially used for sickle cell disease diagnosis and monitoring.

Dielectrophoresis Testing of Nonlinear Viscoelastic Behaviors of Human Red Blood Cells

Dielectrophoresis in microfluidics provides a useful tool to test biomechanics of living cells, regardless of surface charges on cell membranes. We have designed an experimental method to characterize the nonlinear viscoelastic behaviors of single cells using dielectrophoresis in a microfluidic channel. This method uses radio frequency, low voltage excitations through interdigitated microelectrodes, allowing probing multiple cells simultaneously with controllable load levels. Dielectrophoretic force was calibrated using a triaxial ellipsoid model. Using a Kelvin-Voigt model, the nonlinear shear moduli of cell membranes were determined from the steady-state deformations of red blood cells in response to a series of electric field strengths. The nonlinear elastic moduli of cell membranes ranged from 6.05 μN/m to up to 20.85 μN/m, which were identified as a function of extension ratio, rather than the lumped-parameter models as reported in the literature. Value of the characteristic time of the extensional recovery of cell membranes initially deformed to varied extent was found to be about 0.14 s. Shear viscosity of cell membrane was estimated to be 0.8-2.9 (μN/m)·s. This method is particularly valuable for rapid, non-invasive probing of mechanical properties of living cells.

Computational Biomechanics of Human Red Blood Cells in Hematological Disorders

We review recent advances in multiscale modeling of the biomechanical characteristics of red blood cells (RBCs) in hematological diseases, and their relevance to the structure and dynamics of defective RBCs. We highlight examples of successful simulations of blood disorders including malaria and other hereditary disorders, such as sickle-cell anemia, spherocytosis, and elliptocytosis.

Mesoscopic Adaptive Resolution Scheme toward Understanding of Interactions between Sickle Cell Fibers

Understanding of intracellular polymerization of sickle hemoglobin (HbS) and subsequent interaction with the membrane of a red blood cell (RBC) is important to predict the altered morphologies and mechanical properties of sickle RBCs in sickle cell anemia. However, modeling the integrated processes of HbS nucleation, polymerization, HbS fiber interaction, and subsequent distortion of RBCs is challenging as they occur at multispatial scales, ranging from nanometers to micrometers. To make progress toward simulating the integrated processes, we propose a hybrid HbS fiber model, which couples fine-grained and coarse-grained HbS fiber models through a mesoscopic adaptive resolution scheme (MARS). To this end, we apply a microscopic model to capture the dynamic process of polymerization of HbS fibers, while maintaining the mechanical properties of polymerized HbS fibers by the mesoscopic model, thus providing a means of bridging the subcellular and cellular phenomena in sickle cell disease. At the subcellular level, this model can simulate HbS polymerization with preexisting HbS nuclei. At the cellular level, if combined with RBC models, the generated HbS fibers could be applied to study the morphologies and membrane stiffening of sickle RBCs. One important feature of the MARS is that it can be easily employed in other particle-based multiscale simulations where a dynamic coarse-graining and force-blending method is required. As demonstrations, we first apply the hybrid HbS fiber model to simulate the interactions of two growing fibers and find that their final configurations depend on the orientation and interaction distance between two fibers, in good agreement with experimental observations. We also model the formation of fiber bundles and domains so that we explore the mechanism that causes fiber branching.

Modeling of Biomechanics and Biorheology of Red Blood Cells in Type 2 Diabetes Mellitus

Erythrocytes in patients with type-2 diabetes mellitus (T2DM) are associated with reduced cell deformability and elevated blood viscosity, which contribute to impaired blood flow and other pathophysiological aspects of diabetes-related vascular complications. In this study, by using a two-component red blood cell (RBC) model and systematic parameter variation, we perform detailed computational simulations to probe the alteration of the biomechanical, rheological, and dynamic behavior of T2DM RBCs in response to morphological change and membrane stiffening. First, we examine the elastic response of T2DM RBCs subject to static tensile forcing and their viscoelastic relaxation response upon release of the stretching force. Second, we investigate the membrane fluctuations of T2DM RBCs and explore the effect of cell shape on the fluctuation amplitudes. Third, we subject the T2DM RBCs to shear flow and probe the effects of cell shape and effective membrane viscosity on their tank-treading movement. In addition, we model the cell dynamic behavior in a microfluidic channel with constriction and quantify the biorheological properties of individual T2DM RBCs. Finally, we simulate T2DM RBC suspensions under shear and compare the predicted viscosity with experimental measurements. Taken together, these simulation results and their comparison with currently available experimental data are helpful in identifying a specific parametric model-the first of its kind, to our knowledge-that best describes the main hallmarks of T2DM RBCs, which can be used in future simulation studies of hematologic complications of T2DM patients.

OpenRBC: A Fast Simulator of Red Blood Cells at Protein Resolution

We present OpenRBC, a coarse-grained molecular dynamics code, which is capable of performing an unprecedented in silico experiment-simulating an entire mammal red blood cell lipid bilayer and cytoskeleton as modeled by multiple millions of mesoscopic particles-using a single shared memory commodity workstation. To achieve this, we invented an adaptive spatial-searching algorithm to accelerate the computation of short-range pairwise interactions in an extremely sparse three-dimensional space. The algorithm is based on a Voronoi partitioning of the point cloud of coarse-grained particles, and is continuously updated over the course of the simulation. The algorithm enables the construction of the key spatial searching data structure in our code, i.e., a lattice-free cell list, with a time and space cost linearly proportional to the number of particles in the system. The position and the shape of the cells also adapt automatically to the local density and curvature. The code implements OpenMP parallelization and scales to hundreds of hardware threads. It outperforms a legacy simulator by almost an order of magnitude in time-to-solution and >40 times in problem size, thus providing, to our knowledge, a new platform for probing the biomechanics of red blood cells.

Probing the Twisted Structure of Sickle Hemoglobin Fibers via Particle Simulations

Polymerization of sickle hemoglobin (HbS) is the primary pathogenic event of sickle cell disease. For insight into the nature of the HbS polymer fiber formation, we develop a particle model-resembling a coarse-grained molecular model-constructed to match the intermolecular contacts between HbS molecules. We demonstrate that the particle model predicts the formation of HbS polymer fibers by attachment of monomers to rough fiber ends and the growth rate increases linearly with HbS concentration. We show that the characteristic 14-molecule fiber cross section is preserved during growth. We also correlate the asymmetry of the contact sites on the HbS molecular surface with the structure of the polymer fiber composed of seven helically twisted double strands. Finally, we show that the same asymmetry mediates the mechanical and structural properties of the HbS polymer fiber.

Photoplethysmography for the Assessment of Haemorheology

Haemorheology has been long identified as an early biomarker of a wide range of diseases, especially cardiovascular diseases. This study investigates for the first time the suitability of Photoplethysmography (PPG) as a non-invasive diagnostic method for haemorheological changes. The sensitivity of both PPG components (AC and DC) to changes in haemorheology were rigorously investigated in an in vitro experimental setup that mimics the human circulation. A custom-made reflectance PPG sensor, a pressure transducer and an ultrasonic Doppler flowmeter were used to map changes in flow dynamics and optical responses in an arterial model. The study investigated the effect of shear rates by varying fluid pumping frequencies using 4 set-points and the effect of clot formation using a chemical trigger. Both PPG_AC amplitudes and PPG_DC levels showed significant (p < 0.001) changes during the increase in shear rates and an immediate change after thromboplastin activation. The findings highlight that PPG has the potential to be used as a simple non-invasive method for the detection of blood characteristics, including disaggregation, radial migration and cross-linking fibrin formations. Such capability will enable the assessment of the effects of clotting-activators and anticoagulants (including non-pharmacological methods) and might aid in the early non-invasive assessment of cardiovascular pathologies.

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

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

Sickle cell disease is a genetic blood disease that causes vaso-occlusive pain crises. Here, we investigate the individual sickle cell behavior under controlled hypoxic conditions through patient-specific predictive computational simulations that are informed by companion microfluidic experiments. We identify the different dynamic behavior between individual sickle RBCs and normal ones in microfluidic flow, and analyze the hypoxia-induced alteration in individual cell behavior and single-cell capillary obstruction under physiological conditions.

Microfluidic experimental setup for adhesion and recovery measurements of red blood cells in sickle cell disease

Current microfluidic assays, which aim at quantifying mechanical properties of sickle cell red blood cells (SS-RBCs), suffer from a number of drawbacks in functionalization and flow control. Specifically, physical adsorption functionalization techniques produce inconsistent functional surfaces, and common volumetric flow pumps cannot be used to adjust the flow inside microchannels with minimal delay. We have designed an experimental setup that alleviates these complications by implementing aspiration for microchannel assembly that enables the use of most functionalization techniques and a pressure controller that allows instant and precise changes in the microchannel flow. Utilizing this setup, we have quantified SS-RBC adhesion to the integrin αvβ3, a specific adhesion protein expressed on the endothelium, as well as measured the shear modulus and viscosity of the SS-RBC plasma membrane.

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.

MD/DPD Multiscale Framework for Predicting Morphology and Stresses of Red Blood Cells in Health and Disease

Healthy red blood cells (RBCs) have remarkable deformability, squeezing through narrow capillaries as small as 3 microns in diameter without any damage. However, in many hematological disorders the spectrin network and lipid bilayer of diseased RBCs may be significantly altered, leading to impaired functionality including loss of deformability. We employ a two-component whole-cell multiscale model to quantify the biomechanical characteristics of the healthy and diseased RBCs, including Plasmodium falciparum-infected RBCs (Pf-RBCs) and defective RBCs in hereditary disorders, such as spherocytosis and elliptocytosis. In particular, we develop a two-step multiscale framework based on coarse-grained molecular dynamics (CGMD) and dissipative particle dynamics (DPD) to predict the static and dynamic responses of RBCs subject to tensile forcing, using experimental information only on the structural defects in the lipid bilayer, cytoskeleton, and their interaction. We first employ CGMD on a small RBC patch to compute the shear modulus, bending stiffness, and network parameters, which are subsequently used as input to a whole-cell DPD model to predict the RBC shape and corresponding stress field. For Pf-RBCs at trophozoite and schizont stages, the presence of cytoadherent knobs elevates the shear response in the lipid bilayer and stiffens the RBC membrane. For RBCs in spherocytosis and elliptocytosis, the bilayer-cytoskeleton interaction is weakened, resulting in substantial increase of the tensile stress in the lipid bilayer. Furthermore, we investigate the transient behavior of stretching deformation and shape relaxation of the normal and defective RBCs. Different from the normal RBCs possessing high elasticity, our simulations reveal that the defective RBCs respond irreversibly, i.e., they lose their ability to recover the normal biconcave shape in successive loading cycles of stretching and relaxation. Our findings provide fundamental insights into the microstructure and biomechanics of RBCs, and demonstrate that the two-step multiscale framework presented here can be used effectively for in silico studies of hematological disorders based on first principles and patient-specific experimental input at the protein level.

Sub-cellular modeling of platelet transport in blood flow through microchannels with constriction

Platelet transport through arterial constrictions is one of the controlling processes influencing their adhesive functions and the formation of thrombi. We perform high-fidelity mesoscopic simulations of blood flow in microchannels with constriction, resembling arterial stenoses. The wall shear rates inside the constrictions reach levels as high as ≈8000 s(-1), similar to those encountered in moderate atherosclerotic plaques. Both red blood cells and platelets are resolved at sub-cellular resolution using the Dissipative Particle Dynamics (DPD) method. We perform a systematic study on the red blood cell and platelet transport by considering different levels of constriction, blood hematocrit and flow rates. We find that higher levels of constriction and wall shear rates lead to significantly enhanced margination of platelets, which may explain the experimental observations of enhanced post-stenosis platelet aggregation. We also observe similar margination effects for stiff particles of spherical shapes such as leukocytes. To our knowledge, such numerical simulations of dense blood through complex geometries have not been performed before, and our quantitative findings could shed new light on the associated physiological processes such as ATP release, plasma skimming, and thrombus formation.

Integrated multiscale biomaterials experiment and modelling: a perspective

Advances in multiscale models and computational power have enabled a broad toolset to predict how molecules, cells, tissues and organs behave and develop. A key theme in biological systems is the emergence of macroscale behaviour from collective behaviours across a range of length and timescales, and a key element of these models is therefore hierarchical simulation. However, this predictive capacity has far outstripped our ability to validate predictions experimentally, particularly when multiple hierarchical levels are involved. The state of the art represents careful integration of multiscale experiment and modelling, and yields not only validation, but also insights into deformation and relaxation mechanisms across scales. We present here a sampling of key results that highlight both challenges and opportunities for integrated multiscale experiment and modelling in biological systems.