Vesicle dynamics in a confined Poiseuille flow: from steady state to chaos

Dynamics and rheology of vesicles under confined Poiseuille flow

The rheological behavior and dynamics of a vesicle suspension, serving as a simplified model for red blood cells, are explored within a Poiseuille flow under the Stokes limit. Investigating vesicle response has led to the identification of novel solutions that complement previously documented forms like the parachute and slipper shapes. This study has brought to light the existence of alternative configurations, including a fully off-centered form and a multilobe structure. The study unveils the presence of two distinct branches associated with the slipper shape. One branch arises as a consequence of a supercritical bifurcation from the symmetric parachute shape, while the other emerges from a saddle-node bifurcation. Notably, the findings are represented through diagrams that display data collapsing harmoniously based on a combination of independent dimensionless parameters. Delving into the rheological implications, a remarkable observation emerges: the normalized viscosity (i.e. similar to intrinsic viscosity) exhibits a non-monotonic trend as a function of vesicle concentration. Initially, the normalized viscosity diminishes as the concentration increases, followed by a subsequent rise at higher concentrations. Noteworthy is the presence of a minimum value in the normalized viscosity at lower concentrations, aligning well with the concentrations observed in microcirculation scenarios. The intricate behavior of the normalized viscosity can be attributed to a delicate spatial arrangement within the suspension. Importantly, this trend echoes the observations made in a linear shear flow scenario, thereby underscoring the universality of the rheological behavior for confined suspensions.

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.

Stability and Deformation of Vesicles in a Cylindrical Flow

In this work, we used dissipative particle dynamics to study the stability, deformation, and rupture of polymer vesicles confined in cylindrical channels under the flow field. The morphological evolution, elongation, and rupture of vesicles and the corresponding mechanisms were intensively investigated. Bullet-like vesicles, leaking vesicles, spherical micelles, hamburger-like micelles, and bilayers were observed by changing the degree of confinement and dimensionless shear rate. We found that increasing the dimensionless shear rate and the degree of confinement can cause the deformation or rupture of polymeric vesicles. The asphericity parameter was utilized to describe the degree of elongation of vesicles deviating from the sphere in the direction of the flow. The results show that the aggregates are more likely to be spherical when the confinement is weak, while they become elongated bullet-like shapes when the confinement is strong. The investigation of dynamics reveals that the degree of confinement and the dimensionless shear rate can affect the chain stretching and reorganization during the process of vesicle elongation. Furthermore, the rupture time of the vesicle shows a nonlinear decrease with an increase in the dimensionless shear rate, and the confinement also contributes to the rupture. The results are very useful for guiding the application of vesicles in a flow environment.

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 cells in low Reynolds number flow: A vorticity-based characterization of shapes in two dimensions

Studies on the mechanical properties of red blood cells improve the diagnosis of some blood-related diseases. Some existing numerical methods have successfully simulated the coupling between a fluid and red blood cells. This paper introduces an alternative phase-field model formulation of two-dimensional cells that solves the vorticity and stream function that simplifies the numerical implementation. We integrate red blood cell dynamics immersed in a Poiseuille flow and reproduce previously reported morphologies (slippers or parachutes). In the case of flow in a very wide channel, we discover a new metastable shape referred to as 'anti-parachute' that evolves into a horizontal slipper centered on the channel. This sort of metastable morphology may contribute to the dynamical response of the blood.

Red blood cells under flow show maximal ATP release for specific hematocrit

ATP release by red blood cells (RBCs) under shear stress (SS) plays a pivotal role in endothelial biochemical signaling cascades. The aim of this study is to investigate through numerical simulation how RBC spatiotemporal organization depends on flow and geometrical conditions to generate ATP patterns. Numerical simulations were conducted in a straight channel by considering both plasma and explicit presence of RBCs, their shape deformation and cell-cell interaction, and ATP release by RBCs. Two ATP release pathways through cell membrane are taken into account: pannexin 1 channel, sensitive to SS, and cystic fibrosis transmembrane conductance regulator, which responds to cell deformation. Several flow and hematocrit conditions are explored. The problem is solved by the lattice Boltzmann method. Application of SS to the RBC suspension triggers a nontrivial spatial RBC organization and ATP patterns. ATP localizes preferentially in the vicinity of the cell-free layer close to channel wall. Conditions for maximal ATP release per cell are identified, which depend on vessel size and hematocrit Ht. Increasing further Ht beyond optimum enhances the total ATP release but should degrade oxygen transport capacity, a compromise between an efficient ATP release and minimal blood dissipation. Moreover, ATP is boosted in capillaries, suggesting a vasomotor activity coordination throughout the resistance network.

Axial and Nonaxial Migration of Red Blood Cells in a Microtube

