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

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.

A framework of computer vision-enhanced microfluidic approach for automated assessment of the transient sickling kinetics in sickle red blood cells

The occurrence of vaso-occlusive crisis greatly depends on the competition between the sickling delay time and the transit time of individual sickle cells, i.e., red blood cells (RBCs) from sickle cell disease (SCD) patients, while they are traversing the circulatory system. Many drugs for treating SCD work by inhibiting the polymerization of sickle hemoglobin (HbS), effectively delaying the sickling process in sickle cells (SS RBCs). Most previous studies on screening anti-sickling drugs, such as voxelotor, rely on in vitro testing of sickling characteristics, often conducted under prolonged deoxygenation for up to 1 hour. However, since the microcirculation of RBCs typically takes less than 1 minute, the results of these studies may be less accurate and less relevant for in vitro-in vivo correlation. In our current study, we introduce a computer vision-enhanced microfluidic framework designed to automatically capture the transient sickling kinetics of SS RBCs within a 1-min timeframe. Our study has successfully detected differences in the transient sickling kinetics between vehicle control and voxelotor-treated SS RBCs. This approach has the potential for broader applications in screening anti-sickling therapies.

Single-cell characterization of deformation and dynamics of mesenchymal stem cells in microfluidic systems: A computational study

Understanding the homing dynamics of individual mesenchymal stem cells (MSCs) in physiologically relevant microenvironments is crucial for improving the efficacy of MSC-based therapies for therapeutic and targeting purposes. This study investigates the passive homing behavior of individual MSCs in micropores that mimic interendothelial clefts through predictive computational simulations informed by previous microfluidic experiments. Initially, we quantified the size-dependent behavior of MSCs in micropores and elucidated the underlying mechanisms. Subsequently, we analyzed the shape deformation and traversal dynamics of each MSC. In addition, we conducted a systematic investigation to understand how the mechanical properties of MSCs impact their traversal process. We considered geometric and mechanical parameters, such as reduced cell volume, cell-to-nucleus diameter ratio, and cytoskeletal prestress states. Furthermore, we quantified the changes in the MSC traversal process and identified the quantitative limits in their response to variations in micropore length. Taken together, the computational results indicate the complex dynamic behavior of individual MSCs in the confined microflow. This finding offers an objective way to evaluate the homing ability of MSCs in an interendothelial-slit-like microenvironment.

Measuring cell deformation by microfluidics

Microfluidics is an increasingly popular method for studying cell deformation, with various applications in fields such as cell biology, biophysics, and medical research. Characterizing cell deformation offers insights into fundamental cell processes, such as migration, division, and signaling. This review summarizes recent advances in microfluidic techniques for measuring cellular deformation, including the different types of microfluidic devices and methods used to induce cell deformation. Recent applications of microfluidics-based approaches for studying cell deformation are highlighted. Compared to traditional methods, microfluidic chips can control the direction and velocity of cell flow by establishing microfluidic channels and microcolumn arrays, enabling the measurement of cell shape changes. Overall, microfluidics-based approaches provide a powerful platform for studying cell deformation. It is expected that future developments will lead to more intelligent and diverse microfluidic chips, further promoting the application of microfluidics-based methods in biomedical research, providing more effective tools for disease diagnosis, drug screening, and treatment.

Microfluidic methods to advance mechanistic understanding and translational research in sickle cell disease

Sickle cell disease (SCD) is caused by a single point mutation in the β-globin gene of hemoglobin, which produces an altered sickle hemoglobin (HbS). The ability of HbS to polymerize under deoxygenated conditions gives rise to chronic hemolysis, oxidative stress, inflammation, and vaso-occlusion. Herein, we review recent findings using microfluidic technologies that have elucidated mechanisms of oxygen-dependent and -independent induction of HbS polymerization and how these mechanisms elicit the biophysical and inflammatory consequences in SCD pathophysiology. We also discuss how validation and use of microfluidics in SCD provides the opportunity to advance development of numerous therapeutic strategies, including curative gene therapies.

Recent Advances in Computational Modeling of Biomechanics and Biorheology of Red Blood Cells in Diabetes

Diabetes mellitus, a metabolic disease characterized by chronically elevated blood glucose levels, affects about 29 million Americans and more than 422 million adults all over the world. Particularly, type 2 diabetes mellitus (T2DM) accounts for 90-95% of the cases of vascular disease and its prevalence is increasing due to the rising obesity rates in modern societies. Although multiple factors associated with diabetes, such as reduced red blood cell (RBC) deformability, enhanced RBC aggregation and adhesion to the endothelium, as well as elevated blood viscosity are thought to contribute to the hemodynamic impairment and vascular occlusion, clinical or experimental studies cannot directly quantify the contributions of these factors to the abnormal hematology in T2DM. Recently, computational modeling has been employed to dissect the impacts of the aberrant biomechanics of diabetic RBCs and their adverse effects on microcirculation. In this review, we summarize the recent advances in the developments and applications of computational models in investigating the abnormal properties of diabetic blood from the cellular level to the vascular level. We expect that this review will motivate and steer the development of new models in this area and shift the attention of the community from conventional laboratory studies to combined experimental and computational investigations, aiming to provide new inspirations for the development of advanced tools to improve our understanding of the pathogenesis and pathology of T2DM.

