Analysis of red blood cell deformation under fast shear flow for better estimation of hemolysis

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.

A predictive model using risk factor categories for hospital-acquired pneumonia in patients with aneurysmal subarachnoid hemorrhage

Objectives

To identify risk factors for hospital-acquired pneumonia (HAP) in patients with aneurysmal subarachnoid hemorrhage (aSAH) and establish a predictive model to aid evaluation.

Methods

The cohorts of 253 aSAH patients were divided into the HAP group (n = 64) and the non-HAP group (n = 189). Univariate and multivariate logistic regression were performed to identify risk factors. A logistic model (Model-Logit) was established based on the independent risk factors. We used risk factor categories to develop a model (Model-Cat). Receiver operating characteristic curves were generated to determine the cutoff values. Areas under the curves (AUCs) were calculated to assess the accuracy of models and single factors. The Delong test was performed to compare the AUCs.

Results

The multivariate logistic analysis showed that the age [p = 0.012, odds ratio (OR) = 1.059, confidence interval (CI) = 1.013-1.107], blood glucose (BG; >7.22 mmol/L; p = 0.011, OR = 2.781, CI = 1.263-6.119), red blood distribution width standard deviation (RDW-SD; p = 0.024, OR = 1.118, CI = 1.015-1.231), and Glasgow coma scale (GCS; p < 0.001, OR = 0.710, CI = 0.633-0.798) were independent risk factors. The Model-Logit was as follows: Logit(P) = -5.467 + 0.057 * Age + 1.023 * BG (>7.22 mmol/L, yes = 1, no = 0) + 0.111 * RDW-SD-0.342 * GCS. The AUCs values of the Model-Logit, GCS, age, BG (>7.22 mmol/L), and RDW-SD were 0.865, 0.819, 0.634, 0.698, and 0.625, respectively. For clinical use, the Model-Cat was established. In the Model-Cat, the AUCs for GCS, age, BG, and RDW-SD were 0.850, 0.760, 0.700, 0.641, and 0.564, respectively. The AUCs of the Model-Logit were insignificantly higher than the Model-Cat (Delong test, p = 0.157). The total points from -3 to 4 and 5 to 14 were classified as low- and high-risk levels, respectively.

Conclusions

Age, BG (> 7.22 mmol/L), GCS, and RDW-SD were independent risk factors for HAP in aSAH patients. The Model-Cat was convenient for practical evaluation. The aSAH patients with total points from 5 to 14 had a high risk for HAP, suggesting the need for more attention during treatment.

Device-Induced Hemostatic Disorders in Mechanically Assisted Circulation

Mechanically assisted circulation (MAC) sustains the blood circulation in the body of a patients undergoing cardiac surgery with cardiopulmonary bypass (CPB) or on ventricular assistance with a ventricular assist device (VAD) or on extracorporeal membrane oxygenation (ECMO) with a pump-oxygenator system. While MAC provides short-term (days to weeks) support and long-term (months to years) for the heart and/or lungs, the blood is inevitably exposed to non-physiological shear stress (NPSS) due to mechanical pumping action and in contact with artificial surfaces. NPSS is well known to cause blood damage and functional alterations of blood cells. In this review, we discussed shear-induced platelet adhesion, platelet aggregation, platelet receptor shedding, and platelet apoptosis, shear-induced acquired von Willebrand syndrome (AVWS), shear-induced hemolysis and microparticle formation during MAC. These alterations are associated with perioperative bleeding and thrombotic events, morbidity and mortality, and quality of life in MCS patients. Understanding the mechanism of shear-induce hemostatic disorders will help us develop low-shear-stress devices and select more effective treatments for better clinical outcomes.

The ratio of superoxide dismutase to standard deviation of erythrocyte distribution width as a predictor of systemic lupus erythematosus

BACKGROUND

To explore the clinical value of the serum superoxide dismutase-to-standard deviation of erythrocyte distribution width ratio (SRSR) in systemic lupus erythematosus (SLE).

METHODS

A total of 222 SLE patients from the Rheumatology and Immunology Department in the Second Affiliated Hospital of Chongqing Medical University from January 2017 to April 2019 were collected as the experimental group, and a total of 202 healthy physical examiners were extracted as the control group. Neutrophil-to-lymphocyte ratio (NLR), superoxide dismutase-to-standard deviation of erythrocyte distribution width ratio (SRSR), and platelet-to-lymphocyte ratio (PLR) were calculated from the collected data and then compared the level of the above three indexes between the two groups. In addition, we analyzed the association between SRSR and clinically relevant indicators.

RESULTS

We found that the SRSR of SLE patients was significantly lower than healthy control group, by analyzing the receiver operating characteristic (ROC) curve; it revealed that the SRSR had higher specificity and sensitivity than either superoxide dismutase (SOD) or standard deviation of erythrocyte distribution width (RDW-SD) alone. The area under the curve (AUC) for SRSR was significantly larger than either SOD or RDW-SD alone, and the AUC for SRSR was also larger than NLR and PLR. And it was found that SRSR was independently correlated with SLE disease activity through multiple linear regression analysis.

CONCLUSION

SRSR is a useful biomarker for the diagnosis of SLE, and it is of great significance in the clinical application.

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.

Applications of machine learning for simulations of red blood cells in microfluidic devices

Background

For optimization of microfluidic devices for the analysis of blood samples, it is useful to simulate blood cells as elastic objects in flow of blood plasma. In such numerical models, we primarily need to take into consideration the movement and behavior of the dominant component of the blood, the red blood cells. This can be done quite precisely in small channels and within a short timeframe. However, larger volumes or timescales require different approaches. Instead of simplifying the simulation, we use a neural network to predict the movement of the red blood cells.

