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

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.

How fork-length asymmetry affects solvent connectivity and diffusion in grafted polymeric model membranes

The hydrophilic pore morphology and solvent diffusion within model (amphiphilic) polymer membranes are simulated by dissipative particle dynamics (DPD). The polymers are composed of a backbone of 18 covalently bonded A beads to which at regular intervals side chains are attached. The side chains are composed of linear Ap chains (i.e., -A1-A2…Ap) from which two branches, [AsC] and [ArC], split off (s ≤ r). C beads serve as functionalized hydrophilic pendent sites. The branch lengths (s + 1 and r + 1) are varied. Five repeat unit designs (with general formula A3[Ap[AsC][ArC]]) are considered: A2[A3C][A3C] (symmetric branching), A2[A2C][A4C], A2[AC][A5C], A2[C][A6C] (highly asymmetric branching), and A4[AC][A3C]. The distribution of water (W) and W diffusion through nanophase segregated hydrophilic pores is studied. For similar primary length p, an increase in side chain symmetry favors hydrophilic pore connectivity and long-range water transport. C beads located on the longer [ArC] branches reveal the highest C bead mobility and are more strongly associated with water than the C beads on the shorter [AsC] branches. The connectivity of hydrophilic (W and W + C) phases through mapped replica of selected snapshots obtained from Monte Carlo tracer diffusion simulations is in line with trends found from the W bead diffusivities during DPD simulations. The diffusive pathways for protons (H+) in proton exchange membranes and for hydronium (OH-) in anion exchange membranes are the same as for solvents. Therefore, control of the side chain architecture is an interesting design parameter for optimizing membrane conductivities.

Effect of vibration of the vortex mixer on the red blood cells

PURPOSES

Red blood cells (RBCs) are often subject to vibration during processing, transfusion, and transport. Further research is necessary to understand the effects of vibration on human RBCs and to reduce experimental deviations caused by device vibration.

METHODS

Flow cytometry was used in this study to observe the cytokine expression of IgG and IgA and deformation of human red blood cells affected by the vibration of a vortex mixer with varying frequency (750 rpm and 1500 rpm), duration (5 min and 10 min), and container volume (96 well plate and 48 well plate).

RESULTS

The size of RBCs in duration of 10 min is obviously smaller than the duration of 5 min. The 10-minute duration led to visibly smaller RBC sizes compared to the 5-minute duration. There was little effect on the size of RBCs in the 10-minute groups from differences in frequency and container volume. However, decreased RBC size can be observed in the 5-minute groups, where frequency is increased or container volume is decreased. Echinocytes were present in photomicrographs of all 10-minute groups, but microstructure of the RBCs was not impacted by vortex mixer vibration. The elevated frequency or reduced container volume results in an increased cytokine expression of IgG within the 5-minute groupings.

CONCLUSION

It can be inferred that vibration must not be overlooked due to its potential impact on the shape and cytokine expression of RBCs. Hence, the inclusion of vibration must be taken into consideration in experiments and devices pertaining to RBCs.

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.

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

Erythrophagocytosis occurring in the spleen is a critical process for removing senescent and diseased red blood cells (RBCs) from the microcirculation. Although some progress has been made in understanding how the biological signaling pathways mediate the phagocytic processes, the role of the biophysical interaction between RBCs and macrophages, particularly under pathological conditions such as sickle cell disease, has not been adequately studied. Here, we combine computational simulations with microfluidic experiments to quantify RBC-macrophage adhesion dynamics under flow conditions comparable to those in the red pulp of the spleen. We also investigate the RBC-macrophage interaction under normoxic and hypoxic conditions. First, we calibrate key model parameters in the adhesion model using microfluidic experiments for normal and sickle RBCs under normoxia and hypoxia. We then study the adhesion dynamics between the RBC and the macrophage. Our simulation illustrates three typical adhesion states, each characterized by a distinct dynamic motion of the RBCs, namely firm adhesion, flipping adhesion, and no adhesion (either due to no contact with macrophages or detachment from the macrophages). We also track the number of bonds formed when RBCs and macrophages are in contact, as well as the contact area between the two interacting cells, providing mechanistic explanations for the three adhesion states observed in the simulations and microfluidic experiments. Furthermore, we quantify, for the first time to our knowledge, the adhesive forces between RBCs (normal and sickle) and macrophages under different oxygenated conditions. Our results show that the adhesive forces between normal cells and macrophages under normoxia are in the range of 33-58 pN and 53-92 pN for sickle cells under normoxia and 155-170 pN for sickle cells under hypoxia. Taken together, our microfluidic and simulation results improve our understanding of the biophysical interaction between RBCs and macrophages in sickle cell disease and provide a solid foundation for investigating the filtration function of the splenic macrophages under physiological and pathological conditions.

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, there are relatively fewer works on investigating 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 a set of key factors that are 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 five-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.

