In silico and in vitro study of the adhesion dynamics of erythrophagocytosis in sickle cell disease

The dynamics of red blood cells traversing slits of mechanical heart valves under high shear

Hemolysis, including subclinical hemolysis, is a potentially severe complications of mechanical heart valves (MHVs), which leads to shortened red blood cell (RBC) lifespan and hemolytic anemia. Serious hemolysis is usually associated with structural deterioration and regurgitation. However, the shear stress in MHVs' narrow leakage slits is much lower than the shear stress threshold causing hemolysis and the mechanisms in this context remain largely unclear. This study investigated the hemolysis mechanism of RBCs in cell-size slits under high shear rates by establishing in vitro microfluidic devices and a coarse-grained molecular dynamics (CGMD) model, considering both fluid and structural effects simultaneously. Microfluidic experiments and computational simulation revealed six distinct dynamic states of RBC traversal through MHVs' microscale slits under various shear rates and slit sizes. It elucidated that RBC dynamic states were influenced by not only by fluid forces but significantly by the compressive force of slit walls. The variation of the potential energy of the cell membrane indicated its stretching, deformation, and rupture during traversal, corresponding to the six dynamic states. The maximum forces exerted on membrane by water particles and slit walls directly determined membrane rupture, serving as a critical determinant. This analysis helps in understanding the contribution of the slit walls to membrane rupture and identifying the threshold force that leads to membrane rupture. The hemolysis mechanism of traversing microscale slits is revealed to effectively explain the occurrences of hemolysis and subclinical hemolysis.

Heme- and iron-activated macrophages in sickle cell disease: an updated perspective

PURPOSE OF REVIEW

Sickle cell disease (SCD) is a hereditary blood disorder due to a single-point mutation in the β-globin gene. The ensuing hemoglobin has the tendency to polymerize upon deoxygenation, leading to the typical sickle shape of red blood cells. While the primary pathology of sickle cell disease is a direct consequence of altered red blood cells, emerging evidence highlights the central role of macrophages in mediating hemoglobin scavenging, perpetuating oxidative stress and inflammation, and causing endothelial dysfunction and tissue remodeling.

RECENT FINDINGS

Recent research uncovered the impact of heme and iron overload on macrophage polarization and functions in sickle cell disease, and its implication for chronic inflammation and tissue damage in vital organs such as the liver, spleen, lungs and kidneys. By providing a thorough understanding of the dynamic interactions between macrophages and various cellular components within the sickle cell disease milieu, these studies have laid the foundation for the identification of macrophage-related cellular and molecular mechanisms potentially targetable for therapeutic purposes to attenuate sickle complications.

SUMMARY

This review provides a current update about recent discoveries on heme/iron-activated macrophages in SCD, shedding light on their critical role in disease pathophysiology. Ultimately, it proposes avenues for future research aimed at addressing the relevance of these cells for other sickle complications and at targeting them to mitigate disease morbidity and improve patient outcomes.

Splenic filtration of red blood cells in physiology, malaria and sickle cell disease

PURPOSE OF REVIEW

The human spleen clears the blood from circulating microorganisms and red blood cells (RBCs) displaying alterations. This review analyzes how generic mechanisms by which the spleen senses RBC, such pitting, trapping and erythrophagocytosis, impact the pathogenesis of two major spleen-related diseases, malaria and sickle cell disease (SCD).

RECENT FINDINGS

Scintigraphy, functional histology, comparison of circulating and splenic RBC, ex-vivo perfusion of human spleens and in-silico modeling enable relevant exploration of how the spleen retains and processes RBC in health and disease. Iterative cross-validations between medical observations, in-vitro experiments and in-silico modeling point to mechanical sensing of RBC as a central event in both conditions. Spleen congestion is a common pathogenic process explaining anemia and splenomegaly, the latter carrying a risk of severe complications such as acute splenic sequestration crisis and hypersplenism in SCD. Sickling of hemoglobin S-containing RBC may contribute to these complications without necessarily being the trigger.

SUMMARY

Ongoing progress in the exploration and understanding of spleen-related complications in malaria and SCD open the way to optimized prognosis evaluation and therapeutic applications.

Red blood cell passage through deformable interendothelial slits in the spleen: Insights into splenic filtration and hemodynamics