Computational modeling of biomechanics and biorheology of heated red blood cells

Because of their compromised deformability, heat denatured erythrocytes have been used as labeled probes to visualize spleen tissue or to assess the ability of the spleen to retain stiff red blood cells (RBCs) for over three decades, e.g., see Looareesuwan et al. N. Engl. J. Med. (1987). Despite their good accessibility, it is still an open question how heated RBCs compare to certain diseased RBCs in terms of their biomechanical and biorheological responses, which may undermine their effective usage and even lead to misleading experimental observations. To help answering this question, we perform a systematic computational study of the hemorheological properties of heated RBCs with several physiologically relevant static and hemodynamic settings, including optical-tweezers test, relaxation of prestretched RBCs, RBC traversal through a capillary-like channel and a spleen-like slit, and a viscometric rheology test. We show that our in silico RBC models agree well with existing experiments. Moreover, under static tests, heated RBCs exhibit deformability deterioration comparable to certain disease-impaired RBCs such as those in malaria. For RBC traversal under confinement (through microchannel or slit), heated RBCs show prolonged transit time or retention depending on the level of confinement and heating procedure, suggesting that carefully heat-treated RBCs may be useful for studying splenic- or vaso-occlusion in vascular pathologies. For the rheology test, we expand the existing bulk viscosity data of heated RBCs to a wider range of shear rates (1-1000 s^-1) to represent most pathophysiological conditions in macro- or microcirculation. Although heated RBC suspension shows elevated viscosity comparable to certain diseased RBC suspensions under relatively high shear rates (100-1000 s^-1), they underestimate the elevated viscosity (e.g., in sickle cell anemia) at low shear rates (<10 s^-1). Our work provides mechanistic rationale for selective usage of heated RBC as a potentially useful model for studying the abnormal traversal dynamics and hemorheology in certain blood disorders.

Microfluidics in Sickle Cell Disease Research: State of the Art and a Perspective Beyond the Flow Problem

Sickle cell disease (SCD) is the monogenic hemoglobinopathy where mutated sickle hemoglobin molecules polymerize to form long fibers under deoxygenated state and deform red blood cells (RBCs) into predominantly sickle form. Sickled RBCs stick to the vascular bed and obstruct blood flow in extreme conditions, leading to acute painful vaso-occlusion crises (VOCs) – the leading cause of mortality in SCD. Being a blood disorder of deformed RBCs, SCD manifests a wide-range of organ-specific clinical complications of life (in addition to chronic pain) such as stroke, acute chest syndrome (ACS) and pulmonary hypertension in the lung, nephropathy, auto-splenectomy, and splenomegaly, hand-foot syndrome, leg ulcer, stress erythropoiesis, osteonecrosis and osteoporosis. The physiological inception for VOC was initially thought to be only a fluid flow problem in microvascular space originated from increased viscosity due to aggregates of sickled RBCs; however, over the last three decades, multiple molecular and cellular mechanisms have been identified that aid the VOC in vivo. Activation of adhesion molecules in vascular endothelium and on RBC membranes, activated neutrophils and platelets, increased viscosity of the blood, and fluid physics driving sickled and deformed RBCs to the vascular wall (known as margination of flow) – all of these come together to orchestrate VOC. Microfluidic technology in sickle research was primarily adopted to benefit from mimicking the microvascular network to observe RBC flow under low oxygen conditions as models of VOC. However, over the last decade, microfluidics has evolved as a valuable tool to extract biophysical characteristics of sickle red cells, measure deformability of sickle red cells under simulated oxygen gradient and shear, drug testing, in vitro models of intercellular interaction on endothelialized or adhesion molecule-functionalized channels to understand adhesion in sickle microenvironment, characterizing biomechanics and microrheology, biomarker identification, and last but not least, for developing point-of-care diagnostic technologies for low resource setting. Several of these platforms have already demonstrated true potential to be translated from bench to bedside. Emerging microfluidics-based technologies for studying heterotypic cell–cell interactions, organ-on-chip application and drug dosage screening can be employed to sickle research field due to their wide-ranging advantages.

Bibliometric analysis of global sickle cell disease research from 1997 to 2017

Introduction

Sickle cell disease is an autosomal recessive genetic disease caused by a single point mutation in the β-globin chain of the hemoglobin. It has been recognized by the World Health Organization as a public health priority since 2006.

Methods