Modeling Reactive Hyperemia to Better Understand and Assess Microvascular Function: A Review of Techniques

Reactive hyperemia is a well-established technique for the non-invasive evaluation of the peripheral microcirculatory function, measured as the magnitude of limb re-perfusion after a brief period of ischemia. Despite widespread adoption by researchers and clinicians alike, many uncertainties remain surrounding interpretation, compounded by patient-specific confounding factors (such as blood pressure or the metabolic rate of the ischemic limb). Mathematical modeling can accelerate our understanding of the physiology underlying the reactive hyperemia response and guide in the estimation of quantities which are difficult to measure experimentally. In this work, we aim to provide a comprehensive guide for mathematical modeling techniques that can be used for describing the key phenomena involved in the reactive hyperemia response, alongside their limitations and advantages. The reported methodologies can be used for investigating specific reactive hyperemia aspects alone, or can be combined into a computational framework to be used in (pre-)clinical settings.

Circulating cell clusters aggravate the hemorheological abnormalities in COVID-19

Microthrombi and circulating cell clusters (CCCs) are common microscopic findings in patients with 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 CCCs, namely white blood cell (WBC) clusters, platelet clusters and red blood cell (RBC) 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 RBC 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 do 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.

Self-Assembly of Lipid Mixtures in Solutions: Structures, Dynamics Processes and Mechanical Properties

The self-assembly of lipid mixtures in aqueous solution was investigated by dissipative particle dynamics simulation. Two types of lipid molecules were modelled, where three mixed structures, i.e., the membrane, perforated membrane and vesicle, were determined in the self-assembly processes. Phase behaviour was investigated by using the phase diagrams based on the tail chain lengths for the two types of lipids. Several parameters, such as chain number and average radius of gyration, were employed to explore the structural formations of the membrane and perforated membrane in the dynamic processes. Interface tension was used to demonstrate the mechanical properties of the membrane and perforated membrane in the equilibrium state and dynamics processes. Results help us to understand the self-assembly mechanism of the biomolecule mixtures, which has a potential application for designing the lipid molecule-based bio-membranes in solutions.

A particle-based computational model to analyse remodelling of the red blood cell cytoskeleton during malaria infections

Red blood cells can withstand the harsh mechanical conditions in the vasculature only because the bending rigidity of their plasma membrane is complemented by the shear elasticity of the underlying spectrin-actin network. During an infection by the malaria parasite Plasmodium falciparum, the parasite mines host actin from the junctional complexes and establishes a system of adhesive knobs, whose main structural component is the knob-associated histidine rich protein (KAHRP) secreted by the parasite. Here we aim at a mechanistic understanding of this dramatic transformation process. We have developed a particle-based computational model for the cytoskeleton of red blood cells and simulated it with Brownian dynamics to predict the mechanical changes resulting from actin mining and KAHRP-clustering. Our simulations include the three-dimensional conformations of the semi-flexible spectrin chains, the capping of the actin protofilaments and several established binding sites for KAHRP. For the healthy red blood cell, we find that incorporation of actin protofilaments leads to two regimes in the shear response. Actin mining decreases the shear modulus, but knob formation increases it. We show that dynamical changes in KAHRP binding affinities can explain the experimentally observed relocalization of KAHRP from ankyrin to actin complexes and demonstrate good qualitative agreement with experiments by measuring pair cross-correlations both in the computer simulations and in super-resolution imaging experiments.

Malaria is one of the deadliest infectious diseases and its symptoms are related to the blood stage, when the parasite multiplies within red blood cells. In order to avoid clearance by the spleen, the parasite produces specific factors like the adhesion receptor PfEMP1 and the multifunctional protein KAHRP that lead to the formation of adhesive knobs on the surface of the red blood cells and thus increase residence time in the vasculature. We have developed a computational model for the parasite-induced remodelling of the actin-spectrin network to quantitatively predict the dynamical changes in the mechanical properties of the infected red blood cells and the spatial distribution of the different protein components of the membrane skeleton. Our simulations show that KAHRP can relocate to actin junctions due to dynamical changes in binding affinities, in good qualitative agreement with super-resolution imaging experiments. In the future, our simulation framework can be used to gain further mechanistic insight into the way malaria parasites attack the red blood cell cytoskeleton.

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.

How the spleen reshapes and retains young and old red blood cells: A computational investigation

