Shape transitions of fluid vesicles and red blood cells in capillary flows

Numerical insights on the volume phase transition of thermoresponsive hollow microgels

Hollow microgels, consisting of a pNIPAM polymer network with a central cavity, have significant potential due to their tunable softness and encapsulation capabilities. Using molecular dynamics simulations, we thoroughly characterise the swelling behaviour of neutral hollow microgels across the volume phase transition (VPT) upon varying crosslinker concentration, shell thickness, and size. In particular, we examine in detail the onset of cavity filling and its relation to the VPT, detecting the presence of a discontinuity in the radius of gyration of the microgels, if an appropriate balance between shell stiffness and thermoresposiveness is reached. The discontinuity is, however, absent in the behaviour of the hydrodynamic radius, in agreement with experimental observations. We then test our numerical model by direct comparison of form factors with available measurements in the literature and also establish a minimal-size, stable hollow microgel for future computationally feasible bulk investigations. Overall, our findings provide valuable insights into the fundamental swelling properties of hollow microgels that can be useful to control the opening and closing of the cavity for application purposes.

Renormalized mechanics and stochastic thermodynamics of growing model protocells

Uncovering the rules governing the nonequilibrium dynamics of the membranes that define biological cells is of central importance to understanding the physics of living systems. We theoretically and computationally investigate the behavior of model protocells-flexible quasispherical vesicles-that exchange membrane constituents, internal volume, and heat with an external reservoir. The excess chemical potential and osmotic pressure difference imposed by the reservoir act as generalized thermodynamic driving forces that modulate vesicle morphology. We identify an associated nonequilibrium morphological transition between a weakly driven regime, in which growing vesicles remain quasispherical, and a strongly driven regime, in which vesicles accommodate rapid membrane uptake by developing surface wrinkles. This transition emerges due to the renormalization of membrane mechanical properties by nonequilibrium driving. Further, using insights from stochastic thermodynamics we propose a minimal vesicle growth-shape law that remains robust even in strongly driven, far-from-equilibrium regimes.

A Synergistic Overview between Microfluidics and Numerical Research for Vascular Flow and Pathological Investigations

Vascular diseases are widespread, and sometimes such life-threatening medical disorders cause abnormal blood flow, blood particle damage, changes to flow dynamics, restricted blood flow, and other adverse effects. The study of vascular flow is crucial in clinical practice because it can shed light on the causes of stenosis, aneurysm, blood cancer, and many other such diseases, and guide the development of novel treatments and interventions. Microfluidics and computational fluid dynamics (CFDs) are two of the most promising new tools for investigating these phenomena. When compared to conventional experimental methods, microfluidics offers many benefits, including lower costs, smaller sample quantities, and increased control over fluid flow and parameters. In this paper, we address the strengths and weaknesses of computational and experimental approaches utilizing microfluidic devices to investigate the rheological properties of blood, the forces of action causing diseases related to cardiology, provide an overview of the models and methodologies of experiments, and the fabrication of devices utilized in these types of research, and portray the results achieved and their applications. We also discuss how these results can inform clinical practice and where future research should go. Overall, it provides insights into why a combination of both CFDs, and experimental methods can give even more detailed information on disease mechanisms recreated on a microfluidic platform, replicating the original biological system and aiding in developing the device or chip itself.

Delivery Strategy to Enhance the Therapeutic Efficacy of Liver Fibrosis via Nanoparticle Drug Delivery Systems

Liver fibrosis (LF) is a pathological repair reaction caused by a chronic liver injury that affects the health of millions of people worldwide, progressing to life-threatening cirrhosis and liver cancer without timely intervention. Due to the complexity of LF pathology, multiple etiological characteristics, and the deposited extracellular matrix, traditional drugs cannot reach appropriate targets in a time-space matching way, thus decreasing the therapeutic effect. Nanoparticle drug delivery systems (NDDS) enable multidrug co-therapy and develop multifactor delivery strategies targeting pathological processes, showing great potential in LF therapy. Based on the pathogenesis and the current clinical treatment status of LF, we systematically elucidate the targeting mechanism of NDDS used in the treatment of LF. Subsequently, we focus on the progress of drug delivery applications for LF, including combined delivery for the liver fibrotic pathological environment, overcoming biological barriers, precise intracellular regulation, and intelligent responsive delivery for the liver fibrotic microenvironment. We hope that this review will inspire the rational design of NDDS for LF in the future in order to provide ideas and methods for promoting LF regression and cure.

Erythrocyte metabolism

Our aim is to present an updated overview of the erythrocyte metabolism highlighting its richness and complexity. We have manually collected and connected the available biochemical pathways and integrated them into a functional metabolic map. The focus of this map is on the main biochemical pathways consisting of glycolysis, the pentose phosphate pathway, redox metabolism, oxygen metabolism, purine/nucleoside metabolism, and membrane transport. Other recently emerging pathways are also curated, like the methionine salvage pathway, the glyoxalase system, carnitine metabolism, and the lands cycle, as well as remnants of the carboxylic acid metabolism. An additional goal of this review is to present the dynamics of erythrocyte metabolism, providing key numbers used to perform basic quantitative analyses. By synthesizing experimental and computational data, we conclude that glycolysis, pentose phosphate pathway, and redox metabolism are the foundations of erythrocyte metabolism. Additionally, the erythrocyte can sense oxygen levels and oxidative stress adjusting its mechanics, metabolism, and function. In conclusion, fine-tuning of erythrocyte metabolism controls one of the most important biological processes, that is, oxygen loading, transport, and delivery.

Classification of chemically modified red blood cells in microflow using machine learning video analysis

We classify native and chemically modified red blood cells with an AI based video classifier. Using TensorFlow video analysis enables us to capture not only the morphology of the cell but also the trajectories of motion of individual red blood cells and their dynamics. We chemically modify cells in three different ways to model different pathological conditions and obtain classification accuracies for all three classification tasks of more than 90% between native and modified cells. Unlike standard cytometers that are based on immunophenotyping our microfluidic cytometer allows to rapidly categorize cells without any fluorescence labels simply by analysing the shape and flow of red blood cells.

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.

