Computational Biomechanics of Human Red Blood Cells in Hematological Disorders

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.

Piezo1 Regulation Involves Lipid Domains and the Cytoskeleton and Is Favored by the Stomatocyte-Discocyte-Echinocyte Transformation

Piezo1 is a mechanosensitive ion channel required for various biological processes, but its regulation remains poorly understood. Here, we used erythrocytes to address this question since they display Piezo1 clusters, a strong and dynamic cytoskeleton and three types of submicrometric lipid domains, respectively enriched in cholesterol, GM1 ganglioside/cholesterol and sphingomyelin/cholesterol. We revealed that Piezo1 clusters were present in both the rim and the dimple erythrocyte regions. Upon Piezo1 chemical activation by Yoda1, the Piezo1 cluster proportion mainly increased in the dimple area. This increase was accompanied by Ca2+ influx and a rise in echinocytes, in GM1/cholesterol-enriched domains in the dimple and in cholesterol-enriched domains in the rim. Conversely, the effects of Piezo1 activation were abrogated upon membrane cholesterol depletion. Furthermore, upon Piezo1-independent Ca2+ influx, the above changes were not observed. In healthy donors with a high echinocyte proportion, Ca2+ influx, lipid domains and Piezo1 fluorescence were high even at resting state, whereas the cytoskeleton membrane occupancy was lower. Accordingly, upon decreases in cytoskeleton membrane occupancy and stiffness in erythrocytes from patients with hereditary spherocytosis, Piezo1 fluorescence was increased. Altogether, we showed that Piezo1 was differentially controlled by lipid domains and the cytoskeleton and was favored by the stomatocyte-discocyte-echinocyte transformation.

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.

Quantifying Shear-induced Margination and Adhesion of Platelets in Microvascular Blood Flow

Platelet margination and adhesion are two critical and closely related steps in thrombus formation. Using dissipative particle dynamics (DPD) method that seamlessly models blood cells, blood plasma, and vessel walls with functionalized surfaces, we quantify the shear-induced margination and adhesion of platelets in microvascular blood flow. The results show that the occurrence of shear-induced RBC-platelet collisions has a remarkable influence on the degree of platelet margination. We characterize the lateral motion of individual platelets by a mean square displacement analysis of platelet trajectories, and find that the wall-induced lift force and the shear-induced displacement in wall-bounded flow cause the variation in near-wall platelet distribution. We then investigate the platelet adhesive dynamics under different flow conditions, by conducting DPD simulations of blood flow in a microtube with fibrinogen-coated wall surfaces. We find that the platelet adhesion is enhanced with the increase of fibrinogen concentration level but decreased with the increase of shear rate. These results are consistent with available experimental results. In addition, we demonstrate that the adherent platelets have a negative impact on the margination dynamics of the circulating platelets, which is mainly due to the climbing effect induced by the adherent ones. Taken together, these findings provide useful insights into the platelet margination and adhesion dynamics, which may facilitate the understanding of the predominant processes governing the initial stage of thrombus formation.

Effects of membrane viscoelasticity on the red blood cell dynamics in a microcapillary

The mechanical properties of red blood cells (RBCs) play key roles in their biological functions in microcirculation. In particular, RBCs must deform significantly to travel through microcapillaries with sizes comparable with or even smaller than their own. Although the dynamics of RBCs in microcapillaries have received considerable attention, the effect of membrane viscoelasticity has been largely overlooked. In this work, we present a computational study based on the boundary integral method and thin-shell mechanics to examine how membrane viscoelasticity influences the dynamics of RBCs flowing through straight and constricted microcapillaries. Our results reveal that the cell with a viscoelastic membrane undergoes substantially different motion and deformation compared with results based on a purely elastic membrane model. Comparisons with experimental data also suggest the importance of accounting for membrane viscoelasticity to properly capture the transient dynamics of an RBC flowing through a microcapillary. Taken together, these findings demonstrate the significant effects of membrane viscoelasticity on RBC dynamics in different microcapillary environments. The computational framework also lays the groundwork for more accurate quantitative modeling of the mechanical response of RBCs in their mechanotransduction process in subsequent investigations.

Circulating cell clusters aggravate the hemorheological abnormalities in COVID-19

