Dynamic modeling for shear stress induced ATP release from vascular endothelial cells

Fabricating a multi-component microfluidic system for exercise-induced endothelial cell mechanobiology guided by hemodynamic similarity

Generating precise in vivo arterial endothelial hemodynamic microenvironments using microfluidics is essential for exploring endothelial mechanobiology. However, a hemodynamic principle guiding the fabrication of microfluidic systems is still lacking. We propose a hemodynamic similarity principle for quickly obtaining the input impedance of the microfluidic system in vitro derived from that of the arterial system in vivo to precisely generate the desired endothelial hemodynamic microenvironments. First, based on the equivalent of blood pressure (BP) and wall shear stress (WSS) waveforms, we establish a hemodynamic similarity principle to efficiently map the input impedance in vivo to that in vitro, after which the multi-component microfluidic system is designed and fabricated using a lumped parameter hemodynamic model. Second, numerical simulation and experimental studies are carried out to validate the performance of the designed microfluidic system. Finally, the intracellular Ca²⁺ responses after exposure to different intensities of exercise-induced BP and WSS waveforms are measured to improve the reliability of EC mechanobiological studies using the designed microfluidic system. Overall, the proposed hemodynamic similarity principle can guide the fabrication of a multi-component microfluidic system for endothelial cell mechanobiology.

A microfluidic system for precisely reproducing physiological blood pressure and wall shear stress to endothelial cells

To reproduce hemodynamic stress microenvironments of endothelial cells in vitro is of vital significance, by which one could exploit the quantitative impact of hemodynamic stresses on endothelial function and seek innovative approaches to prevent circulatory system diseases. Although microfluidic technology has been regarded as an effective method to create physiological microenvironments, a microfluidic system to precisely reproduce physiological arterial hemodynamic stress microenvironments has not been reported yet. In this paper, a novel microfluidic chip consisting of a cell culture chamber with on-chip afterload components designed by the principle of input impedance to mimic the global hemodynamic behaviors is proposed. An external feedback control system is developed to accurately generate the input pressure waveform. A lumped parameter hemodynamic model (LPHM) is built to represent the input impedance to mimic the on-chip global hemodynamic behaviors. Sensitivity analysis of the model parameters is also elaborated. The performance of reproducing physiological blood pressure and wall shear stress is validated by both numerical characterization and flow experiment. Investigation of intracellular calcium ion dynamics in human umbilical vein endothelial cells is finally conducted to demonstrate the biological applicability of the proposed microfluidic system.

The Story of Wall Shear Stress in Coronary Artery Atherosclerosis: Biochemical Transport and Mechanotransduction

Coronary artery atherosclerosis is a local, multifactorial, complex disease, and the leading cause of death in the US. Complex interactions between biochemical transport and biomechanical forces influence disease growth. Wall shear stress (WSS) affects coronary artery atherosclerosis by inducing endothelial cell mechanotransduction and by controlling the near-wall transport processes involved in atherosclerosis. Each of these processes is controlled by WSS differently and therefore has complicated the interpretation of WSS in atherosclerosis. In this paper, we present a comprehensive theory for WSS in atherosclerosis. First, a short review of shear stress-mediated mechanotransduction in atherosclerosis was presented. Next, subject-specific computational fluid dynamics (CFD) simulations were performed in ten coronary artery models of diseased and healthy subjects. Biochemical-specific mass transport models were developed to study low-density lipoprotein, nitric oxide, adenosine triphosphate, oxygen, monocyte chemoattractant protein-1, and monocyte transport. The transport results were compared with WSS vectors and WSS Lagrangian coherent structures (WSS LCS). High WSS magnitude protected against atherosclerosis by increasing the production or flux of atheroprotective biochemicals and decreasing the near-wall localization of atherogenic biochemicals. Low WSS magnitude promoted atherosclerosis by increasing atherogenic biochemical localization. Finally, the attracting WSS LCS's role was more complex where it promoted or prevented atherosclerosis based on different biochemicals. We present a summary of the different pathways by which WSS influences coronary artery atherosclerosis and compare different mechanotransduction and biotransport mechanisms.

