Search

The generation of fluid flows by autophoretic microswimmers has been proposed as a mechanism to enhance mass transport and mixing at the micro- and nanoscale. Here, we experimentally investigate the ability of model 2D active baths of photocatalytic silica-titania Janus microspheres to enhance the diffusivity of tracer particles at different microswimmer densities below the onset of collective behavior. Inspired by the similarities between our experimental findings and previous results for biological microorganisms, we then model our Janus microswimmers using a general squirmer framework, specifically treating them as neutral squirmers. The numerical simulations faithfully capture our observations, offer an insight into the microscopic mechanism underpinning tracer transport, and allow us to expand the parameter space beyond our experimental system. We find strong evidence that near-field interactions dominate enhancements in tracer diffusivity in active Janus baths, leading to the identification of an operating window for enhanced tracer transport by chemical microswimmers based on scaling arguments. Based on this argumentation, we suggest that for many chemically active colloidal systems, hydrodynamics alone is likely to be insufficient to induce appreciable mixing of passive components with large diffusion coefficients.

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.

Competing Interests: Declaration of interests The authors declare no competing interests.

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.

Competing Interests: Declaration of interests The author declares no competing interests.

Recent advances in micro- and nano-technologies allow the construction of complex active systems from biological and synthetic materials. An interesting example is active vesicles, which consist of a membrane enclosing self-propelled particles, and exhibit several features resembling biological cells. We investigate numerically the behavior of active vesicles, where the enclosed self-propelled particles can adhere to the membrane. A vesicle is represented by a dynamically triangulated membrane, while the adhesive active particles are modelled as active Brownian particles (ABPs) that interact with the membrane via the Lennard-Jones potential. Phase diagrams of dynamic vesicle shapes as a function of ABP activity and particle volume fraction inside the vesicle are constructed for different strengths of adhesive interactions. At low ABP activity, adhesive interactions dominate over the propulsion forces, such that the vesicle attains near static configurations, with protrusions of membrane-wrapped ABPs having ring-like and sheet-like structures. At moderate particle densities and strong enough activities, active vesicles show dynamic highly-branched tethers filled with string-like arrangements of ABPs, which do not occur in the absence of particle adhesion to the membrane. At large volume fractions of ABPs, vesicles fluctuate for moderate particle activities, and elongate and finally split into two vesicles for large ABP propulsion strengths. We also analyze membrane tension, active fluctuations, and ABP characteristics ( e.g. , mobility, clustering), and compare them to the case of active vesicles with non-adhesive ABPs. The adhesion of ABPs to the membrane significantly alters the behavior of active vesicles, and provides an additional parameter for controlling their behavior.

Cells in living organisms are subjected to mechanical strains caused by external forces like overcrowding, resulting in strong deformations that affect cell function. We study the interplay between deformation and crowding of red blood cells (RBCs) in dispersions of nonabsorbing rod-like viruses. We identify a sequence of configurational transitions of RBC doublets, including configurations that can only be induced by long-ranged attraction: highly fluctuating T-shaped and face-to-face configurations at low, and doublets approaching a complete spherical configuration at high, rod concentrations. Complementary simulations are used to explore different energy contributions to deformation as well as the stability of RBC doublet configurations. Our advanced analysis of 3D reconstructed confocal images of RBC doublets quantifies the depletion interaction and the resulting deformation energy. Thus, we introduce a noninvasive, high-throughput platform that is generally applicable to investigate the mechanical response of biological cells to external forces and characterize their mechanical properties., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Clustering of flagellated microswimmers such as sperm is often mediated by hydrodynamic interactions between them. To better understand the interaction of microswimmers in viscoelastic fluids, we perform two-dimensional simulations of two swimming sheets, using a viscoelastic version of the smoothed dissipative particle dynamics method that implements the Oldroyd-B fluid model. Elasticity of sheets (stiff versus soft) defines two qualitatively different regimes of clustering, where stiff sheets exhibit a much more robust clustering than soft sheets. A formed doublet of soft sheets generally swims faster than a single swimmer, while a pair of two stiff sheets normally shows no speed enhancement after clustering. A pair of two identical swimmers is stable for most conditions, while differences in the beating amplitudes and/or frequencies between the two sheets can destroy the doublet stability. Clustering of two distinct swimmers is most stable at Deborah numbers of De = τω ≈ 1 ( τ is the relaxation time of a viscoelastic fluid and ω is the beating frequency), in agreement with experimental observations. Therefore, the clustering of two swimmers depends non-monotonically on De. Our results suggest that the cluster stability is likely a dominant factor which determines the cluster size of collectively moving flagellated swimmers.