Microthrombi and circulating cell clusters are common microscopic findings in patients with coronavirus disease 2019 (COVID-19) at different stages in the disease course, implying that they may function as the primary drivers in disease progression. Inspired by a recent flow imaging cytometry study of the blood samples from patients with COVID-19, we perform computational simulations to investigate the dynamics of different types of circulating cell clusters, namely white blood cell (WBC) clusters, platelet clusters, and red blood cell clusters, over a range of shear flows and quantify their impact on the viscosity of the blood. Our simulation results indicate that the increased level of fibrinogen in patients with COVID-19 can promote the formation of red blood cell clusters at relatively low shear rates, thereby elevating the blood viscosity, a mechanism that also leads to an increase in viscosity in other blood diseases, such as sickle cell disease and type 2 diabetes mellitus. We further discover that the presence of WBC clusters could also aggravate the abnormalities of local blood rheology. In particular, the extent of elevation of the local blood viscosity is enlarged as the size of the WBC clusters grows. On the other hand, the impact of platelet clusters on the local rheology is found to be negligible, which is likely due to the smaller size of the platelets. The difference in the impact of WBC and platelet clusters on local hemorheology provides a compelling explanation for the clinical finding that the number of WBC clusters is significantly correlated with thrombotic events in COVID-19 whereas platelet clusters are not. Overall, our study demonstrates that our computational models based on dissipative particle dynamics can serve as a powerful tool to conduct quantitative investigation of the mechanism causing the pathological alterations of hemorheology and explore their connections to the clinical manifestations in COVID-19.

Interaction of the Blood Components with Ascending Thoracic Aortic Aneurysm Wall: Biomechanical and Fluid Analyses

BACKGROUND

Ascending thoracic aortic aneurysm (ATAA) is an asymptomatic localized dilation of the aorta that is prone to rupture with a high rate of mortality. While diameter is the main risk factor for rupture assessment, it has been shown that the peak wall stress from finite element (FE) simulations may contribute to refinement of clinical decisions. In FE simulations, the intraluminal boundary condition is a single-phase blood flow that interacts with the thoracic aorta (TA). However, the blood is consisted of red blood cells (RBCs), white blood cells (WBCs), and plasma that interacts with the TA wall, so it may affect the resultant stresses and strains in the TA, as well as hemodynamics of the blood.

METHODS

In this study, discrete elements were distributed in the TA lumen to represent the blood components and mechanically coupled using fluid-structure interaction (FSI). Healthy and aneurysmal human TA tissues were subjected to axial and circumferential tensile loadings, and the hyperelastic mechanical properties were assigned to the TA and ATAA FE models.

RESULTS

The ATAA showed larger tensile and shear stresses but smaller fluid velocity compared to the ATA. The blood components experienced smaller shear stress in interaction with the ATAA wall compared to TA. The computational fluid dynamics showed smaller blood velocity and wall shear stress compared to the FSI.

CONCLUSIONS

This study is a first proof of concept, and future investigations will aim at validating the novel methodology to derive a more reliable ATAA rupture risk assessment considering the interaction of the blood components with the TA wall.

Identification of RNA Markers in Red Blood Cells for Doping Control in Autologous Blood Transfusion

The World Anti-Doping Agency (WADA) has prohibited the use of autologous blood transfusion (ABT) as a doping method by athletes. It is difficult to detect this doping method in laboratory tests, and a robust testing method has not yet been established. We conducted an animal experiment and used total RNA sequencing (RNA-Seq) to identify novel RNA markers to detect ABT doping within red blood cells (RBCs) as a pilot study before human trials. This study used whole blood samples from Wistar rats. The whole blood samples were mixed with a citrate–phosphate–dextrose solution with adenine (CPDA) and then stored in a refrigerator at 4 °C for 0 (control), 10, or 20 days. After each storage period, total RNA-Seq and bioinformatics were performed following RNA extraction and the purification of the RBCs. In the results, clear patterns of expression fluctuations were observed depending on the storage period, and it was found that there were large numbers of genes whose expression decreased in the 10- and 20-day periods compared to the control. Moreover, additional bioinformatic analysis identified three significant genes whose expression levels were drastically decreased according to the storage period. These results provide novel insights that may allow future studies to develop a testing method for ABT doping.