The neurovascular unit of capillary blood vessels in the rat nervous system. A rapid-Golgi electron microscopy study

We describe a pericapillary organ in the rat forebrain and cerebellar cortex. It consists of a series of tripartite synapses with synaptic extensions enveloped by astrocytic endfeet that are linked to the capillary wall by synaptic extensions. Reciprocal specializations of the pericyte-capillary blood vessel (CBV) with such specialized synapses suggest a mechanoreceptor role. In Golgi-impregnated and 3D reconstructions of the cerebral cortex and thalamus, a series of TSs appear to be sequentially ordered in a common dendrite, paralleled by synaptic outgrowths termed golf club synaptic extensions (GCE) opposed to a longitudinal crest (LC) from the capillary basal lamina (BL). Our results show that, in the cerebellar cortex, afferent fibers and interneurons display microanatomical structures that strongly suggest an interaction with the capillary wall. Afferent mossy fiber (MF) rosettes and ascending granule cell axons and their dendrites define the pericapillary passage interactions that are entangled by endfeet. The presence of mRNA of the mechanosensitive channel Piezo1 in the MF rosettes, together with the surrounding end-feet and the capillary wall form mechanosensory units. The ubiquity of such units to modulate synaptic transmission is also supported by Piezo1 mRNA expressing pyramidal isocortical and thalamic neurons. This scenario suggests that ascending impulses to the cerebellar and cortical targets are presynaptically modulated by the reciprocal interaction with the mechanosensory pericapillary organ that ultimately modulates the vasomotor response.

Direct Blood Cell Flow Imaging in Microvascular Networks

Blood flow dynamics in microvascular networks are intimately related to the health of tissues and organs. While numerous imaging modalities and techniques have been developed to assess blood flow dynamics for various applications, their utilization has been hampered by limited imaging speed and indirect quantification of blood flow dynamics. Here, direct blood cell flow imaging (DBFI) is demonstrated that provides visualization of individual motions of blood cells over a field of 0.71 mm × 1.42 mm with a time resolution of 0.69 ms (1450 frames s-1 ) without using any exogenous agents. DBFI enables precise dynamic analysis of blood cell flow velocities and fluxes in various vessels over a large field, from capillaries to arteries and veins, with unprecedented time resolution. Three exemplary applications of DBFI, quantification of blood flow dynamics of 3D vascular networks, analysis of heartbeat induced blood flow dynamics, and analysis of blood flow dynamics of neurovascular coupling, illustrate the potential of this new imaging technology.

Shape transitions in a network model of active elastic shells

Morphogenesis involves the transformation of initially simple shapes, such as multicellular spheroids, into more complex 3D shapes. These shape changes are governed by mechanical forces including molecular motor-generated forces as well as hydrostatic fluid pressure, both of which are actively regulated in living matter through mechano-chemical feedback. Inspired by autonomous, biophysical shape change, such as occurring in the model organism hydra, we introduce a minimal, active, elastic model featuring a network of springs in a globe-like spherical shell geometry. In this model there is coupling between activity and the shape of the shell: if the local curvature of a filament represented by a spring falls below a critical value, its elastic constant is actively changed. This results in deformation of the springs that changes the shape of the shell. By combining excitation of springs and pressure regulation, we show that the shell undergoes a transition from spheroidal to either elongated ellipsoidal or a different spheroidal shape, depending on pressure. There exists a critical pressure at which there is an abrupt change from ellipsoids to spheroids, showing that pressure is potentially a sensitive switch for material shape. We thus offer biologically inspired design principles for autonomous shape transitions in active elastic shells.

AI based image analysis of red blood cells in oscillating microchannels

The flow dynamics of red blood cells in vivo in blood capillaries and in vitro in microfluidic channels is complex. Cells can obtain different shapes such as discoid, parachute, slipper-like shapes and various intermediate states depending on flow conditions and their viscoelastic properties. We use artificial intelligence based analysis of red blood cells (RBCs) in an oscillating microchannel to distinguish healthy red blood cells from red blood cells treated with formaldehyde to chemically modify their viscoelastic behavior. We used TensorFlow to train and validate a deep learning model and achieved a testing accuracy of over 97%. This method is a first step to a non-invasive, label-free characterization of diseased red blood cells and will be useful for diagnostic purposes in haematology labs. This method provides quantitative data on the number of affected cells based on single cell classification.

Multi-particle collision dynamics for a coarse-grained model of soft colloids applied to ring polymers

The simulation of polymer solutions often requires the development of methods that accurately include hydrodynamic interactions. Resolution on the atomistic scale is too computationally expensive to cover mesoscopic time and length scales on which the interesting polymer phenomena are observed. Therefore, coarse-graining methods have to be applied. In this work, the solvent is simulated using the well-established multi-particle collision dynamics scheme, and for the polymer, different coarse-graining methods are employed and compared against the monomer resolved Kremer-Grest model by their resulting diffusion coefficients. This research builds on previous work [Ruiz-Franco et al., J. Chem. Phys. 151, 074902 (2019)], in which star polymers and linear chains in a solvent were simulated and two different coarse-graining methods were developed, in order to increase computational efficiency. The present work extends this approach to ring polymers and seeks to refine one of the authors' proposed model: the penetrable soft colloid model. It was found that both proposed models are not well suited to ring polymers; however, the introduction of a factor to the PSC model delivers satisfying results for the diffusion behavior by regulating the interaction intensity with the solvent.

Shape Transitions of Red Blood Cell under Oscillatory Flows in Microchannels

