Simulation of malaria-infected red blood cells in microfluidic channels: Passage and blockage

Lattice Boltzmann simulation of deformable fluid-filled bodies: progress and perspectives

With the rapid development of studies involving droplet microfluidics, drug delivery, cell detection, and microparticle synthesis, among others, many scientists have invested significant efforts to model the flow of these fluid-filled bodies. Motivated by the intricate coupling between hydrodynamics and the interactions of fluid-filled bodies, several methods have been developed. The objective of this review is to present a compact foundation of the methods used in the literature in the context of lattice Boltzmann methods. For hydrodynamics, we focus on the lattice Boltzmann method due to its specific ability to treat time- and spatial-dependent boundary conditions and to incorporate new physical models in a computationally efficient way. We split the existing methods into two groups with regard to the interfacial boundary: fluid-structure and fluid-fluid methods. The fluid-structure methods are characterised by the coupling between fluid dynamics and mechanics of the flowing body, often used in applications involving membranes and similar flexible solid boundaries. We further divide fluid-structure-based methods into two subcategories, those which treat the fluid-structure boundary as a continuum medium and those that treat it as a discrete collection of individual springs and particles. Next, we discuss the fluid-fluid methods, particularly useful for the simulations of fluid-fluid interfaces. We focus on models for immiscible droplets and their interaction in a suspending fluid and describe benchmark tests to validate the models for fluid-filled bodies.

Hydrodynamics of a multicomponent vesicle under strong confinement

We numerically investigate the hydrodynamics and membrane dynamics of a multicomponent vesicle in two strongly confined geometries. This serves as a simplified model for red blood cells undergoing large deformations while traversing narrow constrictions. We propose a new parameterization for the bending modulus that remains positive for all lipid phase parameter values. For a multicomponent vesicle passing through a stenosis, we establish connections between various properties: lipid phase coarsening, size and flow profile of the lubrication layers, excess pressure, and the tank-treading velocity of the membrane. For a multicomponent vesicle passing through a contracting channel, we find that the lipid always phase separates so that the vesicle is stiffer in the front as it passes through the constriction. For both cases of confinement we find that lipid coarsening is arrested under strong confinement, and resumes at a high rate upon relief from extreme confinement. The results may be useful for efficient sorting lipid domains using microfluidic flows by controlled release of vesicles passing through strong confinement.

Numerical investigation on red blood cell flow based on unstructured grid

Prediction of blood cell flow is known as the difficult research by reason of the complexity of blood vessel. In this study, considering the complex structure of blood vessels, a mechanical model for red blood cell (RBC) based on unstructured grid has been established to study the flow characteristics of RBCs in complex blood vessels. In the model, the strain-energy function by Skalak is employed to model the shear elasticity and surface-area conservation of the membrane, and the hinge spring is used to describe the forces originating from local bending of the membrane. The immersed boundary method is utilized to couple the interphase force. Using the model, the stretching test of RBC is compared with the experiment data, and the good agreement verified the validation of the present model. The morphology of red blood cell and the blood viscosity in micro-vessel are studied. RBCs move with a symmetric shape (parachute shape) in small blood vessels, and the buckling instability is observed when the RBC flow slowly through a micro-vessel or a converging-diverging capillary. When the vessel diameter is around 10 μm, the reverse Fahraeus-Lindqvist effect is presented. The blood apparent viscosity shows linear increase with the blood hematocrit. In addition, Malaria infection can make the RBC deformability decreased and the blood apparent viscosity increased.

Characterization of fibronectin properties by integrated micro-fluidic experiments and fluid-structure interaction simulations

Fibronectin (Fn) has been observed to assemble in the extracellular matrix (ECM) of cell culture and stretch in response to the external force. The alteration of molecule domain functions generally follows the extension of Fn. Several researchers have investigated fibronectin extensively in molecular architecture and conformation structure. However, the bulk material behavior of the Fn in the ECM has not been fully depicted at the cell scale, and many studies have ignored physiological conditions. Conversely, microfluidic techniques that explore cellular properties based on cell deformation and adhesion have emerged as a powerful and effective platform to study cell rheological transformation in a physiological environment. However, directly quantifying properties from microfluidic measurements remains a challenge. Therefore, it is an efficient way to combine experimental measurements with a robust and reliable numerical framework to calibrate the mechanical stress distribution in the test sample. In this paper, we present a monolithic Lagrangian fluid-structure interaction (FSI) approach within the Optimal Transportation Meshfree (OTM) framework that enables the investigation of the adherent Red Blood Cell (RBC) interacting with fluid and overcomes the drawbacks of the traditional computational tools such as the mesh entanglement and interface tracking, etc. This study aims to assess the material properties of the RBC and Fn fiber by calibrating the numerical predictions to experimental measurements. Moreover, a physical-based constitutive model will be proposed to describe the bulk behavior of the Fn fiber inflow, and the rate-dependent deformation and separation of the Fn fiber will be discussed.

Persistent red blood cells retain their ability to move in microcapillaries under high levels of oxidative stress

Oxidative stress is one of the key factors that leads to red blood cells (RBCs) aging, and impairs their biomechanics and oxygen delivery. It occurs during numerous pathological processes and causes anaemia, one of the most frequent side effects of cancer chemotherapy. Here, we used microfluidics to simulate the microcirculation of RBCs under oxidative stress induced by tert-Butyl hydroperoxide. Oxidative stress was expected to make RBCs more rigid, which would lead to decrease their transit velocity in microfluidic channels. However, single-cell tracking combined with cytological and AFM studies reveals cell heterogeneity, which increases with the level of oxidative stress. The data indicates that the built-in antioxidant defence system has a limit exceeding which haemoglobin oxidation, membrane, and cytoskeleton transformation occurs. It leads to cell swelling, increased stiffness and adhesion, resulting in a decrease in the transit velocity in microcapillaries. However, even at high levels of oxidative stress, there are persistent cells in the population with an undisturbed biophysical phenotype that retain the ability to move in microcapillaries. Developed microfluidic analysis can be used to determine RBCs' antioxidant capacity for the minimization of anaemia during cancer chemotherapy.