A new membrane formulation for modelling the flow of stomatocyte, discocyte, and echinocyte red blood cells

In this work, a numerical model that enables simulation of the deformation and flow behaviour of differently aged Red Blood Cells (RBCs) is developed. Such cells change shape and decrease in deformability as they age, thus impacting their ability to pass through the narrow capillaries in the body. While the body filters unviable cells from the blood naturally, cell aging poses key challenges for blood stored for transfusions. Therefore, understanding the influence RBC morphology and deformability have on their flow is vital. While several existing models represent young Discocyte RBC shapes well, a limited number of numerical models are developed to model aged RBC morphologies like Stomatocytes and Echinocytes. The existing models are also limited to shear and stretching simulations. Flow characteristics of these morphologies are yet to be investigated. This paper aims to develop a new membrane formulation for the numerical modelling of Stomatocyte, Discocytes and Echinocyte RBC morphologies to investigate their deformation and flow behaviour. The model used represents blood plasma using the Lattice Boltzmann Method (LBM) and the RBC membrane using the discrete element method (DEM). The membrane and the plasma are coupled by the Immersed Boundary Method (IBM). Previous LBM-IBM-DEM formulations represent RBC membrane response based on forces generated from changes in the local area, local length, local bending, and cell volume. In this new model, two new force terms are added: the local area difference force and the local curvature force, which are specially incorporated to model the flow and deformation behaviour of Stomatocytes and Echinocytes. To verify the developed model, the deformation behaviour of the three types of RBC morphologies are compared to well-characterised stretching and shear experiments. The flow modelling capabilities of the method are then demonstrated by modelling the flow of each cell through a narrow capillary. The developed model is found to be as accurate as benchmark Smoothed Particle Hydrodynamics (SPH) approaches while being significantly more computationally efficient.

Multiphysics and multiscale modeling of microthrombosis in COVID-19

Emerging clinical evidence suggests that thrombosis in the microvasculature of patients with Coronavirus disease 2019 (COVID-19) plays an essential role in dictating the disease progression. Because of the infectious nature of SARS-CoV-2, patients’ fresh blood samples are limited to access for in vitro experimental investigations. Herein, we employ a novel multiscale and multiphysics computational framework to perform predictive modeling of the pathological thrombus formation in the microvasculature using data from patients with COVID-19. This framework seamlessly integrates the key components in the process of blood clotting, including hemodynamics, transport of coagulation factors and coagulation kinetics, blood cell mechanics and adhesive dynamics, and thus allows us to quantify the contributions of many prothrombotic factors reported in the literature, such as stasis, the derangement in blood coagulation factor levels and activities, inflammatory responses of endothelial cells and leukocytes to the microthrombus formation in COVID-19. Our simulation results show that among the coagulation factors considered, antithrombin and factor V play more prominent roles in promoting thrombosis. Our simulations also suggest that recruitment of WBCs to the endothelial cells exacerbates thrombogenesis and contributes to the blockage of the blood flow. Additionally, we show that the recent identification of flowing blood cell clusters could be a result of detachment of WBCs from thrombogenic sites, which may serve as a nidus for new clot formation. These findings point to potential targets that should be further evaluated, and prioritized in the anti-thrombotic treatment of patients with COVID-19. Altogether, our computational framework provides a powerful tool for quantitative understanding of the mechanism of pathological thrombus formation and offers insights into new therapeutic approaches for treating COVID-19 associated thrombosis.

Emerging clinical evidence suggests that thrombosis in the microvasculature of patients with Coronavirus disease 2019 (COVID-19) plays an essential role in dictating the disease progression. We employ a novel multiphysics and multiscale computational framework to investigate the underlying mechanism of the pathological formation of microthrombi and circulating cell clusters in COVID-19. We quantify the contributions of many prothrombotic factors reported in the literature, such as stasis, the derangement in blood coagulation factor levels and activities, inflammatory responses of endothelial cells and leukocytes to the microthrombus formation in COVID-19, through which we identify the potential targets that should be further evaluated, and prioritized in the anti-thrombotic treatment of patients with COVID-19.

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 investigation of blood cell transport in retinal microaneurysms

