Computational study of cell adhesion and rolling in flow channel by meshfree method

Tethering and rolling of circulating leukocytes on the surface of endothelium are critical steps during an inflammatory response. A soft solid cell model was proposed to study monocytes tethering and rolling behaviors on substrate surface in shear flow. The interactions between monocytes and micro-channel surface were modeled by a coarse-grained molecular adhesive potential. The computational model was implemented in a Lagrange-type meshfree Galerkin formulation to investigate the monocyte tethering and rolling process with different flow rates. From the simulation results, it was found that the flow rate has profound effects on the rolling velocity, contact area and effective stress of monocytes. As the flow rate increased, the rolling velocity would increase linearly, whereas the contact area and average effective stress in monocyte showed nonlinear increase.     

Influence of sickle hemoglobin polymerization and membrane properties on deformability of sickle erythrocytes in the microcirculation.

The rheological properties of normal erythrocytes appear to be largely determined by those of the red cell membrane. In sickle cell disease, the intracellular polymerization of sickle hemoglobin upon deoxygenation leads to a marked increase in intracellular viscosity and elastic stiffness as well as having indirect effects on the cell membrane. To estimate the components of abnormal cell rheology due to the polymerization process and that due to the membrane abnormalities, we have developed a simple mathematical model of whole cell deformability in narrow vessels. This model uses hydrodynamic lubrication theory to describe the pulsatile flow in the gap between a cell and the vessel wall. The interior of the cell is modeled as a Voigt viscoelastic solid with parameters for the viscous and elastic moduli, while the membrane is assigned an elastic shear modulus. In response to an oscillatory fluid shear stress, the cell--modeled as a cylinder of constant volume and surface area--undergoes a conical deformation which may be calculated. We use published values of normal and sickle cell membrane elastic modulus and of sickle hemoglobin viscous and elastic moduli as a function of oxygen saturation, to estimate normalized tip displacement, d/ho, and relative hydrodynamic resistance, Rr, as a function of polymer fraction of hemoglobin for sickle erythrocytes. These results show the transition from membrane to internal polymer dominance of deformability as oxygen saturation is lowered. More detailed experimental data, including those at other oscillatory frequencies and for cells with higher concentrations of hemoglobin S, are needed to apply fully this approach to understanding the deformability of sickle erythrocytes in the microcirculation. The model should be useful for reconciling the vast and disparate sets of data available on the abnormal properties of sickle cell hemoglobin and sickle erythrocyte membranes, the two main factors that lead to pathology in patients with this disease.     

Investigating the effects of membrane deformability on artificial capsule adhesion to the functionalized surface

Understanding, manipulating and controlling cellular adhesion processes can be critical in developing biomedical technologies. Adhesive mechanisms can be used to the target, pattern and separate cells such as leukocytes from whole blood for biomedical applications. The deformability response of the cell directly affects the rolling and adhesion behavior under viscous linear shear flow conditions. To that end, the primary objective of the present study was to investigate numerically the influence of capsule membrane’s nonlinear material behavior (i.e. elastic-plastic to strain hardening) on the rolling and adhesion behavior of representative artificial capsules. Specifically, spherical capsules with radius of \(3.75\, \upmu \hbox {m}\) were represented using an elastic membrane governed by a Mooney–Rivlin strain energy functions. The surfaces of the capsules were coated with P-selectin glycoprotein-ligand-1 to initiate binding interaction with P-selectin-coated planar surface with density of \(150\,\upmu \hbox {m}^{-2}\) under linear shear flow varying from 100 to \(400\,\hbox {s}^{-1}\). The numerical model is based on the Immersed Boundary Method for rolling of deformable capsule in shear flow coupled with Monte Carlo simulation for receptor/ligand interaction modeled using Bell model. The results reveal that the mechanical properties of the capsule play an important role in the rolling behavior and the binding kinetics between the capsule contact surface and the substrate. The rolling behavior of the strain hardening capsules is relatively smoother and slower compared to the elastic-plastic capsules. The strain hardening capsules exhibits higher contact area at any given shear rate compared to elastic-plastic capsules. The increase in contact area leads to decrease in rolling velocity. The capsule contact surface is not in complete contact with the substrate because of thin lubrication film that is trapped between the capsule and substrate. This creates a concave shape on the bottom surface of the capsule that is referred to as a dimple. In addition, the present study demonstrates that the average total bond force from the capsules lifetime increases by 37 % for the strain hardening capsules compared to elastic-plastic capsules at shear rate of \(400\,\hbox {s}^{-1}\). Finally, the model demonstrates the effect of finite membrane deformation on the coupling between hydrodynamic and receptor/ligand interaction.     