The spleen constantly clears altered red blood cells (RBCs) from the circulation, tuning the balance between RBC formation (erythropoiesis) and removal. The retention and elimination of RBCs occur predominantly in the open circulation of the spleen, where RBCs must cross submicron-wide inter-endothelial slits (IES). Several experimental and computational studies have illustrated the role of IES in filtrating the biomechanically and morphologically altered RBCs based on a rigid wall assumption. However, these studies also reported that when the size of IES is close to the lower end of clinically observed sizes (less than 0.5 μm), an unphysiologically large pressure difference across the IES is required to drive the passage of normal RBCs, sparking debates on the feasibility of the rigid wall assumption. In this work, We propose two deformable IES models, namely the passive model and the active model, aiming to explore the impact of the deformability of IES on the filtration function of the spleen. In the passive model, we implement the worm-like string model to depict the IES's deformation as it interacts with blood plasma and allows RBC to traverse. In contrast, the active model involved regulating the IES deformation based on the local pressure surrounding the slit. To demonstrate the validity of the deformable model, we simulate the filtration of RBCs with varied size and stiffness by IES under three scenarios: (1) a single RBC traversing a single slit; (2) a suspension of RBCs traversing an array of slits, mimicking in vitro spleen-on-a-chip experiments; (3) RBC suspension passing through the 3D spleen filtration unit known as'the splenon'. Our simulation results of RBC passing through a single slit show that the deformable IES model offers more accurate predictions of the critical cell surface area to volume ratio that dictate the removal of aged RBCs from circulation compared to prior rigid-wall models. Our biophysical models of the spleen-on-a-chip indicate a hierarchy of filtration function stringency: rigid model > passive model > active model, providing a possible explanation of the filtration function of IES. We also illustrate that the biophysical model of 'the splenon' enables us to replicate the ex vivo experiments involving spleen filtration of malaria-infected RBCs. Taken together, our simulation findings indicate that the deformable IES model could serve as a mesoscopic representation of spleen filtration function closer to physiological reality, addressing questions beyond the scope of current experimental and computational models and enhancing our understanding of the fundamental flow dynamics and mechanical clearance processes within in the human spleen.

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.

Interplay between endothelial glycocalyx layer and red blood cell in microvascular blood flow: A numerical study

The endothelial glycocalyx layer (EGL) plays a crucial role in regulating blood flow in microvessels. Experimental evidence suggests that there is greater blood flow resistance in vivo compared to in vitro, partially due to the presence of the EGL. However, the complex relationship between EGL deformation and blood cell behavior in shear flow and its quantification remains incompletely understood. To address this gap, we employ a particle-based numerical simulation technique to examine the interaction of the EGL with flowing red blood cells (RBCs) in microtubes. We examine changes in EGL deformation in response to variations in shear rate, EGL graft density, and contour height. Our results indicate that the alterations in EGL height are influenced by the mechanical properties of the EGL, flow conditions, and the RBC-EGL interaction. The flowing RBC compresses the EGL, causing a notable reduction in EGL height near the RBC flow. Additionally, we find that the presence of the EGL in the microtube results in increased RBC deformation and a wider gap between the RBC and tube wall due to spatial occupancy. The significant impact of the EGL on RBC flow is particularly evident in microtubes with diameters ranging from 7 to 10µm, a range consistent with notable differences in vascular flow resistance observed between in vivo and in vitro experiments. The simulation results shed insight on the dynamic interplay between RBC and the EGL in microvascular blood flow.

Two-component macrophage model for active phagocytosis with pseudopod formation

Macrophage phagocytosis is critical for the immune response, homeostasis regulation, and tissue repair. This intricate process involves complex changes in cell morphology, cytoskeletal reorganization, and various receptor-ligand interactions controlled by mechanical constraints. However, there is a lack of comprehensive theoretical and computational models that investigate the mechanical process of phagocytosis in the context of cytoskeletal rearrangement. To address this issue, we propose a novel coarse-grained mesoscopic model that integrates a fluid-like cell membrane and a cytoskeletal network to study the dynamic phagocytosis process. The growth of actin filaments results in the formation of long and thin pseudopods, and the initial cytoskeleton can be disassembled upon target entry and reconstructed after phagocytosis. Through dynamic changes in the cytoskeleton, our macrophage model achieves active phagocytosis by forming a phagocytic cup utilizing pseudopods in two distinct ways. We have developed a new algorithm for modifying membrane area to prevent membrane rupture and ensure sufficient surface area during phagocytosis. In addition, the bending modulus, shear stiffness, and cortical tension of the macrophage model are investigated through computation of the axial force for the tubular structure and micropipette aspiration. With this model, we simulate active phagocytosis at the cytoskeletal level and investigate the mechanical process during the dynamic interplay between macrophage and target particles.

A Simulation of the Mechanical Testing of the Cell Membrane and Cytoskeleton

