Application of multiphase computational fluid dynamics to analyze monocyte adhesion

The numerical study on the effects of cardiac function on the aortic oxygen distribution

Although venoarterial extracorporeal membrane oxygenation (VA-ECMO) was widely used in clinical practice, the effects of cardiac output on the aortic oxygen distribution were still unclear. Hence, the present study aims to evaluate the effect of cardiac function on the aortic oxygen distribution under VA-ECMO support. A novel model, combining computational fluid dynamics, multiphase fluid approach, and oxygen transport theory together, was established. According to the clinical practice, four cardiac output conditions, including 0, 1, 2, and 2.5 L/min, were designed. The results demonstrated that the proposed method could accurately calculate the distribution of oxygen in the aorta. Moreover, the aortic oxygen distribution was significantly regulated by the local blood flow pattern. The deoxygenated blood flow and oxygenated blood flow met at the aortic arch and formed the so-called oxygenshed phenomenon. Along with the cardiac output increase, the oxygenshed was moved from the proximal of the aortic arch to the descending aorta. Meanwhile, the oxygen contents in the brachiocephalic artery and left common carotid artery were reduced along with the increase of cardiac output. The study could provide much useful information on the oxygen distribution in the aorta to surgeons and operators of VA-ECMO. Graphical abstract The results showed the deoxygenated blood and oxygenated blood met at the aortic arch and formed the so-called "oxygenshed" phenomenon. This phenomenon is consistent with the phenomenon called "watershed".

Hemodynamic modeling of leukocyte and erythrocyte transport and interactions in intracranial aneurysms by a multiphase approach

Hemodynamics has been recognized as an important factor in the development, growth, and rupture of cerebral aneurysms, and investigated by computational fluid dynamics techniques using a single phase approach. However, flow-dependent cell transport and interactions are usually ignored in single phase models, in which blood is usually treated as a single phase Newtonian fluid. For getting better insight into the underlying pathology of intracranial aneurysm, cell transport and interactions should be covered in hemodynamic studies. In the present study, a multiphase hemodynamic model incorporating cell transport and interactions was developed, in which blood was modeled as multiphase fluid having a continuous phase (plasma) and two particulate phases (erythrocytes and leukocytes). The model showed good agreement with experimental data and observations in the literature, and was applied to four patient-specific aneurysms in a pulsatile manner. Leukocyte accumulations were predicted at locations with flow disturbance and low wall shear stress. The concentrations of leukocyte at accumulation sites were found to exceed 200 to 500% of normal physiological level on three unstable aneurysms, including two ruptured aneurysms and a growing aneurysm where accumulation was observed near a daughter sac and a secondary aneurysm. This suggested that aneurysms with complex secondary flow patterns could be prone to leukocyte accumulation on the wall. As this is the first study to characterize cell transport and interactions in aneurysm hemodynamics, our model can serve as a foundation for future intracranial aneurysm models.

Circulatory shear flow alters the viability and proliferation of circulating colon cancer cells

During cancer metastasis, circulating tumor cells constantly experience hemodynamic shear stress in the circulation. Cellular responses to shear stress including cell viability and proliferation thus play critical roles in cancer metastasis. Here, we developed a microfluidic approach to establish a circulatory microenvironment and studied circulating human colon cancer HCT116 cells in response to a variety of magnitude of shear stress and circulating time. Our results showed that cell viability decreased with the increase of circulating time, but increased with the magnitude of wall shear stress. Proliferation of cells survived from circulation could be maintained when physiologically relevant wall shear stresses were applied. High wall shear stress (60.5 dyne/cm(2)), however, led to decreased cell proliferation at long circulating time (1 h). We further showed that the expression levels of β-catenin and c-myc, proliferation regulators, were significantly enhanced by increasing wall shear stress. The presented study provides a new insight to the roles of circulatory shear stress in cellular responses of circulating tumor cells in a physiologically relevant model, and thus will be of interest for the study of cancer cell mechanosensing and cancer metastasis.

Computational study of particle size effects on selective binding of nanoparticles in arterial stenosis

In order to elucidate particle size and wall shear effects on the selective binding of nanoparticles to vessel wall, particle binding to the wall of arterial stenosis was computationally analyzed using a transport and reaction model. The attachment rate constant was modeled as a function of shear rate and particle size. The results showed that it had a positive correlation with the shear rate for particles smaller than 600 nm and a negative correlation with the shear rate for particles larger than 800 nm. Small size particles showed high binding selectivity in the stenosis region for the normal and shear-activated wall, whereas large particles showed high binding selectivity in the low and oscillatory zone for the shear-activated wall.