Affinity of red blood cell membrane for particle surfaces measured by the extent of particle encapsulation.

An experimental technique and a simple analysis are presented that can be used to quantitate the affinity of red blood cell membrane for surfaces of small beads or microsomal particles up to 3 micrometers Diam. The technique is demonstrated with an example of dextran-mediated adhesion of small spherical red cell fragments to normal red blood cells. Cells and particles are positioned for contact by manipulation with glass micropipets. The mechanical equilibrium of the adhesive contact is represented by the variational expression that the decrease in interfacial free energy due to a virtual increase in contact area is balanced by the increase in elastic energy of the membrane due to virtual deformation. The surface affinity is the reduction in free energy per unit area of the interface associated with the formation of adhesive contact. From numerical computations of equilibrium configurations, the surface affinity is derived as a function of the fractional extent of particle encapsulation. The range of surface affinities for which the results are applicable is increased over previous techniques to several times the value of the elastic shear modulus. It is shown that bending rigidity of the membrane has little effect on the analytical results for particles 1--3 micrometers Diam and that results are essentially the same for both cup- and disk-shaped red cells. A simple analytical model is shown to give a good approximation for surface affinity (normalized by the elastic shear modulus) as a function of the fractional extent of particle encapsulation. The model predicts that a particle would be almost completely vacuolized for surface affinities greater than or equal to 10 times the elastic shear modulus. Based on an elastic shear modulus of 6.6 x 10(-3) dyn/cm, the range for the red cell-particle surface affinity as measured by this technique is from approximately 7 x 10(-4) to 7 x 10(-2) erg/cm2. Also, an approximate relation is derived for the level of surface affinity necessary to produce particle vacuolization by a phospholipid bilayer surface which possesses bending rigidity and a fixed tension.     

Influence of cell deformation on leukocyte rolling adhesion in shear flow

Blood cell interaction with vascular endothelium is important in microcirculation, where rolling adhesion of circulating leukocytes along the surface of endothelial cells is a prerequisite for leukocyte emigration under flow conditions. HL-60 cell rolling adhesion to surface-immobilized P-selectin in shear flow was investigated using a side-view flow chamber, which permitted measurements of cell deformation and cell-substrate contact length as well as cell rolling velocity. A two-dimensional model was developed based on the assumption that fluid energy input to a rolling cell was essentially distributed into two parts: cytoplasmic viscous dissipation, and energy needed to break adhesion bonds between the rolling cell and its substrate. The flow fields of extracellular fluid and intracellular cytoplasm were solved using finite element methods with a deformable cell membrane represented by an elastic ring. The adhesion energy loss was calculated based on receptor-ligand kinetics equations. It was found that, as a result of shear-flow-induced cell deformation, cell-substrate contact area under high wall shear stresses (20 dyn/cm2) could be as much as twice of that under low stresses (0.5 dyn/cm2). An increase in contact area may cause more energy dissipation to both adhesion bonds and viscous cytoplasm, whereas the fluid energy input may decrease due to the flattened cell shape. Our model predicts that leukocyte rolling velocity will reach a plateau as shear stress increases, which agrees with both in vivo and in vitro experimental observations.     

A 3-D computational model predicts that cell deformation affects selectin-mediated leukocyte rolling