The spleen, the largest secondary lymphoid organ in humans, not only fulfils a broad range of immune functions, but also plays an important role in red blood cell’s (RBC) life cycle. Although much progress has been made to elucidate the critical biological processes involved in the maturation of young RBCs (reticulocytes) as well as removal of senescent RBCs in the spleen, the underlying mechanisms driving these processes are still obscure. Herein, we perform a computational study to simulate the passage of RBCs through interendothelial slits (IES) in the spleen at different stages of their lifespan and investigate the role of the spleen in facilitating the maturation of reticulocytes and in clearing the senescent RBCs. Our simulations reveal that at the beginning of the RBC life cycle, intracellular non-deformable particles in reticulocytes can be biomechanically expelled from the cell upon passage through IES, an insightful explanation of why this peculiar “pitting” process is spleen-specific. Our results also show that immature RBCs shed surface area by releasing vesicles after crossing IES and progressively acquire the biconcave shape of mature RBCs. These findings likely explain why RBCs from splenectomized patients are significantly larger than those from nonsplenectomized subjects. Finally, we show that at the end of their life span, senescent RBCs are not only retained by IES due to reduced deformability but also become susceptible to mechanical lysis under shear stress. This finding supports the recent hypothesis that transformation into a hemolyzed ghost is a prerequisite for phagocytosis of senescent RBCs. Altogether, our computational investigation illustrates critical biological processes in the spleen that cannot be observed in vivo or in vitro and offer insights into the role of the spleen in the RBC physiology.

The spleen, the largest secondary lymphoid organ in humans, not only fulfils a broad range of immune functions, but also plays an important role in red blood cell (RBC) life cycle. In this study, we perform a computational study to simulate the passage of RBCs through interendothelial slits (IES) in the spleen at different stages of their lifespan, a critical biological process that cannot be observed in humans. Our simulation results illustrate a specific role of spleen in shaping young RBCs, which points to a probable missing step in current in vitro RBC culture protocols that fail to generate a majority of typical biconcave RBCs. Our results also reveal that intra-splenic mechanical constraints likely contribute to the final clearance and elimination of aged RBCs. Altogether, we demonstrate that our computational model can provide mechanistic rationales for experimental studies, offer insights into the role of the spleen in the RBC physiology and help the optimization of in vitro RBC culture techniques.

Organic and inorganic equivalent models for analysis of red blood cell mechanical behaviour

Experimental investigation into the mechanical response of red blood cells is presently impeded with the main impediments being the micro dimensions involved and ethical issues associated with in vivo testing. The widely employed alternative approach of computational modelling suffers from its own inherent limitations being reliant on precise constitutive and boundary information. Moreover, and somewhat critically, numerical computational models themselves are required to be validated by means of experimentation and hence suffer similar impediments. An alternative experimental approach is examined in this paper involving large-scale equivalent models manufactured principally from inorganic, and to lesser extent organic, materials. Although there presently exists no known method providing the means to investigate the mechanical response of red blood cells using scaled models simultaneously having different dimensions and materials, the present paper aims to develop a scaled framework based on the new finite-similitude theory that has appeared in the recent open literature. Computational models are employed to test the effectiveness of the proposed method, which in principle can provide experimental solution methods to a wide range of practical applications including the design of red-blood cell nanorobots and drug delivery systems. By means of experimentally validated numerical experiments under impact loading it is revealed that although exact prediction is not achieved good accuracy can nevertheless be obtained. Furthermore, it is demonstrated how the proposed approach for first time provides a means to relate models at different scales founded on different constitutive equations.

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.

Definitions of local density in density-dependent potentials for mixtures

Density-dependent potentials are frequently used in materials simulations because of their approximate description of many-body effects at minimal computational cost. However, in order to apply such models to multicomponent systems, an appropriate definition of total local particle density is required. Here, we discuss two definitions of local density in the context of many-body dissipative particle dynamics. We show that only a potential which combines local densities from all particle types in its argument gives physically meaningful results for all composition ratios. Drawing on the ideas from metal potentials, we redefine local density such that it can accommodate different intertype interactions despite the constraint to keep the main interaction parameter constant, known as Warren's no-go theorem, and generalize the many-body potential to heterogeneous systems. We then show via simulation how liquid-liquid and liquid-solid coexistence can arise just by tuning the interaction parameters.

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.

Identifying the start of a platelet aggregate by the shear rate and the cell-depleted layer