Multi-scale model of clogging in microfluidic devices with grid-like geometries

We propose a coarse-grained theoretical model to capture the ageing of microfluidic devices under different conditions including constant applied flow rate and constant applied pressure gradient. Microfluidic devices that sort cells by their deformability hold significant promise for medical applications. However, clogging in these microfluidic systems causes their properties to change over time and potentially limits their reliability. We compare the results of the coarse-grained model with those of stochastic simulations and with existing theoretical studies. Lastly, we apply the model to experimental data on the clogging of sickle red blood cells and discuss its wider applicability.

Structural Configuration of Blood Cell Membranes Determines Their Nonlinear Deformation Properties

The ability of neutrophils and red blood cells (RBCs) to undergo significant deformations is a key to their normal functioning. Disruptions of these processes can lead to pathologies. This work studied the influence of structural configuration rearrangements of membranes after exposure to external factors on the ability of native membranes of neutrophils and RBCs to undergo deep deformation. The rearrangement of the structural configuration of neutrophil and RBC membranes under the influence of cytological fixatives caused nonlinear deformation phenomena. There were an increase in Young's modulus, a decrease in the depth of homogeneous bending, and a change in the distance between cytoskeletal junctions. Based on the results of the analysis of experimental data, a mathematical model was proposed that describes the process of deep bending of RBСs and neutrophil membranes.

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.

Numerical Simulations of Red-Blood Cells in Fluid Flow: A Discrete Multiphysics Study

In this paper, we present a methodological study of modelling red blood cells (RBCs) in shear-induced flows based on the discrete multiphysics (DMP) approach. The DMP is an alternative approach from traditional multiphysics based on meshless particle-based methods. The proposed technique has been successful in modelling multiphysics and multi-phase problems with large interfacial deformations such as those in biological systems. In this study, we present the proposed method and introduce an accurate geometrical representation of the RBC. The results were validated against available data in the literature. We further illustrate that the proposed method is capable of modelling the rupture of the RBC membrane with minimum computational difficulty.

Quantitative absorption imaging of red blood cells to determine physical and mechanical properties

Red blood cells or erythrocytes, constituting 40 to 45 percent of the total volume of human blood are vesicles filled with hemoglobin with a fluid-like lipid bilayer membrane connected to a 2D spectrin network. The shape, volume, hemoglobin mass, and membrane stiffness of RBCs are important characteristics that influence their ability to circulate through the body and transport oxygen to tissues. In this study, we show that a simple two-LED set up in conjunction with standard microscope imaging can accurately determine the physical and mechanical properties of single RBCs. The Beer–Lambert law and undulatory motion dynamics of the membrane have been used to measure the total volume, hemoglobin mass, membrane tension coefficient, and bending modulus of RBCs. We also show that this method is sensitive enough to distinguish between the mechanical properties of RBCs during morphological changes from a typical discocyte to echinocytes and spherocytes. Measured values of the tension coefficient and bending modulus are 1.27 × 10⁻⁶ J m⁻² and 7.09 × 10⁻²⁰ J for discocytes, 4.80 × 10⁻⁶ J m⁻² and 7.70 × 10⁻²⁰ J for echinocytes, and 9.85 × 10⁻⁶ J m⁻² and 9.69 × 10⁻²⁰ J for spherocytes, respectively. This quantitative light absorption imaging reduces the complexity related to the quantitative imaging of the biophysical and mechanical properties of a single RBC that may lead to enhanced yet simplified point of care devices for analyzing blood cells.

The constant thickness in the microfluidic channel is used for controlled absorption of red and blue light to measure red blood cell hemoglobin and height mapping. High speed recording of the height mapping provides us the membrane fluctuation.

A system for the high-throughput measurement of the shear modulus distribution of human red blood cells

Reduced deformability of red blood cells (RBCs) can affect the hemodynamics of the microcirculation and reduce oxygen transport efficiency. It is also well known that reduced RBC deformability is a signature of various physical disorders, including sepsis, and that the primary determinant of RBC deformability is the membrane shear modulus. To measure the distribution of an individual's RBC shear modulus with high throughput, we a) developed a high-fidelity computational model of RBCs in confined microchannels to inform design decisions; b) created a novel experimental system combining microfluidic flow, imaging, and image analysis; and c) performed automated comparisons between measured quantities and numerical predictions to extract quantitative measures of the RBC shear modulus for each of thousands of cells. We applied our computational simulation platform to construct the appropriate deformability figure(s) of merit to quantify RBC stiffness based on an experimentally measured, steady-state cell shape in flow through a microchannel. In particular, we determined a shape parameter based on the second moment of the cell shape that is sensitive to the changes in the membrane stiffness and cell size. We then conducted microfluidic experiments and developed custom automated image processing codes to identify and track the position and shape of individual RBCs within micro-constrictions. The fabricated microchannels include a square cross-section imaging region (7 by 7 μm) and we applied order 10 kPa pressure differences to induce order 10 mm s^-1 cell velocities. The combination of modeling, microfluidics, and imaging enables, for the first time, quantitative measurement of the shear moduli of thousands of RBCs in human blood samples. We demonstrate the high-throughput features by sensitive quantification of the changes in the distribution of RBC stiffness with aging. This combined measurement and computational platform is ultimately intended to diagnose blood cell disorders in patients.