This paper aims to examine the ability to control Red Blood Cell (RBCs) dynamics and the associated extracellular flow patterns in microfluidic channels via oscillatory flows. Our computational approach employs a hybrid continuum-particle coupling, in which the cell membrane and cytosol fluid are modeled using the Dissipative Particle Dynamics (DPD) method. The blood plasma is modeled as an incompressible fluid via the Immersed Boundary Method (IBM). This coupling is novel because it provides an accurate description of RBC dynamics while the extracellular flow patterns around the RBCs are also captured in detail. Our coupling methodology is validated with available experimental and computational data in the literature and shows excellent agreement. We explore the controlling regimes by varying the shape of the oscillatory flow waveform at the channel inlet. Our simulation results show that a host of RBC morphological dynamics emerges depending on the channel geometry, the incoming flow waveform, and the RBC initial location. Complex dynamics of RBC are induced by the flow waveform. Our results show that the RBC shape is strongly dependent on its initial location. Our results suggest that the controlling of oscillatory flows can be used to induce specific morphological shapes of RBCs and the surrounding fluid patterns in bio-engineering applications.

Biomechanical Investigation of Red Cell Sedimentation Using Blood Shear Stress and Blood Flow Image in a Capillary Chip

Blood image intensity has been used to detect erythrocyte sedimentation rate (ESR). However, it does not give information on the biophysical properties of blood samples under continuous ESR. In this study, to quantify mechanical variations of blood under continuous ESR, blood shear stress and blood image intensity were obtained by analyzing blood flows in the capillary channel. A blood sample is loaded into a driving syringe to demonstrate the proposed method. The blood flow rate is set in a periodic on-off pattern. A blood sample is then supplied into a capillary chip, and microscopic blood images are captured at specific intervals. Blood shear stress is quantified from the interface of the bloodstream in the coflowing channel. τ0 is defined as the maximum shear stress obtained at the first period. Simultaneously, ESRτ is then obtained by analyzing temporal variations of blood shear stress for every on period. AII is evaluated by analyzing the temporal variation of blood image intensity for every off period. According to the experimental results, a shorter period of T = 4 min and no air cavity contributes to the high sensitivity of the two indices (ESRτ and AII). The τ0 exhibits substantial differences with respect to hematocrits (i.e., 30-50%) as well as diluents. The ESRτ and AII showed a reciprocal relationship with each other. Three suggested properties represented substantial differences for suspended blood samples (i.e., hardened red blood cells, different concentrations of dextran solution, and fibrinogen). In conclusion, the present method can detect variations in blood samples under continuous ESR effectively.

Dynamics of fluid bilayer vesicles: Soft meshes and robust curvature energy discretization

Continuum models like the Helfrich Hamiltonian are widely used to describe fluid bilayer vesicles. Here we study the molecular dynamics compatible dynamics of the vertices of two-dimensional meshes representing the bilayer, whose in-plane motion is only weakly constrained. We show (i) that Jülicher's discretization of the curvature energy offers vastly superior robustness for soft meshes compared to the commonly employed expression by Gommper and Kroll and (ii) that for sufficiently soft meshes, the typical behavior of fluid bilayer vesicles can emerge even if the mesh connectivity remains fixed throughout the simulations. In particular, soft meshes can accommodate large shape transformations, and the model can generate the typical ℓ^{-4} signal for the amplitude of surface undulation modes of nearly spherical vesicles all the way up to the longest wavelength modes. Furthermore, we compare results for Newtonian, Langevin, and Brownian dynamics simulations of the mesh vertices to demonstrate that the internal friction of the membrane model is negligible, making it suitable for studying the internal dynamics of vesicles via coupling to hydrodynamic solvers or particle-based solvent models.

Recent advances in micro-sized oxygen carriers inspired by red blood cells

Supplementing sufficient oxygen to cells is always challenging in biomedical engineering fields such as tissue engineering. Originating from the concept of a 'blood substitute', nano-sized artificial oxygen carriers (AOCs) have been studied for a long time for the optimization of the oxygen supplementation and improvement of hypoxia environments in vitro and in vivo. When circulating in our bodies, micro-sized human red blood cells (hRBCs) feature a high oxygen capacity, a unique biconcave shape, biomechanical and rheological properties, and low frictional surfaces, making them efficient natural oxygen carriers. Inspired by hRBCs, recent studies have focused on evolving different AOCs into microparticles more feasibly able to achieve desired architectures and morphologies and to obtain the corresponding advantages. Recent micro-sized AOCs have been developed into additional categories based on their principal oxygen-carrying or oxygen-releasing materials. Various biomaterials such as lipids, proteins, and polymers have also been used to prepare oxygen carriers owing to their rapid oxygen transfer, high oxygen capacity, excellent colloidal stability, biocompatibility, suitable biodegradability, and long storage. In this review, we concentrated on the fabrication techniques, applied biomaterials, and design considerations of micro-sized AOCs to illustrate the advances in their performances. We also compared certain recent micro-sized AOCs with hRBCs where applicable and appropriate. Furthermore, we discussed existing and potential applications of different types of micro-sized AOCs.

A Review of Continuum Mechanics for Mechanical Deformation of Lipid Membranes

Mechanical deformation of lipid membranes plays important roles in various cellular tasks. Curvature deformation and lateral stretching are two major energy contributions to the mechanical deformation of lipid membranes. In this paper, continuum theories for these two major membrane deformation events were reviewed. Theories based on curvature elasticity and lateral surface tension were introduced. Numerical methods as well as biological applications of the theories were discussed.

The stress-free state of human erythrocytes: Data-driven inference of a transferable RBC model

The stress-free state (SFS) of red blood cells (RBCs) is a fundamental reference configuration for the calibration of computational models, yet it remains unknown. Current experimental methods cannot measure the SFS of cells without affecting their mechanical properties, whereas computational postulates are the subject of controversial discussions. Here, we introduce data-driven estimates of the SFS shape and the visco-elastic properties of RBCs. We employ data from single-cell experiments that include measurements of the equilibrium shape of stretched cells and relaxation times of initially stretched RBCs. A hierarchical Bayesian model accounts for these experimental and data heterogeneities. We quantify, for the first time, the SFS of RBCs and use it to introduce a transferable RBC (t-RBC) model. The effectiveness of the proposed model is shown on predictions of unseen experimental conditions during the inference, including the critical stress of transitions between tumbling and tank-treading cells in shear flow. Our findings demonstrate that the proposed t-RBC model provides predictions of blood flows with unprecedented accuracy and quantified uncertainties.

Dynamic response of red blood cells in health and disease