Cell models play a crucial role in analyzing the mechanical response of cells and quantifying cellular damage incurred during micromanipulation. While traditional models can capture the overall mechanical behavior of cells, they often lack the ability to discern among distinct cellular components. Consequently, by employing dissipative particle dynamics, this study constructed a triangular network-like representation of the cell membrane along with cross-linked cytoskeletal chains. The mechanical properties of both the membrane and cytoskeleton were then analyzed through a series of simulated mechanical tests, validated against real-world experiments. The investigation utilized particle-tracking rheology to monitor changes in the mean square displacements of membrane particles over time, facilitating the analysis of the membrane's storage and loss moduli. Additionally, the cytoskeletal network's storage and loss moduli were examined via a double-plate oscillatory shear experiment. The simulation results revealed that both the membrane and cytoskeleton exhibit viscoelastic behavior, as evidenced by the power-law dependency of their storage and loss moduli on frequency. Furthermore, indentation and microinjection simulations were conducted to examine the overall mechanical properties of cells. In the indentation experiments, an increase in the shear modulus of the membrane's WLCs correlated with a higher Young's modulus for the entire cell. Regarding the microinjection experiment, augmenting the microinjection speed resulted in reduced deformation of the cell at the point of membrane rupture and a lower percentage of high strain.

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.

Signaling-biophysical modeling unravels mechanistic control of red blood cell phagocytosis by macrophages in sickle cell disease

Red blood cell (RBC) aging manifests through progressive changes in cell morphology, rigidity, and expression of membrane proteins. To maintain the quality of circulating blood, splenic macrophages detect the biochemical signals and biophysical changes of RBCs and selectively clear them through erythrophagocytosis. In sickle cell disease (SCD), RBCs display alterations affecting their interaction with macrophages, leading to aberrant phagocytosis that may cause life-threatening spleen sequestration crises. To illuminate the mechanistic control of RBC engulfment by macrophages in SCD, we integrate a system biology model of RBC-macrophage signaling interactions with a biophysical model of macrophage engulfment, as well as in vitro phagocytosis experiments using the spleen-on-a-chip technology. Our modeling framework accurately predicts the phagocytosis dynamics of RBCs under different disease conditions, reveals patterns distinguishing normal and sickle RBCs, and identifies molecular targets including Src homology 2 domain-containing protein tyrosine phosphatase-1 (SHP1) and cluster of differentiation 47 (CD47)/signal regulatory protein α (SIRPα) as therapeutic targets to facilitate the controlled clearance of sickle RBCs in the spleen.

A combined computational and experimental investigation of the filtration function of splenic macrophages in sickle cell disease

Being the largest lymphatic organ in the body, the spleen also constantly controls the quality of red blood cells (RBCs) in circulation through its two major filtration components, namely interendothelial slits (IES) and red pulp macrophages. In contrast to the extensive studies in understanding the filtration function of IES, fewer works investigate how the splenic macrophages retain the aged and diseased RBCs, i.e., RBCs in sickle cell disease (SCD). Herein, we perform a computational study informed by companion experiments to quantify the dynamics of RBCs captured and retained by the macrophages. We first calibrate the parameters in the computational model based on microfluidic experimental measurements for sickle RBCs under normoxia and hypoxia, as those parameters are not available in the literature. Next, we quantify the impact of key factors expected to dictate the RBC retention by the macrophages in the spleen, namely, blood flow conditions, RBC aggregation, hematocrit, RBC morphology, and oxygen levels. Our simulation results show that hypoxic conditions could enhance the adhesion between the sickle RBCs and macrophages. This, in turn, increases the retention of RBCs by as much as four-fold, which could be a possible cause of RBC congestion in the spleen of patients with SCD. Our study on the impact of RBC aggregation illustrates a 'clustering effect', where multiple RBCs in one aggregate can make contact and adhere to the macrophages, leading to a higher retention rate than that resulting from RBC-macrophage pair interactions. Our simulations of sickle RBCs flowing past macrophages for a range of blood flow velocities indicate that the increased blood velocity could quickly attenuate the function of the red pulp macrophages on detaining aged or diseased RBCs, thereby providing a possible rationale for the slow blood flow in the open circulation of the spleen. Furthermore, we quantify the impact of RBC morphology on their tendency to be retained by the macrophages. We find that the sickle and granular-shaped RBCs are more likely to be filtered by macrophages in the spleen. This finding is consistent with the observation of low percentages of these two forms of sickle RBCs in the blood smear of SCD patients. Taken together, our experimental and simulation results aid in our quantitative understanding of the function of splenic macrophages in retaining the diseased RBCs and provide an opportunity to combine such knowledge with the current knowledge of the interaction between IES and traversing RBCs to apprehend the complete filtration function of the spleen in SCD.

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.