Tank-treading dynamics of red blood cells in shear flow: On the membrane viscosity rheology

In this article, extensive three-dimensional simulations are conducted for tank-treading (TT) red blood cells (RBCs) in shear flow with different cell viscous properties and flow conditions. Apart from recent numerical studies on TT RBCs, this research considers the uncertainty in cytoplasm viscosity, covers a more complete range of shear flow situations of available experiments, and examines the TT behaviors in more details. Key TT characteristics, including the rotation frequency, deformation index, and inclination angle, are compared with available experimental results of similar shear flow conditions. Fairly good simulation-experiment agreements for these parameters can be obtained by adjusting the membrane viscosity values; however, different rheological relationships between the membrane viscosity and the flow shear rate are noted for these comparisons: shear thinning from the TT frequency, Newtonian from the inclination angle, and shear thickening from the cell deformation. Previous studies claimed a shear-thinning membrane viscosity model based on the TT frequency results; however, such a conclusion seems premature from our results and more carefully designed and better controlled investigations are required for the RBC membrane rheology. In addition, our simulation results reveal complicate RBC TT features and such information could be helpful for a better understanding of in vivo and in vitro RBC dynamics.     

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.     

Red blood cell membrane water permeability increases with length of ex vivo storage

Water transport across the red blood cell (RBC) membrane is an essential cell function that needs to be preserved during ex vivo storage. Progressive biochemical depletion during storage can result in significant conformational and compositional changes to the membrane. Characterizing the changes to RBC water permeability can help in evaluating the quality of stored blood products and aid in the development of improved methods for the cryopreservation of red blood cells. This study aimed to characterize the water permeability (L_p), osmotically inactive fraction (b), and Arrhenius activation energy (E_a) at defined storage time-points throughout storage and to correlate the observed results with other in vitro RBC quality parameters. RBCs were collected from age- and sex-matched blood donors. A stopped flow spectrophotometer was used to determine L_p and b by monitoring changes in hemoglobin autofluorescence when RBCs were exposed to anisotonic solutions. Experimental values of L_p were characterized at three different temperatures (4, 20 and 37 °C) to determine the E_a. Results showed that L_p, b, and E_a of stored RBCs significantly increase by day 21 of storage. Degradation of the RBC membrane with length of storage was seen as an increase in hemolysis and supernatant potassium, and a decrease in deformability, mean corpuscular hemoglobin concentration and supernatant sodium. RBC osmotic characteristics were shown to change with storage and correlate with changes in RBC membrane quality metrics. Monitoring water parameters is a predictor of membrane damage and loss of membrane integrity in ex vivo stored RBCs.     

Computational modeling of biomechanics and biorheology of heated red blood cells

Because of their compromised deformability, heat denatured erythrocytes have been used as labeled probes to visualize spleen tissue or to assess the ability of the spleen to retain stiff red blood cells (RBCs) for over three decades, e.g., see Looareesuwan et al. N. Engl. J. Med. (1987). Despite their good accessibility, it is still an open question how heated RBCs compare to certain diseased RBCs in terms of their biomechanical and biorheological responses, which may undermine their effective usage and even lead to misleading experimental observations. To help answering this question, we perform a systematic computational study of the hemorheological properties of heated RBCs with several physiologically relevant static and hemodynamic settings, including optical-tweezers test, relaxation of prestretched RBCs, RBC traversal through a capillary-like channel and a spleen-like slit, and a viscometric rheology test. We show that our in silico RBC models agree well with existing experiments. Moreover, under static tests, heated RBCs exhibit deformability deterioration comparable to certain disease-impaired RBCs such as those in malaria. For RBC traversal under confinement (through microchannel or slit), heated RBCs show prolonged transit time or retention depending on the level of confinement and heating procedure, suggesting that carefully heat-treated RBCs may be useful for studying splenic- or vaso-occlusion in vascular pathologies. For the rheology test, we expand the existing bulk viscosity data of heated RBCs to a wider range of shear rates (1–1000 s⁻¹) to represent most pathophysiological conditions in macro- or microcirculation. Although heated RBC suspension shows elevated viscosity comparable to certain diseased RBC suspensions under relatively high shear rates (100–1000 s⁻¹), they underestimate the elevated viscosity (e.g., in sickle cell anemia) at low shear rates (<10 s⁻¹). Our work provides mechanistic rationale for selective usage of heated RBC as a potentially useful model for studying the abnormal traversal dynamics and hemorheology in certain blood disorders.     

