Multiscale modeling of blood flow: from single cells to blood rheology

Collective behavior of squirmers in thin films

Bacteria in biofilms form complex structures and can collectively migrate within mobile aggregates, which is referred to as swarming. This behavior is influenced by a combination of various factors, including morphological characteristics and propulsive forces of swimmers, their volume fraction within a confined environment, and hydrodynamic and steric interactions between them. In our study, we employ the squirmer model for microswimmers and the dissipative particle dynamics method for fluid modeling to investigate the collective motion of swimmers in thin films. The film thickness permits a free orientation of non-spherical squirmers, but constraints them to form a two-layered structure at maximum. Structural and dynamic properties of squirmer suspensions confined within the slit are analyzed for different volume fractions of swimmers, motility types (e.g., pusher, neutral squirmer, puller), and the presence of a rotlet dipolar flow field, which mimics the counter-rotating flow generated by flagellated bacteria. Different states are characterized, including a gas-like phase, swarming, and motility-induced phase separation, as a function of increasing volume fraction. Our study highlights the importance of an anisotropic swimmer shape, hydrodynamic interactions between squirmers, and their interaction with the walls for the emergence of different collective behaviors. Interestingly, the formation of collective structures may not be symmetric with respect to the two walls. Furthermore, the presence of a rotlet dipole significantly mitigates differences in the collective behavior between various swimmer types. These results contribute to a better understanding of the formation of bacterial biofilms and the emergence of collective states in confined active matter.

A method to prevent clogging and clustering in microfluidic systems using microbubble streaming

This paper presents an innovative strategy to address the issues of clogging and cluster-related challenges in microchannels within microfluidic devices. Leveraging three-dimensional (3D) microbubble streaming as a dynamic solution, our approach involves the controlled activation of microbubbles near channel constrictions, inducing microstreaming with distinctive features. This microstreaming, characterized by a high non-uniform 3D gradient and significant shear stress, effectively inhibits arch formation at constrictions and disintegrates particle clusters, demonstrating real-time prevention of clogging incidents and blockages. This study includes experimental validation of the anti-clogging technique, a detailed examination of microstreaming phenomena, and their effects on clogging and clustering issues. It also incorporates statistical analyses performed in various scenarios to verify the method's effectiveness and adaptability. Moreover, a versatile control system has been designed that operates in event-triggered, continuous, or periodic modes, which suits different lab-on-a-chip applications and improves the overall functionality of microfluidic systems.

Dynamic framework for large-scale modeling of membranes and peripheral proteins

In this chapter, we present a novel computational framework to study the dynamic behavior of extensive membrane systems, potentially in interaction with peripheral proteins, as an alternative to conventional simulation methods. The framework effectively describes the complex dynamics in protein-membrane systems in a mesoscopic particle-based setup. Furthermore, leveraging the hydrodynamic coupling between the membrane and its surrounding solvent, the coarse-grained model grounds its dynamics in macroscopic kinetic properties such as viscosity and diffusion coefficients, marrying the advantages of continuum- and particle-based approaches. We introduce the theoretical background and the parameter-space optimization method in a step-by-step fashion, present the hydrodynamic coupling method in detail, and demonstrate the application of the model at each stage through illuminating examples. We believe this modeling framework to hold great potential for simulating membrane and protein systems at biological spatiotemporal scales, and offer substantial flexibility for further development and parametrization.

Coronary hemodynamic simulation study

In this paper, a two-way fluid-structure coupling model is developed to simulate and analyze the hemodynamic process based on dynamic coronary angiography, and examine the influence of different hemodynamic parameters on coronary arteries in typical coronary stenosis lesions. Using the measured FFR pressure data of a patient, the pressure-time function curve is fitted to ensure the accuracy of the boundary conditions. The average error of the simulation pressure results compared to the test data is 6.74%. In addition, the results related to blood flow, pressure contour and wall shear stress contour in a typical cardiac cycle are obtained by simulation analysis. These results are found to be in good agreement with the laws of the real cardiac cycle, which verifies the rationality of the simulation. In conclusion, based on the modeling and hemodynamic simulation analysis process of dynamic coronary angiography, this paper proposes a method to assist the analysis and evaluation of coronary hemodynamic and functional parameters, which has certain practical significance.

Motion of microswimmers in cylindrical microchannels

Biological and artificial microswimmers often have to propel through a variety of environments, ranging from heterogeneous suspending media to strong geometrical confinement. Under confinement, local flow fields generated by microswimmers, and steric and hydrodynamic interactions with their environment determine the locomotion. We propose a squirmer-like model to describe the motion of microswimmers in cylindrical microchannels, where propulsion is generated by a fixed surface slip velocity. The model is studied using an approximate analytical solution for cylindrical swimmer shapes, and by numerical hydrodynamics simulations for spherical and spheroidal shapes. For the numerical simulations, we employ the dissipative particle dynamics method for modelling fluid flow. Both the analytical model and simulations show that the propulsion force increases with increasing confinement. However, the swimming velocity under confinement remains lower than the swimmer speed without confinement for all investigated conditions. In simulations, different swimming modes (i.e. pusher, neutral, puller) are investigated, and found to play a significant role in the generation of propulsion force when a swimmer approaches a dead end of a capillary tube. Propulsion generation in confined systems is local, such that the generated flow field generally vanishes beyond the characteristic size of the swimmer. These results contribute to a better understanding of microswimmer force generation and propulsion under strong confinement, including the motion in porous media and in narrow channels.

Confinement effect on the microcapillary flow and shape of red blood cells

The ability to change shape is essential for the proper functioning of red blood cells (RBCs) within the microvasculature. The shape of RBCs significantly influences blood flow and has been employed in microfluidic lab-on-a-chip devices, serving as a diagnostic biomarker for specific pathologies and enabling the assessment of RBC deformability. While external flow conditions, such as the vessel size and the flow velocity, are known to impact microscale RBC flow, our comprehensive understanding of how their shape-adapting ability is influenced by channel confinement in biomedical applications remains incomplete. This study explores the impact of various rectangular and square channels, each with different confinement and aspect ratios, on the in vitro RBC flow behavior and characteristic shapes. We demonstrate that rectangular microchannels, with a height similar to the RBC diameter in combination with a confinement ratio exceeding 0.9, are required to generate distinctive well-defined croissant and slipper-like RBC shapes. These shapes are characterized by their equilibrium positions in the channel cross section, and we observe a strong elongation of both stable shapes in response to the shear rate across the different channels. Less confined channel configurations lead to the emergence of unstable other shape types that display rich shape dynamics. Our work establishes an experimental framework to understand the influence of channel size on the single-cell flow behavior of RBCs, providing valuable insights for the design of biomicrofluidic single-cell analysis applications.