Microaneurysms (MAs) are one of the earliest clinically visible signs of diabetic retinopathy (DR). MA leakage or rupture may precipitate local pathology in the surrounding neural retina that impacts visual function. Thrombosis in MAs may affect their turnover time, an indicator associated with visual and anatomic outcomes in the diabetic eyes. In this work, we perform computational modeling of blood flow in microchannels containing various MAs to investigate the pathologies of MAs in DR. The particle-based model employed in this study can explicitly represent red blood cells (RBCs) and platelets as well as their interaction in the blood flow, a process that is very difficult to observe in vivo. Our simulations illustrate that while the main blood flow from the parent vessels can perfuse the entire lumen of MAs with small body-to-neck ratio (BNR), it can only perfuse part of the lumen in MAs with large BNR, particularly at a low hematocrit level, leading to possible hypoxic conditions inside MAs. We also quantify the impacts of the size of MAs, blood flow velocity, hematocrit and RBC stiffness and adhesion on the likelihood of platelets entering MAs as well as their residence time inside, two factors that are thought to be associated with thrombus formation in MAs. Our results show that enlarged MA size, increased blood velocity and hematocrit in the parent vessel of MAs as well as the RBC-RBC adhesion promote the migration of platelets into MAs and also prolong their residence time, thereby increasing the propensity of thrombosis within MAs. Overall, our work suggests that computational simulations using particle-based models can help to understand the microvascular pathology pertaining to MAs in DR and provide insights to stimulate and steer new experimental and computational studies in this area.

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.

Modeling Red Blood Cell Viscosity Contrast Using Inner Soft Particle Suspension

The inner viscosity of a biological red blood cell is about five times larger than the viscosity of the blood plasma. In this work, we use dissipative particles to enable the proper viscosity contrast in a mesh-based red blood cell model. Each soft particle represents a coarse-grained virtual cluster of hemoglobin proteins contained in the cytosol of the red blood cell. The particle interactions are governed by conservative and dissipative forces. The conservative forces have purely repulsive character, whereas the dissipative forces depend on the relative velocity between the particles. We design two computational experiments that mimic the classical viscometers. With these experiments we study the effects of particle suspension parameters on the inner cell viscosity and provide parameter sets that result in the correct viscosity contrast. The results are validated with both static and dynamic biological experiment, showing an improvement in the accuracy of the original model without major increase in computational complexity.

Digital blood in massively parallel CPU/GPU systems for the study of platelet transport

We propose a highly versatile computational framework for the simulation of cellular blood flow focusing on extreme performance without compromising accuracy or complexity. The tool couples the lattice Boltzmann solver Palabos for the simulation of blood plasma, a novel finite-element method (FEM) solver for the resolution of deformable blood cells, and an immersed boundary method for the coupling of the two phases. The design of the tool supports hybrid CPU-GPU executions (fluid, fluid-solid interaction on CPUs, deformable bodies on GPUs), and is non-intrusive, as each of the three components can be replaced in a modular way. The FEM-based kernel for solid dynamics outperforms other FEM solvers and its performance is comparable to state-of-the-art mass-spring systems. We perform an exhaustive performance analysis on Piz Daint at the Swiss National Supercomputing Centre and provide case studies focused on platelet transport, implicitly validating the accuracy of our tool. The tests show that this versatile framework combines unprecedented accuracy with massive performance, rendering it suitable for upcoming exascale architectures.

An Introductory Overview of Image-Based Computational Modeling in Personalized Cardiovascular Medicine

Cardiovascular diseases account for the number one cause of deaths in the world. Part of the reason for such grim statistics is our limited understanding of the underlying mechanisms causing these devastating pathologies, which is made difficult by the invasiveness of the procedures associated with their diagnosis (e.g., inserting catheters into the coronal artery to measure blood flow to the heart). Likewise, it is also difficult to design and test assistive devices without implanting them in vivo. However, with the recent advancements made in biomedical scanning technologies and computer simulations, image-based modeling (IBM) has arisen as the next logical step in the evolution of non-invasive patient-specific cardiovascular medicine. Yet, due to its novelty, it is still relatively unknown outside of the niche field. Therefore, the goal of this manuscript is to review the current state-of-the-art and the limitations of the methods used in this area of research, as well as their applications to personalized cardiovascular investigations and treatments. Specifically, the modeling of three different physics – electrophysiology, biomechanics and hemodynamics – used in the cardiovascular IBM is discussed in the context of the physiology that each one of them describes and the mechanisms of the underlying cardiac diseases that they can provide insight into. Only the “bare-bones” of the modeling approaches are discussed in order to make this introductory material more accessible to an outside observer. Additionally, the imaging methods, the aspects of the unique cardiac anatomy derived from them, and their relation to the modeling algorithms are reviewed. Finally, conclusions are drawn about the future evolution of these methods and their potential toward revolutionizing the non-invasive diagnosis, virtual design of treatments/assistive devices, and increasing our understanding of these lethal cardiovascular diseases.