Microfluidic deformability-activated sorting of single particles

Mechanical properties have emerged as a significant label-free marker for characterizing deformable particles such as cells. Here, we demonstrated the first single-particle-resolved, cytometry-like deformability-activated sorting in the continuous flow on a microfluidic chip. Compared with existing deformability-based sorting techniques, the microfluidic device presented in this work measures the deformability and immediately sorts the particles one-by-one in real time. It integrates the transit-time-based deformability measurement and active hydrodynamic sorting onto a single chip. We identified the critical factors that affect the sorting dynamics by modeling and experimental approaches. We found that the device throughput is determined by the summation of the sensing, buffering, and sorting time. A total time of ~100 ms is used for analyzing and sorting a single particle, leading to a throughput of 600 particles/min. We synthesized poly(ethylene glycol) diacrylate (PEGDA) hydrogel beads as the deformability model for device validation and performance evaluation. A deformability-activated sorting purity of 88% and an average efficiency of 73% were achieved. We anticipate that the ability to actively measure and sort individual particles one-by-one in a continuous flow would find applications in cell-mechanotyping studies such as correlational studies of the cell mechanical phenotype and molecular mechanism.

Methods to Investigate the Deformability of RBC During Malaria

Despite a 30% decline in mortality since 2000, malaria still affected 219 million subjects and caused 435,000 deaths in 2017. Red blood cells (RBC) host Plasmodium parasites that cause malaria, of which Plasmodium falciparum is the most pathogenic. The deformability of RBC is markedly modified by invasion and development of P. falciparum. Surface membrane area is potentially impacted by parasite entry and development, the cytoskeleton is modified by parasite proteins and cytosol viscosity is altered by parasite metabolism. RBC hosting mature parasites (second half of the asexual erythrocytic cycle) are abnormally stiff but the main reason for their absence from the circulation is their adherence to endothelial cells, mediated by parasite proteins exposed at the infected-RBC surface. By contrast, the circulation of non-adherent rings and gametocytes, depends predominantly on deformability. Altered deformability of rings and of uninfected-RBC altered by malaria infection is an important determinant of malaria pathogenesis. It also impacts the response to antimalarial therapy. Unlike conventional antimalarials that target mature stages, currently recommended first-line artemisinin derivatives and the emerging spiroindolones act on circulating rings. Methods to investigate the deformability of RBC are therefore critical to understand the clearance of infected- and uninfected-RBC in malaria. Herein, we review the main methods to assess the deformability of P. falciparum infected-RBC, and their contribution to the understanding of how P. falciparum infection causes disease, how the parasite is transmitted and how antimalarial drugs induce parasite clearance.

Numerical simulation of optical control for a soft particle in a microchannel

Technologies that use optical force to actively control particles in microchannels are a significant area of research interest in various fields. An optical force is generated by the momentum change caused by the refraction and reflection of light, which changes the particle surface as a function of the angle of incidence of light and which in turn feeds back and modifies the force on the particle. Simulating this phenomenon is a complex task. The deformation of a particle, the interaction between the surrounding fluid and the particle, and the reflection and refraction of light should be analyzed simultaneously. Herein, a deformable particle in a microchannel subjected to optical interactions is simulated using the three-dimensional lattice Boltzmann immersed-boundary method. The laser from the optical source is analyzed by dividing it into individual rays. To calculate the optical forces exerted on the particle, the intensity, momentum, and ray direction are calculated. The optical-separator problem with one optical source is analyzed by measuring the distance traveled because of the optical force. The optical-stretcher problem with two optical sources is then studied by analyzing the relation between the intensity of the optical source and particle deformation. This simulation will help the design of sorting and measuring by optical force.

Shear stress in the microvasculature: influence of red blood cell morphology and endothelial wall undulation

The effect of red blood cells and the undulation of the endothelium on the shear stress in the microvasculature is studied numerically using the lattice Boltzmann-immersed boundary method. The results demonstrate a significant effect of both the undulation of the endothelium and red blood cells on wall shear stress. Our results also reveal that morphological alterations of red blood cells, as occur in certain pathologies, can significantly affect the values of wall shear stress. The resulting fluctuations in wall shear stress greatly exceed the nominal values, emphasizing the importance of the particulate nature of blood as well as a more realistic description of vessel wall geometry in the study of hemodynamic forces. We find that within the channel widths investigated, which correspond to those found in the microvasculature, the inverse minimum distance normalized to the channel width between the red blood cell and the wall is predictive of the maximum wall shear stress observed in straight channels with a flowing red blood cell. We find that the maximum wall shear stress varies several factors more over a range of capillary numbers (dimensionless number relating strength of flow to membrane elasticity) and reduced areas (measure of deflation of the red blood cell) than the minimum wall shear stress. We see that waviness reduces variation in minimum and maximum shear stresses among different capillary and reduced areas.

A monolithic fluid-structure interaction framework applied to red blood cells