Modeling of Endothelial Calcium Responses within a Microfluidic Generator of Spatio-Temporal ATP and Shear Stress Signals

Intracellular calcium dynamics play essential roles in the proper functioning of cellular activities. It is a well known important chemosensing and mechanosensing process regulated by the spatio-temporal microenvironment. Nevertheless, how spatio-temporal biochemical and biomechanical stimuli affect calcium dynamics is not fully understood and the underlying regulation mechanism remains missing. Herein, based on a developed microfluidic generator of biochemical and biomechanical signals, we theoretically analyzed the generation of spatio-temporal ATP and shear stress signals within the microfluidic platform and investigated the effect of spatial combination of ATP and shear stress stimuli on the intracellular calcium dynamics. The simulation results demonstrate the capacity and flexibility of the microfluidic system in generating spatio-temporal ATP and shear stress. Along the transverse direction of the microchannel, dynamic ATP signals of distinct amplitudes coupled with identical shear stress are created, which induce the spatio-temporal diversity in calcium responses. Interestingly, to the multiple combinations of stimuli, the intracellular calcium dynamics reveal two main modes: unimodal and oscillatory modes, showing significant dependence on the features of the spatio-temporal ATP and shear stress stimuli. The present study provides essential information for controlling calcium dynamics by regulating spatio-temporal biochemical and biomechanical stimuli, which shows the potential in directing cellular activities and understanding the occurrence and development of disease.

The intracellular calcium dynamics in a single vascular endothelial cell being squeezed through a narrow microfluidic channel

Revealing the mechanisms underlying the intracellular calcium responses in vascular endothelial cells (VECs) induced by mechanical stimuli contributes to a better understanding for vascular diseases, including hypertension, atherosclerosis, and aneurysm. Combining with experimental measurement and Computational Fluid Dynamics simulation, we developed a mechanobiological model to investigate the intracellular [Ca²⁺] response in a single VEC being squeezed through narrow microfluidic channel. The time-dependent cellular surface tension dynamics was quantified throughout the squeezing process. In our model, the various Ca²⁺ signaling pathways activated by mechanical stimulation is fully considered. The simulation results of our model exhibited well agreement with our experimental results. By using the model, we theoretically explored the mechanism of the two-peak intracellular [Ca²⁺] response in single VEC being squeezed through narrow channel and made some testable predictions for guiding experiment in the future.

Effect of spatial heterogeneity and colocalization of eNOS and capacitative calcium entry channels on shear stress-induced NO production by endothelial cells: A modeling approach

Introduction

Colocalization of endothelial nitric oxide synthase (eNOS) and capacitative Ca²⁺ entry (CCE) channels in microdomains such as cavaeolae in endothelial cells (ECs) has been shown to significantly affect intracellular Ca²⁺ dynamics and NO production, but the effect has not been well quantified.

Methods

We developed a two-dimensional continuum model of an EC integrating shear stress-mediated ATP production, intracellular Ca²⁺ mobilization, and eNOS activation to investigate the effects of spatial colocalization of plasma membrane eNOS and CCE channels on Ca²⁺ dynamics and NO production in response to flow-induced shear stress. Our model examines the hypothesis that subcellular colocalization of cellular components can be critical for optimal coupling of NO production to blood flow.

Results

Our simulations predict that heterogeneity of CCE can result in formation of microdomains with significantly higher Ca²⁺ compared to the average cytosolic Ca²⁺. Ca²⁺ buffers with lower or no mobility further enhanced Ca²⁺ gradients relative to mobile buffers. Colocalization of eNOS to CCE channels significantly increased NO production.

Conclusions

Our results provide quantitative understanding for the role of spatial heterogeneity and the compartmentalization of signals in regulation of shear stress-induced NO production.

A microfluidic device with spatiotemporal wall shear stress and ATP signals to investigate the intracellular calcium dynamics in vascular endothelial cells