Improved Model of Fluorescence Recovery Expands the Application of Multiphoton Fluorescence Recovery after Photobleaching in Vivo

Multiphoton fluorescence recovery after photobleaching is a well-established microscopy technique used to measure the diffusion of macromolecules in biological systems. We have developed an improved model of the fluorescence recovery that includes the effects of convective flows within a system. We demonstrate the validity of this two-component diffusion-convection model through in vitro experimentation in systems with known diffusion coefficients and known flow speeds, and show that the diffusion-convection model broadens the applicability of the multiphoton fluorescence recovery after photobleaching technique by enabling accurate determination of the diffusion coefficient, even when significant flows are present. Additionally, we find that this model allows for simultaneous measurement of the flow speed in certain regimes. Finally, we demonstrate the effectiveness of the diffusion-convection model in vivo by measuring the diffusion coefficient and flow speed within tumor vessels of 4T1 murine mammary adenocarcinomas implanted in the dorsal skinfold chamber.     

Role of Macromolecular Depletion in Red Blood Cell Adhesion

The adhesion of red blood cells (RBCs) to cells or surfaces is of current basic science and clinical interest yet the details governing this process are still being explored. In this study, the effects of nonadsorbing polymers on the adhesion of RBC to albumin-coated glass have been investigated employing interference reflection microscopy. Our experimental results indicate that adhesion can be induced in the presence of dextran with a molecular mass ≥70 kDa and that the induced forces are strong enough to significantly suppress membrane undulations. The overall dependence of the adhesion energies on the polymer concentration is consistent with the assumption that macromolecular depletion induces this attractive interaction. In conclusion, our results indicate that depletion interaction might play a significant role in RBC adhesion via initiating close contacts, and thus suggest the importance of depletion forces for RBC interactions and its relevance to a wide variety of in vitro and in vivo cell-cell and cell-surface interactions.     

Lung Ischemia: A Model for Endothelial Mechanotransduction

Endothelial cells in vivo are constantly exposed to shear associated with blood flow and altered shear stress elicits cellular responses (mechanotransduction). This review describes the role of shear sensors and signal transducers in these events. The major focus is the response to removal of shear as occurs when blood flow is compromised (i.e., ischemia). Pulmonary ischemia studied with the isolated murine lung or flow adapted pulmonary microvascular endothelial cells in vitro results in endothelial generation of reactive oxygen species (ROS) and NO. The response requires caveolae and is initiated by endothelial cell depolarization via K_ATP channel closure followed by activation of NADPH oxidase (NOX2) and NO synthase (eNOS), signaling through MAP kinases, and endothelial cell proliferation. These physiological mediators can promote vasodilation and angiogenesis as compensation for decreased tissue perfusion.     

Quantification of Blood Flow and Topology in Developing Vascular Networks

Since fluid dynamics plays a critical role in vascular remodeling, quantification of the hemodynamics is crucial to gain more insight into this complex process. Better understanding of vascular development can improve prediction of the process, and may eventually even be used to influence the vascular structure. In this study, a methodology to quantify hemodynamics and network structure of developing vascular networks is described. The hemodynamic parameters and topology are derived from detailed local blood flow velocities, obtained by in vivo micro-PIV measurements. The use of such detailed flow measurements is shown to be essential, as blood vessels with a similar diameter can have a large variation in flow rate. Measurements are performed in the yolk sacs of seven chicken embryos at two developmental stages between HH 13+ and 17+. A large range of flow velocities (1 µm/s to 1 mm/s) is measured in blood vessels with diameters in the range of 25–500 µm. The quality of the data sets is investigated by verifying the flow balances in the branching points. This shows that the quality of the data sets of the seven embryos is comparable for all stages observed, and the data is suitable for further analysis with known accuracy. When comparing two subsequently characterized networks of the same embryo, vascular remodeling is observed in all seven networks. However, the character of remodeling in the seven embryos differs and can be non-intuitive, which confirms the necessity of quantification. To illustrate the potential of the data, we present a preliminary quantitative study of key network topology parameters and we compare these with theoretical design rules.     

Loss of strumpellin in the melanocytic lineage impairs the WASH Complex but does not affect coat colour