A parallel fully coupled (monolithic) fluid-structure interaction (FSI) algorithm has been applied to the deformation of red blood cells (RBCs) in capillaries, where cell deformability has significant effects on blood rheology. In the present FSI algorithm, fluid domain is discretized using the side-centered unstructured finite volume method based on the Arbitrary Lagrangian-Eulerian (ALE) formulation; meanwhile, solid domain is discretized with the classical Galerkin finite element formulation for the Saint Venant-Kirchhoff material in a Lagrangian frame. In addition, the compatible kinematic boundary condition is enforced at the fluid-solid interface in order to conserve the mass of cytoplasmic fluid within the red cell at machine precision. In order to solve the resulting large-scale algebraic linear systems in a fully coupled manner, a new matrix factorization is introduced similar to that of the projection method, and the parallel algebraic multigrid solver BoomerAMG is used for the scaled discrete Laplacian provided by the HYPRE library, which we access through the PETSc library. Three important physical parameters for the blood flow are simulated and analyzed: (1) the effect of capillary diameter, (2) the effect of red cell membrane thickness, and (3) the effect of red cell spacing (hematocrit). The numerical calculations initially indicate a shape deformation in which biconcave discoid shape changes to a parachute-like shape. Furthermore, the parachute-like cell shape in small capillaries undergoes a cupcake-shaped buckling instability, which has not been observed in the literature. The instability forms thin riblike features, and the red cell deformation is not axisymmetric but three-dimensional.

Simultaneous measurement of blood pressure and RBC aggregation by monitoring on–off blood flows supplied from a disposable air-compressed pump

A simple method for simultaneously measuring RBC aggregation and blood pressure is demonstrated by analyzing blood flows supplied from a disposable air-compressed pump.

Modeling Cell Adhesion and Extravasation in Microvascular System

The blood flow behaviors in the microvessels determine the transport modes and further affect the metastasis of circulating tumor cells (CTCs). Much biochemical and biological efforts have been made on CTC metastasis; however, precise experimental measurement and accurate theoretical prediction on its mechanical mechanism are limited. To complement these, numerical modeling of a CTC extravasation from the blood circulation, including the steps of adhesion and transmigration, is discussed in this chapter. The results demonstrate that CTCs prefer to adhere at the positive curvature of curved microvessels, which is attributed to the positive wall shear stress/gradient. Then, the effects of particulate nature of blood on CTC adhesion are investigated and are found to be significant in the microvessels. Furthermore, the presence of red blood cell (RBC) aggregates is also found to promote the CTC adhesion by providing an additional wall-directed force. Finally, a single cell passing through a narrow slit, mimicking CTC transmigration, was examined under the effects of cell deformability. It showed that the cell shape and surface area increase play a more important role than the cell elasticity in cell transit across the narrow slit.

The dynamics of a healthy and infected red blood cell in flow through constricted channels: A DPD simulation

Understanding the dynamics of red blood cell (RBC) motion under in silico conditions is central to the development of cost-effective diagnostic tools. Specifically, unraveling the relationship between the rheological properties and the nature of shape change in the RBC (healthy or infected) can be extremely useful. In case of malarial infection, RBC progressively loses its deformability and tends to occlude the microvessel. In the present study, detailed mesoscopic simulations are performed to investigate the deformation dynamics of an RBC in flow through a constricted channel. Specifically, the manifestation of viscous forces (through flow rates) on the passage and blockage characteristics of a healthy red blood cell (hRBC) vis-á-vis an infected red blood cell (iRBC) are investigated. A finite-sized dissipative particle dynamics framework is used to model plasma in conjunction with a discrete model for the RBC. Instantaneous wall boundary method was used to model no-slip wall boundary conditions with a good control on the near-wall density fluctuations and compressibility effects. To investigate the microvascular occlusion, the RBC motion through 2 types of constricted channels, viz, (1) a tapered microchannel and (2) a stenosed-type microchannel, were simulated. It was observed that the deformation of an infected cell was much less compared with a healthy cell, with an attendant increase in the passage time. Apart from the qualitative features, deformation indices were obtained. The deformation of hRBC was sudden, while the iRBC deformed slowly as it traversed through the constriction. For higher flow rates, both hRBC and iRBC were found to undergo severe deformation. Even under low flow rates, hRBC could easily traverse past the constricted channel. However, for sufficiently slow flow rates (eg, capillary flows), the microchannel was found to be completely blocked by the iRBC.

Mechanical property characterization of hundreds of single nuclei based on microfluidic constriction channel

As label-free biomarkers, the mechanical properties of nuclei are widely treated as promising biomechanical markers for cell type classification and cellular status evaluation. However, previously reported mechanical parameters were derived from only around 10 nuclei, lacking statistical significances due to low sample numbers. To address this issue, nuclei were first isolated from SW620 and A549 cells, respectively, using a chemical treatment method. This was followed by aspirating them through two types of microfluidic constriction channels for mechanical property characterization. In this study, hundreds of nuclei were characterized, producing passage times of 0.5 ± 1.2 s for SW620 nuclei in type I constriction channel (n = 153), 0.045 ± 0.047 s for SW620 nuclei in type II constriction channel (n = 215) and 0.50 ± 0.86 s for A549 nuclei in type II constriction channel. In addition, neural network based pattern recognition was used to classify the nuclei isolated from SW620 and A549 cells, producing successful classification rates of 87.2% for diameters of nuclei, 85.5% for passage times of nuclei and 89.3% for both passage times and diameters of nuclei. These results indicate that the characterization of the mechanical properties of nuclei may contribute to the classification of different tumor cells.

Local Hematocrit Fluctuation Induced by Malaria-Infected Red Blood Cells and Its Effect on Microflow