The viscoelastic response of the red blood cells (RBCs) affected by hematological disorders become severely impaired by the altered biophysical and morphological properties. These include traits like reduced deformability, increased membrane viscosity, and change in cell shape, causing substantial changes in the overall hemodynamics. RBCs, by virtue of their highly elastic membrane and low bending rigidity, exhibit complex dynamics when exposed to cyclic, transient forces in the microcirculation. Here, we employ mesoscopic numerical simulations based on the dissipative particle dynamics (DPD) framework to explore the dynamics of healthy, schizont stage malaria-infected and type 2 diabetes mellitus affected RBCs subjected to external time-dependent loads. The paper focuses on the imposition and cessation of external forcing on the cells of two different typologies, saw-tooth cyclic wave loading and sudden loads in the form of creep and relaxation phenomena. The effects of varying the rate of stress and the applied stress magnitude were investigated. Our simulations disclosed unique shape transitions of the hysteresis curves at varied loading rates. A careful analysis reveals a critical threshold of half cycle time of the from wherein the deformation of all cells observed, healthy or otherwise, falls under the nearly reversible deformation regime displaying minimal energy dissipation. Finally, we also examined the individual effects of the different constitutive and geometric characteristics attributed to the pathological cells and observed interesting recovery dynamics of spherocytes and cells having high shear moduli. The distinguished deformation behaviour of healthy and diseased cells could establish external force as a valuable initial biomarker.

Impact of interfacial rheology on finger tip splitting

Fluid-fluid interfaces, laden with polymers, surfactants, lipid bilayers, proteins, solid particles, or other surface-active agents, often exhibit a rheologically complex response to deformations. Despite its academic and practical relevance to fluid dynamics and various other fields of research, the role of interfacial rheology in viscous fingering remains fairly underexplored. A noteworthy exception is the work by Li and Manikantan [Phys. Rev. Fluids 6, 074001 (2021)2469-990X10.1103/PhysRevFluids.6.074001], who used linear stability analysis to show that surface rheological stresses act to stabilize the development of radial viscous fingering at the linear regime. In this paper, we perform a perturbative, second-order mode-coupling analysis of the system and investigate the influence of interfacial rheology on the morphology of the fingering structures at early nonlinear stages of the dynamics. In particular, we focus on understanding how interfacial rheology impacts the emblematic finger tip-widening and finger tip-splitting phenomena that take place in radial viscous fingering in Hele-Shaw cells. We describe the viscous Newtonian fluid-fluid interface by using a Boussinesq-Scriven model, and derive a generalized Young-Laplace pressure jump condition at the fluid-fluid interface. In this framing, we go beyond the purely linear description and use Darcy's law to obtain a perturbative mode-coupling differential equation which describes the time evolution of the perturbation amplitudes, accurate to second order. Our early nonlinear mode-coupling results indicate that regardless of their stabilizing action at the linear regime, interfacial rheology effects favor finger tip widening, leading to the occurrence of enhanced finger tip-splitting events.

Circulating cell clusters aggravate the hemorheological abnormalities in COVID-19

Microthrombi and circulating cell clusters are common microscopic findings in patients with coronavirus disease 2019 (COVID-19) at different stages in the disease course, implying that they may function as the primary drivers in disease progression. Inspired by a recent flow imaging cytometry study of the blood samples from patients with COVID-19, we perform computational simulations to investigate the dynamics of different types of circulating cell clusters, namely white blood cell (WBC) clusters, platelet clusters, and red blood cell clusters, over a range of shear flows and quantify their impact on the viscosity of the blood. Our simulation results indicate that the increased level of fibrinogen in patients with COVID-19 can promote the formation of red blood cell clusters at relatively low shear rates, thereby elevating the blood viscosity, a mechanism that also leads to an increase in viscosity in other blood diseases, such as sickle cell disease and type 2 diabetes mellitus. We further discover that the presence of WBC clusters could also aggravate the abnormalities of local blood rheology. In particular, the extent of elevation of the local blood viscosity is enlarged as the size of the WBC clusters grows. On the other hand, the impact of platelet clusters on the local rheology is found to be negligible, which is likely due to the smaller size of the platelets. The difference in the impact of WBC and platelet clusters on local hemorheology provides a compelling explanation for the clinical finding that the number of WBC clusters is significantly correlated with thrombotic events in COVID-19 whereas platelet clusters are not. Overall, our study demonstrates that our computational models based on dissipative particle dynamics can serve as a powerful tool to conduct quantitative investigation of the mechanism causing the pathological alterations of hemorheology and explore their connections to the clinical manifestations in COVID-19.

Mem3DG: Modeling membrane mechanochemical dynamics in 3D using discrete differential geometry

Biomembranes adopt varying morphologies that are vital to cellular functions. Many studies use computational modeling to understand how various mechanochemical factors contribute to membrane shape transformations. Compared with approximation-based methods (e.g., finite element method [FEM]), the class of discrete mesh models offers greater flexibility to simulate complex physics and shapes in three dimensions; its formulation produces an efficient algorithm while maintaining coordinate-free geometric descriptions. However, ambiguities in geometric definitions in the discrete context have led to a lack of consensus on which discrete mesh model is theoretically and numerically optimal; a bijective relationship between the terms contributing to both the energy and forces from the discrete and smooth geometric theories remains to be established. We address this and present an extensible framework, Mem3DG, for modeling 3D mechanochemical dynamics of membranes based on discrete differential geometry (DDG) on triangulated meshes. The formalism of DDG resolves the inconsistency and provides a unifying perspective on how to relate the smooth and discrete energy and forces. To demonstrate, Mem3DG is used to model a sequence of examples with increasing mechanochemical complexity: recovering classical shape transformations such as 1) biconcave disk, dumbbell, and unduloid; and 2) spherical bud on spherical, flat-patch membrane; investigating how the coupling of membrane mechanics with protein mobility jointly affects phase and shape transformation. As high-resolution 3D imaging of membrane ultrastructure becomes more readily available, we envision Mem3DG to be applied as an end-to-end tool to simulate realistic cell geometry under user-specified mechanochemical conditions.