The ENUF method-Ewald summation based on nonuniform fast Fourier transform: Implementation, parallelization, and application

Computer simulations of model systems are widely used to explore striking phenomena in promising applications spanning from physics, chemistry, biology, to materials science and engineering. The long range electrostatic interactions between charged particles constitute a prominent factor in determining structures and states of model systems. How to efficiently calculate electrostatic interactions in simulation systems subjected to partial or full periodic boundary conditions has been a grand challenging task. In the past decades, a large variety of computational schemes has been proposed, among which the Ewald summation method is the most reliable route to accurately deal with electrostatic interactions between charged particles in simulation systems. In addition, extensive efforts have been done to improve computational efficiencies of the Ewald summation based methods. Representative examples are approaches based on cutoffs, reaction fields, multi-poles, multi-grids, and particle-mesh schemes. We sketched an ENUF method, an abbreviation for the Ewald summation method based on the nonuniform fast Fourier transform technique, and have implemented this method in particle-based simulation packages to calculate electrostatic energies and forces at micro- and mesoscopic levels. Extensive computational studies of conformational properties of polyelectrolytes, dendrimer-membrane complexes, and ionic fluids demonstrated that the ENUF method and its derivatives conserve both energy and momentum to floating point accuracy, and exhibit a computational complexity of O N log N with optimal physical parameters. These ENUF based methods are attractive alternatives in molecular simulations where high accuracy and efficiency of simulation methods are needed to accelerate calculations of electrostatic interactions at extended spatiotemporal scales.

Predictive modelling of thrombus formation in diabetic retinal microaneurysms

Microaneurysms (MAs) are one of the earliest clinically visible signs of diabetic retinopathy (DR). Vision can be reduced at any stage of DR by MAs, which may enlarge, rupture and leak fluid into the neural retina. Recent advances in ophthalmic imaging techniques enable reconstruction of the geometries of MAs and quantification of the corresponding haemodynamic metrics, such as shear rate and wall shear stress, but there is lack of computational models that can predict thrombus formation in individual MAs. In this study, we couple a particle model to a continuum model to simulate the platelet aggregation in MAs with different shapes. Our simulation results show that under a physiologically relevant blood flow rate, thrombosis is more pronounced in saccular-shaped MAs than fusiform-shaped MAs, in agreement with recent clinical findings. Our model predictions of the size and shape of the thrombi in MAs are consistent with experimental observations, suggesting that our model is capable of predicting the formation of thrombus for newly detected MAs. This is the first quantitative study of thrombosis in MAs through simulating platelet aggregation, and our results suggest that computational models can be used to predict initiation and development of intraluminal thrombus in MAs as well as provide insights into their role in the pathophysiology of DR.

Valid Presumption of Shiga Toxin-Mediated Damage of Developing Erythrocytes in EHEC-Associated Hemolytic Uremic Syndrome

The global emergence of clinical diseases caused by enterohemorrhagic Escherichia coli (EHEC) is an issue of great concern. EHEC release Shiga toxins (Stxs) as their key virulence factors, and investigations on the cell-damaging mechanisms toward target cells are inevitable for the development of novel mitigation strategies. Stx-mediated hemolytic uremic syndrome (HUS), characterized by the triad of microangiopathic hemolytic anemia, thrombocytopenia, and acute renal injury, is the most severe outcome of an EHEC infection. Hemolytic anemia during HUS is defined as the loss of erythrocytes by mechanical disruption when passing through narrowed microvessels. The formation of thrombi in the microvasculature is considered an indirect effect of Stx-mediated injury mainly of the renal microvascular endothelial cells, resulting in obstructions of vessels. In this review, we summarize and discuss recent data providing evidence that HUS-associated hemolytic anemia may arise not only from intravascular rupture of erythrocytes, but also from the extravascular impairment of erythropoiesis, the development of red blood cells in the bone marrow, via direct Stx-mediated damage of maturing erythrocytes, leading to “non-hemolytic” anemia.