Results

The neural network uses data from the numerical simulation for learning, however, the simulation needs only be run once. Alternatively, the data could come from video processing of a recording of a biological experiment. Afterwards, the network is able to predict the movement of the red blood cells because it is a system of bases that gives an approximate cell velocity at each point of the simulation channel as a linear combination of bases.In a simple box geometry, the neural network gives results comparable to predictions using fluid streamlines, however in a channel with obstacles forming slits, the neural network is about five times more accurate.The network can also be used as a discriminator between different situations. We observe about two-fold increase in mean relative error when a network trained on one geometry is used to predict trajectories in a modified geometry. Even larger increase was observed when it was used to predict trajectories of cells with different elastic properties.

Conclusions

While for uncomplicated box channels there is no advantage in using a system of bases instead of a simple prediction using fluid streamlines, in a more complicated geometry, the neural network is significantly more accurate. Another application of this system of bases is using it as a comparison tool for different modeled situations. This has a significant future potential when applied to processing data from videos of microfluidic flows.

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.

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.

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.

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.

Expediting red blood cell transfusions by syringing causes significant hemolysis

BACKGROUND

Techniques commonly used to expedite blood transfusions include pneumatically pressurizing red blood cell (RBC) bags or manual syringing its contents. We compared these techniques on RBC hemolysis using a simulated transfusion model.

STUDY DESIGN AND METHODS

Fifteen warmed RBC units that were 12.3 ± 4.3 (95% confidence interval [CI], 10.1-14.5) days old were each subjected to two experimental rapid transfusion techniques. RBCs from each technique were directed through 18- and 22-gauge cannulas attached to blood administration sets. One technique involved RBC bag pressurization to 300 mmHg. The other employed a 20-mL syringe to effect forceful, manual aspiration from the RBC bag followed by forceful, manual RBC injection. The control group was gravity driven without cannulas. Free hemoglobin (Hb) concentrations were measured and percent hemolysis was calculated.

RESULTS

Free Hb concentrations and percent hemolysis (median [95% CI]) were similar in the control (0.05 [0.03-0.08] g/dL and 0.13% [0.09%-0.17%], respectively) and pressurized experiments (0.06 [0.05-0.09] g/dL; 0.14% [0.12%-0.22%]), respectively. Syringing resulted in 10-fold higher free Hb concentrations (0.55 [0.38-0.92] g/dL) and percent hemolysis (1.28% [1.03%-2.15%]) than when employing the control (p < 0.0001) or pressurization (p < 0.0001) techniques. Cannula sizes studied did not affect hemolysis.

CONCLUSION

Forceful manual syringing caused significant hemolysis and high free Hb concentrations. Pressurizing RBC bags induced no more hemolysis than after gravity-facilitated transfusions. Syringing to expedite RBC transfusions should be avoided in favor of pneumatic RBC bag pressurization.

Extending the Power-Law Hemolysis Model to Complex Flows

Hemolysis (damage to red blood cells) is a long-standing problem in blood contacting devices, and its prediction has been the goal of considerable research. The most popular model relating hemolysis to fluid stresses is the power-law model, which was developed from experiments in pure shear only. In the absence of better data, this model has been extended to more complex flows by replacing the shear stress in the power-law equation with a von Mises-like scalar stress. While the validity of the scalar stress also remains to be confirmed, inconsistencies exist in its application, in particular, two forms that vary by a factor of 2 have been used. This article will clarify the proper extension of the power law to complex flows in a way that maintains correct results in the limit of pure shear.

Characterization of Nanoparticle Dispersion in Red Blood Cell Suspension by the Lattice Boltzmann-Immersed Boundary Method

Nanodrug-carrier delivery in the blood stream is strongly influenced by nanoparticle (NP) dispersion. This paper presents a numerical study on NP transport and dispersion in red blood cell (RBC) suspensions under shear and channel flow conditions, utilizing an immersed boundary fluid-structure interaction model with a lattice Boltzmann fluid solver, an elastic cell membrane model and a particle motion model driven by both hydrodynamic loading and Brownian dynamics. The model can capture the multiphase features of the blood flow. Simulations were performed to obtain an empirical formula to predict NP dispersion rate for a range of shear rates and cell concentrations. NP dispersion rate predictions from the formula were then compared to observations from previous experimental and numerical studies. The proposed formula is shown to accurately predict the NP dispersion rate. The simulation results also confirm previous findings that the NP dispersion rate is strongly influenced by local disturbances in the flow due to RBC motion and deformation. The proposed formula provides an efficient method for estimating the NP dispersion rate in modeling NP transport in large-scale vascular networks without explicit RBC and NP models.

Multi-scale biological and physical modelling of the tumour micro-environment

Paced by advances in high performance computing, and algorithms for multi-physics and multi-scale simulation, a number of groups have recently established numerical models of flowing blood systems, where cell-scale interactions are explicitly resolved. To be biologically representative, these models account for some or all of: (1) fluid dynamics of the carrier flow, (2) structural dynamics of the cells and vessel walls, (3) interaction and transport biochemistry, and, (4) methods for scaling to physiologically representative numbers of cells. In this article, our interest is the modelling of the tumour micro-environment. We review the broader area of cell-scale resolving blood flow modelling, while focusing on the particular interactions of tumour cells and white blood cells, known to play an important role in metastasis.