Leukocyte recruitment to sites of inflammation is initiated by their tethering and rolling on the activated endothelium under flow. Even though the fast kinetics and high tensile strength of selectin-ligand bonds are primarily responsible for leukocyte rolling, experimental evidence suggests that cellular properties such as cell deformability and microvillus elasticity actively modulate leukocyte rolling behavior. Previous theoretical models either assumed cells as rigid spheres or were limited to two-dimensional representations of deformable cells with deterministic receptor-ligand kinetics, thereby failing to accurately predict leukocyte rolling. We therefore developed a three-dimensional computational model based on the immersed boundary method to predict receptor-mediated rolling of deformable cells in shear flow coupled to a Monte Carlo method simulating the stochastic receptor-ligand interactions. Our model predicts for the first time that the rolling of more compliant cells is relatively smoother and slower compared to cells with stiffer membranes, due to increased cell-substrate contact area. At the molecular level, we show that the average number of bonds per cell as well as per single microvillus decreases with increasing membrane stiffness. Moreover, the average bond lifetime decreases with increasing shear rate and with increasing membrane stiffness, due to higher hydrodynamic force experienced by the cell. Taken together, our model captures the effect of cellular properties on the coupling between hydrodynamic and receptor-ligand bond forces, and successfully explains the stable leukocyte rolling at a wide range of shear rates over that of rigid microspheres.     

A numerical analysis of forces exerted by laminar flow on spreading cells in a parallel plate flow chamber assay

Exposure of spreading anchorage-dependent cells to laminar flow is a common technique to measure the strength of cell adhesion. Since cells protrude into the flow stream, the force exerted by the fluid on the cells is a function of cell shape. To assess the relationship between cell shape and the hydrodynamic force on adherent cells, we obtained numerical solutions of the velocity and stress fields around bovine aortic endothelial cells during various stages of spreading and calculated the force required to detach the cells. Morphometric parameters were obtained from light and scanning electron microscopy measurements. Cells were assumed to have a constant volume, but the surface area increased during spreading until the membrane was stretched taut. Two-dimensional models of steady flow were generated using the software packages ANSYS (mesh generation) and FIDAP (problem solution). The validity of the numerical results was tested by comparison with published results for a semicircle in contact with the surface. The drag force and torque were greatest for round cells making initial contact with the surface. During spreading, the drag force and torque declined by factors of 2 and 20, respectively. The calculated forces and moments were used in adhesion models to predict the wall shear stress at which the cells detached. Based upon published values for the bond force and receptor number, round cells should detach at shear stresses between 2.5 and 6 dyn/cm(2), whereas substantially higher stresses are needed to detach spreading and fully spread cells. Results from the simulations indicate that (1) the drag force varies little with cell shape whereas the torque is very sensitive to cell shape, and (2) the increase in the strength of adhesion during spreading is due to increased contact area and receptor densities within the contact area. (c) 1993 John Wiley & Sons, Inc.     

Rolling adhesion of leukocytes on soft substrates: Does substrate stiffness matter?

Cell rolling on vascular endothelium under hydrodynamic blood flow is critical for realization of many physiological and pathological processes, such as inflammatory response and tumor metastasis. The blood-borne cells are in direct contact with the inner layer of endothelium, formed by a highly compliant layer of endothelial cells. The effect of endothelial stiffness on the adhesion and motion of rolling cells is poorly understood. Inspired by recent in vitro studies, here we implemented a computational method to model the specific adhesion of a rolling cell onto a soft substrate, subjected to a creeping shear flow. The substrate is modeled as an elastic half-space, coated with P- and E-selectin receptors with specific affinity for the complementary ligands located on the moving cell. Of particular importance is to predict the effect of substrate stiffness on cell adhesion and its kinematics and kinetics of motion. Simulation results show that the effect of substrate compliance is minimal when coated with P-selectin. Conversely, the trajectory of rolling cells on E-selectin coated substrates is sensitive to the substrate compliance. This is attributed to the moderation of binding forces applied by the soft substrate which leads to a higher average translational velocity of cells.     