Venous Thromboembolism: Review of Clinical Challenges, Biology, Assessment, Treatment, and Modeling

Venous thromboembolism (VTE) is a massive clinical challenge, annually affecting millions of patients globally. VTE is a particularly consequential pathology, as incidence is correlated with extremely common risk factors, and a large cohort of patients experience recurrent VTE after initial intervention. Altered hemodynamics, hypercoagulability, and damaged vascular tissue cause deep-vein thrombosis and pulmonary embolism, the two permutations of VTE. Venous valves have been identified as likely locations for initial blood clot formation, but the exact pathway by which thrombosis occurs in this environment is not entirely clear. Several risk factors are known to increase the likelihood of VTE, particularly those that increase inflammation and coagulability, increase venous resistance, and damage the endothelial lining. While these risk factors are useful as predictive tools, VTE diagnosis prior to presentation of outward symptoms is difficult, chiefly due to challenges in successfully imaging deep-vein thrombi. Clinically, VTE can be managed by anticoagulants or mechanical intervention. Recently, direct oral anticoagulants and catheter-directed thrombolysis have emerged as leading tools in resolution of venous thrombosis. While a satisfactory VTE model has yet to be developed, recent strides have been made in advancing in silico models of venous hemodynamics, hemorheology, fluid-structure interaction, and clot growth. These models are often guided by imaging-informed boundary conditions or inspired by benchtop animal models. These gaps in knowledge are critical targets to address necessary improvements in prediction and diagnosis, clinical management, and VTE experimental and computational models.

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

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

Energy spectrum analysis on a red blood cell model

It is important to understand the dynamics of red blood cells (RBCs) in blood flow. This requires the formulation of coarse-grained RBC models that reproduce the hydrodynamic properties of blood accurately. One of the models that successfully reproduces the rheology and morphology of blood has been proposed by Fedosov et al. [Comput. Methods Appl. Mech. Eng. 199, 1937-1948 (2010)]. The proposed RBC model contains several parameters whose values are determined by either various experiments or physical requirements. In this study, we developed a new method of determining parameter values precisely from the fluctuations of the RBC membrane. Specifically, we studied the relationship between the spectra of the fluctuations and model parameters. Characteristic peaks were observed in the spectra, whose peak frequencies were dependent on the parameter values. In addition, we investigated the spectra of the radius of gyration. We identified the peaks originating from the spring potential and the volume-conserving potential appearing in the spectra. These results lead to the precise experimental determination of the parameters used in the RBC model.

Nanoparticle-mediated cancer cell therapy: basic science to clinical applications

Every sixth person in the world dies due to cancer, making it the second leading severe cause of death after cardiovascular diseases. According to WHO, cancer claimed nearly 10 million deaths in 2020. The most common types of cancers reported have been breast (lung, colon and rectum, prostate cases), skin (non-melanoma) and stomach. In addition to surgery, the most widely used traditional types of anti-cancer treatment are radio- and chemotherapy. However, these do not distinguish between normal and malignant cells. Additional treatment methods have evolved over time for early detection and targeted therapy of cancer. However, each method has its limitations and the associated treatment costs are quite high with adverse effects on the quality of life of patients. Use of individual atoms or a cluster of atoms (nanoparticles) can cause a paradigm shift by virtue of providing point of sight sensing and diagnosis of cancer. Nanoparticles (1-100 nm in size) are 1000 times smaller in size than the human cell and endowed with safer relocation capability to attack mechanically and chemically at a precise location which is one avenue that can be used to destroy cancer cells precisely. This review summarises the extant understanding and the work done in this area to pave the way for physicians to accelerate the use of hybrid mode of treatments by leveraging the use of various nanoparticles.

Assessing Differential Particle Deformability under Microfluidic Flow Conditions

Assessing the mechanical behavior of nano- and micron-scale particles with complex shapes is fundamental in drug delivery. Although different techniques are available to quantify the bulk stiffness in static conditions, there is still uncertainty in assessing particle deformability in dynamic conditions. Here, a microfluidic chip is designed, engineered, and validated as a platform to assess the mechanical behavior of fluid-borne particles. Specifically, potassium hydroxide (KOH) wet etching was used to realize a channel incorporating a series of micropillars (filtering modules) with different geometries and openings, acting as microfilters in the direction of the flow. These filtering modules were designed with progressively decreasing openings, ranging in size from about 5 down to 1 μm. Discoidal polymeric nanoconstructs (DPNs), with a diameter of 5.5 μm and a height of 400 nm, were realized with different poly(lactic-co-glycolic acid) (PLGA) and poly(ethylene glycol) (PEG) ratios (PLGA/PEG), namely, 5:1 and 1:0, resulting in soft and rigid particles, respectively. Given the peculiar geometry of DPNs, the channel height was kept to 5 μm to limit particle tumbling or flipping along the flow. After thorough physicochemical and morphological characterization, DPNs were tested within the microfluidic chip to investigate their behavior under flow. As expected, most rigid DPNs were trapped in the first series of pillars, whereas soft DPNs were observed to cross multiple filtering modules and reach the micropillars with the smallest opening (1 μm). This experimental evidence was also supported by computational tools, where DPNs were modeled as a network of springs and beads immersed in a Newtonian fluid using the smoothed particle hydrodynamics (SPH) method. This preliminary study presents a combined experimental-computational framework to quantify, compare, and analyze the characteristics of particles having complex geometrical and mechanical attributes under flow conditions.

Effect of Cell Age and Membrane Rigidity on Red Blood Cell Shape in Capillary Flow