The five‐subunit WASH complex generates actin networks that participate in endocytic trafficking, migration and invasion in various cell types. Loss of one of the two subunits WASH or strumpellin in mice is lethal, but little is known about their role in mammals in vivo. We explored the role of strumpellin, which has previously been linked to hereditary spastic paraplegia, in the mouse melanocytic lineage. Strumpellin knockout in melanocytes revealed abnormal endocytic vesicle morphology but no impairment of migration in vitro or in vivo and no change in coat colour. Unexpectedly, WASH and filamentous actin could still localize to vesicles in the absence of strumpellin, although the shape and size of vesicles was altered. Blue native PAGE revealed the presence of two distinct WASH complexes, even in strumpellin knockout cells, revealing that the WASH complex can assemble and localize to endocytic compartments in cells in the absence of strumpellin.     

Bioenergetic modulators hamper cancer cell viability and enhance response to chemotherapy

Gliomas are characterized by a marked glycolytic metabolism with a consequent production of massive amounts of lactate, even in the presence of normal levels of oxygen, associated to increased invasion capacity and to higher resistance to conventional treatment. This work aimed to understand how the metabolic modulation can influence tumour aggressive features and its potential to be used as complementary therapy. We assessed the effect of bioenergetic modulators (BMs) targeting different metabolic pathways in glioma cell characteristics. The in vivo effect of BMs was evaluated using the chicken chorioallantoic membrane model. Additionally, the effect of pre‐treatment with BMs in the response to the antitumour drug temozolomide (TMZ) was analysed in vitro. Cell treatment with the BMs induced a decrease in cell viability and in migratory/invasion abilities, as well as modifications in metabolic parameters (glucose, lactate and ATP) and increased the cytotoxicity of the conventional drug TMZ. Furthermore, all BMs decreased the tumour growth and the number of blood vessels in an in vivo model. Our results demonstrate that metabolic modulation has the potential to be used as therapy to decrease the aggressiveness of the tumours or to be combined with conventional drugs used in glioma treatment.     

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

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

Calibration of Optical Tweezers for In Vivo Force Measurements: How do Different Approaches Compare?

There is significant interest in quantifying force production inside cells, but since conditions in vivo are less well controlled than those in vitro, in vivo measurements are challenging. In particular, the in vivo environment may vary locally as far as its optical properties, and the organelles manipulated by the optical trap frequently vary in size and shape. Several methods have been proposed to overcome these difficulties. We evaluate the relative merits of these methods and directly compare two of them, a refractive index matching method, and a light-momentum-change method. Since in vivo forces are frequently relatively high (e.g., can exceed 15 pN for lipid droplets), a high-power laser is employed. We discover that this high-powered trap induces local temperature changes, and we develop an approach to compensate for uncertainties in the magnitude of applied force due to such temperature variations.     

Simulation of Platelet,Thrombus and Erythrocyte Hydrodynamic Interactions in a 3D Arteriole with In Vivo Comparison

Cylindrical blood vessels, ellipsoid platelets and biconcave-shaped deformable erythrocytes (RBCs) are important participants in hemostasis and thrombosis. However, due to the challenge of combining these components in simulation tools, few simulation studies have included all of them in realistic three-dimensional models. In the present study, we apply a recently developed simulation model to incorporate these components and analyze the flow in a thrombotic tubular arteriole, particularly the detailed hydrodynamic interactions between the thrombus shape, RBCs and platelets. It was found that at certain azimuth positions, the velocity drops in the proximity of both the upstream and downstream edge of the thrombus, which is accompanied by a rapid velocity increase in the narrowed region. The RBCs alter the flow profiles significantly from the typical low Reynolds (Re) number flow, and also enhance the deposition of free flowing platelets onto the thrombus. By evaluating the platelet-thrombus interaction and platelet-RBC interaction together, several mechanisms of platelet deposition augmentation are identified. With in vivo data comparison, our model illustrates the potential of future thrombosis studies that incorporate detailed receptor-ligand adhesion modules.     

The circulatory dynamics of human red blood cell homeostasis: Oxy-deoxy and PIEZO1-triggered changes