Human red blood cells (RBCs) are subjected to high viscous shear stress, especially during microcirculation, resulting in stable deformed shapes such as parachute or slipper shape. Those unique deformed RBC shapes, accompanied with axial or nonaxial migration, cannot be fully described according to traditional knowledge about lateral movement of deformable spherical particles. Although several experimental and numerical studies have investigated RBC behavior in microchannels with similar diameters as RBCs, the detailed mechanical characteristics of RBC lateral movement—in particular, regarding the relationship between stable deformed shapes, equilibrium radial RBC position, and membrane load—has not yet been fully described. Thus, we numerically investigated the behavior of single RBCs with radii of 4 μm in a circular microchannel with diameters of 15 μm. Flow was assumed to be almost inertialess. The problem was characterized by the capillary number, which is the ratio between fluid viscous force and membrane elastic force. The power (or energy dissipation) associated with membrane deformations was introduced to quantify the state of membrane loads. Simulations were performed with different capillary numbers, viscosity ratios of the internal to external fluids of RBCs, and initial RBC centroid positions. Our numerical results demonstrated that axial or nonaxial migration of RBC depended on the stable deformed RBC shapes, and the equilibrium radial position of the RBC centroid correlated well with energy expenditure associated with membrane deformations.

Hydrodynamic Shape Changes Underpin Nuclear Rerouting in Branched Hyphae of an Oomycete Pathogen

Multinucleate fungi and oomycetes are phylogenetically distant but structurally similar. To address whether they share similar nuclear dynamics, we carried out time-lapse imaging of fluorescently labeled Phytophthora palmivora nuclei. Nuclei underwent coordinated bidirectional movements during plant infection. Within hyphal networks growing in planta or in axenic culture, nuclei either are dragged passively with the cytoplasm or actively become rerouted toward nucleus-depleted hyphal sections and often display a very stretched shape. Benomyl-induced depolymerization of microtubules reduced active movements and the occurrence of stretched nuclei. A centrosome protein localized at the leading end of stretched nuclei, suggesting that, as in fungi, astral microtubule-guided movements contribute to nuclear distribution within oomycete hyphae. The remarkable hydrodynamic shape adaptations of Phytophthora nuclei contrast with those in fungi and likely enable them to migrate over longer distances. Therefore, our work summarizes mechanisms which enable a near-equal nuclear distribution in an oomycete. We provide a basis for computational modeling of hydrodynamic nuclear deformation within branched tubular networks.IMPORTANCE Despite their fungal morphology, oomycetes constitute a distinct group of protists related to brown algae and diatoms. Many oomycetes are pathogens and cause diseases of plants, insects, mammals, and humans. Extensive efforts have been made to understand the molecular basis of oomycete infection, but durable protection against these pathogens is yet to be achieved. We use a plant-pathogenic oomycete to decipher a key physiological aspect of oomycete growth and infection. We show that oomycete nuclei travel actively and over long distances within hyphae and during infection. Such movements require microtubules anchored on the centrosome. Nuclei hydrodynamically adapt their shape to travel in or against the flow. In contrast, fungi lack a centrosome and have much less flexible nuclei. Our findings provide a basis for modeling of flexible nuclear shapes in branched hyphal networks and may help in finding hard-to-evade targets to develop specific antioomycete strategies and achieve durable crop disease protection.

Migration velocity of red blood cells in microchannels

The lateral migration of red blood cells (RBCs) in confined channel flows is an important ingredient of microcirculatory hydrodynamics and is involved in the development of a cell free layer near vessel walls and influences the distribution of RBCs in networks. It is also relevant to a number of lab-on-chip applications. This migration is a consequence of their deformability and is due to the combined effects of hydrodynamic wall repulsion and the curvature of the fluid velocity profile. We performed microfluidic experiments with dilute suspensions of RBCs in which the trajectories and migration away from the channel wall are analyzed to extract the mean behavior, from which we propose a generic scaling law for the transverse migration velocity valid in a whole range of parameters relevant to microcirculatory and practical situations. Experiments with RBCs of different mechanical properties (separated by density gradient sedimentation or fixed with glutaraldehyde) show the influence of this parameter which can induce significant dispersion of the trajectories.

ATP Release by Red Blood Cells under Flow: Model and Simulations

ATP is a major player as a signaling molecule in blood microcirculation. It is released by red blood cells (RBCs) when they are subjected to shear stresses large enough to induce a sufficient shape deformation. This prominent feature of chemical response to shear stress and RBC deformation constitutes an important link between vessel geometry, flow conditions, and the mechanical properties of RBCs, which are all contributing factors affecting the chemical signals in the process of vasomotor modulation of the precapillary vessel networks. Several in vitro experiments have reported on ATP release by RBCs due to mechanical stress. These studies have considered both intact RBCs as well as cells within which suspected pathways of ATP release have been inhibited. This has provided profound insights to help elucidate the basic governing key elements, yet how the ATP release process takes place in the (intermediate) microcirculation zone is not well understood. We propose here an analytical model of ATP release. The ATP concentration is coupled in a consistent way to RBC dynamics. The release of ATP, or the lack thereof, is assumed to depend on both the local shear stress and the shape change of the membrane. The full chemo-mechanical coupling problem is written in a lattice-Boltzmann formulation and solved numerically in different geometries (straight channels and bifurcations mimicking vessel networks) and under two kinds of imposed flows (shear and Poiseuille flows). Our model remarkably reproduces existing experimental results. It also pinpoints the major contribution of ATP release when cells traverse network bifurcations. This study may aid in further identifying the interplay between mechanical properties and chemical signaling processes involved in blood microcirculation.