Intracellular calcium dynamics plays an important role in the regulation of vascular endothelial cellular functions. In order to probe the intracellular calcium dynamic response under synergistic effect of wall shear stress (WSS) and adenosine triphosphate (ATP) signals, a novel microfluidic device, which provides the adherent vascular endothelial cells (VECs) on the bottom of microchannel with WSS signal alone, ATP signal alone, and different combinations of WSS and ATP signals, is proposed based upon the principles of fluid mechanics and mass transfer. The spatiotemporal profiles of extracellular ATP signals from numerical simulation and experiment studies validate the implementation of our design. The intracellular calcium dynamics of VECs in response to either WSS signal or ATP signal alone, and different combinations of WSS and ATP signals have been investigated. It is found that the synergistic effect of the WSS and ATP signals plays a more significant role in the signal transduction of VECs rather than that from either WSS signal or ATP signal alone. In particular, under the combined stimuli of WSS and ATP signals with different amplitudes and frequencies, the amplitudes and frequencies of the intracellular Ca²⁺ dynamic signals are observed to be closely related to the amplitudes and frequencies of WSS or ATP signals.

A multi-component parallel-plate flow chamber system for studying the effect of exercise-induced wall shear stress on endothelial cells

BACKGROUND

In vivo studies have demonstrated that reasonable exercise training can improve endothelial function. To confirm the key role of wall shear stress induced by exercise on endothelial cells, and to understand how wall shear stress affects the structure and the function of endothelial cells, it is crucial to design and fabricate an in vitro multi-component parallel-plate flow chamber system which can closely replicate exercise-induced wall shear stress waveforms in artery.

METHODS

The in vivo wall shear stress waveforms from the common carotid artery of a healthy volunteer in resting and immediately after 30 min acute aerobic cycling exercise were first calculated by measuring the inner diameter and the center-line blood flow velocity with a color Doppler ultrasound. According to the above in vivo wall shear stress waveforms, we designed and fabricated a parallel-plate flow chamber system with appropriate components based on a lumped parameter hemodynamics model. To validate the feasibility of this system, human umbilical vein endothelial cells (HUVECs) line were cultured within the parallel-plate flow chamber under abovementioned two types of wall shear stress waveforms and the intracellular actin microfilaments and nitric oxide (NO) production level were evaluated using fluorescence microscope.

RESULTS

Our results show that the trends of resting and exercise-induced wall shear stress waveforms, especially the maximal, minimal and mean wall shear stress as well as oscillatory shear index, generated by the parallel-plate flow chamber system are similar to those acquired from the common carotid artery. In addition, the cellular experiments demonstrate that the actin microfilaments and the production of NO within cells exposed to the two different wall shear stress waveforms exhibit different dynamic behaviors; there are larger numbers of actin microfilaments and higher level NO in cells exposed in exercise-induced wall shear stress condition than resting wall shear stress condition.

CONCLUSION

The parallel-plate flow chamber system can well reproduce wall shear stress waveforms acquired from the common carotid artery in resting and immediately after exercise states. Furthermore, it can be used for studying the endothelial cells responses under resting and exercise-induced wall shear stress environments in vitro.

Mathematical model for shear stress dependent NO and adenine nucleotide production from endothelial cells

We developed a mass transport model for a parallel-plate flow chamber apparatus to predict the concentrations of nitric oxide (NO) and adenine nucleotides (ATP, ADP) produced by cultured endothelial cells (ECs) and investigated how the net rates of production, degradation, and mass transport for these three chemical species vary with changes in wall shear stress (τw). These simulations provide an improved understanding of experimental results obtained with parallel-plate flow chambers and allows quantitative analysis of the relationship between τw, adenine nucleotide concentrations, and NO produced by ECs. Experimental data obtained after altering ATP and ADP concentrations with apyrase were analyzed to quantify changes in the rate of NO production (RNO). The effects of different isoforms of apyrase on ATP and ADP concentrations and nucleotide-dependent changes in RNO could be predicted with the model. A decrease in ATP was predicted with apyrase, but an increase in ADP was simulated due to degradation of ATP. We found that a simple proportional relationship relating a component of RNO to the sum of ATP and ADP provided a close match to the fitted curve for experimentally measured changes in RNO with apyrase. Estimates for the proportionality constant ranged from 0.0067 to 0.0321 μM/s increase in RNO per nM nucleotide concentration, depending on which isoform of apyrase was modeled, with the largest effect of nucleotides on RNO at low τw (<6 dyn/cm(2)).