Active vesicles, constructed through the confinement of self-propelled particles (SPPs) inside a lipid membrane shell, exhibit a large variety of non-equilibrium shapes, ranging from the formation of local tethers and dendritic conformations, to prolate and bola-like structures. To better understand the behavior of active vesicles, we perform simulations of membranes modelled as dynamically triangulated surfaces enclosing active Brownian particles. A systematic analysis of membrane deformations and SPP clustering, as a function of SPP activity and volume fraction inside the vesicle is carried out. Distributions of membrane local curvature, and the clustering and mobility of SPPs obtained from simulations of active vesicles are analysed. There exists a feedback mechanism between the enhancement of membrane curvature, the formation of clusters of active particles, and local or global changes in vesicle shape. The emergence of active tension due to the activity of SPPs can well be captured by the Young-Laplace equation. Furthermore, a simple numerical method for tether detection is presented and used to determine correlations between the number of tethers, their length, and local curvature. We also provide several geometrical arguments to explain different tether characteristics for various conditions. These results contribute to the future development of steerable active vesicles or soft micro-robots whose behaviour can be controlled and used for potential applications.

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.

The erythrocyte (or red blood cell) sedimentation rate (ESR) is commonly interpreted as a measure of cell aggregation and as a biomarker of inflammation. It is well known that an increase of fibrinogen concentration, an aggregation-inducing protein for erythrocytes, leads to an increase of the sedimentation rate of erythrocytes, which is generally explained through the formation and faster settling of large disjoint aggregates. However, many aspects of erythrocyte sedimentation conform well with the collapse of a particle gel rather than with the sedimentation of disjoint aggregates. Using experiments and cell-level numerical simulations, we systematically investigate the dependence of ESR on fibrinogen concentration and its relation to the microstructure of the gel-like erythrocyte suspension. We show that for physiological aggregation interactions, an increase in the attraction strength between cells results in a cell network with larger void spaces. This geometrical change in the network structure occurs due to anisotropic shape and deformability of erythrocytes and leads to an increased gel permeability and faster sedimentation. Our results provide a comprehensive relation between the ESR and the cell-level structure of erythrocyte suspensions and support the gel hypothesis in the interpretation of blood sedimentation.

During the blood stage of malaria pathogenesis, parasites invade healthy red blood cells (RBC) to multiply inside the host and evade the immune response. When attached to RBC, the parasite first has to align its apex with the membrane for a successful invasion. Since the parasite's apex sits at the pointed end of an oval (egg-like) shape with a large local curvature, apical alignment is in general an energetically unfavorable process. Previously, using coarse-grained mesoscopic simulations, we have shown that optimal alignment time is achieved due to RBC membrane deformation and the stochastic nature of bond-based interactions between the parasite and RBC membrane (Hillringhaus et al., 2020). Here, we demonstrate that the parasite's shape has a prominent effect on the alignment process. The alignment times of spherical parasites for intermediate and large bond off-rates (or weak membrane-parasite interactions) are found to be close to those of an egg-like shape. However, for small bond off-rates (or strong adhesion and large membrane deformations), the alignment time for a spherical shape increases drastically. Parasite shapes with large aspect ratios such as oblate and long prolate ellipsoids are found to exhibit very long alignment times in comparison to the egg-like shape. At a stiffened RBC, a spherical parasite aligns faster than any other investigated shape. This study shows that the original egg-like shape performs not worse for parasite alignment than other considered shapes but is more robust with respect to different adhesion interactions and RBC membrane rigidities., Competing Interests: AD, SH, GG, DF No competing interests declared, (© 2021, Dasanna et al.)

Background: Chorea-acanthocytosis (ChAc) is a rare hereditary neurodegenerative disease with deformed red blood cells (RBCs), so-called acanthocytes, as a typical marker of the disease. Erythrocyte sedimentation rate (ESR) was recently proposed as a diagnostic biomarker. To date, there is no treatment option for affected patients, but promising therapy candidates, such as dasatinib, a Lyn-kinase inhibitor, have been identified., Methods: RBCs of two ChAc patients during and after dasatinib treatment were characterized by the ESR, clinical hematology parameters and the 3D shape classification in stasis based on an artificial neural network. Furthermore, mathematical modeling was performed to understand the contribution of cell morphology and cell rigidity to the ESR. Microfluidic measurements were used to compare the RBC rigidity between ChAc patients and healthy controls., Results: The mechano-morphological characterization of RBCs from two ChAc patients in an off-label treatment with dasatinib revealed differences in the ESR and the acanthocyte count during and after the treatment period, which could not directly be related to each other. Clinical hematology parameters were in the normal range. Mathematical modeling indicated that RBC rigidity is more important for delayed ESR than cell shape. Microfluidic experiments confirmed a higher rigidity in the normocytes of ChAc patients compared to healthy controls., Conclusions: The results increase our understanding of the role of acanthocytes and their associated properties in the ESR, but the data are too sparse to answer the question of whether the ESR is a suitable biomarker for treatment success, whereas a correlation between hematological and neuronal phenotype is still subject to verification.