We investigate numerically the microscale blood flow in which red blood cells (RBCs) are partially infected by Plasmodium falciparum, the malaria parasite. The infected RBCs are modeled as more rigid cells with less deformability than healthy ones. Our study illustrates that, in a 10 μm microvessel in low-hematocrit conditions (18% and 27%), the Plasmodium falciparum-infected red blood cells (Pf-IRBCs) and healthy ones first form a train of cells. Because of the slow moving of the Pf-IRBCs, the local hematocrit (H_ct) near the Pf-IRBCs is then increased, to approximately 40% or even higher values. This increase of the local hematocrit is temporary and is kept for a longer length of time because of the long RBC train formed in 27%-H_ct condition. Similar hematocrit elevation at the downstream region with 45%-H_ct in the same 10 μm microvessel is also observed with the cells randomly located. In 20 μm microvessels with 45%-H_ct, the Pf-IRBCs slow down the velocity of the healthy red blood cells (HRBCs) around them and then locally elevate the volume fraction and result in the accumulation of the RBCs at the center of the vessels, thus leaving a thicker cell free layer (CFL) near the vessel wall than normal. Variation of wall shear stress (WSS) is caused by the fluctuation of local H_ct and the distance between the wall and the RBCs. Moreover, in high-hematocrit condition (45%), malaria-infected cells have a tendency to migrate to the edge of the aggregates which is due to the uninterrupted hydrodynamic interaction between the HRBCs and Pf-IRBC. Our results suggest that the existence of Pf-IRBCs is a nonnegligible factor for the fluctuation of hematocrit and WSS and also contributes to the increase of CFL of pathological blood flow in microvessels. The numerical approach presented has the potential to be utilized to RBC disorders and other hematologic diseases.

Investigation of red blood cell mechanical properties using AFM indentation and coarse-grained particle method

Background

Red blood cells (RBCs) deform significantly and repeatedly when passing through narrow capillaries and delivering dioxygen throughout the body. Deformability of RBCs is a key characteristic, largely governed by the mechanical properties of the cell membrane. This study investigated RBC mechanical properties using atomic force microscopy (AFM) with the aim to develop a coarse-grained particle method model to study for the first time RBC indentation in both 2D and 3D. This new model has the potential to be applied to further investigate the local deformability of RBCs, with accurate control over adhesion, probe geometry and position of applied force.

Results

The model considers the linear stretch capacity of the cytoskeleton, bending resistance and areal incompressibility of the bilayer, and volumetric incompressibility of the internal fluid. The model’s performance was validated against force–deformation experiments performed on RBCs under spherical AFM indentation. The model was then used to investigate the mechanisms which absorbed energy through the indentation stroke, and the impact of varying stiffness coefficients on the measured deformability. This study found the membrane’s bending stiffness was most influential in controlling RBC physical behaviour for indentations of up to 200 nm.

Conclusions

As the bilayer provides bending resistance, this infers that structural changes within the bilayer are responsible for the deformability changes experienced by deteriorating RBCs. The numerical model presented here established a foundation for future investigations into changes within the membrane that cause differences in stiffness between healthy and deteriorating RBCs, which have already been measured experimentally with AFM.

Electronic supplementary material

The online version of this article (10.1186/s12938-017-0429-5) contains supplementary material, which is available to authorized users.

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.

Off-center motion of a trapped elastic capsule in a microfluidic channel with a narrow constriction

Owing to their significance in capsule-related engineering and biomedical applications, a number of studies have considered the dynamics of elastic capsules flowing in constricted microchannels. However, these studies have focused on capsules moving along the channel centerline. In the present study, we numerically investigate the transient motion of an elastic capsule in a microfluidic channel with a rectangular constriction, which is initially trapped at the constriction inlet while off the channel centerline (i.e., on the channel bottom-wall). Under the push of the surrounding flow, the capsule can squeeze into the constriction, but only if the capsule deformability or the constriction size is sufficiently large. We find that the critical capillary number leading to the penetration of the capsule into the constriction is larger for off-centerline capsules compared to centered capsules. The centered capsule is stationary at the steady state when it remains stuck at the constriction; in contrast, the off-centerline capsule is not stationary but exhibits a tank-treading motion, i.e., its overall shape maintains a nonspherical shape with a protrusion into the constriction while its membrane exhibits a continuous rotation. Further, we examine the dependence of the capsule motion type, capsule deformation degree and membrane tension distribution on the capillary number (measuring the effects of flow strength and membrane mechanics) and constriction geometries (including the constriction height and width). Finally, we discuss the mechanism governing the capsule motion by analyzing the hydrodynamic forces acting on the capsule. The shear force acting on the capsule top owing to the fluid flow in the gap between the capsule top and the channel top-wall is the main source inducing the membrane tank-treading rotation.

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.

Hybrid smoothed dissipative particle dynamics and immersed boundary method for simulation of red blood cells in flows

In biofluid flow systems, often the flow problems of fluids of complex structures, such as the flow of red blood cells (RBCs) through complex capillary vessels, need to be considered. The smoothed dissipative particle dynamics (SDPD), a particle-based method, is one of the easy and flexible methods to model such complex structure fluids. It couples the best features of the smoothed particle hydrodynamics (SPH) and dissipative particle dynamics (DPD), with parameters having specific physical meaning (coming from SPH discretization of the Navier-Stokes equations), combined with thermal fluctuations in a mesoscale simulation, in a similar manner to the DPD. On the other hand, the immersed boundary method (IBM), a preferred method for handling fluid-structure interaction problems, has also been widely used to handle the fluid-RBC interaction in RBC simulations. In this paper, we aim to couple SDPD and IBM together to carry out the simulations of RBCs in complex flow problems. First, we develop the SDPD-IBM model in details, including the SDPD model for the evolving fluid flow, the RBC model for calculating RBC deformation force, the IBM for treating fluid-RBC interaction, and the solid boundary treatment model as well. We then conduct the verification and validation of the combined SDPD-IBM method. Finally, we demonstrate the capability of the SDPD-IBM method by simulating the flows of RBCs in rectangular, cylinder, curved, bifurcated, and constricted tubes, respectively.