Continuum modeling of non-conservative fluid membrane for simulating long-term cell dynamics

Living cells actively deform and move by their force generations in three-dimensional (3D) space. These 3D cell dynamics occur over a long-term time scale, ranging from tens of minutes to days. On such a time scale, turnover of cell membrane constituents due to endocytosis and exocytosis cannot be ignored, i.e., the surface membrane dynamically deforms without mass conservation. Although membrane turnover is essential for large deformation of cells, there is no computational framework yet to simulate long-term cell dynamics with a non-conservative fluidic membrane. In this paper, we proposed a computational framework for simulating the long-term dynamics of a cell membrane in 3D space. For this purpose, in the proposed framework, the cell surface membrane is treated as a viscous fluid membrane without mass conservation. Cell shape is discretized by a triangular mesh, and its dynamics are expressed by effective energy and dissipation function. The mesh structure, distorted by membrane motion, is dynamically optimized by introducing a modified dynamic remeshing method. To validate the proposed framework, numerical simulations were performed, showing that the membrane flow is reproduced in a physically consistent manner and that the artificial effects of the remeshing method were negligible. To further demonstrate the applicability of the proposed framework, numerical simulations of cell migration induced by a mechanism similar to the Marangoni effect, i.e., the polarized surface tension actively generated by the cell, were performed. The observed cell behaviors agreed with existing analytical solutions, indicating that the proposed computational framework can quantitatively reproduce long-term active cell dynamics with membrane turnover. Based on the simple description of cell membrane dynamics, this framework provides a useful basis for analyzing various cell shaping and movement.

Development and validation of quantitative optical index of skin blood content

SIGNIFICANCE

We present an approach to estimate with simple instrumentation the amount of red blood cells in the skin microvasculature, designated as parameter LRBC. Variations of parameter LRBC are shown to reflect local changes in the quantity of skin red blood cells during a venous occlusion challenge.

AIM

To validate a simple algebraic model of light transport in skin using the Monte Carlo method and to develop a measure of the red blood cell content in skin microvessels using the Monte Carlo predictions; to guide the development of an instrument to measure experimentally variations of the amount of red blood cells in the skin.

APPROACH

Monte Carlo simulations were carried out in a multilayer model of the skin to compute remitted light intensities as a function of distance from the illumination locus for different values of the skin blood content. The simulation results were used to compute parameter LRBC and its variations with local skin blood content. An experimental setup was developed to measure parameter LRBC in human volunteers in whom skin blood content of the forearm increased during temporary interruption of the venous outflow.

RESULTS

In the simulations, parameter LRBC was ∼16  μm in baseline conditions, and it increased in near proportion with the blood content of the skin layers. Measuring the diffusely reflected light intensity 0.5 to 1.2 mm away from the illumination locus was optimal to detect appreciable changes of the reflected light intensity as skin blood content was altered. Parameter LRBC measured experimentally on the human forearm was 17  ±  2  μm in baseline conditions it increased at a rate of 4  ±  2  μm  /  min when venous outflow was temporarily interrupted.

CONCLUSION

Parameter LRBC derived experimentally with a two-wavelength diffuse reflectometer can be used to measure local variations of the amount of red blood cells in skin microvessels.

Shape transformations of red blood cells in the capillary and their possible connections to oxygen transportation

In this work, a series of numerical simulations have been performed to obtain the steady shapes of red blood cells under a shear force field in the capillary. Two possible classes of steady shapes, the axisymmetric parachute and the non-axisymmetric parachute, are found. If we assume that oxygen diffusion across the red cell membrane is mediated by membrane curvature, it is found that the non-axisymmetric parachute will be more favorable due to its special shape which enables it to have a larger portion of membrane patch capable of releasing oxygen to tissues.

The effect of stiffened diabetic red blood cells on wall shear stress in a reconstructed 3D microaneurysm

Blood flow within the vasculature of the retina has been found to influence the progression of diabetic retinopathy. In this research cell resolved blood flow simulations are used to study the pulsatile flow of whole blood through a segmented retinal microaneurysm. Images were collected using adaptive optics optical coherence tomography of the retina of a patient with diabetic retinopathy, and a sidewall (sacciform) microaneurysm was segmented from the volumetric data. The original microaneurysm neck width was varied to produce two additional aneurysm geometries in order to probe the influence of neck width on the transport of red blood cells and platelets into the aneurysm. Red blood cell membrane stiffness was also increased to resolve the impact of rigid red blood cells, as a result of diabetes, in blood flow. Wall shear stress and wall shear stress gradients were calculated throughout the aneurysm domains, and the quantification of the influence of the red blood cells is presented. Average wall shear stress and wall shear stress gradients increased due to the increase of red blood cell membrane stiffness. Stiffened red blood cells were also found to induce higher local wall shear stress and wall shear stress gradients as they passed through the leading and draining parental vessels. Stiffened red blood cells were found to penetrate the aneurysm sac more than healthy red blood cells, as well as decreasing the margination of platelets to the vessel walls of the parental vessel, which caused a decrease in platelet penetration into the aneurysm sac.

Erythrocyte sedimentation: Effect of aggregation energy on gel structure during collapse

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.

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.

Nonlinear Transient and Steady State Stretching of Deflated Vesicles in Flow

Membrane-bound vesicles and organelles exhibit a wide array of nonspherical shapes at equilibrium, including biconcave and tubular morphologies. Despite recent progress, the stretching dynamics of deflated vesicles is not fully understood, particularly far from equilibrium where complex nonspherical shapes undergo large deformations in flow. Here, we directly observe the transient and steady-state nonlinear stretching dynamics of deflated vesicles in extensional flow using a Stokes trap. Automated flow control is used to observe vesicle dynamics over a wide range of flow rates, shape anisotropy, and viscosity contrast. Our results show that deflated vesicle membranes stretch into highly deformed shapes in flow above a critical capillary number Ca_c1. We further identify a second critical capillary number Ca_c2, above which vesicle stretch diverges in flow. Vesicles are robust to multiple nonlinear stretch-relax cycles, evidenced by relaxation of dumbbell-shaped vesicles containing thin lipid tethers following flow cessation. An analytical model is developed for vesicle deformation in flow, which enables comparison of nonlinear steady-state stretching results with theories for different reduced volumes. Our results show that the model captures the steady-state stretching of moderately deflated vesicles; however, it underpredicts the steady-state nonlinear stretching of highly deflated vesicles. Overall, these results provide a new understanding of the nonlinear stretching dynamics and membrane mechanics of deflated vesicles in flow.