The vital function of red blood cells (RBCs) is to mediate the transport of oxygen from lungs to tissues and of CO₂ from tissues to lungs. The gas exchanges occur during capillary transits within fractions of a second. Each oxygenation-deoxygenation and deoxygenation-reoxygenation transition on hemoglobin triggers sharp changes in RBC pH, leading to downstream changes in ion fluxes, membrane potential, and cell volume. The dynamics of these changes during the variable periods between capillary transits in vivo remains a mystery inaccessible to study by current methodologies, a knowledge gap on a fundamental physiological process that is the focus of the present study. The use of a computational model of human RBC homeostasis of tested accreditation enabled a detailed investigation of the expected RBC changes during intercapillary transits, with results advancing novel insights and predictions. The predicted rates of relative RBC volume change on oxygenation-deoxygenation (oxy-deoxy) and deoxygenation-reoxygenation transitions were about 1.5%/min and ?0.9%/min, respectively, far too slow to allow the cells to reach steady states in the intervals between capillary transits. The amplitude of the oxy-deoxy-reoxygenation volume fluctuations varied in proportion with the duration of the intercapillary transit intervals. Upon capillary entry, oxy-deoxy-induced changes occur concurrently with deformation-induced PIEZO1 channel activation, both processes affecting cell pH, membrane potential, and cell volume during intertransit periods. The model showed that the effects were strictly additive as expected from processes operating independently on the cell’s homeostatic fabric. Analysis of the mechanisms behind these predictions revealed, for the first time, the complex interactions between oxy-deoxy and ion transport processes that ensure the long-term homeostatic stability of RBCs for optimal gas transport in physiological conditions and how these may become altered in diseased states. Possible designs of microfluidic devices to test the model predictions are discussed.     

Visualization study of motion and deformation of red blood cells in a microchannel with straight,divergent and convergent sections

The size of red blood cells (RBC) is on the same order as the diameter of microvascular vessels. Therefore, blood should be regarded as a two-phase flow system of RBCs suspended in plasma rather than a continuous medium of microcirculation. It is of great physiological and pathological significance to investigate the effects of deformation and aggregation of RBCs on microcirculation. In this study, a visualization experiment was conducted to study the microcirculatory behavior of RBCs in suspension. Motion and deformation of RBCs in a microfluidic chip with straight, divergent, and convergent microchannel sections have been captured by microscope and high-speed camera. Meanwhile, deformation and movement of RBCs were investigated under different viscosity, hematocrit, and flow rate in this system. For low velocity and viscosity, RBCs behaved in their normal biconcave disc shape and their motion was found as a flipping motion: they not only deformed their shapes along the flow direction, but also rolled and rotated themselves. RBCs were also found to aggregate, forming rouleaux at very low flow rate and viscosity. However, for high velocity and viscosity, RBCs deformed obviously under the shear stress. They elongated along the flow direction and performed a tank-treading motion.     

Atomic force microscopy study of red blood cell membrane nanostructure during oxidation‐reduction processes

The morphology and functional state of red blood cells (RBCs) mainly depends on the configuration of the spectrin network, which can be broken under the influence of intoxication because of oxidation processes in the cells. Measurement of these processes is a complex problem. The most suitable and prospective method that resolves this problem is atomic force microscopy (AFM). We used AFM to study the changes in the spectrin matrix and RBC morphology during oxidation processes caused by ultraviolet (UV) irradiation in RBC suspension. The number of discocytes decreased from 98% (in control) to 12%. We obtained AFM images of the spectrin matrix in RBC ghosts. Atomic force microscopy allows for the direct observation and quantitative measurement of the disturbances in the structure of the spectrin matrix during oxidation processes in RBCs. The typical section size of the spectrin network changed from approximately 80 to 200 nm (in control) to 600 nm and even to 1000 nm after UV irradiation. An AFM study showed that incubation of RBCs with Cytoflavin® after UV irradiation preserved the forms of RBCs almost at control levels; 89% of the cells remained as discocytes. To quantify the intensity of the oxidation‐reduction processes, the percentage of haemoglobin derivatives was measured. The content of methaemoglobin varied in the range of 1% to 70% during the experiments. These evidence‐based studies are important for the fundamental research of interactions during redox processes in RBCs at the molecular level.     

Role of erythrocytes in leukocyte-endothelial interactions: mathematical model and experimental validation.