The Scopus database was used in this study with the search descriptors: “sickle cell” and “sickle cell disease”. We applied common bibliometric indicators to evaluate the trend in scientific literature in sickle cell disease research.

Results

We retrieved a total of 19,921 pieces of scientific literature in the repertoire from 1997 to 2017. The Price law was fulfilled in the trend of production of scientific literature on SCD as the growth of scientific literature was more exponential (r = 0.9751; r² = 0.9509) than linear (r = 0.9721; r² = 0.9449). We observed a duplication time of 4.52 years. The Bradford core was made up of 69 journals with Blood at the top, publishing the greatest number of articles. The most productive institutions were mostly United States agencies and hospitals. The United States was the most productive country. The National Institute of Health was the most productive institution and also had the highest number of citations. Vichinsky E was the most productive author, while the most cited article was published by Circulation.

Conclusion

The growth of scientific literature in Sickle cell disease was found to be high. However, the exponential growth trend shows a “yet-to-be-explored” area of research. This study will be useful for physicians, researchers, research funders and policy-cum-decision makers.

Quantifying Fibrinogen-Dependent Aggregation of Red Blood Cells in Type 2 Diabetes Mellitus

Fibrinogen is regarded as the main glycoprotein in the aggregation of red blood cells (RBCs), a normally occurring phenomenon that has a major impact on blood rheology and hemodynamics, especially under pathological conditions, including type 2 diabetes mellitus (T2DM). In this study, we investigate the fibrinogen-dependent aggregation dynamics of T2DM RBCs through patient-specific predictive computational simulations that invoke key parameters derived from microfluidic experiments. We first calibrate our model parameters at the doublet (a rouleau consisting of two aggregated RBCs) level for healthy blood samples by matching the detaching force required to fully separate RBC doublets with measurements using atomic force microscopy and optical tweezers. Using results from companion microfluidic experiments that also provide in vitro quantitative information on cell-cell adhesive dynamics, we then quantify the rouleau dissociation dynamics at the doublet and multiplet (a rouleau consisting of three or more aggregated RBCs) levels for obese patients with or without T2DM. Specifically, we examine the rouleau breakup rate when it passes through microgates at doublet level and investigate the effect of rouleau alignment in altering its breakup pattern at multiplet level. This study seamlessly integrates in vitro experiments and simulations and consequently enhances our understanding of the complex cell-cell interaction, highlighting the importance of the aggregation and disaggregation dynamics of RBCs in patients at increased risk of microvascular complications.

Gene expression of HIF-1α and VEGF in response to hypoxia in sickle cell anaemia: Influence of hydroxycarbamide

Hypoxia and hemoglobin S polymerization are two triggers responsible for initiating erythrocyte sickling and the consequent clinical sickle cell anemia (SCA) events. The objective of this study was to investigate the expression of hypoxia-responsive genes in SCA, testing for correlation with the clinical-laboratorial characteristics of the patient and hydroxyurea therapy. Our results showed, for the first time, a significantly increased expression of HIF-1α and VEGF genes in patients with SCA and an inverse dose-response relationship with hydroxyurea therapy. These results suggest that hypoxic stress may be involved in both the severity of SCA and its response to treatment.

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.

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.

Red Blood Cell Adhesion to Heme-Activated Endothelial Cells Reflects Clinical Phenotype in Sickle Cell Disease

In sickle cell disease (SCD), 'disease severity' associates with increased RBC adhesion to quiescent endothelium, but the impact on activated endothelium is not known. Increased concentrations of free heme result from intravascular hemolysis in SCD. Heme is essential for aerobic metabolism, and plays an important role in numerous biological processes. Excess free heme induces reactive oxygen species generation and endothelial activation, which are associated with cardiovascular disorders including atherosclerosis, hypertension, and thrombosis. Here, we utilized an endothelialized microfluidic platform (Endothelium-on-a-chip) to assess adhesion of sickle hemoglobin-containing red blood cells (HbS RBCs), from adults with homozygous SCD, to heme-activated human endothelial cells (EC) in vitro. Confluent EC monolayers in microchannels were treated with pathophysiologically relevant levels of heme in order to simulate the highly hemolytic intravascular milieu seen in SCD. RBC adhesion to heme-activated ECs varied from subject to subject, and was associated with plasma markers of hemolysis (LDH) and reticulocytosis, thereby linking those RBCs that are most likely to adhere with those that are most likely to hemolyze. These results re-emphasize the critical contribution made by heterogeneous adhesive HbS RBCs to the pathophysiology of SCD. We found that adhesion of HbS RBCs to heme-activated ECs varied amongst individuals in the study population, and associated with biomarkers of hemolysis and inflammation, age, and a recent history of transfusion. Importantly, the microfluidic approach described herein holds promise as a clinically feasible Endothelium-on-a-chip platform with which to study complex heterocellular adhesive interactions in SCD. This article is protected by copyright. All rights reserved.

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.