A Two-Dimensional Numerical Investigation of Transport of Malaria-Infected Red Blood Cells in Stenotic Microchannels

The malaria-infected red blood cells experience a significant decrease in cell deformability and increase in cell membrane adhesion. Blood hemodynamics in microvessels is significantly affected by the alteration of the mechanical property as well as the aggregation of parasitized red blood cells. In this study, we aim to numerically study the connection between cell-level mechanobiological properties of human red blood cells and related malaria disease state by investigating the transport of multiple red blood cell aggregates passing through microchannels with symmetric stenosis. Effects of stenosis magnitude, aggregation strength, and cell deformability on cell rheology and flow characteristics were studied by a two-dimensional model using the fictitious domain-immersed boundary method. The results indicated that the motion and dissociation of red blood cell aggregates were influenced by these factors and the flow resistance increases with the increase of aggregating strength and cell stiffness. Further, the roughness of the velocity profile was enhanced by cell aggregation, which considerably affected the blood flow characteristics. The study may assist us in understanding cellular-level mechanisms in disease development.

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

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

Modeling the Mechanosensitivity of Neutrophils Passing through a Narrow Channel

Recent experiments have found that neutrophils may be activated after passing through microfluidic channels and filters. Mechanical deformation causes disassembly of the cytoskeleton and a sudden drop of the elastic modulus of the neutrophil. This fluidization is followed by either activation of the neutrophil with protrusion of pseudopods or a uniform recovery of the cytoskeleton network with no pseudopod. The former occurs if the neutrophil traverses the narrow channel at a slower rate. We propose a chemo-mechanical model for the fluidization and activation processes. Fluidization is treated as mechanical destruction of the cytoskeleton by sufficiently rapid bending. Loss of the cytoskeleton removes a pathway by which cortical tension inhibits the Rac protein. As a result, Rac rises and polarizes through a wave-pinning mechanism if the chemical reaction rate is fast enough. This leads to recovery and reinforcement of the cytoskeleton at the front of the neutrophil, and hence protrusion and activation. Otherwise the Rac signal returns to a uniform pre-deformation state and no activation occurs. Thus, mechanically induced neutrophil activation is understood as the competition between two timescales: that of chemical reaction and that of mechanical deformation. The model captures the main features of the experimental observation.

The development of malaria diagnostic techniques: a review of the approaches with focus on dielectrophoretic and magnetophoretic methods

The large number of deaths caused by malaria each year has increased interest in the development of effective malaria diagnoses. At the early-stage of infection, patients show non-specific symptoms or are asymptomatic, which makes it difficult for clinical diagnosis, especially in non-endemic areas. Alternative diagnostic methods that are timely and effective are required to identify infections, particularly in field settings. This article reviews conventional malaria diagnostic methods together with recently developed techniques for both malaria detection and infected erythrocyte separation. Although many alternative techniques have recently been proposed and studied, dielectrophoretic and magnetophoretic approaches are among the promising new techniques due to their high specificity for malaria parasite-infected red blood cells. The two approaches are discussed in detail, including their principles, types, applications and limitations. In addition, other recently developed techniques, such as cell deformability and morphology, are also overviewed in this article.

A numerical study on the elastic modulus of volume and area dilation for a deformable cell in a microchannel

A red blood cell (RBC) in a microfluidic channel is highly interesting for scientists in various fields of research on biological systems. This system has been studied extensively by empirical, analytical, and numerical methods. Nonetheless, research of predicting the behavior of an RBC in a microchannel is still an interesting area. The complications arise from deformation of an RBC and interactions among the surrounding fluid, wall, and RBCs. In this study, a pressure-driven RBC in a microchannel was simulated with a three-dimensional lattice Boltzmann method of an immersed boundary. First, the effect of boundary thickness on the interaction between the wall and cell was analyzed by measuring the time of passage through the narrow channel. Second, the effect of volume conservation stiffness was studied. Finally, the effect of global area stiffness was analyzed.

Microfluidic converging/diverging channels optimised for homogeneous extensional deformation

In this work, we optimise microfluidic converging/diverging geometries in order to produce constant strain-rates along the centreline of the flow, for performing studies under homogeneous extension. The design is examined for both two-dimensional and three-dimensional flows where the effects of aspect ratio and dimensionless contraction length are investigated. Initially, pressure driven flows of Newtonian fluids under creeping flow conditions are considered, which is a reasonable approximation in microfluidics, and the limits of the applicability of the design in terms of Reynolds numbers are investigated. The optimised geometry is then used for studying the flow of viscoelastic fluids and the practical limitations in terms of Weissenberg number are reported. Furthermore, the optimisation strategy is also applied for electro-osmotic driven flows, where the development of a plug-like velocity profile allows for a wider region of homogeneous extensional deformation in the flow field.

Dynamic and rheological properties of soft biological cell suspensions

Quantifying dynamic and rheological properties of suspensions of soft biological particles such as vesicles, capsules, and red blood cells (RBCs) is fundamentally important in computational biology and biomedical engineering. In this review, recent studies on dynamic and rheological behavior of soft biological cell suspensions by computer simulations are presented, considering both unbounded and confined shear flow. Furthermore, the hemodynamic and hemorheological characteristics of RBCs in diseases such as malaria and sickle cell anemia are highlighted.