Computer simulations were performed to study the transport of red blood cells and platelets in high shear flows, mimicking earlier published in vitro experiments in microfluidic devices with high affinity for platelet aggregate formation. The goal is to understand and predict where thrombus formation starts. Additionally, the need of cell-based modelling in these microfluidic devices is demonstrated by comparing our results with macroscopic models, wherein blood is modelled as a continuous fluid. Hemocell, a cell-based blood flow simulation framework is used to investigate the transport physics in the microfluidic devices. The simulations show an enlarged cell-depleted layer at the site where a platelet aggregate forms in the experiments. In this enlarged cell-depleted layer, the probability to find a platelet is higher than in the rest of the microfluidic device. In addition, the shear rates are sufficiently high to allow for the von Willebrand factor to elongate in this region. We hypothesize that the enlarged cell-depleted layer combined with a sufficiently large platelet flux and sufficiently high shear rates result in an haemodynamic environment that is a preferred location for initial platelet aggregation.

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.

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.

Penetration of nanoparticles across a lipid bilayer: effects of particle stiffness and surface hydrophobicity

The cellular uptake of nanoparticles (NPs) has drawn significant attention due to their great importance and potential in drug delivery, bioimaging, and specific targeting. Here, we conduct a computational study on the translocation process of soft nanoparticles with different elasticities and surface hydrophobicities through a lipid bilayer membrane. It is shown that the translocation abilities of hydrophilic NPs can be enhanced by increasing their stiffness, while the penetrability of hydrophobic NPs is weakened by increasing the particle stiffness. The free energy analysis indicates that rigid hydrophilic NPs and soft hydrophobic NPs encounter lower energy barriers during penetration. In direct translocation, different deformation modes are observed for NPs with different surface hydrophobicities during cellular internalization. Further, deformation analysis demonstrates that hydrophilic NPs are flattened in the membrane plane, while hydrophobic NPs are elongated along the membrane norm during penetration. We conclude that the elasticity of NPs has an obvious impact on their ability to penetrate across the lipid bilayer membrane through different morphological responses of hydrophilic and hydrophobic NPs. These results shed light on the coupled effects of particle elasticity and surface hydrophobicity on the cellular uptake of elastic NPs, which may provide useful guidelines for designing effective nanocarrier systems for drug delivery.

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.

Adaptive resolution simulations coupling atomistic water to dissipative particle dynamics

Multiscale methods are the most efficient way to address the interlinked spatiotemporal scales encountered in soft matter and molecular liquids. In the literature reported hybrid approaches span from quantum to atomistic, coarse-grained, and continuum length scales. In this article, we present the hybrid coupling of the molecular dynamics (MD) and dissipative particle dynamics (DPD) methods, bridging the micro- and mesoscopic descriptions. The interfacing is performed within the adaptive resolution scheme (AdResS), which is a linear momentum conserving coupling technique. Our methodology is hence suitable to simulate fluids on the micro/mesoscopic scale, where hydrodynamics plays an important role. The presented approach is showcased for water at ambient conditions. The supramolecular coupling is enabled by a recently developed clustering algorithm SWINGER that assembles, disassembles, and reassembles clusters as needed during the course of the simulation. This allows for a seamless coupling between standard atomistic MD and DPD models. The developed framework can be readily applied to various applications in the fields of materials and life sciences, e.g., simulations of phospholipids and polymer melts, or to study the red blood cells behavior in normal and disease states.

Computational Biomechanics of Human Red Blood Cells in Hematological Disorders

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

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

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

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

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

Cellular Level In-silico Modeling of Blood Rheology with An Improved Material Model for Red Blood Cells

Many of the intriguing properties of blood originate from its cellular nature. Therefore, accurate modeling of blood flow related phenomena requires a description of the dynamics at the level of individual cells. This, however, presents several computational challenges that can only be addressed by high performance computing. We present Hemocell, a parallel computing framework which implements validated mechanical models for red blood cells and is capable of reproducing the emergent transport characteristics of such a complex cellular system. It is computationally capable of handling large domain sizes, thus it is able to bridge the cell-based micro-scale and macroscopic domains. We introduce a new material model for resolving the mechanical responses of red blood cell membranes under various flow conditions and compare it with a well established model. Our new constitutive model has similar accuracy under relaxed flow conditions, however, it performs better for shear rates over 1,500 s^-1. We also introduce a new method to generate randomized initial conditions for dense mixtures of different cell types free of initial positioning artifacts.

Biomechanics and biorheology of red blood cells in sickle cell anemia

Sickle cell anemia (SCA) is an inherited blood disorder that causes painful crises due to vaso-occlusion of small blood vessels. The primary cause of the clinical phenotype of SCA is the intracellular polymerization of sickle hemoglobin resulting in sickling of red blood cells (RBCs) in deoxygenated conditions. In this review, we discuss the biomechanical and biorheological characteristics of sickle RBCs and sickle blood as well as their implications toward a better understanding of the pathophysiology and pathogenesis of SCA. Additionally, we highlight the adhesive heterogeneity of RBCs in SCA and their specific contribution to vaso-occlusive crisis.