Numerical-experimental observation of shape bistability of red blood cells flowing in a microchannel

Red blood cells flowing through capillaries assume a wide variety of different shapes owing to their high deformability. Predicting the realized shapes is a complex field as they are determined by the intricate interplay between the flow conditions and the membrane mechanics. In this work we construct the shape phase diagram of a single red blood cell with a physiological viscosity ratio flowing in a microchannel. We use both experimental in vitro measurements as well as 3D numerical simulations to complement the respective other one. Numerically, we have easy control over the initial starting configuration and natural access to the full 3D shape. With this information we obtain the phase diagram as a function of initial position, starting shape and cell velocity. Experimentally, we measure the occurrence frequency of the different shapes as a function of the cell velocity to construct the experimental diagram which is in good agreement with the numerical observations. Two different major shapes are found, namely croissants and slippers. Notably, both shapes show coexistence at low (<1 mm s^-1) and high velocities (>3 mm s^-1) while in-between only croissants are stable. This pronounced bistability indicates that RBC shapes are not only determined by system parameters such as flow velocity or channel size, but also strongly depend on the initial conditions.

The buckling instability of aggregating red blood cells

Plasma proteins such as fibrinogen induce the aggregation of red blood cells (RBC) into rouleaux, which are responsible for the pronounced shear thinning behavior of blood, control the erythrocyte sedimentation rate (ESR) – a common hematological test – and are involved in many situations of physiological relevance such as structuration of blood in the microcirculation or clot formation in pathological situations. Confocal microscopy is used to characterize the shape of RBCs within rouleaux at equilibrium as a function of macromolecular concentration, revealing the diversity of contact zone morphology. Three different configurations that have only been partly predicted before are identified, namely parachute, male-female and sigmoid shapes, and quantitatively recovered by numerical simulations. A detailed experimental and theoretical analysis of clusters of two cells shows that the deformation increases nonlinearly with the interaction energy. Models indicate a forward bifurcation in which the contacting membrane undergoes a buckling instability from a flat to a deformed contact zone at a critical value of the interaction energy. These results are not only relevant for the understanding of the morphology and stability of RBC aggregates, but also for a whole class of interacting soft deformable objects such as vesicles, capsules or cells in tissues.

Theory and algorithms to compute Helfrich bending forces: a review

Cell membranes are vital to shield a cell's interior from the environment. At the same time they determine to a large extent the cell's mechanical resistance to external forces. In recent years there has been considerable interest in the accurate computational modeling of such membranes, driven mainly by the amazing variety of shapes that red blood cells and model systems such as vesicles can assume in external flows. Given that the typical height of a membrane is only a few nanometers while the surface of the cell extends over many micrometers, physical modeling approaches mostly consider the interface as a two-dimensional elastic continuum. Here we review recent modeling efforts focusing on one of the computationally most intricate components, namely the membrane's bending resistance. We start with a short background on the most widely used bending model due to Helfrich. While the Helfrich bending energy by itself is an extremely simple model equation, the computation of the resulting forces is far from trivial. At the heart of these difficulties lies the fact that the forces involve second order derivatives of the local surface curvature which by itself is the second derivative of the membrane geometry. We systematically derive and compare the different routes to obtain bending forces from the Helfrich energy, namely the variational approach and the thin-shell theory. While both routes lead to mathematically identical expressions, so-called linear bending models are shown to reproduce only the leading order term while higher orders differ. The main part of the review contains a description of various computational strategies which we classify into three categories: the force, the strong and the weak formulation. We finally give some examples for the application of these strategies in actual simulations.

Clusters of red blood cells in microcapillary flow: hydrodynamic versus macromolecule induced interaction

We present experiments on RBCs that flow through micro-capillaries under physiological conditions. The strong flow-shape coupling of these deformable objects leads to a rich variety of cluster formation. We show that the RBC clusters form as a subtle imbrication between hydrodynamic interactions and adhesion forces because of plasma proteins, mimicked by the polymer dextran. Clusters form along the capillaries and macromolecule-induced adhesion contributes to their stability. However, at high yet physiological flow velocities, shear stresses overcome part of the adhesion forces, and cluster stabilization due to hydrodynamics becomes stronger. For the case of pure hydrodynamic interaction, cell-to-cell distances have a pronounced bimodal distribution. Our 2D-numerical simulations on vesicles capture the transition between adhesive and non-adhesive clusters at different flow velocities.

Dynamics of flexible fibers in shear flow

Dynamics of flexible non-Brownian fibers in shear flow at low-Reynolds-number are analyzed numerically for a wide range of the ratios A of the fiber bending force to the viscous drag force. Initially, the fibers are aligned with the flow, and later they move in the plane perpendicular to the flow vorticity. A surprisingly rich spectrum of different modes is observed when the value of A is systematically changed, with sharp transitions between coiled and straightening out modes, period-doubling bifurcations from periodic to migrating solutions, irregular dynamics, and chaos.