Blood flow in the microcirculatory system is crucially affected by intrinsic red blood cell (RBC) properties, such as their deformability. In the smallest vessels of this network, RBCs adapt their shapes to the flow conditions. Although it is known that the age of RBCs modifies their physical properties, such as increased cytosol viscosity and altered viscoelastic membrane properties, the evolution of their shape-adapting abilities during senescence remains unclear. In this study, we investigated the effect of RBC properties on the microcapillary in vitro flow behavior and their characteristic shapes in microfluidic channels. For this, we fractioned RBCs from healthy donors according to their age. Moreover, the membranes of fresh RBCs were chemically rigidified using diamide to study the effect of isolated graded-membrane rigidity. Our results show that a fraction of stable, asymmetric, off-centered slipper-like cells at high velocities decreases with increasing age or diamide concentration. However, while old cells form an enhanced number of stable symmetric croissants at the channel centerline, this shape class is suppressed for purely rigidified cells with diamide. Our study provides further knowledge about the distinct effects of age-related changes of intrinsic cell properties on the single-cell flow behavior of RBCs in confined flows due to inter-cellular age-related cell heterogeneity.

A computational model for the transit of a cancer cell through a constricted microchannel

We propose a three-dimensional computational model to simulate the transient deformation of suspended cancer cells flowing through a constricted microchannel. We model the cell as a liquid droplet enclosed by a viscoelastic membrane, and its nucleus as a smaller stiffer capsule. The cell deformation and its interaction with the suspending fluid are solved through a well-tested immersed boundary lattice Boltzmann method. To identify a minimal mechanical model that can quantitatively predict the transient cell deformation in a constricted channel, we conduct extensive parametric studies of the effects of the rheology of the cell membrane, cytoplasm and nucleus and compare the results with a recent experiment conducted on human leukaemia cells. We find that excellent agreement with the experiment can be achieved by employing a viscoelastic cell membrane model with the membrane viscosity depending on its mode of deformation (shear versus elongation). The cell nucleus limits the overall deformation of the whole cell, and its effect increases with the nucleus size. The present computational model may be used to guide the design of microfluidic devices to sort cancer cells, or to inversely infer cell mechanical properties from their flow-induced deformation.

Measuring the red blood cell shape in capillary flow using spectrally encoded flow cytometry

Red blood cells in small capillaries exhibit a wide variety of deformations that reflect their true physiological conditions at these important locations. By applying a technique for the high-speed microscopy of flowing cells, termed spectrally encoded flow cytometry (SEFC), we image the light reflected from the red blood cells in human capillaries, and propose an analytical slipper-like model for the cell morphology that can reproduce the experimental in vivo images. The results of this work would be useful for studying the unique flow conditions in these vessels, and for extracting useful clinical parameters that reflect the true physiology of the blood cells in situ.

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.

Patient-Specific Inverse Modeling of In Vivo Cardiovascular Mechanics with Medical Image-Derived Kinematics as Input Data: Concepts, Methods, and Applications

Inverse modeling approaches in cardiovascular medicine are a collection of methodologies that can provide non-invasive patient-specific estimations of tissue properties, mechanical loads, and other mechanics-based risk factors using medical imaging as inputs. Its incorporation into clinical practice has the potential to improve diagnosis and treatment planning with low associated risks and costs. These methods have become available for medical applications mainly due to the continuing development of image-based kinematic techniques, the maturity of the associated theories describing cardiovascular function, and recent progress in computer science, modeling, and simulation engineering. Inverse method applications are multidisciplinary, requiring tailored solutions to the available clinical data, pathology of interest, and available computational resources. Herein, we review biomechanical modeling and simulation principles, methods of solving inverse problems, and techniques for image-based kinematic analysis. In the final section, the major advances in inverse modeling of human cardiovascular mechanics since its early development in the early 2000s are reviewed with emphasis on method-specific descriptions, results, and conclusions. We draw selected studies on healthy and diseased hearts, aortas, and pulmonary arteries achieved through the incorporation of tissue mechanics, hemodynamics, and fluid-structure interaction methods paired with patient-specific data acquired with medical imaging in inverse modeling approaches.

Of vascular defense, hemostasis, cancer, and platelet biology: an evolutionary perspective

We have established considerable expertise in studying the role of platelets in cancer biology. From this expertise, we were keen to recognize the numerous venous-, arterial-, microvascular-, and macrovascular thrombotic events and immunologic disorders are caused by severe, acute-respiratory-syndrome coronavirus 2 (SARS-CoV-2) infections. With this offering, we explore the evolutionary connections that place platelets at the center of hemostasis, immunity, and adaptive phylogeny. Coevolutionary changes have also occurred in vertebrate viruses and their vertebrate hosts that reflect their respective evolutionary interactions. As mammals adapted from aquatic to terrestrial life and the heavy blood loss associated with placentalization-based live birth, platelets evolved phylogenetically from thrombocytes toward higher megakaryocyte-blebbing-based production rates and the lack of nuclei. With no nuclei and robust RNA synthesis, this adaptation may have influenced viral replication to become less efficient after virus particles are engulfed. Human platelets express numerous receptors that bind viral particles, which developed from archetypal origins to initiate aggregation and exocytic-release of thrombo-, immuno-, angiogenic-, growth-, and repair-stimulatory granule contents. Whether by direct, evolutionary, selective pressure, or not, these responses may help to contain virus spread, attract immune cells for eradication, and stimulate angiogenesis, growth, and wound repair after viral damage. Because mammalian and marsupial platelets became smaller and more plate-like their biophysical properties improved in function, which facilitated distribution near vessel walls in fluid-shear fields. This adaptation increased the probability that platelets could then interact with and engulf shedding virus particles. Platelets also generate circulating microvesicles that increase membrane surface-area encounters and mark viral targets. In order to match virus-production rates, billions of platelets are generated and turned over per day to continually provide active defenses and adaptation to suppress the spectrum of evolving threats like SARS-CoV-2.

Erythrocyte Sedimentation: Collapse of a High-Volume-Fraction Soft-Particle Gel

The erythrocyte sedimentation rate is one of the oldest medical diagnostic methods whose physical mechanisms remain debatable today. Using both light microscopy and mesoscale cell-level simulations, we show that erythrocytes form a soft-particle gel. Furthermore, the high volume fraction of erythrocytes, their deformability, and weak attraction lead to unusual properties of this gel. A theoretical model for the gravitational collapse is developed, whose predictions are in agreement with detailed macroscopic measurements of the interface velocity.

Blood Viscosity in Subjects With Type 2 Diabetes Mellitus: Roles of Hyperglycemia and Elevated Plasma Fibrinogen