Modelling of Red Blood Cell Morphological and Deformability Changes during In-Vitro Storage

Storage lesion is a critical issue facing transfusion treatments, and it adversely affects the quality and viability of stored red blood cells (RBCs). RBC deformability is a key indicator of cell health. Deformability measurements of each RBC unit are a key challenge in transfusion medicine research and clinical haematology. In this paper, a numerical study, inspired from the previous research for RBC deformability and morphology predictions, is conducted for the first time, to investigate the deformability and morphology characteristics of RBCs undergoing storage lesion. This study investigates the evolution of the cell shape factor, elongation index and membrane spicule details, where applicable, of discocyte, echinocyte I, echinocyte II, echinocyte III and sphero-echinocyte morphologies during 42 days of in-vitro storage at 4 °C in saline-adenine-glucose-mannitol (SAGM). Computer simulations were performed to investigate the influence of storage lesion-induced membrane structural defects on cell deformability and its recoverability during optical tweezers stretching deformations. The predicted morphology and deformability indicate decreasing quality and viability of stored RBCs undergoing storage lesion. The loss of membrane structural integrity due to the storage lesion further degrades the cell deformability and recoverability during mechanical deformations. This numerical approach provides a potential framework to study the RBC deformation characteristics under varying pathophysiological conditions for better diagnostics and treatments.

Deformation behaviour of stomatocyte, discocyte and echinocyte red blood cell morphologies during optical tweezers stretching

The red blood cell (RBC) deformability is a critical aspect, and assessing the cell deformation characteristics is essential for better diagnostics of healthy and deteriorating RBCs. There is a need to explore the connection between the cell deformation characteristics, cell morphology, disease states, storage lesion and cell shape-transformation conditions for better diagnostics and treatments. A numerical approach inspired from the previous research for RBC morphology predictions and for analysis of RBC deformations is proposed for the first time, to investigate the deformation characteristics of different RBC morphologies. The present study investigates the deformability characteristics of stomatocyte, discocyte and echinocyte morphologies during optical tweezers stretching and provides the opportunity to study the combined contribution of cytoskeletal spectrin network and the lipid-bilayer during RBC deformation. The proposed numerical approach predicts agreeable deformation characteristics of the healthy discocyte with the analogous experimental observations and is extended to further investigate the deformation characteristics of stomatocyte and echinocyte morphologies. In particular, the computer simulations are performed to investigate the influence of direct stretching forces on different equilibrium cell morphologies on cell spectrin link extensions and cell elongation index, along with a parametric analysis on membrane shear modulus, spectrin link extensibility, bending modulus and RBC membrane-bead contact diameter. The results agree with the experimentally observed stiffer nature of stomatocyte and echinocyte with respect to a healthy discocyte at experimentally determined membrane characteristics and suggest the preservation of relevant morphological characteristics, changes in spectrin link densities and the primary contribution of cytoskeletal spectrin network on deformation behaviour of stomatocyte, discocyte and echinocyte morphologies during optical tweezers stretching deformation. The numerical approach presented here forms the foundation for investigations into deformation characteristics and recoverability of RBCs undergoing storage lesion.

Cytoskeleton Remodeling Induces Membrane Stiffness and Stability Changes of Maturing Reticulocytes