(1) Background: Chorea-acanthocytosis and McLeod syndrome are the core diseases among the group of rare neurodegenerative disorders called neuroacanthocytosis syndromes (NASs). NAS patients have a variable number of irregularly spiky erythrocytes, so-called acanthocytes. Their detection is a crucial but error-prone parameter in the diagnosis of NASs, often leading to misdiagnoses. (2) Methods: We measured the standard Westergren erythrocyte sedimentation rate (ESR) of various blood samples from NAS patients and healthy controls. Furthermore, we manipulated the ESR by swapping the erythrocytes and plasma of different individuals, as well as replacing plasma with dextran. These measurements were complemented by clinical laboratory data and single-cell adhesion force measurements. Additionally, we followed theoretical modeling approaches. (3) Results: We show that the acanthocyte sedimentation rate (ASR) with a two-hour read-out is significantly prolonged in chorea-acanthocytosis and McLeod syndrome without overlap compared to the ESR of the controls. Mechanistically, through modern colloidal physics, we show that acanthocyte aggregation and plasma fibrinogen levels slow down the sedimentation. Moreover, the inverse of ASR correlates with the number of acanthocytes (R2=0.61, p=0.004). (4) Conclusions: The ASR/ESR is a clear, robust and easily obtainable diagnostic marker. Independently of NASs, we also regard this study as a hallmark of the physical view of erythrocyte sedimentation by describing anticoagulated blood in stasis as a percolating gel, allowing the application of colloidal physics theory.

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., (© 2020 The Authors. Microcirculation published by John Wiley & Sons Ltd.)

Biological cells generate intricate structures by sculpting their membrane from within to actively sense and respond to external stimuli or to explore their environment 1-4 . Several pathogenic bacteria also provide examples of how localized forces strongly deform cell membranes from inside, leading to the invasion of neighbouring healthy mammalian cells 5 . Giant unilamellar vesicles have been successfully used as a minimal model system with which to mimic biological cells 6-11 , but the realization of a minimal system with localized active internal forces that can strongly deform lipid membranes from within and lead to dramatic shape changes remains challenging. Here we present a combined experimental and simulation study that demonstrates how self-propelled particles enclosed in giant unilamellar vesicles can induce a plethora of non-equilibrium shapes and active membrane fluctuations. Using confocal microscopy, in the experiments we explore the membrane response to local forces exerted by self-phoretic Janus microswimmers. To quantify dynamic membrane changes, we perform Langevin dynamics simulations of active Brownian particles enclosed in thin membrane shells modelled by dynamically triangulated surfaces. The most pronounced shape changes are observed at low and moderate particle loadings, with the formation of tether-like protrusions and highly branched, dendritic structures, whereas at high volume fractions globally deformed vesicle shapes are observed. The resulting state diagram predicts the conditions under which local internal forces generate various membrane shapes. A controlled realization of such distorted vesicle morphologies could improve the design of artificial systems such as small-scale soft robots and synthetic cells.