Cell/Bacteria-Based Bioactive Materials for Cancer Immune Modulation and Precision Therapy

Numerous clinical trials for cancer precision medicine research are limited due to the drug resistance, side effects, and low efficacy. Unsatisfactory outcomes are often caused by complex physiologic barriers and abnormal immune events in tumors, such as tumor target alterations and immunosuppression. Cell/bacteria-derived materials with unique bioactive properties have emerged as attractive tools for personalized therapy in cancer. Naturally derived bioactive materials, such as cell and bacterial therapeutic agents with native tropism or good biocompatibility, can precisely target tumors and effectively modulate immune microenvironments to inhibit tumors. Here, the recent advances in the development of cell/bacteria-based bioactive materials for immune modulation and precision therapy in cancer are summarized. Cell/bacterial constituents, including cell membranes, bacterial vesicles, and other active substances have inherited their unique targeting properties and antitumor capabilities. Strategies for engineering living cell/bacteria to overcome complex biological barriers and immunosuppression to promote antitumor efficacy are also summarized. Moreover, past and ongoing trials involving personalized bioactive materials and promising agents such as cell/bacteria-based micro/nano-biorobotics are further discussed, which may become another powerful tool for treatment in the near future.

EFFECT OF ARTERIOLAR DISTENSIBILITY ON THE LATERAL MIGRATION OF LIQUID-FILLED MICROPARTICLES FLOWING IN A HUMAN ARTERIOLE

A promising advance of bioengineering consists in the development of micro-nanoparticles as drug delivery vehicles injected intravenously or intraarterialy for targeted treatment. Proficient functioning of drug carries is conditioned by a reliable prediction of pharmacokinetics in human as well as their dynamical behavior once injected in blood stream. In this study, we aim to provide a reliable numerical prediction of dynamical behavior of microparticles in human arteriole focusing on the crucial mechanism of lateral migration. The dynamical response of the microparticle upon blood flow and arteriolar distensibility is investigated by varying main controlling parameters: viscosity ratio, confinement and capillary number. The influence of the hyperelastic arteriolar wall is highlighted through comparison with an infinitely rigid arteriolar wall. The hydrodynamic interaction in a microparticle train is examined. Fluid–structure interaction is solved by the Arbitrary Lagrangian–Eulerian method using the COMSOL Multiphysics software.

Multi-particle collision dynamics with a non-ideal equation of state. II. Collective dynamics of elongated squirmer rods

Simulations of flow fields around microscopic objects typically require methods that both solve the Navier-Stokes equations and also include thermal fluctuations. One such method popular in the field of soft-matter physics is the particle-based simulation method of multi-particle collision dynamics (MPCD). However, in contrast to the typically incompressible real fluid, the fluid of the traditional MPCD methods obeys the ideal-gas equation of state. This can be problematic because most fluid properties strongly depend on the fluid density. In a recent article, we proposed an extended MPCD algorithm and derived its non-ideal equation of state and an expression for the viscosity. In the present work, we demonstrate its accuracy and efficiency for the simulations of the flow fields of single squirmers and of the collective dynamics of squirmer rods. We use two exemplary squirmer-rod systems for which we compare the outcome of the extended MPCD method to the well-established MPCD version with an Andersen thermostat. First, we explicitly demonstrate the reduced compressibility of the MPCD fluid in a cluster of squirmer rods. Second, for shorter rods, we show the interesting result that in simulations with the extended MPCD method, dynamic swarms are more pronounced and have a higher polar order. Finally, we present a thorough study of the state diagram of squirmer rods moving in the center plane of a Hele-Shaw geometry. From a small to large aspect ratio and density, we observe a disordered state, dynamic swarms, a single swarm, and a jammed cluster, which we characterize accordingly.

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.

Energetic Lifestyle Drives Size and Shape of Avian Erythrocytes

The size and shape of red blood cells (erythrocytes) is determined by key life history strategies in vertebrates. They have a fundamental role to deliver oxygen to tissues, and their ability to do so is shaped by the tissue's need and their shape. Despite considerable interest in how other components of blood are shaped by ecology and life history, few studies have considered erythrocytes themselves. We tested how erythrocyte size and shape varied in relation to energetically demanding activities using a dataset of 631 bird species. We found that in general, birds undergoing greater activities such as long distance migration had smaller and more elongated cells, while those with greater male-male competition had smaller and rounder cells. Smaller, more elongated erythrocytes allow more rapid oxygenation/deoxygenation and support greater aerobic activity. The rounder erythrocytes found in species with strong male–male competition may stem from younger erythrocytes deriving from androgen-induced erythropoiesis rates. Finally, diving species of bird had larger erythrocytes, indicating that erythrocytes are acting as a vital oxygen store. In summary, erythrocyte size and shape in birds are driven by the need to deliver oxygen during energetically costly activities.

Hydrodynamic coupling for particle-based solvent-free membrane models

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

Hydrodynamics of immiscible binary fluids with viscosity contrast: a multiparticle collision dynamics approach

We present a multiparticle collision dynamics (MPC) implementation of layered immiscible fluids A and B of different shear viscosities separated by planar interfaces. The simulated flow profile for imposed steady shear motion and the time-dependent shear stress functions are in excellent agreement with our continuum hydrodynamics results for the composite fluid. The wave-vector dependent transverse velocity auto-correlation functions (TVAF) in the bulk-fluid regions of the layers decay exponentially, and agree with those of single-phase isotropic MPC fluids. In addition, we determine the hydrodynamic mobilities of an embedded colloidal sphere moving steadily parallel or transverse to a fluid-fluid interface, as functions of the distance from the interface. The obtained mobilities are in good agreement with hydrodynamic force multipoles calculations, for a no-slip sphere moving under creeping flow conditions near a clean, ideally flat interface. The proposed MPC fluid-layer model can be straightforwardly implemented, and it is computationally very efficient. Yet, owing to the spatial discretization inherent to the MPC method, the model can not reproduce all hydrodynamic features of an ideally flat interface between immiscible fluids.