Local oxidative and nitrosative stress increases in the microcirculation during leukocytes-endothelial cell interactions

Leukocyte-endothelial cell interactions and leukocyte activation are important factors for vascular diseases including nephropathy, retinopathy and angiopathy. In addition, endothelial cell dysfunction is reported in vascular disease condition. Endothelial dysfunction is characterized by increased superoxide (O₂^•−) production from endothelium and reduction in NO bioavailability. Experimental studies have suggested a possible role for leukocyte-endothelial cell interaction in the vessel NO and peroxynitrite levels and their role in vascular disorders in the arterial side of microcirculation. However, anti-adhesion therapies for preventing leukocyte-endothelial cell interaction related vascular disorders showed limited success. The endothelial dysfunction related changes in vessel NO and peroxynitrite levels, leukocyte-endothelial cell interaction and leukocyte activation are not completely understood in vascular disorders. The objective of this study was to investigate the role of endothelial dysfunction extent, leukocyte-endothelial interaction, leukocyte activation and superoxide dismutase therapy on the transport and interactions of NO, O₂^•− and peroxynitrite in the microcirculation. We developed a biotransport model of NO, O₂^•− and peroxynitrite in the arteriolar microcirculation and incorporated leukocytes-endothelial cell interactions. The concentration profiles of NO, O₂^•− and peroxynitrite within blood vessel and leukocytes are presented at multiple levels of endothelial oxidative stress with leukocyte activation and increased superoxide dismutase accounted for in certain cases. The results showed that the maximum concentrations of NO decreased ∼0.6 fold, O₂^•− increased ∼27 fold and peroxynitrite increased ∼30 fold in the endothelial and smooth muscle region in severe oxidative stress condition as compared to that of normal physiologic conditions. The results show that the onset of endothelial oxidative stress can cause an increase in O₂^•− and peroxynitrite concentration in the lumen. The increased O₂^•− and peroxynitrite can cause leukocytes priming through peroxynitrite and leukocytes activation through secondary stimuli of O₂^•− in bloodstream without endothelial interaction. This finding supports that leukocyte rolling/adhesion and activation are independent events.

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.

Computational analysis of nanoparticle adhesion to endothelium: effects of kinetic rate constants and wall shear rates

Various nanoparticles have been developed as imaging probes and drug carriers, and their selectivity in binding to target cells determines the efficacy of these functionalized nanoparticles. Since target cells in different arterial segments experience different hemodynamic environments, we study the effects of wall shear rate waveforms on particle binding. We also explore the effects of the kinetic rate constant, which is determined by particle design parameters, on particle binding. A transport and reaction model is used to evaluate nanoparticle binding to the substrate in a laminar flow chamber. Flow and particle concentration fields are solved by using a computational fluid dynamics. The particle binding rate increases as the mean value of wall shear increases, and the amplitudes of sinusoidal shear waveform do not affect the bound particle density profiles significantly. Particle binding rates increase with the rate constant of attachment (k(A)), and are more sensitively affected by low k(A) values and less by k(A) values higher than 1 × 10⁻⁶ m s⁻¹. Since binding selectivity is affected by k(A) and the wall shear rate, the results of this study can be used for designing functionalized nanoparticles targeting for the specific cells that experience a specific shear rate.

Application of Population Dynamics to Study Heterotypic Cell Aggregations in the Near-Wall Region of a Shear Flow

Our research focused on the polymorphonuclear neutrophils (PMNs) tethering to the vascular endothelial cells (EC) and the subsequent melanoma cell emboli formation in a shear flow, an important process of tumor cell extravasation from the circulation during metastasis. We applied population balance model based on Smoluchowski coagulation equation to study the heterotypic aggregation between PMNs and melanoma cells in the near-wall region of an in vitro parallel-plate flow chamber, which simulates in vivo cell-substrate adhesion from the vasculatures by combining mathematical modeling and numerical simulations with experimental observations. To the best of our knowledge, a multiscale near-wall aggregation model was developed, for the first time, which incorporated the effects of both cell deformation and general ratios of heterotypic cells on the cell aggregation process. Quantitative agreement was found between numerical predictions and in vitro experiments. The effects of factors, including: intrinsic binding molecule properties, near-wall heterotypic cell concentrations, and cell deformations on the coagulation process, are discussed. Several parameter identification approaches are proposed and validated which, in turn, demonstrate the importance of the reaction coefficient and the critical bond number on the aggregation process.