Modeling of TRPV₄-C₁ -mediated calcium signaling in vascular endothelial cells induced by fluid shear stress and ATP

The calcium signaling plays a vital role in flow-dependent vascular endothelial cell (VEC) physiology. Variations in fluid shear stress and ATP concentration in blood vessels can activate dynamic responses of cytosolic-free [Formula: see text] through various calcium channels on the plasma membrane. In this paper, a novel dynamic model has been proposed for transient receptor potential vanilloid 4 [Formula: see text]-mediated intracellular calcium dynamics in VECs induced by fluid shear stress and ATP. Our model includes [Formula: see text] signaling pathways through P2Y receptors and [Formula: see text] channels (indirect mechanism) and captures the roles of the [Formula: see text] compound channels in VEC [Formula: see text] signaling in response to fluid shear stress (direct mechanism). In particular, it takes into account that the [Formula: see text] compound channels are regulated by intracellular [Formula: see text] and [Formula: see text] concentrations. The simulation studies have demonstrated that the dynamic responses of calcium concentration produced by the proposed model correlate well with the existing experimental observations. We also conclude from the simulation studies that endogenously released ATP may play an insignificant role in the process of intracellular [Formula: see text] response to shear stress.

Modeling of progesterone-induced intracellular calcium signaling in human spermatozoa

Calcium ion is a secondary messenger of mammalian spermatozoa. The dynamic change of its concentration plays a vital role in the process of sperm motility, capacitation, acrosome and fertilization. Progesterone released by the cumulus cells, as a potent stimulator of fertilization, can activate the calcium channels on the plasma membrane, which in turn triggers the dynamic change of intracellular calcium concentration. In this paper, a mathematical model of calcium dynamic response in mammalian spermatozoa induced by progesterone is proposed and numerical simulation of the dynamic model is conducted. The results show that the dynamic response of calcium concentration predicted by the model is in accordance with experimental evidence. The proposed dynamic model can be used to explain the phenomena observed in the experiments and predict new phenomena to be revealed by experimental investigations, which will provide the basis to quantitatively investigate the fluid mechanics and biochemistry for the sperm motility induced by progesterone.

Flow-induced ATP release in patient-specific arterial geometries--a comparative study of computational models

The importance of the endothelium in the local regulation of blood flow is reflected by its influence on vascular tone by means of vasodilatory responses to many physiological stimuli. Regulatory pathways are affected by mass transport and wall shear stress (WSS), via mechanotransduction mechanisms. In the present work, we review the most relevant computational models that have been proposed to date, and introduce a general framework for modelling the responses of the endothelium to alteration in the flow, with a view to understanding the biomechanical processes involved in the pathways to endothelial dysfunction. Simulations are performed on two different patient-specific stenosed carotid artery geometries to investigate the influence of WSS and mass transport phenomena upon the agonist coupling response at the endothelium. In particular, results presented for two different models of WSS-dependent adenosine-5'-triphosphate (ATP) release reveal that existing paradigms may not account for the conditions encountered in vivo and may therefore not be adequate to model the kinetics of ATP at the endothelium.

Local detection of mechanically induced ATP release from bone cells with ATP microbiosensors

The mechanically induced release of adenosine-5'-triphosphate (ATP) from osteoblastic cells (MC3T3-E1) was measured in real time. A stretching device integrated into scanning electrochemical microscopy was developed to apply controlled mechanical strain to MC3T3-E1 cells. For ATP secretion, a stepwise yet uniform mechanical stress was imposed onto MC3T3-E1 cells. The ATP biosensors were positioned at a distance of approximately 30-40 μm above the cell surface. Calibration functions were recorded prior to the cell measurements and revealed a linear response up to 40 μM with a sensitivity of 1-5pA/μM ATP. Stretching MC3T3-E1 cells up to 21% resulted in a concentration of 30.57±4.82 μM of extracellular ATP (N=12) detected above the cell surface. As a control experiment, nifedipine, a L-type voltage sensitive calcium channel (L-VSCC) inhibitor was applied, which blocks Ca(2+)entry from the outer medium into the cell. Inhibition resulted in a significantly smaller amount of released ATP, i.e., 7.08±1.93 μM ATP (N=10). Further control experiments with glucose microbiosensors did not yield significant changes of the baseline current (N=8).