Critical stresses for cancer cell detachment in microchannels

We present experiments involving cancer cells adhering to microchannels, subjected to increasing shear stresses (0.1–30 Pa). Morphological studies were carried out at different shear stresses. Cells exhibit spreading patterns similar to those observed under static conditions, as long as the shear stress is not too high. At critical wall shear stresses (around 2−5 Pa), cell-substrate contact area decreases until detachment at the larger stresses. Critical shear stresses are found to be lower for higher confinements (i.e. smaller cell height to channel height ratio). Fluorescent techniques were used to locate focal adhesions (typically 1 μm² in size) under various shearing conditions, showing that cells increase the number of focal contacts in the region facing the flow. To analyze such data, we propose a model to determine the critical stress, resulting from the competition between hydrodynamic forces and the adhesive cell resistance. With this model, typical adhesive stresses exerted at each focal contact can be determined and are in agreement with previous works.     

Design of a side-view particle imaging velocimetry flow system for cell-substrate adhesion studies

Experimental models that mimic the flow conditions in microcapillaries have suggested that the local shear stresses and shear rates can mediate tumor cell and leukocyte arrest on the endothelium and subsequent sustained adhesion. However, further investigation has been limited by the lack of experimental models that allow quantitative measurement of the hydrodynamic environment over adherent cells. The purpose of this study was to develop a system capable of acquiring quantitative flow profiles over adherent cells. By combining the techniques of side-view imaging and particle image velocimetry (PIV), an in vitro model was constructed that is capable of obtaining quantitative flow data over cells adhering to the endothelium. The velocity over an adherent leukocyte was measured and the shear rate was calculated under low and high upstream wall shear. The microcapillary channel was modeled using computational fluid dynamics (CFD) and the calculated velocity profiles over cells under the low and high shear rates were compared to experimental results. The drag force applied to each cell by the fluid was then computed. This system provides a means for future study of the forces underlying adhesion by permitting characterization of the local hydrodynamic conditions over adherent cells.     

A Multiscale Red Blood Cell Model with Accurate Mechanics, Rheology, and Dynamics

Red blood cells (RBCs) have highly deformable viscoelastic membranes exhibiting complex rheological response and rich hydrodynamic behavior governed by special elastic and bending properties and by the external/internal fluid and membrane viscosities. We present a multiscale RBC model that is able to predict RBC mechanics, rheology, and dynamics in agreement with experiments. Based on an analytic theory, the modeled membrane properties can be uniquely related to the experimentally established RBC macroscopic properties without any adjustment of parameters. The RBC linear and nonlinear elastic deformations match those obtained in optical-tweezers experiments. The rheological properties of the membrane are compared with those obtained in optical magnetic twisting cytometry, membrane thermal fluctuations, and creep followed by cell recovery. The dynamics of RBCs in shear and Poiseuille flows is tested against experiments and theoretical predictions, and the applicability of the latter is discussed. Our findings clearly indicate that a purely elastic model for the membrane cannot accurately represent the RBC's rheological properties and its dynamics, and therefore accurate modeling of a viscoelastic membrane is necessary.     

Thermoelasticity of red blood cell membrane.