Reticulocytes, the precursors of erythrocytes, undergo drastic alterations in cell size, shape, and deformability during maturation. Experimental evidence suggests that young reticulocytes are stiffer and less stable than their mature counterparts; however, the underlying mechanism is yet to be fully understood. Here, we develop a coarse-grained molecular-dynamics reticulocyte membrane model to elucidate how the membrane structure of reticulocytes contributes to their particular biomechanical properties and pathogenesis in blood diseases. First, we show that the extended cytoskeleton in the reticulocyte membrane is responsible for its increased shear modulus. Subsequently, we quantify the effect of weakened cytoskeleton on the stiffness and stability of reticulocytes, via which we demonstrate that the extended cytoskeleton along with reduced cytoskeleton connectivity leads to the seeming paradox that reticulocytes are stiffer and less stable than the mature erythrocytes. Our simulation results also suggest that membrane budding and the consequent vesiculation of reticulocytes can occur independently of the endocytosis-exocytosis pathway, and thus, it may serve as an additional means of removing unwanted membrane proteins from reticulocytes. Finally, we find that membrane budding is exacerbated when the cohesion between the lipid bilayer and the cytoskeleton is compromised, which is in accord with the clinical observations that erythrocytes start shedding membrane surface at the reticulocyte stage in hereditary spherocytosis. Taken together, our results quantify the stiffness and stability change of reticulocytes during their maturation and provide, to our knowledge, new insights into the pathogenesis of hereditary spherocytosis and malaria.

Modeling of the axon plasma membrane structure and its effects on protein diffusion

The axon plasma membrane consists of the membrane skeleton, which comprises ring-like actin filaments connected to each other by spectrin tetramers, and the lipid bilayer, which is tethered to the skeleton via, at least, ankyrin. Currently it is unknown whether this unique axon plasma membrane skeleton (APMS) sets the diffusion rules of lipids and proteins in the axon. To answer this question, we developed a coarse-grain molecular dynamics model for the axon that includes the APMS, the phospholipid bilayer, transmembrane proteins (TMPs), and integral monotopic proteins (IMPs) in both the inner and outer lipid layers. We first showed that actin rings limit the longitudinal diffusion of TMPs and the IMPs of the inner leaflet but not of the IMPs of the outer leaflet. To reconcile the experimental observations, which show restricted diffusion of IMPs of the outer leaflet, with our simulations, we conjectured the existence of actin-anchored proteins that form a fence which restricts the longitudinal diffusion of IMPs of the outer leaflet. We also showed that spectrin filaments could modify transverse diffusion of TMPs and IMPs of the inner leaflet, depending on the strength of the association between lipids and spectrin. For instance, in areas where spectrin binds to the lipid bilayer, spectrin filaments would restrict diffusion of proteins within the skeleton corrals. In contrast, in areas where spectrin and lipids are not associated, spectrin modifies the diffusion of TMPs and IMPs of the inner leaflet from normal to confined-hop diffusion. Overall, we showed that diffusion of axon plasma membrane proteins is deeply anisotropic, as longitudinal diffusion is of different type than transverse diffusion. Finally, we investigated how accumulation of TMPs affects diffusion of TMPs and IMPs of both the inner and outer leaflets by changing the density of TMPs. We showed that the APMS structure acts as a fence that restricts the diffusion of TMPs and IMPs of the inner leaflet within the membrane skeleton corrals. Our findings provide insight into how the axon skeleton acts as diffusion barrier and maintains neuronal polarity.

The axon plasma membrane skeleton consists of repeated periodic actin ring-like structures along its length connected via spectrin tetramers and anchored to the lipid bilayer at least via ankyrin. However, it is currently unclear whether this structure controls diffusion of lipids and proteins in the axon. Here, we developed a coarse-grain molecular dynamics computational model for the axon plasma membrane that comprises minimal representations for the APMS and the lipid bilayer. In a departure from current models, we found that actin rings limit diffusion of proteins only in the inner membrane leaflet. Then, we showed that actin anchored proteins likely act as “fences” confining diffusion of proteins in the outer leaflet. Our simulations, unexpectedly, also revealed that spectrin filaments could impede transverse diffusion in the inner leaflet of the axon and in some conditions modify diffusion from normal to abnormal. We predicted that diffusion of axon plasma membrane proteins is anisotropic as longitudinal diffusion is of different type than transverse (azimuthal) diffusion. We conclude that the periodic structure of the axon plays a critical role in controlling diffusion of proteins and lipids in the axon plasma membrane.

A coarse-grained red blood cell membrane model to study stomatocyte-discocyte-echinocyte morphologies