Flow-dependent mass transfer may trigger endothelial signaling cascades

It is well known that fluid mechanical forces directly impact endothelial signaling pathways. But while this general observation is clear, less apparent are the underlying mechanisms that initiate these critical signaling processes. This is because fluid mechanical forces can offer a direct mechanical input to possible mechanotransducers as well as alter critical mass transport characteristics (i.e., concentration gradients) of a host of chemical stimuli present in the blood stream. However, it has recently been accepted that mechanotransduction (direct mechanical force input), and not mass transfer, is the fundamental mechanism for many hemodynamic force-modulated endothelial signaling pathways and their downstream gene products. This conclusion has been largely based, indirectly, on accepted criteria that correlate signaling behavior and shear rate and shear stress, relative to changes in viscosity. However, in this work, we investigate the negative control for these criteria. Here we computationally and experimentally subject mass-transfer limited systems, independent of mechanotransduction, to the purported criteria. The results showed that the negative control (mass-transfer limited system) produced the same trends that have been used to identify mechanotransduction-dominant systems. Thus, the widely used viscosity-related shear stress and shear rate criteria are insufficient in determining mechanotransduction-dominant systems. Thus, research should continue to consider the importance of mass transfer in triggering signaling cascades.

Fluid flow induced calcium response in osteoblasts: mathematical modeling

Fluid flow in the bone lacuno-canalicular network can induce dynamic fluctuation of intracellular calcium concentration ([Ca(2+)](i)) in osteoblasts, which plays an important role in bone remodeling. There has been limited progress in the mathematical modeling of this process probably due to its complexity, which is controlled by various factors such as Ca(2+) channels and extracellular messengers. In this study we developed a mathematical model to describe [Ca(2+)](i) response induced by fluid shear stress (SS) by integrating the major factors involved and analyzed the effects of different experimental setups (e.g. [Ca(2+)](i) baseline, pretreatment with ATP). In this model we considered the ATP release process and the activities of multiple ion channels and purinergic receptors. The model was further verified quantitatively by comparing the simulation results with experimental data reported in literature. The results showed that: (i) extracellular ATP concentration has more significant effect on [Ca(2+)](i) baseline (73% increase in [Ca(2+)](i) with extracellular ATP concentration varying between 0 and 10 μM), as compared to that induced by SS (25% variation in [Ca(2+)](i) with SS varying from 0 to 3.5 Pa); (ii) Pretreatment with ATP-medium results in different [Ca(2+)](i) response as compared to the control group (ATP-free medium) under SS; (iii) Relative [Ca(2+)](i) fluctuation over baseline is more reliable to show the [Ca(2+)](i) response process than the absolute [Ca(2+)](i) response peak. The developed model may improve the experimental design and facilitate our understanding of the mechanotransduction process in osteoblasts.

Extracellular ATP may contribute to tissue repair by rapidly stimulating purinergic receptor X7-dependent vascular endothelial growth factor release from primary human monocytes

Extracellular ATP has been proposed to act as a danger signal to alert the immune system of cell damage. Release of high local concentrations of ATP activates the nucleotide receptor, purinergic receptor X7 (P2RX7), on monocytic cells, which promotes the processing/release of proinflammatory mediators. Although the proinflammatory actions of P2RX7 are well recognized, little is known regarding the potential function of P2RX7 in repair responses. Because the resolution of inflammation is characterized by monocytic cell-dependent production of proangiogenic factors, we evaluated the contribution of P2RX7 to this process. We observed that both short-term and long-term P2RX7 activation promotes the robust release of vascular endothelial growth factor from primary human monocytes. This vascular endothelial growth factor release is calcium dependent and associated with reactive oxygen species production. This previously unrecognized action of P2RX7 suggests that it may not only participate in inflammation and cell death, but that it is also likely to be important in the control of angiogenesis and wound repair.