The advent of microfluidics in the 1990s promised a revolution in multiple industries from healthcare to chemical processing. Deterministic lateral displacement (DLD) is a continuous-flow microfluidic particle separation method discovered in 2004 that has been applied successfully and widely to the separation of blood cells, yeast, spores, bacteria, viruses, DNA, droplets, and more. Deterministic lateral displacement is conceptually simple and can deliver consistent performance over a wide range of flow rates and particle concentrations. Despite wide use and in-depth study, DLD has not yet been fully elucidated or optimized, with different approaches to the same problem yielding varying results. We endeavor here to provide up-to-date expert opinion on the state-of-art and current fundamental, practical, and commercial challenges with DLD as well as describe experimental and modeling opportunities. Because these challenges and opportunities arise from constraints on hydrodynamics, fabrication, and operation at the micro- and nanoscale, we expect this Perspective to serve as a guide for the broader micro- and nanofluidic community to identify and to address open questions in the field.

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., Competing Interests: SH, AD, GG, DF No competing interests declared, (© 2020, Hillringhaus et al.)

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.

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., (Copyright © 2019 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

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.

Red blood cells in shear flow show a variety of different shapes due to the complex interplay between hydrodynamics and membrane elasticity. Malaria-infected red blood cells become generally adhesive and less deformable. Adhesion to a substrate leads to a reduction in shape variability and to a flipping motion of the non-spherical shapes during the mid-stage of infection. Here, we present a complete state diagram for wall adhesion of red blood cells in shear flow obtained by simulations, using a particle-based mesoscale hydrodynamics approach, multiparticle collision dynamics. We find that cell flipping at a substrate is replaced by crawling beyond a critical shear rate, which increases with both membrane stiffness and viscosity contrast between the cytosol and suspending medium. This change in cell dynamics resembles the transition between tumbling and tank-treading for red blood cells in free shear flow. In the context of malaria infections, the flipping-crawling transition would strongly increase the adhesive interactions with the vascular endothelium, but might be suppressed by the combined effect of increased elasticity and viscosity contrast.

The motion of red blood cells (RBCs) in microchannels is important for microvascular blood flow and biomedical applications such as blood analysis in microfluidics. The current understanding of the complexity of RBC shapes and dynamics in microchannels is mainly based on several simulation studies, but there are a few systematic experimental investigations. Here, we present a combined study that systematically characterizes RBC behavior for a wide range of flow rates and channel sizes. Even though simulations and experiments generally show good agreement, experimental observations demonstrate that there is no single well-defined RBC state for fixed flow conditions but rather a broad distribution of states. This result can be attributed to the inherent variability in RBC mechanical properties, which is confirmed by a model that takes the variation in RBC shear elasticity into account. This represents a significant step toward a quantitative connection between RBC behavior in microfluidic devices and their mechanical properties, which is essential for a high-throughput characterization of diseased cells., (Copyright © 2019 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

A recent study of red blood cells (RBCs) in shear flow [Lanotte et al., Proc. Natl. Acad. Sci. U.S.A. 113, 13289 (2016)PNASA60027-842410.1073/pnas.1608074113] has demonstrated that RBCs first tumble, then roll, transit to a rolling and tumbling stomatocyte, and finally attain polylobed shapes with increasing shear rate, when the viscosity contrast between cytosol and blood plasma is large enough. Using two different simulation techniques, we construct a state diagram of RBC shapes and dynamics in shear flow as a function of shear rate and viscosity contrast, which is also supported by microfluidic experiments. Furthermore, we illustrate the importance of RBC shear elasticity for its dynamics in flow and show that two different kinds of membrane buckling trigger the transition between subsequent RBC states.

Intravascular drug delivery technologies majorly utilize spherical nanoparticles as carrier vehicles. Their targets are often at the blood vessel wall or in the tissue beyond the wall, such that vehicle localization towards the wall (margination) becomes a pre-requisite for their function. To this end, some studies have indicated that under flow environment, micro-particles have a higher propensity than nano-particles to marginate to the wall. Also, non-spherical particles theoretically have a higher area of surface-adhesive interactions than spherical particles. However, detailed systematic studies that integrate various particle size and shape parameters across nano-to-micro scale to explore their wall-localization behavior in RBC-rich blood flow, have not been reported. We address this gap by carrying out computational and experimental studies utilizing particles of four distinct shapes (spherical, oblate, prolate, rod) spanning nano- to-micro scale sizes. Computational studies were performed using the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) package, with Dissipative Particle Dynamics (DPD). For experimental studies, model particles were made from neutrally buoyant fluorescent polystyrene spheres, that were thermo-stretched into non-spherical shapes and all particles were surface-coated with biotin. Using microfluidic setup, the biotin-coated particles were flowed over avidin-coated surfaces in absence versus presence of RBCs, and particle adhesion and retention at the surface was assessed by inverted fluorescence microscopy. Our computational and experimental studies provide a simultaneous analysis of different particle sizes and shapes for their retention in blood flow and indicate that in presence of RBCs, micro-scale non-spherical particles undergo enhanced 'margination + adhesion' compared to nano-scale spherical particles, resulting in their higher binding. These results provide important insight regarding improved design of vascularly targeted drug delivery systems.

Red blood cell (RBC) aggregates play an important role in determining blood rheology. RBCs in plasma or polymer solution interact attractively to form various shapes of RBC doublets, where the attractive interactions can be varied by changing the solution conditions. A systematic numerical study on RBC doublet formation is performed, which takes into account the shear elasticity of the RBC membrane due to the spectrin cytoskeleton, in addition to the membrane bending rigidity. RBC membranes are modeled by two-dimensional triangular networks of linked vertices, which represent three-dimensional cell shapes. The phase space of RBC doublet shapes in a wide range of adhesion strengths, reduced volumes, and shear elasticities is obtained. The shear elasticity of the RBC membrane changes the doublet phases significantly. Experimental images of RBC doublets in different solutions show similar configurations. Furthermore, we show that rouleau formation is affected by the doublet structure.

Adhesion of polymers and proteins to substrates plays a crucial role in many technological applications and biological processes. A prominent example is the von Willebrand factor (VWF) protein, which is essential in blood clotting as it mediates adhesion of blood platelets to the site of injury at high shear rates. VWF is activated by flow and is able to bind efficiently to damaged vessel walls even under extreme flow-stress conditions; however, its adhesion is reversible when the flow strength is significantly reduced or the flow is ceased. Motivated by the properties and behavior of VWF in flow, we investigate adhesion of shear-activated polymers to a planar wall in flow and whether the adhesion is reversible under flow stasis. The main ingredients of the polymer model are cohesive inter-monomer interactions, a catch bond with the adhesive surface, and the shear activation/deactivation of polymer adhesion correlated with its stretching in flow. The cohesive interactions within the polymer maintain a globular conformation under low shear stresses and allow polymer stretching if a critical shear rate is exceeded, which is directly associated with its activation for adhesion. Our results show that polymer adhesion at high shear rates is significantly stabilized by catch bonds, while at the same time they also permit polymer dissociation from a surface at low or no flow stresses. In addition, the activation/deactivation mechanism for adhesion plays a crucial role in the reversibility of its adhesion. These observations help us better understand the adhesive behavior of VWF in flow and interpret its adhesion malfunctioning in VWF-related diseases.

The protein von Willebrand factor (VWF) is essential in primary hemostasis, as it mediates platelet adhesion to vessel walls. VWF retains its compact (globule-like) shape in equilibrium due to internal molecular associations, but is able to stretch when a high enough shear stress is applied. Even though the shear-flow sensitivity of VWF conformation is well accepted, the behavior of VWF under realistic blood flow conditions remains poorly understood. We perform mesoscopic numerical simulations together with microfluidic experiments in order to characterize VWF behavior in blood flow for a wide range of flow-rate and hematocrit conditions. In particular, our results demonstrate that the compact shape of VWF is important for its migration (or margination) toward vessel walls and that VWF stretches primarily in a near-wall region in blood flow making its adhesion possible. Our results show that VWF is a highly optimized protein in terms of its size and internal associations which are necessary to achieve its vital function. A better understanding of the relevant mechanisms for VWF behavior in microcirculation provides a further step toward the elucidation of the role of mutations in various VWF-related diseases.

Von Willebrand factor (VWF) is a key protein in hemostasis as it mediates adhesion of blood platelets to a site of vascular injury. A proper distribution of VWF lengths is important for normal functioning of hemostatic processes, because a diminished number of long VWF chains may significantly limit blood clotting and lead to bleeding, while an abundant number of long VWFs may result in undesired thrombotic events. VWF size distribution is controlled by ADAMTS13 protease, which can cleave VWF chains beyond a critical shear rate when the chains are stretched enough such that cleavage sites become accessible. To better understand the cleavage process, we model VWF cleavage in shear flow using mesoscopic hydrodynamic simulations. Two cleavage models are proposed, a geometrical model based on the degree of local stretching of VWF, and a tension-force model based on instantaneous tension force within VWF bonds. Both models capture the susceptibility of VWF to cleavage at high shear rates; however, the geometrical model appears to be much more robust than the force model. Our simulations show that VWF susceptibility to cleavage in shear flow becomes a universal function of shear rate, independent of VWF length for long enough chains. Furthermore, VWF is cleaved with a higher probability close to its ends in comparison to cleaving in the middle, which results into longer circulation lifetimes of VWF multimers. Simulations of dynamic cleavage of VWF show an exponential distribution of chain lengths, consistently with available in vitro experiments. The proposed cleavage models can be used in realistic simulations of hemostatic processes in blood flow., (Copyright © 2017 IPEM. Published by Elsevier Ltd. All rights reserved.)

Light scattering is a well-established experimental technique, which gains more and more popularity in the biological field because it offers the means for non-invasive imaging and detection. However, the interpretation of light-scattering signals remains challenging due to the complexity of most biological systems. Here, we investigate static and dynamic scattering properties of red blood cells (RBCs) using two mesoscopic hydrodynamics simulation methods-multi-particle collision dynamics and dissipative particle dynamics. Light scattering is studied for various membrane shear elasticities, bending rigidities, and RBC shapes (e.g., biconcave and stomatocyte). Simulation results from the two simulation methods show good agreement, and demonstrate that the static light scattering of a diffusing RBC is not very sensitive to the changes in membrane properties and moderate alterations in cell shapes. We also compute dynamic light scattering of a diffusing RBC, from which dynamic properties of RBCs such as diffusion coefficients can be accessed. In contrast to static light scattering, the dynamic measurements can be employed to differentiate between the biconcave and stomatocytic RBC shapes and generally allow the differentiation based on the membrane properties. Our simulation results can be used for better understanding of light scattering by RBCs and the development of new non-invasive methods for blood-flow monitoring.

Blood viscosity decreases with shear stress, a property essential for an efficient perfusion of the vascular tree. Shear thinning is intimately related to the dynamics and mutual interactions of RBCs, the major component of blood. Because of the lack of knowledge about the behavior of RBCs under physiological conditions, the link between RBC dynamics and blood rheology remains unsettled. We performed experiments and simulations in microcirculatory flow conditions of viscosity, shear rates, and volume fractions, and our study reveals rich RBC dynamics that govern shear thinning. In contrast to the current paradigm, which assumes that RBCs align steadily around the flow direction while their membranes and cytoplasm circulate, we show that RBCs successively tumble, roll, deform into rolling stomatocytes, and, finally, adopt highly deformed polylobed shapes for increasing shear stresses, even for semidilute volume fractions of the microcirculation. Our results suggest that any pathological change in plasma composition, RBC cytosol viscosity, or membrane mechanical properties will affect the onset of these morphological transitions and should play a central role in pathological blood rheology and flow behavior., Competing Interests: The authors declare no conflict of interest.

Recent advances in cell sorting aim at the development of novel methods that are sensitive to various mechanical properties of cells. Microfluidic technologies have a great potential for cell sorting; however, the design of many micro-devices is based on theories developed for rigid spherical particles with size as a separation parameter. Clearly, most bioparticles are non-spherical and deformable and therefore exhibit a much more intricate behavior in fluid flow than rigid spheres. Here, we demonstrate the use of cells' mechanical and dynamical properties as biomarkers for separation by employing a combination of mesoscale hydrodynamic simulations and microfluidic experiments. The dynamic behavior of red blood cells (RBCs) within deterministic lateral displacement (DLD) devices is investigated for different device geometries and viscosity contrasts between the intra-cellular fluid and suspending medium. We find that the viscosity contrast and associated cell dynamics clearly determine the RBC trajectory through a DLD device. Simulation results compare well to experiments and provide new insights into the physical mechanisms which govern the sorting of non-spherical and deformable cells in DLD devices. Finally, we discuss the implications of cell dynamics for sorting schemes based on properties other than cell size, such as mechanics and morphology.

Microvascular blood flow determines a number of important physiological processes of an organism in health and disease. Therefore, a detailed understanding of microvascular blood flow would significantly advance biophysical and biomedical research and its applications. Current developments in modeling of microcirculatory blood flow already allow to go beyond available experimental measurements and have a large potential to elucidate blood flow behavior in normal and diseased microvascular networks. There exist detailed models of blood flow on a single cell level as well as simplified models of the flow through microcirculatory networks, which are reviewed and discussed here. The combination of these models provides promising prospects for better understanding of blood flow behavior and transport properties locally as well as globally within large microvascular networks., (© 2015 Wiley Periodicals, Inc.)

Targeted delivery of drugs and imaging agents is very promising to develop new strategies for the treatment of various diseases such as cancer. For an efficient targeted adhesion, the particles have to migrate toward the walls in blood flow - a process referred to as margination. Due to a huge diversity of available carriers, a good understanding of their margination properties in blood flow depending on various flow conditions and particle properties is required. We employ a particle-based mesoscopic hydrodynamic simulation approach to investigate the margination of different carriers for a wide range of hematocrits (volume fraction of red blood cells) and flow rates. Our results show that margination strongly depends on the thickness of the available free space close to the wall, the so-called red blood cell-free layer (RBC-FL), in comparison to the carrier size. The carriers with a few micrometers in size are comparable with the RBC-FL thickness and marginate better than their sub-micrometer counterparts. Deformable carriers, in general, show worse margination properties than rigid particles. Particle margination is also found to be most pronounced in small channels with a characteristic size comparable to blood capillaries. Finally, different margination mechanisms are discussed., (Copyright © 2015 IPEM. Published by Elsevier Ltd. All rights reserved.)

Deterministic lateral displacement (DLD) devices have great potential for the separation and sorting of various suspended particles based on their size, shape, deformability, and other intrinsic properties. Currently, the basic idea for the separation mechanism is that the structure and geometry of DLDs uniquely determine the flow field, which in turn defines a critical particle size and the particle lateral displacement within a device. We employ numerical simulations using coarse-grained mesoscopic methods and two-dimensional models to elucidate the dynamics of both rigid spherical particles and deformable red blood cells (RBCs) in different DLD geometries. Several shapes of pillars, including circular, diamond, square, and triangular structures, and a few particle sizes are considered. The simulation results show that a critical particle size can be well defined for rigid spherical particles and depends on the details of the DLD structure and the corresponding flow field within the device. However, non-isotropic and deformable particles such as RBCs exhibit much more complex dynamics within a DLD device, which cannot properly be described by a single parameter such as the critical size. The dynamics and deformation of soft particles within a DLD device become also important, indicating that not only size sorting, but additional sorting targets (e.g., shape, deformability, internal viscosity) are possible.

Janus colloids propelled by light, e.g., thermophoretic particles, offer promising prospects as artificial microswimmers. However, their swimming behavior and its dependence on fluid properties and fluid-colloid interactions remain poorly understood. Here, we investigate the behavior of a thermophoretic Janus colloid in its own temperature gradient using numerical simulations. The dissipative particle dynamics method with energy conservation is used to investigate the behavior in non-ideal and ideal-gas like fluids for different fluid-colloid interactions, boundary conditions, and temperature-controlling strategies. The fluid-colloid interactions appear to have a strong effect on the colloid behavior, since they directly affect heat exchange between the colloid surface and the fluid. The simulation results show that a reduction of the heat exchange at the fluid-colloid interface leads to an enhancement of colloid's thermophoretic mobility. The colloid behavior is found to be different in non-ideal and ideal fluids, suggesting that fluid compressibility plays a significant role. The flow field around the colloid surface is found to be dominated by a source-dipole, in agreement with the recent theoretical and simulation predictions. Finally, different temperature-control strategies do not appear to have a strong effect on the colloid's swimming velocity.

Microvascular blood flow resistance has a strong impact on cardiovascular function and tissue perfusion. The flow resistance in microcirculation is governed by flow behavior of blood through a complex network of vessels, where the distribution of red blood cells across vessel cross-sections may be significantly distorted at vessel bifurcations and junctions. In this paper, the development of blood flow and its resistance starting from a dispersed configuration of red blood cells is investigated in simulations for different hematocrit levels, flow rates, vessel diameters, and aggregation interactions between red blood cells. Initially dispersed red blood cells migrate toward the vessel center leading to the formation of a cell-free layer near the wall and to a decrease of the flow resistance. The development of cell-free layer appears to be nearly universal when scaled with a characteristic shear rate of the flow. The universality allows an estimation of the length of a vessel required for full flow development, lc ≲ 25D, for vessel diameters in the range 10 μm < D < 100 μm. Thus, the potential effect of red blood cell dispersion at vessel bifurcations and junctions on the flow resistance may be significant in vessels which are shorter or comparable to the length lc. Aggregation interactions between red blood cells generally lead to a reduction of blood flow resistance. The simulations are performed using the same viscosity for both external and internal fluids and the RBC membrane viscosity is not considered; however, we discuss how the viscosity contrast may affect the results. Finally, we develop a simple theoretical model which is able to describe the converged cell-free-layer thickness at steady-state flow with respect to flow rate. The model is based on the balance between a lift force on red blood cells due to cell-wall hydrodynamic interactions and shear-induced effective pressure due to cell-cell interactions in flow. We expect that these results can also be used to better understand the flow behavior of other suspensions of deformable particles such as vesicles, capsules, and cells., (Copyright © 2015 Elsevier Inc. All rights reserved.)

The motion of red blood cells (RBCs) in microcirculation plays an important role in blood flow resistance and in the cell partitioning within a microvascular network. Different shapes and dynamics of RBCs in microvessels have been previously observed experimentally including the parachute and slipper shapes. We employ mesoscale hydrodynamic simulations to predict the phase diagram of shapes and dynamics of RBCs in cylindrical microchannels, which serve as idealized microvessels, for a wide range of channel confinements and flow rates. A rich dynamical behavior is found, with snaking and tumbling discocytes, slippers performing a swinging motion, and stationary parachutes. We discuss the effects of different RBC states on the flow resistance, and the influence of RBC properties, characterized by the Föppl-von Kármán number, on the shape diagram. The simulations are performed using the same viscosity for both external and internal fluids surrounding a RBC; however, we discuss how the viscosity contrast would affect the shape diagram.

Proper functioning of white blood cells is not possible without their ability to adhere to vascular endothelium, which may occur only if they are close enough to vessel walls. To facilitate the adhesion, white blood cells migrate toward the vessel walls in blood flow through a process called margination. The margination of white cells depends on a number of conditions including local hematocrit, flow rate, red blood cell aggregation, and the deformability of both red and white cells. To better understand the margination process of white blood cells, we employ mesoscopic hydrodynamic simulations of a three-dimensional model of blood flow, which has been previously shown to capture quantitatively realistic blood flow properties and rheology. The margination properties of white blood cells are studied for a wide range of hematocrit values and flow conditions. Efficient white blood cell margination is found in an intermediate range of hematocrit values of Ht ≈ 0.2-0.4 and at relatively low flow rates, characteristic of the venular part of microcirculation. In addition, aggregation interactions between red blood cells lead to enhanced white-blood-cell margination. This simulation study provides a quantitative description of the margination of white blood cells, and is also highly relevant for the margination of particles or cells of similar size such as circulating tumor cells.

Drug delivery by micro- and nano-carriers enables controlled transport of pharmaceuticals to targeted sites. Even though carrier fabrication has made much progress recently, the delivery including controlled particle distribution and adhesion within the body remains a great challenge. The adhesion of carriers is strongly affected by their margination properties (migration toward walls) in the microvasculature. To investigate margination characteristics of carriers of different shapes and sizes and to elucidate the relevant physical mechanisms, we employ mesoscopic hydrodynamic simulations of blood flow. Particle margination is studied for a wide range of hematocrit values, vessel sizes, and flow rates, using two- and three-dimensional models. The simulations show that the margination properties of particles improve with increasing carrier size. Spherical particles yield slightly better margination than ellipsoidal carriers; however, ellipsoidal particles exhibit a slower rotational dynamics near a wall favoring their adhesion. In conclusion, micron-sized ellipsoidal particles are favorable for drug delivery in comparison with sub-micron spherical particles.

Mesoscale simulations of blood flow, where the red blood cells are described as deformable closed shells with a membrane characterized by bending rigidity and stretching elasticity, have made much progress in recent years to predict the flow behavior of blood cells and other components in various flows. To numerically investigate blood flow and blood-related processes in complex geometries, a highly efficient simulation technique for the plasma and solutes is essential. In this review, we focus on the behavior of single and several cells in shear and microcapillary flows, the shear-thinning behavior of blood and its relation to the blood cell structure and interactions, margination of white blood cells and platelets, and modeling hematologic diseases and disorders. Comparisons of the simulation predictions with existing experimental results are made whenever possible, and generally very satisfactory agreement is obtained.

Hematologic disorders arising from infectious diseases, hereditary factors and environmental influences can lead to, and can be influenced by, significant changes in the shape, mechanical and physical properties of red blood cells (RBCs), and the biorheology of blood flow. Hence, modeling of hematologic disorders should take into account the multiphase nature of blood flow, especially in arterioles and capillaries. We present here an overview of a general computational framework based on dissipative particle dynamics (DPD) which has broad applicability in cell biophysics with implications for diagnostics, therapeutics and drug efficacy assessments for a wide variety of human diseases. This computational approach, validated by independent experimental results, is capable of modeling the biorheology of whole blood and its individual components during blood flow so as to investigate cell mechanistic processes in health and disease. DPD is a Lagrangian method that can be derived from systematic coarse-graining of molecular dynamics but can scale efficiently up to arterioles and can also be used to model RBCs down to the spectrin level. We start from experimental measurements of a single RBC to extract the relevant biophysical parameters, using single-cell measurements involving such methods as optical tweezers, atomic force microscopy and micropipette aspiration, and cell-population experiments involving microfluidic devices. We then use these validated RBC models to predict the biorheological behavior of whole blood in healthy or pathological states, and compare the simulations with experimental results involving apparent viscosity and other relevant parameters. While the approach discussed here is sufficiently general to address a broad spectrum of hematologic disorders including certain types of cancer, this paper specifically deals with results obtained using this computational framework for blood flow in malaria and sickle cell anemia.

Cardiovascular pathologies, such as a brain aneurysm, are affected by the global blood circulation as well as by the local microrheology. Hence, developing computational models for such cases requires the coupling of disparate spatial and temporal scales often governed by diverse mathematical descriptions, e.g., by partial differential equations (continuum) and ordinary differential equations for discrete particles (atomistic). However, interfacing atomistic-based with continuum-based domain discretizations is a challenging problem that requires both mathematical and computational advances. We present here a hybrid methodology that enabled us to perform the first multi-scale simulations of platelet depositions on the wall of a brain aneurysm. The large scale flow features in the intracranial network are accurately resolved by using the high-order spectral element Navier-Stokes solver εκ αr . The blood rheology inside the aneurysm is modeled using a coarse-grained stochastic molecular dynamics approach (the dissipative particle dynamics method) implemented in the parallel code LAMMPS. The continuum and atomistic domains overlap with interface conditions provided by effective forces computed adaptively to ensure continuity of states across the interface boundary. A two-way interaction is allowed with the time-evolving boundary of the (deposited) platelet clusters tracked by an immersed boundary method. The corresponding heterogeneous solvers ( εκ αr and LAMMPS) are linked together by a computational multilevel message passing interface that facilitates modularity and high parallel efficiency. Results of multiscale simulations of clot formation inside the aneurysm in a patient-specific arterial tree are presented. We also discuss the computational challenges involved and present scalability results of our coupled solver on up to 300K computer processors. Validation of such coupled atomistic-continuum models is a main open issue that has to be addressed in future work.

In small vessels blood is usually treated as a Newtonian fluid down to diameters of ~200 μm. We investigate the flow of red blood cell (RBC) suspensions driven through small tubes (diameters 10-150 μm) in the range marking the transition from arterioles and venules to the largest capillary vessels. The results of the simulations combined with previous simulations of uniform shear flow and experimental data show that for diameters less than ~100 μm the suspension's stress cannot be described as a continuum, even a heterogeneous one. We employ the dissipative particle dynamics (DPD) model, which has been successfully used to predict human blood bulk viscosity in homogeneous shear flow. In tube flow the cross-stream stress gradient induces an inhomogeneous distribution of RBCs featuring a centreline cell density peak, and a cell-free layer (CFL) next to the wall. For a neutrally buoyant suspension the imposed linear shear-stress distribution together with the differentiable velocity distribution allow the calculation of the local viscosity across the tube section. The viscosity across the section as a function of the strain rate is found to be essentially independent of tube size for the larger diameters and is determined by the local haematocrit ( H ) and shear rate. Other RBC properties such as asphericity, deformation, and cell-flow orientation exhibit similar dependence for the larger tube diameters. As the tube size decreases below ~100 μm in diameter, the viscosity in the central region departs from the large-tube similarity function of the shear rate, since H increases significantly towards the centreline. The dependence of shear stress on tube size, in addition to the expected local shear rate and local haematocrit, implies that blood flow in small tubes cannot be described as a heterogeneous continuum. Based on the analysis of the DPD simulations and on available experimental results, we propose a simple velocity-slip model that can be used in conjunction with continuum-based simulations.

We investigate structural and dynamical properties of ultra-soft colloids in dilute and semi-dilute solutions by hybrid mesoscale simulations under linear shear flow. In particular, the influence of functionality on these properties is addressed. Our study combines molecular dynamics simulations for the solute with the multiparticle collision dynamics approach for the coarse-grained solvent. The star polymers exhibit large conformational and orientational changes in shear flow, which we characterize by the radius of gyration tensor and the alignment angle. These quantities show a universal dependence on a concentration- and functionality-dependent Weissenberg number with slight deviations at high shear rates. Moreover, the star polymers display a rotational dynamics with a shear-rate- and concentration-dependent rotation frequency. We attribute the concentration dependence to the screening of hydrodynamic interactions in semi-dilute star-polymer solutions.

Margination of white blood cells (WBCs) towards vessel walls is an essential precondition for their efficient adhesion to the vascular endothelium. We perform numerical simulations with a two-dimensional blood flow model to investigate the dependence of WBC margination on hydrodynamic interactions of blood cells with the vessel walls, as well as on their collective behavior and deformability. We find WBC margination to be optimal in intermediate ranges of red blood cell (RBC) volume fractions and flow rates, while, beyond these ranges, it is substantially attenuated. RBC aggregation enhances WBC margination, while WBC deformability reduces it. These results are combined in state diagrams, which identify WBC margination for a wide range of flow and cell suspension conditions.