An improved red blood cell (RBC) membrane model is developed based on the bilayer coupling model (BCM) to accurately predict the complete sequence of stomatocyte-discocyte-echinocyte (SDE) transformation of a RBC. The coarse-grained (CG)-RBC membrane model is proposed to predict the minimum energy configuration of the RBC from the competition between lipid-bilayer bending resistance and cytoskeletal shear resistance under given reference constraints. In addition to the conventional membrane surface area, cell volume and bilayer-leaflet-area-difference constraints, a new constraint: total-membrane-curvature is proposed in the model to better predict RBC shapes in agreement with experimental observations. A quantitative evaluation of several cellular measurements including length, thickness and shape factor, is performed for the first time, between CG-RBC model predicted and three-dimensional (3D) confocal microscopy imaging generated RBC shapes at equivalent reference constraints. The validated CG-RBC membrane model is then employed to investigate the effect of reduced cell volume and elastic length scale on SDE transformation, to evaluate the RBC deformability during SDE transformation, and to identify the most probable RBC cytoskeletal reference state. The CG-RBC membrane model can predict the SDE shape behaviour under diverse shape-transforming scenarios, in-vitro RBC storage, microvascular circulation and flow through microfluidic devices.

Testosterone Therapy: An Assessment of the Clinical Consequences of Changes in Hematocrit and Blood Flow Characteristics

INTRODUCTION

Clinical guidelines indicate that hematocrit should be monitored during testosterone replacement therapy (TTh), with action taken if a level of 0.54 is exceeded.

AIM

To consider the extent of changes in hematocrit and putative effects on viscosity, blood flow, and mortality rates after TTh.

METHODS

We focused on literature describing benefits and possible pitfalls of TTh, including increased hematocrit. We used data from the BLAST RCT to determine change in hematocrit after 30 weeks of TTh and describe a clinical case showing the need for monitoring. We consider the validity of the current hematocrit cutoff value at which TTh may be modified. Ways in which hematocrit alters blood flow in the micro- and macro-vasculature are also considered.

MAIN OUTCOME MEASURES

The following measures were assessed: (i) change in hematocrit, (ii) corresponding actions taken in clinical practice, and (iii) possible blood flow changes following change in hematocrit.

RESULTS

Analysis of data from the BLAST RCT showed a significant increase in mean hematocrit of 0.01, the increase greater in men with lower baseline values. Although 0 of 61 men given TTh breached the suggested cutoff of 0.54 after 30 weeks, a clinical case demonstrates the need to monitor hematocrit. An association between hematocrit and morbidity and mortality appears likely but not proven and may be evident only in patient subgroups. The consequences of an increased hematocrit may be mediated by alterations in blood viscosity, oxygen delivery, and flow. Their relative impact may vary in different vascular beds.

CONCLUSIONS

TTh can effect an increased hematocrit via poorly understood mechanisms and may have harmful effects on blood flow that differ in patient subgroups. At present, there appears no scientific basis for using a hematocrit of 0.54 to modify TTh; other values may be more appropriate in particular patient groups. König CS, Balabani S, Hackett GI, et al. Testosterone Therapy: An Assessment of the Clinical Consequences of Changes in Hematocrit and Blood Flow Characteristics. Sex Med Rev 2019;7:650-660.

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.

Nano- and microscale mechanical properties of erythrocytes in hereditary spherocytosis

Hereditary spherocytosis (HS), an erythrocyte membranopathy, is a heterogeneous disease, even at the level of the erythrocyte population. The paper aims at studying the mechanical properties (the Young's modulus, median and RMS roughness of friction force maps; fractal dimension, lacunarity and spatial distribution parameters of lateral force maps) of the cell surface layer of the erythrocytes of two different morphologies (discocytes and spherocytes) in HS using atomic force microscopy. The results of spatial-spectral and fractal analysis showed that the mechanical property maps of the HS spherocyte surface were more structurally homogeneous compared to the maps of HS discocytes. HS spherocytes also had a reduced RMS roughness and lacunarity of the mechanical property maps. The Young's modulus and averaged friction forces over the microscale HS spherocyte surface regions were approximately 20% higher than that of HS discocytes. The revealed significant difference at the nano- and microscales in the structural and mechanical properties of main (discoidal and spheroidal) morphological types of HS erythrocytes can potentially cause blood flow disturbance in the vascular system in HS.

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.

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.

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.