Numerical simulation of a single cell passing through a narrow slit

The narrow slit between endothelial cells that line the microvessel wall is the principal pathway for tumor cell extravasation to the surrounding tissue. To understand this crucial step for tumor hematogenous metastasis, we used dissipative particle dynamics method to investigate an individual cell passing through a narrow slit numerically. The cell membrane was simulated by a spring-based network model which can separate the internal cytoplasm and surrounding fluid. The effects of the cell elasticity, cell shape, nucleus and slit size on the cell transmigration through the slit were investigated. Under a fixed driving force, the cell with higher elasticity can be elongated more and pass faster through the slit. When the slit width decreases to 2/3 of the cell diameter, the spherical cell becomes jammed despite reducing its elasticity modulus by 10 times. However, transforming the cell from a spherical to ellipsoidal shape and increasing the cell surface area by merely 9.3 % can enable the cell to pass through the narrow slit. Therefore, the cell shape and surface area increase play a more important role than the cell elasticity in cell passing through the narrow slit. In addition, the simulation results indicate that the cell migration velocity decreases during entrance but increases during exit of the slit, which is qualitatively in agreement with the experimental observation.

Deformability measurement of red blood cells using a microfluidic channel array and an air cavity in a driving syringe with high throughput and precise detection of subpopulations

Red blood cell (RBC) deformability has been considered a potential biomarker for monitoring pathological disorders. High throughput and detection of subpopulations in RBCs are essential in the measurement of RBC deformability. In this paper, we propose a new method to measure RBC deformability by evaluating temporal variations in the average velocity of blood flow and image intensity of successively clogged RBCs in the microfluidic channel array for specific time durations. In addition, to effectively detect differences in subpopulations of RBCs, an air compliance effect is employed by adding an air cavity into a disposable syringe. The syringe was equally filled with a blood sample (V(blood) = 0.3 mL, hematocrit = 50%) and air (V(air) = 0.3 mL). Owing to the air compliance effect, blood flow in the microfluidic device behaved transiently depending on the fluidic resistance in the microfluidic device. Based on the transient behaviors of blood flows, the deformability of RBCs is quantified by evaluating three representative parameters, namely, minimum value of the average velocity of blood flow, clogging index, and delivered blood volume. The proposed method was applied to measure the deformability of blood samples consisting of homogeneous RBCs fixed with four different concentrations of glutaraldehyde solution (0%-0.23%). The proposed method was also employed to evaluate the deformability of blood samples partially mixed with normal RBCs and hardened RBCs. Thereafter, the deformability of RBCs infected by human malaria parasite Plasmodium falciparum was measured. As a result, the three parameters significantly varied, depending on the degree of deformability. In addition, the deformability measurement of blood samples was successfully completed in a short time (∼10 min). Therefore, the proposed method has significant potential in deformability measurement of blood samples containing hematological diseases with high throughput and precise detection of subpopulations in RBCs.

Inflow/Outflow Boundary Conditions for Particle-Based Blood Flow Simulations: Application to Arterial Bifurcations and Trees

When blood flows through a bifurcation, red blood cells (RBCs) travel into side branches at different hematocrit levels, and it is even possible that all RBCs enter into one branch only, leading to a complete separation of plasma and RBCs. To quantify this phenomenon via particle-based mesoscopic simulations, we developed a general framework for open boundary conditions in multiphase flows that is effective even for high hematocrit levels. The inflow at the inlet is duplicated from a fully developed flow generated in a pilot simulation with periodic boundary conditions. The outflow is controlled by adaptive forces to maintain the flow rate and velocity gradient at fixed values, while the particles leaving the arteriole at the outlet are removed from the system. Upon validation of this approach, we performed systematic 3D simulations to study plasma skimming in arterioles of diameters 20 to 32 microns. For a flow rate ratio 6:1 at the branches, we observed the “all-or-nothing” phenomenon with plasma only entering the low flow rate branch. We then simulated blood-plasma separation in arteriolar bifurcations with different bifurcation angles and same diameter of the daughter branches. Our simulations predict a significant increase in RBC flux through the main daughter branch as the bifurcation angle is increased. Finally, we demonstrated the effectiveness of the new methodology in simulations of blood flow in vessels with multiple inlets and outlets, constructed using an angiogenesis model.

Blood tests, which provide a wealth of information on the state of human health, are often performed on cell-free samples. Therefore, blood-plasma separation needs to be achieved. A simple but effective solution for isolating plasma from blood utilizes capillary bifurcations. In a particle-based simulation study of plasma skimming in capillary bifurcations, the blood flow properties such as velocity and pressure fields differ drastically at the inlet and outlet regions. Therefore, a new open (non-periodic) boundary is required. In this paper, we have developed and validated a general parallel framework for open boundary conditions. This is a non-trivial enabling technology that could be used in all open boundary systems and all particle-based Lagrangian simulations. We performed systematic 3D simulations of blood flow in arteriolar bifurcations and elucidated the biophysical mechanism of blood-plasma separation as well as quantified the effects of branch size and bifurcation angle on cell separation efficiency, which have not been addressed before. We also demonstrated the applicability of the methodology in arterial trees with multiple inlets and outlets.

Biomechanical properties of red blood cells in health and disease towards microfluidics