The viscosity of blood is an indicator in the understanding and treatment of disease. An elevated blood viscosity has been demonstrated in patients with Type 2 Diabetes Mellitus (T2DM), which might represent a risk factor for cardiovascular complications. However, the roles of glycated hemoglobin (HbA1c) and plasma fibrinogen levels on the elevated blood viscosity in subjects with T2DM at different chronic glycemic conditions are still not clear. Here, we evaluate the relationship between the blood viscosity and HbA1c as well as plasma fibrinogen levels in patients with T2DM. The experimental data show that the mean values of the T2DM blood viscosity are higher in groups with higher HbA1c levels, but the correlation between the T2DM blood viscosity and the HbA1c level is not obvious. Instead, when we investigate the influence of plasma fibrinogen level on the blood viscosity in T2DM subjects, we find that the T2DM blood viscosity is significantly and positively correlated with the plasma fibrinogen level. Further, to probe the combined effects of multiple factors (including the HbA1c and plasma fibrinogen levels) on the altered blood viscosity in T2DM, we regroup the experimental data based on the T2DM blood viscosity values at both the low and high shear rates, and our results suggest that the influence of the elevated HbA1c level on blood viscosity is quite limited, although it is an important indicator of glycemic control in T2DM patients. Instead, the elevated blood hematocrit, the enhanced red blood cell (RBC) aggregation induced by the increased plasma fibrinogen level, and the reduced RBC deformation play key roles in the determination of blood viscosity in T2DM. Together, these experimental results are helpful in identifying the key determinants for the altered T2DM blood viscosity, which can be used in future studies of the hemorheological disturbances of T2DM patients.

Normalization of Blood Viscosity According to the Hematocrit and the Shear Rate

The rheological properties of blood depend highly on the properties of its red blood cells: concentration, membrane elasticity, and aggregation. These properties affect the viscosity of blood as well as its shear thinning behavior. Using an experimental analysis of the interface advancement of blood in a microchannel, we determine the viscosity of different samples of blood. In this work, we present two methods that successfully normalize the viscosity of blood for a single and for different donors, first according to the concentration of erythrocytes and second according to the shear rate. The proposed methodology is able to predict the health conditions of the blood samples by introducing a non-dimensional coefficient that accounts for the response to shear rate of the different donors blood samples. By means of these normalization methods, we were able to determine the differences between the red blood cells of the samples and define a range where healthy blood samples can be described by a single behavior.

Microfluidics Approach to the Mechanical Properties of Red Blood Cell Membrane and Their Effect on Blood Rheology

In this article, we describe the general features of red blood cell membranes and their effect on blood flow and blood rheology. We first present a basic description of membranes and move forward to red blood cell membranes' characteristics and modeling. We later review the specific properties of red blood cells, presenting recent numerical and experimental microfluidics studies that elucidate the effect of the elastic properties of the red blood cell membrane on blood flow and hemorheology. Finally, we describe specific hemorheological pathologies directly related to the mechanical properties of red blood cells and their effect on microcirculation, reviewing microfluidic applications for the diagnosis and treatment of these diseases.

The effect of rigid cells on blood viscosity: linking rheology and sickle cell anemia

Sickle cell anemia (SCA) is a disease that affects red blood cells (RBCs). Healthy RBCs are highly deformable objects that under flow can penetrate blood capillaries smaller than their typical size. In SCA there is an impaired deformability of some cells, which are much stiffer and with a different shape than healthy cells, and thereby affect regular blood flow. It is known that blood from patients with SCA has a higher viscosity than normal blood. However, it is unclear how the rigidity of cells is related to the viscosity of blood, in part because SCA patients are often treated with transfusions of variable amounts of normal RBCs and only a fraction of cells will be stiff. Here, we report systematic experimental measurements of the viscosity of a suspension varying the fraction of rigid particles within a suspension of healthy cells. We also perform systematic numerical simulations of a similar mixed suspension of soft RBCs, rigid particles, and their hydrodynamic interactions. Our results show that there is a rheological signature within blood viscosity to clearly identify the fraction of rigidified cells among healthy deformable cells down to a 5% volume fraction of rigidified cells. Although aggregation of RBCs is known to affect blood rheology at low shear rates, and our simulations mimic this effect via an adhesion potential, we show that such adhesion, or aggregation, is unlikely to provide a physical rationalization for the viscosity increase observed in the experiments at moderate shear rates due to rigidified cells. Through numerical simulations, we also highlight that most of the viscosity increase of the suspension is due to the rigidity of the particles rather than their sickled or spherical shape. Our results are relevant to better characterize SCA, provide useful insights relevant to rheological consequences of blood transfusions, and, more generally, extend to the rheology of mixed suspensions having particles with different rigidities, as well as offering possibilities for developments in the field of soft material composites.

Red blood cell shape transitions and dynamics in time-dependent capillary flows

The dynamics of single red blood cells (RBCs) determine microvascular blood flow by adapting their shape to the flow conditions in the narrow vessels. In this study, we explore the dynamics and shape transitions of RBCs on the cellular scale under confined and unsteady flow conditions using a combination of microfluidic experiments and numerical simulations. Tracking RBCs in a comoving frame in time-dependent flows reveals that the mean transition time from the symmetric croissant to the off-centered, nonsymmetric slipper shape is significantly faster than the opposite shape transition, which exhibits pronounced cell rotations. Complementary simulations indicate that these dynamics depend on the orientation of the RBC membrane in the channel during the time-dependent flow. Moreover, we show how the tank-treading movement of slipper-shaped RBCs in combination with the narrow channel leads to oscillations of the cell's center of mass. The frequency of these oscillations depends on the cell velocity, the viscosity of the surrounding fluid, and the cytosol viscosity. These results provide a potential framework to identify and study pathological changes in RBC properties.

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.

Recent advances in blood rheology: a review

Due to the potential impact on the diagnosis and treatment of various cardiovascular diseases, work on the rheology of blood has significantly expanded in the last decade, both experimentally and theoretically. Experimentally, blood has been confirmed to demonstrate a variety of non-Newtonian rheological characteristics, including pseudoplasticity, viscoelasticity, and thixotropy. New rheological experiments and the development of more controlled experimental protocols on more extensive, broadly physiologically characterized, human blood samples demonstrate the sensitivity of aspects of hemorheology to several physiological factors. For example, at high shear rates the red blood cells elastically deform, imparting viscoelasticity, while at low shear rates, they form "rouleaux" structures that impart additional, thixotropic behavior. In addition to the advances in experimental methods and validated data sets, significant advances have also been made in both microscopic simulations and macroscopic, continuum, modeling, as well as novel, multiscale approaches. We outline and evaluate the most promising of these recent developments. Although we primarily focus on human blood rheology, we also discuss recent observations on variations observed across some animal species that provide some indication on evolutionary effects.