Functionalized supported membranes for quantifying adhesion of P. falciparum-infected erythrocytes

The pathology of Plasmodium falciparum malaria is largely defined by the cytoadhesion of infected erythrocytes to the microvascular endothelial lining. The complexity of the endothelial surface and the large range of interactions available for the infected erythrocyte via parasite-encoded adhesins make analysis of critical contributions during cytoadherence challenging to define. Here, we have explored supported membranes functionalized with two important adhesion receptors, ICAM1 or CD36, as a quantitative biomimetic surface to help understand the processes involved in cytoadherence. Parasitized erythrocytes bound to the receptor-functionalized membranes with high efficiency and selectivity under both static and flow conditions, with infected wild-type erythrocytes displaying a higher binding capacity than do parasitized heterozygous sickle cells. We further show that the binding efficiency decreased with increasing intermolecular receptor distance and that the cell-surface contacts were highly dynamic and increased with rising wall shear stress as the cell underwent a shape transition. Computer simulations using a deformable cell model explained the wall-shear-stress-induced dynamic changes in cell shape and contact area via the specific physical properties of erythrocytes, the density of adhesins presenting knobs, and the lateral movement of receptors in the supported membrane.

Microfluidic Obstacle Arrays Induce Large Reversible Shape Change in Red Blood Cells

Red blood cell (RBC) shape change under static and dynamic shear stress has been a source of interest for at least 50 years. High-speed time-lapse microscopy was used to observe the rate of deformation and relaxation when RBCs are subjected to periodic shear stress and deformation forces as they pass through an obstacle. We show that red blood cells are reversibly deformed and take on characteristic shapes not previously seen in physiological buffers when the maximum shear stress was between 2.2 and 25 Pa (strain rate 2200 to 25,000 s⁻¹). We quantify the rates of RBC deformation and recovery using Kaplan–Meier survival analysis. The time to deformation decreased from 320 to 23 milliseconds with increasing flow rates, but the distance traveled before deformation changed little. Shape recovery, a measure of degree of deformation, takes tens of milliseconds at the lowest flow rates and reached saturation at 2.4 s at a shear stress of 11.2 Pa indicating a maximum degree of deformation was reached. The rates and types of deformation have relevance in red blood cell disorders and in blood cell behavior in microfluidic devices.

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

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

Numerical Investigation of the Dynamical Behavior of a Fluid-Filled Microparticle Suspended in Human Arteriole

The study of artificial microparticles (capsules and vesicles) has gained a growing interest with the emergence of bio-engineering. One of their promoting applications is their use as therapeutic vectors for drug delivery, when capsules and vesicles release their capacity in a targeted environment. The dynamic behavior of capsules and vesicles in confined or unbounded flows was widely studied in the literature and their mechanical response was truthfully described using constitutive laws with good agreement with experiences. However, in a context of biological application, to our knowledge, none of published studies investigating the mechanical response of deformable microparticle took into account the real physiological conditions: the rheological properties of blood such as carrying fluid and the mechanical properties of blood vessels. In this paper, we consider a hyperelastic microparticle suspended in human arteriole. We investigate the deformation of the microparticle resulting from its interaction with blood flow and the arteriolar wall using various capillary numbers and respecting physiological properties of blood and arterial wall. The influence of the blood viscosity model (Newtonian versus shear thinning) is investigated and a comparison with a rigid microchannel and a muscle-embedded arteriole is carried out. The fluid structure interaction (FSI) problem is solved using arbitrary Lagrangian Eulerian (ALE) method. Our simulations have revealed that the arteriolar wall distensibility deeply influences both the deformation and velocity of the microparticle: the deformation strongly increases while the velocity decreases in comparison to an infinitely rigid wall. In the context of therapeutic procedure of targeted drug-delivery, a particular attention should be addressed to these observations, in particular for their implication in the burst mechanism.

Symmetry properties of nonlinear hydrodynamic interactions between responsive particles

Two identical particles driven by the same steady force through a viscous fluid may move relative to one another due to hydrodynamic interactions. The presence or absence of this relative translation has a profound effect on the dynamics of a driven suspension consisting of many particles. We consider a pair of particles which, to linear order in the force, do not interact hydrodynamically. If the system possesses an intrinsic property (such as the shape of the particles, their position with respect to a boundary, or the shape of the boundary) which is affected by the external forcing, hydrodynamic interactions that depend nonlinearly on the force may emerge. We study the general properties of such nonlinear response. Analysis of the symmetries under particle exchange and under force reversal leads to general conclusions concerning the appearance of relative translation and the motion's time reversibility. We demonstrate the applicability of the conclusions in three specific examples: (a) two spheres driven parallel to a wall; (b) two deformable objects driven parallel to their connecting line; and (c) two spheres driven along a curved path. The breaking of time reversibility suggests a possible use of nonlinear hydrodynamic interactions to disperse or assemble particles by an alternating force.

Vortical flow structures induced by red blood cells in capillaries

OBJECTIVE

Knowledge about the flow field of the plasma around the red blood cells in capillary flow is important for a physical understanding of blood flow and the transport of micro- and nanoparticles and molecules in the flowing plasma. We conducted an experimental study on the flow field around red blood cells in capillary flow that is complemented by simulations of vortical flow between red blood cells.

METHODS

Red blood cells were injected in a 10 × 12 µm rectangular microchannel at a low hematocrit, and the flow field around one or two cells was captured by a high-speed camera that tracked 250 nm nanoparticles in the flow field, acting as tracers.

RESULTS