The binding of circulating cells to the vascular wall is a central process in inflammation, metastasis, and therapeutic cell delivery. Previous in vitro studies have identified the adhesion molecules on various circulating cells and the endothelium that govern the process under static conditions. Other studies have attempted to simulate in vivo conditions by subjecting adherent cells to shear stress as they interact with the endothelial cells in vitro. These experiments are generally performed with the cells suspended in Newtonian solutions. However, in vivo conditions are more complex because of the non-Newtonian flow of blood, which is a suspension consisting of 20-40% erythrocytes by volume. The forces imparted by the erythrocytes in the flow can contribute to the process of cell adhesion. A number of experimental and theoretical studies have suggested that the rheology of blood can influence the binding of circulating leukocytes by increasing the normal and axial forces on leukocytes or the frequency of their collision with the vessel wall, but there have been no systematic investigations of these phenomena to date. The present study quantifies the contribution of red blood cells (RBCs) in cell capture and adhesion to endothelial monolayers using a combination of mathematical modeling and in vitro studies. Mathematical modeling of the flow experiments suggested a physical mechanism involving RBC-induced leukocyte dispersion and/or increased normal adhesive contact. Flow chamber studies performed with and without RBCs in the suspending medium showed increases in wall collision and binding frequencies, and a decrease in rolling velocity in the presence of erythrocytes. Increased fluid viscosity alone did not influence the binding frequency, and the differences could not be attributed to large near-wall excesses of the lymphocytes. The results indicate that RBCs aid in the transport and initial engagement of lymphocytes to the vascular wall, modifying the existing paradigm for immune cell surveillance of the vascular endothelium by adding the erythrocyte as an essential contributor to this process.     

Influence of microflow on hepatic sinusoid blood flow and red blood cell deformation

Hepatic sinusoids present complex anatomical structures such as the endothelial sieve pores and the Disse space, which govern the microscopic blood flow in the sinusoids and are associated with structural variations in liver fibrosis and cirrhosis. However, the contributions of the permeability of endothelial and collagen layers and the roughness of hepatocyte microvilli to the features of this microflow remain largely unknown. Here, an immersed boundary method coupled with a lattice Boltzmann method was adopted in an in vitro hepatic sinusoidal model, and flow field and erythrocyte deformation analyses were conducted by introducing three new source terms including permeability of the endothelial layer, resistance of hepatocyte microvilli and collagen layers, and deformation of red blood cells (RBCs). Numerical calculations indicated that alterations in endothelial permeability could significantly affect the flow velocity and flow rate distributions in hepatic sinusoids. Interestingly, a biphasic regulating pattern of shear stress occurred simultaneously on the surface of hepatocytes and the lower side of endothelium, i.e., the shear stress increased with increased thickness of hepatocyte microvilli and collagen layer when the endothelial permeability was high but decreased with the increase of the thickness at low endothelial permeability. Additionally, this specified microflow manipulates typical RBC deformation inside the sinusoid, yielding one-third of the variation of deformable index with varied endothelial permeability. These simulations not only are consistent with experimental measurements using in vitro liver sinusoidal chip but also elaborate the contributions of endothelial and collagen layer permeability and wall roughness. Thus, our results provide a basis for further characterizing this microflow and understanding its effects on cellular migration and deformation in the hepatic sinusoids.     

In situ solid-state NMR study of antimicrobial peptide interactions with erythrocyte membranes

Antimicrobial peptides are promising therapeutic agents to mitigate the global rise of antibiotic resistance. They generally act by perturbing the bacterial cell membrane and are thus less likely to induce resistance. Because they are membrane-active molecules, it is critical to verify and understand their potential action toward eukaryotic cells to help design effective and safe drugs. In this work, we studied the interaction of two antimicrobial peptides, aurein 1.2 and caerin 1.1, with red blood cell (RBC) membranes using in situ ³¹P and ²H solid-state NMR (SS-NMR). We established a protocol to integrate up to 25% of deuterated fatty acids in the membranes of ghosts, which are obtained when hemoglobin is removed from RBCs. Fatty acid incorporation and the integrity of the lipid bilayer were confirmed by SS-NMR and fluorescence confocal microscopy. Leakage assays were performed to assess the lytic power of the antimicrobial peptides. The in situ perturbation of the ghost membranes by aurein 1.2 and caerin 1.1 revealed by ³¹P and ²H SS-NMR is consistent with membrane perturbation through a carpet mechanism for aurein 1.2, whereas caerin 1.1 acts on RBCs via pore formation. These results are compatible with fluorescence microscopy images of the ghosts. The peptides interact with eukaryotic membranes following similar mechanisms that take place in bacteria, highlighting the importance of hydrophobicity when determining such interactions. Our work bridges model membranes and in vitro studies and provides an analytical toolbox to assess drug toxicity toward eukaryotic cells.