Lattice Boltzmann simulations on the tumbling to tank-treading transition: effects of membrane viscosity

The tumbling to tank-treading (TB-TT) transition for red blood cells (RBCs) has been widely investigated, with a main focus on the effects of the viscosity ratio [Formula: see text] (i.e., the ratio between the viscosities of the fluids inside and outside the membrane) and the shear rate [Formula: see text] applied to the RBC. However, the membrane viscosity [Formula: see text] plays a major role in a realistic description of RBC dynamics, and only a few works have systematically focused on its effects on the TB-TT transition. In this work, we provide a parametric investigation on the effect of membrane viscosity [Formula: see text] on the TB-TT transition for a single RBC. It is found that, at fixed viscosity ratios [Formula: see text], larger values of [Formula: see text] lead to an increased range of values of capillary number at which the TB-TT transition occurs; moreover, we found that increasing [Formula: see text] or increasing [Formula: see text] results in a qualitatively but not quantitatively similar behaviour. All results are obtained by means of mesoscale numerical simulations based on the lattice Boltzmann models. This article is part of the theme issue 'Progress in mesoscale methods for fluid dynamics simulation'.

Hydrodynamic coupling for particle-based solvent-free membrane models

The great challenge with biological membrane systems is the wide range of scales involved, from nanometers and picoseconds for individual lipids to the micrometers and beyond millisecond for cellular signaling processes. While solvent-free coarse-grained membrane models are convenient for large-scale simulations and promising to provide insight into slow processes involving membranes, these models usually have unrealistic kinetics. One major obstacle is the lack of an equally convenient way of introducing hydrodynamic coupling without significantly increasing the computational cost of the model. To address this, we introduce a framework based on anisotropic Langevin dynamics, for which major in-plane and out-of-plane hydrodynamic effects are modeled via friction and diffusion tensors from analytical or semi-analytical solutions to Stokes hydrodynamic equations. Using this framework, in conjunction with our recently developed membrane model, we obtain accurate dispersion relations for planar membrane patches, both free-standing and in the vicinity of a wall. We briefly discuss how non-equilibrium dynamics is affected by hydrodynamic interactions. We also measure the surface viscosity of the model membrane and discuss the affecting dissipative mechanisms.

The Provocative Roles of Platelets in Liver Disease and Cancer

Both platelets and the liver play important roles in the processes of coagulation and innate immunity. Platelet responses at the site of an injury are rapid; their immediate activation and structural changes minimize the loss of blood. The majority of coagulation proteins are produced by the liver—a multifunctional organ that also plays a critical role in many processes: removal of toxins and metabolism of fats, proteins, carbohydrates, and drugs. Chronic inflammation, trauma, or other causes of irreversible damage to the liver can dysregulate these pathways leading to organ and systemic abnormalities. In some cases, platelet-to-lymphocyte ratios can also be a predictor of disease outcome. An example is cirrhosis, which increases the risk of bleeding and prothrombotic events followed by activation of platelets. Along with a triggered coagulation cascade, the platelets increase the risk of pro-thrombotic events and contribute to cancer progression and metastasis. This progression and the resulting tissue destruction is physiologically comparable to a persistent, chronic wound. Various cancers, including colorectal cancer, have been associated with increased thrombocytosis, platelet activation, platelet-storage granule release, and thrombosis; anti-platelet agents can reduce cancer risk and progression. However, in cancer patients with pre-existing liver disease who are undergoing chemotherapy, the risk of thrombotic events becomes challenging to manage due to their inherent risk for bleeding. Chemotherapy, also known to induce damage to the liver, further increases the frequency of thrombotic events. Depending on individual patient risks, these factors acting together can disrupt the fragile balance between pro- and anti-coagulant processes, heightening liver thrombogenesis, and possibly providing a niche for circulating tumor cells to adhere to—thus promoting both liver metastasis and cancer-cell survival following treatment (that is, with minimal residual disease in the liver).

Loading and relaxation dynamics of a red blood cell

We use mesoscale numerical simulations to investigate the unsteady dynamics of a single red blood cell (RBC) subjected to an external mechanical load. We carry out a detailed comparison between the loading (L) dynamics, following the imposition of the mechanical load on the RBC at rest, and the relaxation (R) dynamics, allowing the RBC to relax to its original shape after the sudden arrest of the mechanical load. Such a comparison is carried out by analyzing the characteristic times of the two corresponding dynamics, i.e., t_L and t_R. When the intensity of the mechanical load is small enough, the two kinds of dynamics are symmetrical (t_L≈t_R) and independent of the typology of mechanical load (intrinsic dynamics); otherwise, in marked contrast, an asymmetry is found, wherein the loading dynamics is typically faster than the relaxation one. This asymmetry manifests itself with non-universal characteristics, e.g., dependency on the applied load and/or on the viscoelastic properties of the RBC membrane. To deepen such a non-universal behaviour, we consider the viscosity of the erythrocyte membrane as a variable parameter and focus on three different typologies of mechanical load (mechanical stretching, shear flow, elongational flow): this allows to clarify how non-universality builds up in terms of the deformation and rotational contributions induced by the mechanical load on the membrane. Finally, we also investigate the effect of the elastic shear modulus on the characteristic times t_L and t_R. Our results provide crucial and quantitative information on the unsteady dynamics of RBC and its membrane response to the imposition/cessation of external mechanical loads.

In-Flow dynamics of an area-difference-energy spring-particle red blood cell model on non-uniform grids

In this paper the area-difference-energy spring-particle (ADE-SP) red blood cell (RBC) structural model developed by Chen and Boyle is coupled with a lattice Boltzmann flux solver to simulate RBC dynamics. The novel ADE-SP model accounts for bending resistance due to the membrane area difference of RBCs while the lattice Boltzmann flux solver offers reduced computational runtimes through GPU parallelisation and enabling the employment of non-uniform meshes. This coupled model is used to simulate RBC dynamics and predictions are compared with existing experimental measurements. The simulations successfully predict tumbling, tank-treading, swinging and intermittent behaviour of an RBC in shear flow, and demonstrate the capability of the model in capturing in-flow RBC behaviours.