The elastic properties of the human red blood cell membrane have been measured as functions of temperature. The area compressibility modulus and the elastic shear modulus, which together characterize the surface elastic behavior of the membrane, have been measured over the temperature range of 2-50 degrees C with micropipette aspiration of flaccid and osmotically swollen red cells. In addition, the fractional increase in membrane surface area from 2-50 degrees C has been measured to give a value for the thermal area expansivity. The value of the elastic shear modulus at 25 degrees C was measured to be 6.6 X 10(-3) dyne/cm. The change in the elastic shear modulus with temperature was -6 X 10(-5) dyne/cm degrees C. Fractional forces were shown to be only on the order of 10-15%. The area compressibility modulus at 25 degrees C was measured to be 450 dyne/cm. The change in the area compressibility modulus with temperature was -6 dyne/cm degrees C. The thermal area expansivity for red cell membrane was measured to be 1.2 X 10(-3)/degrees C. With this data and thermoelastic relations the heat of expansion is determined to be 110-200 ergs/cm2; the heat of extension is 2 X 10(-2) ergs/cm2 for unit extension of the red cell membrane. The heat of expansion is of the order anticipated for a lipid bilayer idealized as twice the behavior of a monolayer at an oil-water interface. The observation that the heat of extension is positive demonstrates that the entropy of the material increases with extension, and that the dominant mechanism of elastic energy storage is energetic. Assuming that the red cell membrane shear rigidity is associated with "spectrin," unit extension of the membrane increases the configurational entropy of spectrin by 500 cal/mol.     

Study of local hydrodynamic environment in cell-substrate adhesion using side-view μPIV technology

Tumor cell adhesion to the endothelium under shear flow conditions is a critical step that results in circulation-mediated tumor metastasis. This study presents experimental and computational techniques for studying the local hydrodynamic environment around adherent cells and how local shear conditions affect cell-cell interactions on the endothelium in tumor cell adhesion. To study the local hydrodynamic profile around heterotypic adherent cells, a side-view flow chamber assay coupled with micro particle imaging velocimetry (μPIV) technique was developed, where interactions between leukocytes and tumor cells in the near-endothelial wall region and the local shear flow environment were characterized. Computational fluid dynamics (CFD) simulations were also used to obtain quantitative flow properties around those adherent cells. Results showed that cell dimension and relative cell-cell positions had strong influence on local shear rates. The velocity profile above leukocytes and tumor cells displayed very different patterns. Larger cell deformations led to less disturbance to the flow. Local shear rates above smaller cells were observed to be more affected by relative positions between two cells.     

Effects of flowing RBCs on adhesion of a circulating tumor cell in microvessels

Adhesion of circulating tumor cells (CTCs) to the microvessel wall largely depends on the blood hydrodynamic conditions, one of which is the blood viscosity. Since blood is a non-Newtonian fluid, whose viscosity increases with hematocrit, in the microvessels at low shear rate. In this study, the effects of hematocrit, vessel size, flow rate and red blood cell (RBC) aggregation on adhesion of a CTC in the microvessels were numerically investigated using dissipative particle dynamics. The membrane of cells was represented by a spring-based network connected by elastic springs to characterize its deformation. RBC aggregation was modeled by a Morse potential function based on depletion-mediated assumption, and the adhesion of the CTC to the vessel wall was achieved by the interactions between receptors and ligands at the CTC and those at the endothelial cells forming the vessel wall. The results demonstrated that in the microvessel of \(15\,\upmu \hbox {m}\) diameter, the CTC has an increasing probability of adhesion with the hematocrit due to a growing wall-directed force, resulting in a larger number of receptor–ligand bonds formed on the cell surface. However, with the increase in microvessel size, an enhanced lift force at higher hematocrit detaches the initial adherent CTC quickly. If the microvessel is comparable to the CTC in diameter, CTC adhesion is independent of Hct. In addition, the velocity of CTC is larger than the average blood flow velocity in smaller microvessels and the relative velocity of CTC decreases with the increase in microvessel size. An increased blood flow resistance in the presence of CTC was also found. Moreover, it was found that the large deformation induced by high flow rate and the presence of aggregation promote the adhesion of CTC.     

Dynamics of the primary cilium in shear flow

