Heterogeneous ATP patterns in microvascular networks.

ATP is not only an energy carrier but also serves as an important signalling molecule in many physiological processes. Abnormal ATP level in blood vessel is known to be related to several pathologies, such as inflammation, hypoxia and atherosclerosis. Using advanced numerical methods, we analysed ATP released by red blood cells (RBCs) and its degradation by endothelial cells (ECs) in a cat mesentery-inspired vascular network, accounting for RBC mutual interaction and interactions with vascular walls. Our analysis revealed a heterogeneous ATP distribution in the network, with higher concentrations in the cell-free layer, concentration peaks around bifurcations and heterogeneity among vessels of the same level. These patterns arise from the spatio-temporal organization of RBCs induced by the network geometry. It is further shown that an alteration of hematocrit and flow strength significantly affects ATP level as well as heterogeneity in the network. These findings constitute a first building block to elucidate the intricate nature of ATP patterns in vascular networks and the far reaching consequences for other biochemical signalling, such as calcium, by ECs.

Mathematical modeling of intracellular calcium in presence of receptor: a homeostatic model for endothelial cell.

Calcium is a ubiquitous molecule and second messenger that regulates many cellular functions ranging from exocytosis to cell proliferation at different time scales. In the vasculature, a constant adenosine triphosphate (ATP) concentration is maintained because of ATP released by red blood cells (RBCs). These ATP molecules continuously react with purinergic receptors on the surface of endothelial cells (ECs). Consequently, a cascade of chemical reactions are triggered that result in a transient cytoplasmic calcium (Ca[Formula: see text]), followed by return to its basal concentration. The mathematical models proposed in the literature are able to reproduce the transient peak. However, the trailing concentration is always higher than the basal cytoplasmic Ca[Formula: see text] concentrations, and the Ca[Formula: see text] concentration in endoplasmic reticulum (ER) remains lower than its initial concentration. This means that the intracellular homeostasis is not recovered. We propose, herein, a minimal model of calcium kinetics. We find that the desensitization of EC surface receptors due to phosphorylation and recycling plays a vital role in maintaining calcium homeostasis in the presence of a constant stimulus (ATP). The model is able to capture several experimental observations such as refilling of Ca[Formula: see text] in the ER, variation of cytoplasmic Ca[Formula: see text] transient peak in ECs, the resting cytoplasmic Ca[Formula: see text] concentration, the effect of removing ATP from the plasma on Ca[Formula: see text] homeostasis, and the saturation of cytoplasmic Ca[Formula: see text] transient peak with increase in ATP concentration. Direct confrontation with several experimental results is conducted. This work paves the way for systematic studies on coupling between blood flow and chemical signaling, and should contribute to a better understanding of the relation between (patho)physiological conditions and Ca[Formula: see text] kinetics.

Flow driven vesicle unbinding under mechanosensitive adhesion.

Ligand receptor based adhesion is the primary mode of interaction of cellular blood constituents with the endothelium. These adhered entities also experience shear flow imposed by the blood which may lead to their detachment due to the viscous lift forces. Here, we have studied the role of the ligand-receptor bond kinetics in the detachment of an adhered vesicle (a simplified cell model) under shear flow. Using boundary integral formulation we performed numerical simulation of a two dimensional vesicle under shear flow for different values of applied shear rates and time scale of bond kinetics. We observe that the vesicle demonstrates three steady state configurations - adhered, pinned and detached for fast enough ligand-receptor kinetics (akin to Lennard-Jones adhesion). However, for slow bond kinetics the pinned state is not observed. We present scaling laws for the critical shear rates corresponding to the transitions among these three states. These results can help with identifying the processes of cell adhesion/detachment in the blood stream, prevalent features during the immune response and cancer metastasis.

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

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

Amoeboid Swimming Is Propelled by Molecular Paddling in Lymphocytes.

Mammalian cells developed two main migration modes. The slow mesenchymatous mode, like crawling of fibroblasts, relies on maturation of adhesion complexes and actin fiber traction, whereas the fast amoeboid mode, observed exclusively for leukocytes and cancer cells, is characterized by weak adhesion, highly dynamic cell shapes, and ubiquitous motility on two-dimensional and in three-dimensional solid matrix. In both cases, interactions with the substrate by adhesion or friction are widely accepted as a prerequisite for mammalian cell motility, which precludes swimming. We show here experimental and computational evidence that leukocytes do swim, and that efficient propulsion is not fueled by waves of cell deformation but by a rearward and inhomogeneous treadmilling of the cell external membrane. Our model consists of a molecular paddling by transmembrane proteins linked to and advected by the actin cortex, whereas freely diffusing transmembrane proteins hinder swimming. Furthermore, continuous paddling is enabled by a combination of external treadmilling and selective recycling by internal vesicular transport of cortex-bound transmembrane proteins. This mechanism explains observations that swimming is five times slower than the retrograde flow of cortex and also that lymphocytes are motile in nonadherent confined environments. Resultantly, the ubiquitous ability of mammalian amoeboid cells to migrate in two dimensions or three dimensions and with or without adhesion can be explained for lymphocytes by a single machinery of heterogeneous membrane treadmilling.

Amoeboid swimming in a compliant channel.

Several prokaryotes and eukaryotic cells swim in the presence of deformable and rigid surfaces that form confinement. The most commonly observed examples from biological systems are motility of leukocytes and pathogens present within the blood suspension through a microvascular network, and locomotion of eukaryotic cells such as immune system cells and cancerous cells through interstices between soft interstitial cells and the extracellular matrix within the interstitial tissue. This motivated us to investigate numerically the flow dynamics of amoeboid swimming in a flexible channel. The effects of wall stiffness and channel confinement on the flow dynamics and swimmer motion are studied. The swimmer motion through the flexible channel is substantially decelerated compared to the rigid channel. The strong confinement in the amply flexible channel imprisons the swimmer by severely restricting its forward motion. The swimmer velocity in a stiff channel displays nonmonotonic variation with the confinement while it shows monotonic reduction in a highly flexible channel. The physical rationale behind such distinct velocity behaviour in flexible and rigid channels is illustrated using an instantaneous flow field and flow history displayed by the swimmer. This behavior follows from a subtle interplay between the shape changes exhibited by the swimmer and the wall compliance. This study may aid in understanding the influence of elasticity of the surrounding environment on cell motility in immunological surveillance and invasiveness of cancer cells.

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

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