In silico biophysics and hemorheology of blood hyperviscosity syndrome

Hyperviscosity syndrome (HVS) is characterized by an increase of the blood viscosity by up to seven times the normal blood viscosity, resulting in disturbances to the circulation in the vasculature system. HVS is commonly associated with an increase of large plasma proteins and abnormalities in the properties of red blood cells, such as cell interactions, cell stiffness, and increased hematocrit. Here, we perform a systematic study of the effect of each biophysical factor on the viscosity of blood by employing the dissipative particle dynamic method. Our in silico platform enables manipulation of each parameter in isolation, providing a unique scheme to quantify and accurately investigate the role of each factor in increasing the blood viscosity. To study the effect of these four factors independently, each factor was elevated more than its values for a healthy blood while the other factors remained constant, and viscosity measurement was performed for different hematocrits and flow rates. Although all four factors were found to increase the overall blood viscosity, these increases were highly dependent on the hematocrit and the flow rates imposed. The effect of cell aggregation and cell concentration on blood viscosity were predominantly observed at low shear rates, in contrast to the more magnified role of cell rigidity and plasma viscosity at high shear rates. Additionally, cell-related factors increase the whole blood viscosity at high hematocrits compared with the relative role of plasma-related factors at lower hematocrits. Our results, mapped onto the flow rates and hematocrits along the circulatory system, provide a correlation to underpinning mechanisms for HVS findings in different blood vessels.

Equivalent Scalar Stress Formulation Taking into Account Non-Resolved Turbulent Scales

PURPOSE

Cardiovascular engineering includes flows with fluid-dynamical stresses as a parameter of interest. Mechanical stresses are high-risk factors for blood damage and can be assessed by computational fluid dynamics. By now, it is not described how to calculate an adequate scalar stress out of turbulent flow regimes when the whole share of turbulence is not resolved by the simulation method and how this impacts the stress calculation.

METHODS

We conducted direct numerical simulations (DNS) of test cases (a turbulent channel flow and the FDA nozzle) in order to access all scales of flow movement. After validation of both DNS with literature und experimental data using magnetic resonance imaging, the mechanical stress is calculated as a baseline. Afterwards, same flows are calculated using state-of-the-art turbulence models. The stresses are computed for every result using our definition of an equivalent scalar stress, which includes the influence from respective turbulence model, by using the parameter dissipation. Afterwards, the results are compared with the baseline data.

RESULTS

The results show a good agreement regarding the computed stress. Even when no turbulence is resolved by the simulation method, the results agree well with DNS data. When the influence of non-resolved motion is neglected in the stress calculation, it is underpredicted in all cases.

CONCLUSION

With the used scalar stress formulation, it is possible to include information about the turbulence of the flow into the mechanical stress calculation even when the used simulation method does not resolve any turbulence.

Effect of cytosol viscosity on the flow behavior of red blood cell suspensions in microvessels

OBJECTIVE

The flow behavior of blood is strongly affected by red blood cell (RBC) properties, such as the viscosity ratio C between cytosol and suspending medium, which can significantly be altered in several pathologies (e.g. sickle-cell disease, malaria). The main objective of this study is to understand the effect of C on macroscopic blood flow properties such as flow resistance in microvessels, and to link it to the deformation and dynamics of single RBCs.

METHODS

We employ mesoscopic hydrodynamic simulations to investigate flow properties of RBC suspensions with different cytosol viscosities for various flow conditions in cylindrical microchannels.

RESULTS

Starting from a dispersed cell configuration which approximates RBC dispersion at vessel bifurcations in the microvasculature, we find that the flow convergence and development of RBC-free layer (RBC-FL) depend only weakly on C, and require a convergence length in the range of 25D-50D, where D is channel diameter. In vessels with D ≤ 20 μ m , the final resistance of developed flow is nearly the same for C = 5 and C = 1, while for D = 40 μ m , the flow resistance for C = 5 is about 10% larger than for C = 1. The similarities and differences in flow resistance can be explained by viscosity-dependent RBC-FL thicknesses, which are associated with the viscosity-dependent dynamics of single RBCs.

CONCLUSIONS

The weak effect on the flow resistance and RBC-FL explains why RBCs can contain a high concentration of hemoglobin for efficient oxygen delivery, without a pronounced increase in the flow resistance. Furthermore, our results suggest that significant alterations in microvascular flow in various pathologies are likely not due to mere changes in cytosolic viscosity.

Computational fluid dynamics: a primer for congenital heart disease clinicians

Computational fluid dynamics has become an important tool for studying blood flow dynamics. As an in-silico collection of methods, computational fluid dynamics is noninvasive and provides numerical values for the most important parameters of blood flow, such as velocity and pressure that are crucial in hemodynamic studies. In this primer, we briefly explain the basic theory and workflow of the two most commonly applied computational fluid dynamics techniques used in the congenital heart disease literature: the finite element method and the finite volume method. We define important terminology and include specific examples of how using these methods can answer important clinical questions in congenital cardiac surgery planning and perioperative patient management.

On the effects of membrane viscosity on transient red blood cell dynamics

Computational Fluid Dynamics (CFD) is currently used to design and improve the hydraulic properties of biomedical devices, wherein the large scale blood circulation needs to be simulated by accounting for the mechanical response of red blood cells (RBCs) at the mesoscale. In many practical instances, biomedical devices work on time-scales comparable to the intrinsic relaxation time of RBCs: thus, a systematic understanding of the time-dependent response of erythrocyte membranes is crucial for the effective design of such devices. So far, this information has been deduced from experimental data, which do not necessarily adapt to the broad variety of fluid dynamic conditions that can be encountered in practice. This work explores the novel possibility of studying the time-dependent response of an erythrocyte membrane to external mechanical loads via mesoscale numerical simulations, with a primary focus on the detailed characterisation of the RBC relaxation time tc following the arrest of the external mechanical load. The adopted mesoscale model exploits a hybrid Immersed Boundary-Lattice Boltzmann Method (IB-LBM), coupled with the Standard Linear Solid (SLS) model to account for the RBC membrane viscosity. We underscore the key importance of the 2D membrane viscosity μm to correctly reproduce the relaxation time of the RBC membrane. A detailed assessment of the dependencies on the typology and strength of the applied mechanical loads is also provided. Overall, our findings open interesting future perspectives for the study of the non-linear response of RBCs immersed in time-dependent strain fields.