In this work, the equilibrium shape and dynamics of a primary cilium under flow are investigated by using both theoretical modeling and experiment. The cilium is modeled as an elastic beam that may undergo large deflection due to the hydrodynamic load. Equilibrium results show that the anchoring effects of the basal body on the cilium axoneme behave as a nonlinear rotational spring. Details of the rotational spring are elucidated by coupling the elastic beam with an elastic shell. We further study the dynamics of cilium under shear flow with the cilium base angle determined from the nonlinear rotational spring, and obtain good agreement in cilium bending and relaxing dynamics when comparing between modeling and experimental results. These results potentially shed light on the physics underlying the mechanosensitive ion channel transport through the ciliary membrane.     

Impact of aeration strategy on CHO cell performance during antibody production

Stirred tank bioreactors using suspension adapted mammalian cells are typically used for the production of complex therapeutic proteins. The hydrodynamic conditions experienced by cells within this environment have been shown to directly impact growth, productivity, and product quality and therefore an improved understanding of the cellular response is critical. Here we investigate the sub‐lethal effects of different aeration strategies on Chinese hamster ovary cells during monoclonal antibody production. Two gas delivery systems were employed to study the presence and absence of the air–liquid interface: bubbled direct gas sparging and a non‐bubbled diffusive silicone membrane system. Additionally, the effect of higher gas flow rate in the sparged bioreactor was examined. Both aeration systems were run using chemically defined media with and without the shear protectant Pluronic F‐68 (PF‐68). Cells were unable to grow with direct gas sparging without PF‐68; however, when a silicone membrane aeration system was implemented growth was comparable to the sparged bioreactor with PF‐68, indicating the necessity of shear protectants in the presence of bubbles. The cultures exposed to increased hydrodynamic stress were shown by flow cytometry to have decreased F‐actin intensity within the cytoskeleton and enter apoptosis earlier. This indicates that these conditions elicit a sub‐lethal physiological change in cells that would not be detected by the at‐line assays which are normally implemented during cell culture. These physiological changes only result in a difference in continuous centrifugation performance under high flow rate conditions. Product quality was more strongly affected by culture age than the hydrodynamic conditions tested. © 2012 American Institute of Chemical Engineers Biotechnol. Prog., 2013.     

Crossflow microfiltration of animal cells

Laminar shear is the primary mechanism of cell damage, limiting flow rate (and hence flux) in crossflow microfiltration of animal cells. Sensitivity to hydrodynamic and interfacial stress is reduced by the addition of 0.1% Pluronic polyol. A critical average wall shear rate of 3000 s(-1) (above which damage occurs) is found for several cell types, including mammalian and insect cells. Hydrodynamic stress also limits the maximum tip speed in a rotary lobe pump to less than 350 cm/s. Turbulent flow in the recirculation loop piping at Reynolds numbers of up to 71,000 does not cause cell damage. Maximum sustainable flux decreases with cell concentration and increases with cell size (in qualitative agreement with the hydrodynamic lift model). A flux of 30 to 75 L/m(2) h (depending on cell size) can be sustained during 20-fold concentration from 2.5 x 10(6) cells/ml, while maintaining high cell viability.     

Aggregation and disaggregation of red blood cells

The aggregation of red blood cells may be analyzed as an interaction of an adhesive surface energy and the elastic stored energy which results from deformation of the cell. The adhesive surface energy is the work required to separate a unit adhered area and is the resultant of adhesive forces due to the bridging molecules that bind the cells together and the electrostatic repulsion due to surface charge. The elastic strain energy in the case of the red blood is associated with the membrane elasticity only since the interior of the cell is liquid. The membrane elasticity is due both to bending stiffness and shear. The area expansion is small and may be neglected. These assumptions allow realistic computation of red cell shapes in rouleaux. The disaggregation of rouleaux requires an external force which must overcome the adhesive energy and also supply additional elastic energy of deformation. Depending on the geometry, the initial effect of elastic energy may tend to aid disaggregation. In a shear flow, the stresses on a suspended rouleau alternately tend to compress and to disaggregate the cells if they are free to rotate. This introduces a time dependence so that viscous effects due to the viscosity of the cell membrane, the cell cytoplasm and the external fluid may play a role in determining whether disaggregation proceeds to completion or not.