Red blood cells (RBCs) possess a unique capacity for undergoing cellular deformation to navigate across various human microcirculation vessels, enabling them to pass through capillaries that are smaller than their diameter and to carry out their role as gas carriers between blood and tissues. Since there is growing evidence that red blood cell deformability is impaired in some pathological conditions, measurement of RBC deformability has been the focus of numerous studies over the past decades. Nevertheless, reports on healthy and pathological RBCs are currently limited and, in many cases, are not expressed in terms of well-defined cell membrane parameters such as elasticity and viscosity. Hence, it is often difficult to integrate these results into the basic understanding of RBC behaviour, as well as into clinical applications. The aim of this review is to summarize currently available reports on RBC deformability and to highlight its association with various human diseases such as hereditary disorders (e.g., spherocytosis, elliptocytosis, ovalocytosis, and stomatocytosis), metabolic disorders (e.g., diabetes, hypercholesterolemia, obesity), adenosine triphosphate-induced membrane changes, oxidative stress, and paroxysmal nocturnal hemoglobinuria. Microfluidic techniques have been identified as the key to develop state-of-the-art dynamic experimental models for elucidating the significance of RBC membrane alterations in pathological conditions and the role that such alterations play in the microvasculature flow dynamics.

Microfluidic-based measurement of erythrocyte sedimentation rate for biophysical assessment of blood in an in vivo malaria-infected mouse

This study suggests a new erythrocyte sedimentation rate (ESR) measurement method for the biophysical assessment of blood by using a microfluidic device. For an effective ESR measurement, a disposable syringe filled with blood is turned upside down and aligned at 180° with respect to gravitational direction. When the blood sample is delivered into the microfluidic device from the top position of the syringe, the hematocrit of blood flowing in the microfluidic channel decreases because the red blood cell-depleted region is increased from the top region of the syringe. The variation of hematocrit is evaluated by consecutively capturing images and conducting digital image processing technique for 10 min. The dynamic variation of ESR is quantitatively evaluated using two representative parameters, namely, time constant (λ) and ESR-area (AESR). To check the performance of the proposed method, blood samples with various ESR values are prepared by adding different concentrations of dextran solution. λ and AESR are quantitatively evaluated by using the proposed method and a conventional method, respectively. The proposed method can be used to measure ESR with superior reliability, compared with the conventional method. The proposed method can also be used to quantify ESR of blood collected from malaria-infected mouse under in vivo condition. To indirectly compare with the results obtained by the proposed method, the viscosity and velocity of the blood are measured using the microfluidic device. As a result, the biophysical properties, including ESR and viscosity of blood, are significantly influenced by the parasitemia level. These experimental demonstrations support the notion that the proposed method is capable of effectively monitoring the biophysical properties of blood.

Continuous size-based separation of microparticles in a microchannel with symmetric sharp corner structures

A new microchannel with a series of symmetric sharp corner structures is reported for passive size-dependent particle separation. Micro particles of different sizes can be completely separated based on the combination of the inertial lift force and the centrifugal force induced by the sharp corner structures in the microchannel. At appropriate flow rate and Reynolds number, the centrifugal force effect on large particles, induced by the sharp corner structures, is stronger than that on small particles; hence after passing a series of symmetric sharp corner structures, large particles are focused to the center of the microchannel, while small particles are focused at two particle streams near the two side walls of the microchannel. Particles of different sizes can then be completely separated. Particle separation with this device was demonstrated using 7.32 μm and 15.5 μm micro particles. Experiments show that in comparison with the prior multi-orifice flow fractionation microchannel and multistage-multiorifice flow fractionation microchannel, this device can completely separate two-size particles with narrower particle stream band and larger separation distance between particle streams. In addition, it requires no sheath flow and complex multi-stage separation structures, avoiding the dilution of analyte sample and complex operations. The device has potentials to be used for continuous, complete particle separation in a variety of lab-on-a-chip and biomedical applications.

Red blood cell dynamics in polymer brush-coated microcapillaries: A model of endothelial glycocalyx in vitro

The confined flow of red blood cells (RBCs) in microvasculature is essential for oxygen delivery to body tissues and has been extensively investigated in the literature, both in vivo and in vitro. One of the main problems still open in microcirculation is that flow resistance in microcapillaries in vivo is higher than that in vitro. This discrepancy has been attributed to the glycocalyx, a macromolecular layer lining the inner walls of vessels in vivo, but no direct experimental evidence of this hypothesis has been provided so far. Here, we investigate the flow behavior of RBCs in glass microcapillaries coated with a polymer brush (referred to as "hairy" microcapillaries as opposed to "bare" ones with no coating), an experimental model system of the glycocalyx. By high-speed microscopy imaging and image analysis, a velocity reduction of RBCs flowing in hairy microcapillaries as compared to bare ones is indeed found at the same pressure drop. Interestingly, such slowing down is larger than expected from lumen reduction due to the polymer brush and displays an on-off trend with a threshold around 70 nm of polymer brush dry thickness. Above this threshold, the presence of the polymer brush is associated with an increased RBC deformation, and RBC velocity is independent on polymer brush thickness (at the same pressure drop). In conclusion, this work provides direct support to the hypothesis that the glycocalyx is the main factor responsible of the higher flow resistance found in microcapillaries in vivo.

Real-time visualization of 3-D dynamic microscopic objects using optical diffraction tomography

3-D refractive index (RI) distribution is an intrinsic bio-marker for the chemical and structural information about biological cells. Here we develop an optical diffraction tomography technique for the real-time reconstruction of 3-D RI distribution, employing sparse angle illumination and a graphic processing unit (GPU) implementation. The execution time for the tomographic reconstruction is 0.21 s for 96(3) voxels, which is 17 times faster than that of a conventional approach. We demonstrated the real-time visualization capability with imaging the dynamics of Brownian motion of an anisotropic colloidal dimer and the dynamic shape change in a red blood cell upon shear flow.