Stochastic bond dynamics facilitates alignment of malaria parasite at erythrocyte membrane upon invasion

Malaria parasites invade healthy red blood cells (RBCs) during the blood stage of the disease. Even though parasites initially adhere to RBCs with a random orientation, they need to align their apex toward the membrane in order to start the invasion process. Using hydrodynamic simulations of a RBC and parasite, where both interact through discrete stochastic bonds, we show that parasite alignment is governed by the combination of RBC membrane deformability and dynamics of adhesion bonds. The stochastic nature of bond-based interactions facilitates a diffusive-like re-orientation of the parasite at the RBC membrane, while RBC deformation aids in the establishment of apex-membrane contact through partial parasite wrapping by the membrane. This bond-based model for parasite adhesion quantitatively captures alignment times measured experimentally and demonstrates that alignment times increase drastically with increasing rigidity of the RBC membrane. Our results suggest that the alignment process is mediated simply by passive parasite adhesion.

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.

Multiple Particle Manipulation under Dielectrophoresis Effect: Modeling and Experiments

The dissipative particle dynamics (DPD) technique was employed to design multiple microfluidic devices for investigating the motion of bioparticles at low Reynolds numbers. A DPD in-house FORTRAN code was developed to simulate the trajectories of two microparticles in the presence of hydrodynamic and transverse deflecting force fields via considering interparticle interaction forces. The particle-particle interactions were described by using a simplified version of the Morse potential. The transverse deflecting force considered in this microfluidic application was the dielectrophoresis (DEP) force. Multiple microfluidic devices with different configurations of microelectrodes were numerically designed to investigate the dielectrophoretic behavior of bioparticles for their trajectories and the focusing of bioparticles into a single stream in the middle of the microchannel. The DPD simulation results were verified and validated against previously reported numerical and experimental works in the literature. The computationally designed microdevices were fabricated by employing standard lithographic techniques, and experiments were conducted via taking red blood cells as the representative bioparticles. The experimental results for the trajectories and focusing showed good agreement with the numerical results.

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.

Dissipative particle dynamics with energy conservation: Isoenergetic integration and transport properties

Simulations of nano- to micro-meter scale fluidic systems under thermal gradients require consistent mesoscopic methods accounting for both hydrodynamic interactions and proper transport of energy. One such method is dissipative particle dynamics with energy conservation (DPDE), which has been used for various fluid systems with non-uniform temperature distributions. We propose an easily parallelizable modification of the velocity-Verlet algorithm based on local energy redistribution for each DPDE particle such that the total energy in a simulated system is conserved up to machine precision. Furthermore, transport properties of a DPDE fluid are analyzed in detail. In particular, an analytical approximation for the thermal conductivity coefficient is derived, which allows its a priori estimation for a given parameter set. Finally, we provide approximate expressions for the dimensionless Prandtl and Schmidt numbers, which characterize fluid transport properties and can be adjusted independently by a proper selection of model parameters. In conclusion, our results strengthen the DPDE method as a very robust approach for the investigation of mesoscopic systems with temperature inhomogeneities.

Effect of Gravity on the Scale of Compliant Shells

Thin shells are found across scales ranging from biological blood cells to engineered large-span roof structures. The engineering design of thin shells used as mechanisms has occasionally been inspired by biomimetic concept generators. The research goal of this paper is to establish the physical limits of scalability of shells. Sixty-four instances of shells across length scales have been organized into five categories: engineering stiff and compliant, plant compliant, avian egg stiff, and micro-scale compliant shells. Based on their thickness and characteristic dimensions, the mechanical behavior of these 64 shells can be characterized as 3D solids, thick or thin shells, or membranes. Two non-dimensional indicators, the Föppl–von Kármán number and a novel indicator, namely the gravity impact number, are adopted to establish the scalability limits of these five categories. The results show that these shells exhibit similar mechanical behavior across scales. As a result, micro-scale shell geometries found in biology, can be upscaled to engineered shell geometries. However, as the characteristic shell dimension increases, gravity (and its associated loading) becomes a hindrance to the adoption of thin shells as compliant mechanisms at the larger scales-the physical limit of compliance in the scaling of thin shells is found to be around 0.1 m.

Importance of Erythrocyte Deformability for the Alignment of Malaria Parasite upon Invasion

Invasion of erythrocytes by merozoites is an essential step for the survival and progression of malaria parasites. To invade red blood cells (RBCs), apicomplexan parasites have to adhere with their apex to the RBC membrane. This necessary apex-membrane contact (or alignment) is not immediately established because the orientation of a free merozoite with respect to the RBC membrane is random when an adhesion contact first occurs. Therefore, it has been suggested that after the initial adhesion, merozoites facilitate their proper alignment by inducing considerable membrane deformations, frequently observed before the invasion process. This proposition is based on a positive correlation between RBC membrane deformation and successful invasion; however, the role of RBC mechanics and its deformation in the alignment process remains elusive. Using a mechanically realistic model of a deformable RBC, we investigate numerically the importance of RBC deformability for merozoite alignment. Adhesion between the parasite and RBC membrane is modeled by an attractive potential that might be inhomogeneous, mimicking possible adhesion gradients at the surface of a parasite. Our results show that RBC membrane deformations are crucial for successful merozoite alignment and require interaction strengths comparable to adhesion forces measured experimentally. Adhesion gradients along the parasite body further improve its alignment. Finally, an increased membrane rigidity is found to result in poor merozoite alignment, which can be a possible reason for a reduction in the invasion susceptibility of RBCs in several blood diseases associated with membrane stiffening.

Near-Newtonian Blood Behavior - Is It Good to Be a Camel?