While the flow field around a steady "croissant" shape is found to be similar to that of a rigid sphere, the flow field around a "slipper" shape exhibits a small vortex at the rear of the red blood cell. Even more pronounced are vortex-like structures observed in the central region between two neighboring croissants.

CONCLUSIONS

The rotation frequency of the vortices is to a good approximation, inversely proportional to the distance between the cells. Our experimental data are complemented by numerical simulations.

Interplay between cooperativity of intercellular receptor-ligand binding and coalescence of nanoscale lipid clusters in adhering membranes

Adhesion of biological cells is mediated by the specific binding of receptors and ligands which are typically large proteins spanning through the plasma membranes of the contacting cells. The receptors and ligands can exhibit affinity for nanoscale lipid clusters that form within the plasma membrane. A central question is how these nanoscale lipid clusters physically affect and respond to the receptor-ligand binding during cell adhesion. Within the framework of classical statistical mechanics we find that the receptor-ligand binding reduces the threshold energy for lipid clusters to coalesce into mesoscale domains by up to ∼50%, and that the formation of these domains induces significant cooperativity of the receptor-ligand binding. The interplay between the receptor-ligand binding cooperativity and the lipid domain formation manifests acute sensitivity of the membrane system to changes in control parameters. This sensitivity can be crucial in cell signaling and immune responses.

Buckling transitions and soft-phase invasion of two-component icosahedral shells

What is the optimal distribution of two types of crystalline phases on the surface of icosahedral shells, such as of many viral capsids? We here investigate the distribution of a thin layer of soft material on a crystalline convex icosahedral shell. We demonstrate how the shapes of spherical viruses can be understood from the perspective of elasticity theory of thin two-component shells. We develop a theory of shape transformations of an icosahedral shell upon addition of a softer, but still crystalline, material onto its surface. We show how the soft component "invades" the regions with the highest elastic energy and stress imposed by the 12 topological defects on the surface. We explore the phase diagram as a function of the surface fraction of the soft material, the shell size, and the incommensurability of the elastic moduli of the rigid and soft phases. We find that, as expected, progressive filling of the rigid shell by the soft phase starts from the most deformed regions of the icosahedron. With a progressively increasing soft-phase coverage, the spherical segments of domes are filled first (12 vertices of the shell), then the cylindrical segments connecting the domes (30 edges) are invaded, and, ultimately, the 20 flat faces of the icosahedral shell tend to be occupied by the soft material. We present a detailed theoretical investigation of the first two stages of this invasion process and develop a model of morphological changes of the cone structure that permits noncircular cross sections. In conclusion, we discuss the biological relevance of some structures predicted from our calculations, in particular for the shape of viral capsids.

Multi-particle collision dynamics with a non-ideal equation of state. I

The method of multi-particle collision dynamics (MPCD) and its different implementations are commonly used in the field of soft matter physics to simulate fluid flow at the micron scale. Typically, the coarse-grained fluid particles are described by the equation of state of an ideal gas, and the fluid is rather compressible. This is in contrast to conventional fluids, which are incompressible for velocities much below the speed of sound, and can cause inhomogeneities in density. We propose an algorithm for MPCD with a modified collision rule that results in a non-ideal equation of state and a significantly decreased compressibility. It allows simulations at less computational costs compared to conventional MPCD algorithms. We derive analytic expressions for the equation of state and the corresponding compressibility as well as shear viscosity. They show overall very good agreement with simulations, where we determine the pressure by simulating a quiet bulk fluid and the shear viscosity by simulating a linear shear flow and a Poiseuille flow.

Modeling of Gas Exchange in the Lungs

This overview presents the recent progress in our understanding of gas transfer by the lungs during the respiratory cycle and during breath holding. Different phenomena intervene in gas transfer, convection and diffusion in the gas, dissolution, diffusion across the alveolar-capillary membrane, diffusion across blood plasma, and finally diffusion and reaction with hemoglobin inside blood cells. The different gases, O₂ , CO, and NO, have very different reaction times with hemoglobin ranging from a few microseconds to tens of milliseconds. This is leading to different outcomes. For O₂ , the solutions to the coupled nonlinear gas and blood equations are obtained at the acinus level. They include the fact that the acinar internal ventilation is strongly heterogeneous due to the arborescent structure. Also, in the dynamic calculation, one takes care of the delay between the start of inhalation and arrival of fresh air in the acinus. This "dead" time is the dynamic equivalent of the dead space ventilation. The question of the dependence of Vo₂ on ventilation and perfusion takes a different form. The results show that Vo₂ is not only a function of the ventilation/perfusion ratio but also depends on the variables: acinar ventilation VE_ac and perfusion Q_ac . The ratio VE_ac /Q_ac roughly determines arterial O₂ saturation and arterial and alveolar O₂ partial pressure. The classic Roughton-Forster interpretation of DLCO (separation between independent membrane and blood resistance) was a mathematical conjecture. It was shown recently that this conjecture was violated. This article presents an alternative interpretation that uses time concepts instead of resistance. © 2021 American Physiological Society. Compr Physiol 11:1289-1314, 2021.

Targeting spectrin redox switches to regulate the mechanoproperties of red blood cells

The mechanical properties of red blood cells (RBCs) are fundamental for their physiological role as gas transporters. RBC flexibility and elasticity allow them to survive the hemodynamic changes in the different regions of the vascular tree, to dynamically contribute to the flow thereby decreasing vascular resistance, and to deform during the passage through narrower vessels. RBC mechanoproperties are conferred mainly by the structural characteristics of their cytoskeleton, which consists predominantly of a spectrin scaffold connected to the membrane via nodes of actin, ankyrin and adducin. Changes in redox state and treatment with thiol-targeting molecules decrease the deformability of RBCs and affect the structure and stability of the spectrin cytoskeleton, indicating that the spectrin cytoskeleton may contain redox switches. In this perspective review, we revise current knowledge about the structural and functional characterization of spectrin cysteine redox switches and discuss the current lines of research aiming to understand the role of redox regulation on RBC mechanical properties. These studies may provide novel functional targets to modulate RBC function, blood viscosity and flow, and tissue perfusion in disease conditions.

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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