From a certain level of exercise-intensity onward, hematocrit increases in horses, which brings more oxygen carriers into the bloodstream. Camels, however, when used in competitive racing could be even in need of iron supplementation and blood transfusions due to a severe reduction of their available hematocrit compared to their resting hematocrit. Since the extrinsic and intrinsic mechanical properties of camel erythrocytes (RBC) are so different compared to RBCs of other mammals, the question arises whether this observation might be a response to endurance exercise aiming at keeping the RBC count low. Rheometry indicated dromedary camel blood to behave almost Newtonian, which is unique amongst mammals. Shear thinning did increase with the hematocrit, but remained marginal compared to horses. As a result, camel whole blood viscosity (WBV) exceeded horse WBV at high shear rates, an effect, which was significantly augmented when the packed cell volume (PCV) was increased. Therefore, in camels any infusion of RBCs into the bloodstream can increase the cardiac work and the energy input into the endothelium more effectively, which should generate vascular remodeling in the long term. Yielding, however, was completely absent in camel blood, confirming low cohesion between its components at quasi-static flow. Camel blood remained a viscous liquid without a threshold even at unphysiologically high PCVs. This can help to washout lactate when camels start to dehydrate and might contribute to the sustained working ability of these animals. The subtle pseudoplastic behavior and the high viscosity contrast across the RBC membrane point to weak coupling between blood flow and red cell behavior. We predict that RBCs flow as separate entities and can show various types of motion, which can lead to friction instead of being collectively aligned to the flow direction. In comparison to horses, this behavior will become relevant at higher RBC counts in front of flow obstacles and possibly cause vascular remodeling if the PCV rises during strenuous exercise, a matter that should be avoided.

Deformation and dynamics of erythrocytes govern their traversal through microfluidic devices with a deterministic lateral displacement architecture

Deterministic lateral displacement (DLD) microfluidic devices promise versatile and precise processing of biological samples. However, this prospect has been realized so far only for rigid spherical particles and remains limited for biological cells due to the complexity of cell dynamics and deformation in microfluidic flow. We employ mesoscopic hydrodynamics simulations of red blood cells (RBCs) in DLD devices with circular posts to better understand the interplay between cell behavior in complex microfluidic flow and sorting capabilities of such devices. We construct a mode diagram of RBC behavior (e.g., displacement, zig-zagging, and intermediate modes) and identify several regimes of RBC dynamics (e.g., tumbling, tank-treading, and trilobe motion). Furthermore, we link the complex interaction dynamics of RBCs with the post to their effective cell size and discuss relevant physical mechanisms governing the dynamic cell states. In conclusion, sorting of RBCs in DLD devices based on their shear elasticity is, in general, possible but requires fine-tuning of flow conditions to targeted mechanical properties of the RBCs.

Bioactive lipid metabolism in platelet "first responder" and cancer biology

Platelets can serve as "first responders" in cancer and metastasis. This is partly due to bioactive lipid metabolism that drives both platelet and cancer biology. The two primary eicosanoid metabolites that maintain platelet rapid response homeostasis are prostacyclin made by endothelial cells that inhibits platelet function, which is counterbalanced by thromboxane produced by platelets during activation, aggregation, and platelet recruitment. Both of these arachidonic acid metabolites are inherently unstable due to their chemical structure. Tumor cells by contrast predominantly make more chemically stable prostaglandin E₂, which is the primary bioactive lipid associated with inflammation and oncogenesis. Pharmacological, clinical, and epidemiologic studies demonstrate that non-steroidal anti-inflammatory drugs (NSAIDs), which target cyclooxygenases, can help prevent cancer. Much of the molecular and biological impact of these drugs is generally accepted in the field. Cyclooxygenases catalyze the rate-limiting production of substrate used by all synthase molecules, including those that produce prostaglandins along with prostacyclin and thromboxane. Additional eicosanoid metabolites include lipoxygenases, leukotrienes, and resolvins that can also influence platelets, inflammation, and carcinogenesis. Our knowledge base and technology are now progressing toward identifying newer molecular and cellular interactions that are leading to revealing additional targets. This review endeavors to summarize new developments in the field.

L-cysteine improves blood fluidity impaired by acetaldehyde: In vitro evaluation

Blood fluidity is reportedly influenced by the volume and function of blood cells and plasma and is a predictor of primary cardiovascular events in patients with traditional cardiovascular risk factors. Heavy alcohol consumption was shown to be associated with a higher risk for cardiovascular diseases. Acetaldehyde (ACD), an oxidizing substance formed from ethanol, reportedly stimulates monocyte adhesion, causes abnormalities in the red blood cell (RBC) membrane, and decreases RBC deformability. In addition, it was reported that blood ACD levels are reduced in mice pretreated with L-cysteine. However, there are no studies on the effect of ACD and/or L-cysteine on blood fluidity. In the present study, we evaluated whether ACD impairs blood fluidity. In addition, the effect of L-cysteine on blood fluidity impaired by ACD was examined. Blood samples were obtained from 10 healthy, non-smoking, male volunteers (age: 23.4 ± 1.2 years, body mass index: 21.8 ± 2.6 kg/m²). ACD or ACD and L-cysteine were added to the blood samples before each experiment. We measured the passage time of 100 μL blood and RBC suspension using Kikuchi’s microchannel method. Percentage of microchannel obstruction and the number of adherent white blood cells (WBCs) on microchannel terrace were counted. The blood passage time, percentage of microchannel obstruction, and numbers of adherent WBCs on the microchannel terrace increased after adding ACD in a concentration-dependent manner, whereas they decreased after adding ACD and L-cysteine in a L-cysteine concentration-dependent manner. No significant effects were observed in passage time for 100 μL RBC suspension after adding ACD and L-cysteine. This study suggested that blood fluidity impaired by ACD might improve after adding L-cysteine.

Numerical design of a microfluidic chip for probing mechanical properties of cells

Microfluidic chips have been widely used to probe the mechanical properties of cells, which are recognized as a promising label-free biomarker for some diseases. In our previous work (Ye et al., 2018), we have studied the relationships between the transit time and the mechanical properties of a cell flowing through a microchannel with a single constriction, which potentially forms a basis for a microfluidic chip to measure cell's mechanical properties. Here, we investigate this microfluidic chip design and examine its potential in performances. We first develop the simultaneous dependence of the transit time on both the shear and bending moduli of a cell, and then examine the chip sensitivity with respect to the cell mechanical properties while serializing a single constriction along the flow direction. After that, we study the effect of the flow velocity on the transit time, and also test the chip's ability to identify heterogeneous cells with different mechanical properties. The results show that the microfluidic chip designed is capable of identifying heterogeneous cells, even when only one unhealthy cell is included. The serialization of chip can greatly increase the chip sensitivity with respect to the mechanical properties of cells. The flow with a higher velocity helps in not only promoting the chip throughput, but also in providing more accurate transit time measurements, because the cell prefers a symmetric deformation under a high velocity.