Focal adhesion clustering drives endothelial cell morphology on patterned surfaces.

In many cell types, shape and function are intertwined. In vivo, vascular endothelial cells (ECs) are typically elongated and aligned in the direction of blood flow; however, near branches and bifurcations where atherosclerosis develops, ECs are often cuboidal and have no preferred orientation. Thus, understanding the factors that regulate EC shape and alignment is important. In vitro, EC morphology and orientation are exquisitely sensitive to the composition and topography of the substrate on which the cells are cultured; however, the underlying mechanisms remain poorly understood. Different strategies of substrate patterning for regulating EC shape and orientation have been reported including adhesive motifs on planar surfaces and micro- or nano-scale gratings that provide substrate topography. Here, we explore how ECs perceive planar bio-adhesive versus microgrooved topographic surfaces having identical feature dimensions. We show that while the two types of patterned surfaces are equally effective in guiding and directing EC orientation, the cells are considerably more elongated on the planar patterned surfaces than on the microgrooved surfaces. We also demonstrate that the key factor that regulates cellular morphology is focal adhesion clustering which subsequently drives cytoskeletal organization. The present results promise to inform design strategies of novel surfaces for the improved performance of implantable cardiovascular devices.

A Mathematical Model for the Sounds Produced by Knuckle Cracking.

The articular release of the metacarpophalangeal joint produces a typical cracking sound, resulting in what is commonly referred to as the cracking of knuckles. Despite over sixty years of research, the source of the knuckle cracking sound continues to be debated due to inconclusive experimental evidence as a result of limitations in the temporal resolution of non-invasive physiological imaging techniques. To support the available experimental data and shed light onto the source of the cracking sound, we have developed a mathematical model of the events leading to the generation of the sound. The model resolves the dynamics of a collapsing cavitation bubble in the synovial fluid inside a metacarpophalangeal joint during an articular release. The acoustic signature from the resulting bubble dynamics is shown to be consistent in both magnitude and dominant frequency with experimental measurements in the literature and with our own experiments, thus lending support for cavitation bubble collapse as the source of the cracking sound. Finally, the model also shows that only a partial collapse of the bubble is needed to replicate the experimentally observed acoustic spectra, thus allowing for bubbles to persist following the generation of sound as has been reported in recent experiments.

Dynamic monitoring of cell mechanical properties using profile microindentation.

We have developed a simple and relatively inexpensive system to visualize adherent cells in profile while measuring their mechanical properties using microindentation. The setup allows simultaneous control of cell microenvironment by introducing a micropipette for the delivery of soluble factors or other cell types. We validate this technique against atomic force microscopy measurements and, as a proof of concept, measure the viscoelastic properties of vascular endothelial cells in terms of an apparent stiffness and a dimensionless parameter that describes stress relaxation. Furthermore, we use this technique to monitor the time evolution of these mechanical properties as the cells' actin is depolymerized using cytochalasin-D.

Computational Fluid Dynamic Simulations of Maternal Circulation: Wall Shear Stress in the Human Placenta and Its Biological Implications.

INTRODUCTION: In the human placenta the maternal blood circulates in the intervillous space (IVS). The syncytiotrophoblast (STB) is in direct contact with maternal blood. The wall shear stress (WSS) exerted by the maternal blood flow on the STB has not been evaluated. Our objective was to determine the physiological WSS exerted on the surface of the STB during the third trimester of pregnancy. MATERIAL AND METHODS: To gain insight into the shear stress levels that the STB is expected to experience in vivo, we have formulated three different computational models of varying levels of complexity that reflect different physical representations of the IVS. Computations of the flow fields in all models were performed using the CFD module of the finite element code COMSOL Multiphysics 4.4. The mean velocity of maternal blood in the IVS during the third trimester was measured in vivo with dynamic MRI (0.94±0.14 mm.s-1). To investigate if the in silico results are consistent with physiological observations, we studied the cytoadhesion of human parasitized (Plasmodium falciparum) erythrocytes to primary human STB cultures, in flow conditions with different WSS values. RESULTS: The WSS applied to the STB is highly heterogeneous in the IVS. The estimated average values are relatively low (0.5±0.2 to 2.3±1.1 dyn.cm-2). The increase of WSS from 0.15 to 5 dyn.cm-2 was associated with a significant decrease of infected erythrocyte cytoadhesion. No cytoadhesion of infected erythrocytes was observed above 5 dyn.cm-2 applied for one hour. CONCLUSION: Our study provides for the first time a WSS estimation in the maternal placental circulation. In spite of high maternal blood flow rates, the average WSS applied at the surface of the chorionic villi is low (<5 dyn.cm-2). These results provide the basis for future physiologically-relevant in vitro studies of the biological effects of WSS on the STB.

Modulation of ATP/ADP concentration at the endothelial surface by shear stress: effect of flow disturbance.

The adenine nucleotides ATP and ADP modulate the release of endothelial-derived relaxing factors and hence play an important role in flow-mediated arterial vasoregulation. Adenine nucleotide concentration at the endothelial cell (EC) surface within an artery is determined by a balance of convective-diffusive delivery of blood-borne nucleotides to the EC surface, hydrolysis of these nucleotides at the cell surface, and flow-induced ATP release from ECs. Previous numerical simulations in a parallel plate flow chamber had demonstrated that flow-induced ATP release has a profound effect on nucleotide concentration under both steady and pulsatile flow conditions. In the present study, we have extended this analysis to probe the impact of disturbed flow downstream of a backward facing step on adenine nucleotide concentration at the EC surface. The results have demonstrated that over a wide range of applied wall shear stress, the ATP concentration at the EC surface drops abruptly within the disturbed flow zone due to increased nucleotide residence time within this region. The concentration is intricately sensitive to the kinetics of flow-induced ATP release, and this sensitivity is more pronounced at lower levels of wall shear stress.

Unsteady and three-dimensional simulation of blood flow in the human aortic arch.

A three-dimensional and pulsatile blood flow in a human aortic arch and its three major branches has been studied numerically for a peak Reynolds number of 2500 and a frequency (or Womersley) parameter of 10. The simulation geometry was derived from the three-dimensional reconstruction of a series of two-dimensional slices obtained in vivo using CAT scan imaging on a human aorta. The numerical simulations were obtained using a projection method, and a finite-volume formulation of the Navier-Stokes equations was used on a system of overset grids. Our results demonstrate that the primary flow velocity is skewed towards the inner aortic wall in the ascending aorta, but this skewness shifts to the outer wall in the descending thoracic aorta. Within the arch branches, the flow velocities were skewed to the distal walls with flow reversal along the proximal walls. Extensive secondary flow motion was observed in the aorta, and the structure of these secondary flows was influenced considerably by the presence of the branches. Within the aorta, wall shear stresses were highly dynamic, but were generally high along the outer wall in the vicinity of the branches and low along the inner wall, particularly in the descending thoracic aorta. Within the branches, the shear stresses were considerably higher along the distal walls than along the proximal walls. Wall pressure was low along the inner aortic wall and high around the branches and along the outer wall in the ascending thoracic aorta. Comparison of our numerical results with the localization of early atherosclerotic lesions broadly suggests preferential development of these lesions in regions of extrema (either maxima or minima) in wall shear stress and pressure.

Modulation of ATP/ADP concentration at the endothelial surface by shear stress: effect of flow-induced ATP release.

The adenine nucleotides ATP and ADP induce the production of vasoactive compounds in vascular endothelial cells (ECs). Therefore, knowledge of how flow affects the concentration of ATP and ADP at the EC surface may be important for understanding shear stress-mediated vasoregulation. The concentration of ATP and ADP is determined by convective and diffusive transport as well as by hydrolysis of these nucleotides by ectonucleotidases at the EC surface. Previous mathematical modeling has demonstrated that for steady flow in a parallel plate flow chamber, the combined ATP+ADP concentration does not change considerably over a wide range of shear stress. This finding has been used to argue that the effect of flow on adenine nucleotide transport could not account for the dependence of endothelial responses to ATP on the magnitude of applied shear stress. The present study extends the previous modeling to include pulsatile flow as well as flow-induced endothelial ATP release. Our results demonstrate that flow-induced ATP release has a pronounced effect on nucleotide concentration under both steady and pulsatile flow conditions. While the combined ATP+ADP concentration at the EC surface in the absence of flow-induced ATP release changes by only approximately 10% over the wall shear stress range 0.1-10 dyne cm(-2), inclusion of this release leads to a concentration change of approximately 34%-106% over the same shear stress range, depending on how ATP release is modeled. These results suggest that the dependence of various endothelial responses to shear stress on the magnitude of the applied shear stress may be partially attributable to flow-induced changes in cell-surface adenine nucleotide concentration.

Endothelial cellular response to altered shear stress.

Endothelial cells are normally exposed constantly to mechanical forces that significantly influence their phenotype. This symposium presented recent information concerning endothelial cell responses to shear stress associated with blood flow. Endothelial cell shear stress mechanosensors that have been proposed include membrane receptor kinases, integrins, G proteins, ion channels, intercellular junction proteins, membrane lipids (e.g., those associated with caveolae), and the cytoskeleton. These sensors are linked to signaling cascades that interact with or result in generation of reactive oxygen species, nitric oxide, and various transcription factors among other responses. Endothelial cells adapt to sustained shear stress, and either an increase or decrease from normal shear leads to signaling events. In vitro models for the study of endothelial cell responses must consider the pattern of shear stress (e.g., steady vs. oscillatory flow), the scaffold for cell growth (e.g., basement membrane or other cell types such as smooth muscle cells), and the extent of flow adaptation. These cellular responses have major relevance for understanding the pathophysiological effects of increased shear stress associated with hypertension or decreased shear stress associated with thrombotic occlusion.

A model for shear stress-induced deformation of a flow sensor on the surface of vascular endothelial cells.

Fluid mechanical shear stress elicits humoral, metabolic, and structural responses in vascular endothelial cells (ECs); however, the mechanisms involved in shear stress sensing and transduction remain incompletely understood. Beyond being responsive to shear stress, ECs distinguish among and respond differently to different types of shear stress. Recent observations suggest that endothelial shear stress sensing may occur through direct interaction of the flow with cell-surface structures that act as primary flow sensors. This paper presents a mathematical model for the shear stress-induced deformation of a flow sensor on the EC surface. The sensor is modeled as a cytoskeleton-coupled viscoelastic structure exhibiting standard linear solid behavior. Since ECs respond differently to different types of flow, the deformation and resulting velocity of the sensor in response to steady, non-reversing pulsatile, and oscillatory flow have been studied. Furthermore, the sensitivity of the results to changes in various model parameters including the magnitude of applied shear stress, the constants that characterize the viscoelastic behavior, and the pulsatile flow frequency (f) has been investigated. The results have demonstrated that in response to a suddenly applied shear stress, the sensor exhibits a level of instantaneous deformation followed by gradual creeping to the long-term response. The peak deformation increases linearly with the magnitude of the applied shear stress and decreases for viscoelastic constants that correspond to stiffer sensors. While the sensor deformation depends on f for low f values, the deformation becomes f -independent above a critical threshold frequency. Finally, the peak sensor deformation is considerably larger for steady and non-reversing pulsatile flow than for oscillatory flow. If the extent of sensor deformation correlates with the intensity of flow-mediated endothelial signaling, then our results suggest possible mechanisms by which ECs distinguish among steady, non-reversing pulsatile, and oscillatory shear stress.

Flow-induced expression of endothelial Na-K-Cl cotransport: dependence on K(+) and Cl(-) channels.

Steady laminar shear stress has been shown previously to markedly increase Na-K-Cl cotransporter mRNA and protein in human umbilical vein endothelial cells and also to rapidly increase endothelial K(+) and Cl(-) channel conductances. The present study was done to evaluate the effects of shear stress on Na-K-Cl cotransporter activity and protein expression in bovine aortic endothelial cells (BAEC) and to determine whether changes in cotransporter expression may be dependent on early changes in K(+) and Cl(-) channel conductances. Confluent BAEC monolayers were exposed in a parallel-plate flow chamber to either steady shear stress (19 dyn/cm(2)) or purely oscillatory shear stress (0 +/- 19 dyn/cm(2)) for 6-48 h. After shearing, BAEC monolayers were assessed for Na-K-Cl cotransporter activity or were subjected to Western blot analysis of cotransporter protein. Steady shear stress led to a 2- to 4-fold increase in BAEC cotransporter protein levels and a 1.5- to 1.8-fold increase in cotransporter activity, increases that were sustained over the longest time periods studied. Oscillatory flow, in contrast, had no effect on cotransporter protein levels. In the presence of flow-sensitive K(+) and Cl(-) channel pharmacological blockers, the steady shear stress-induced increase in cotransporter protein was virtually abolished. These results suggest that shear stress modulates the expression of the BAEC Na-K-Cl cotransporter by mechanisms that are dependent on flow-activated ion channels.

Influence of different forms of fluid shear stress on vascular endothelial TGF-beta1 mRNA expression.

This study investigated the effects of long-term exposure of steady (19 dyne/cm2), 1-Hz non-reversing pulsatile (19+/-6 dyne/cm2) and 1-Hz purely oscillatory (0+/-19 dyne/cm2) shear stress on endothelial transforming growth factor-beta1 (TGF-beta1) mRNA expression. Cultured bovine aortic endothelial cells (BAECs) were systematically exposed to the three flow conditions for periods of 2, 6, 12 and 24 h, and relative differences in TGF-beta1 mRNA levels were measured by semi-quantitative RT-PCR. In response to steady shear stress, TGF-beta1 mRNA levels normalized to no-flow controls were 1.24, 1.42, 1.30 and 1.47 at the 2-, 6-, 12- and 24-h time points, respectively. In response to non-reversing pulsatile flow, these levels were 1.49, 1.64, 1.64 and 1.73, while the respective transcript levels for oscillatory flow were 1.33, 1.12, 1.12 and 1. 93. These results indicate that BAEC TGF-beta1 mRNA was up-regulated with the kinetics of the up-regulation faster for steady and non-reversing pulsatile flow than for oscillatory flow. Given the preferential localization of early atherosclerotic lesions in arterial regions exposed to low and/or oscillatory shear stress and the implication of TGF-beta1 as an athero-protective gene, these results are consistent with the notion that regions transiently exposed to oscillatory flow may be particularly prone to atherosclerosis.

A flow-activated chloride-selective membrane current in vascular endothelial cells.

Shear stress-induced activation of endothelial ion channels, one of the earliest responses to flow, is implicated in mechano-signal transduction that results in the regulation of vascular tone. The effects of laminar flow on endothelial membrane potential were studied in vitro using both fluorescent potentiometric dye measurements and whole-cell patch-clamp recordings. The application of flow stimulated membrane hyperpolarization, which was reversed to depolarization within 35 to 160 seconds. The depolarization was caused by a Cl(-)-selective membrane current activated by flow independently of the K(+) channel-mediated hyperpolarization. Thus, flow activated both K(+) and Cl(-) currents, with the net membrane potential being determined by the balance of the responses. Membrane potential sensitivity to flow was unchanged by flow preconditioning that elongated and aligned the cells.

Responsiveness of vascular endothelium to shear stress: potential role of ion channels and cellular cytoskeleton (review).

The frictional forces associated with blood flow expose vascular endothelium in arteries to a complex and highly dynamic shear stress distribution. The ability of endothelial cells to respond to shear stress is essential for arterial vasoregulation in response to acute hemodynamic changes and for vascular wall remodeling following chronic changes in blood flow. Furthermore, endothelial responsiveness to shear stress may play a role in the localization of early atherosclerotic lesions. Shear stress elicits a wide range of humoral, metabolic, and structural responses in endothelial cells. These include activation of ion channels and of G proteins, induction of oscillations in intracellular calcium concentration, alterations in the expression of various important genes, and extensive cytoskeletal reorganization. Mechanisms of shear stress sensing and transmission in endothelium are discussed in light of the complex shear stress distribution to which endothelial cells are exposed in vivo and with particular emphasis on the potentially central role of flow-sensitive ion channels and the cellular cytoskeleton. Finally, the ability of endothelial cells to distinguish among and to respond differentially to different types of shear stress is highlighted.

Computational study of the effect of geometric and flow parameters on the steady flow field at the rabbit aorto-celiac bifurcation.

Arterial hemodynamic forces may play a role in the localization of early atherosclerotic lesions. We have been developing numerical techniques based on overset or "Chimera" type formulations to solve the Navier-Stokes equations in complex geometries simulating arterial bifurcations. This paper presents three-dimensional steady flow computations in a model of the rabbit aorto-celiac bifurcation. The computational methods were validated by comparing the numerical results to previously-obtained flow visualization data. Once validated, the numerical algorithms were used to investigate the sensitivity of the computed flow field and resulting wall shear stress distribution to various geometric and hemodynamic parameters. The results demonstrated that a decrease in the extent of aortic taper downstream of the celiac artery induced looping fluid motion along the lateral walls of the aorta and shifted the peak wall shear stress from downstream of the celiac artery to upstream. Increasing the flow Reynolds number led to a sharp increase in spatial gradients of wall shear stress. The flow field was highly sensitive to the flow division ratio, i.e., the fraction of total flow rate that enters the celiac artery, with larger values of this ratio leading to the occurrence of flow separation along the dorsal wall of the aorta. Finally, skewness of the inlet velocity profile had a profound impact on the wall shear stress distribution near the celiac artery. While not physiological due to the assumption of steady flow, these results provide valuable insight into the fluid physics at geometries simulating arterial bifurcations.

Microcinematographic studies of flow patterns in the excised rabbit aorta and its major branches.

Arterial fluid mechanics may play a role as a localizing factor for early atherosclerosis. Flow patterns in natural rabbit aortas rendered transparent were studied using a microcinematographic visualization technique. The aortic arch exhibited a single cell of clockwise-rotating helical secondary flow along the ventral and inner walls. Flow separation occurred proximal to the two arch branches with flow reversal proximal to the brachiocephalic artery. Sinusoidal flow rendered the helical motion more pronounced in systole, while the reverse flow zone periodically expanded and contracted. Steady flow in the abdominal aorta revealed streamlines which follow slow looping trajectories lateral to ostia before tracing helical paths into the branches. Flow separation was present along the dorsal wall of the aorta opposite the superior mesenteric artery. With the exception of the left renal artery, steady flow wall shear stresses were higher distal to ostia than proximal. Spatial gradients of wall shear stress were larger around branches than elsewhere. Similar to observed flow patterns, sites of enhanced macromolecular permeability, as observed previously in the normal rabbit aorta, follow a clockwise helical pattern in the arch and exhibit a distribution around ostia that correlates to some degree with regions of elevated shear stress gradients.

Measurement of flow rates through aortic branches in the anesthetized rabbit.

The volumetric flow rates, mean and pulsatile, in the aorta and its major branches were measured in nonfed, anesthetized rabbits, using a transit time Doppler ultrasonic flowmeter. Anesthesia was maintained with isoflurane, and a vasodilator was applied topically during the measurements to avoid introducing additional flow resistance due to vasoconstriction. The cranial mesenteric and celiac arteries received the bulk of the aortic flow, (mean +/- SD) 29.5 +/- 6.6% and 23.3 +/- 5.8%, respectively, for mean flow. The brachiocephalic artery received as much as 14.7 +/- 3.2%, while each of the other branches received a considerably smaller fraction: 7.1 +/- 2.5% for the left subclavian artery, 6.2 +/- 2.6% and 5.1 +/- 2.2%, respectively, for the right and left renal arteries, and 6.0 +/- 2.5% for each of the two iliac arteries. Flow divisions were nearly the same in paired vessels. Peak pulsatile flow divisions were similar to their steady flow counterparts in the brachiocephalic, left subclavian, celiac, and cranial mesenteric arteries, but were smaller in the renal and iliac arteries, although the difference was not statistically significant. Reverse flow from one or more of the branches back into the aorta occurred in diastole in seven of eight rabbits studied.

Spatial relationships in early signaling events of flow-mediated endothelial mechanotransduction.

Blood flow interactions with the vascular endothelium represent a specialized example of mechanical regulation of cell function that has important physiological and pathological cardiovascular consequences. The endothelial monolayer in vivo acts as a signal transduction interface for forces associated with flowing blood (hemodynamic forces) in the acute regulation of artery tone and chronic structural remodeling of arteries, including the pathology of atherosclerosis. Mechanisms related to spatial relationships at the cell surfaces and throughout the cell that influence flow-mediated endothelial mechanotransduction are discussed. In particular, flow-mediated ion channel activation and cytoskeletal dynamics are considered in relation to topographic analyses of the luminal and abluminal surfaces of living endothelial cells.

Topographical mapping of sites of enhanced HRP permeability in the normal rabbit aorta.

The spatial distribution of sites of enhanced permeability to the macromolecule horseradish peroxidase (HRP) in the normal rabbit aorta after one min circulation was studied using image analysis. These sites, referred to as "HRP spots," exhibit a nonuniform distribution that is qualitatively similar in all rabbits studied. The density of HRP spots is highest in the aortic arch, decreases distally, reaches a minimum in the lower descending thoracic aorta, and then increases again in the abdominal aorta. The region of highest spot density follows a clockwise helical pattern in the aortic arch and outside the arch occurs in streaks largely oriented in the bulk flow direction. The streaks in the abdominal aorta localize along the anatomical right lateral wall and occasionally along the left lateral wall proximal to the celiac artery and along the ventral wall between the celiac and superior mesenteric arteries. The density of spots is high in the immediate vicinity of aortic ostia with the most elevated density being distal to ostia in most cases. At a short distance from the ostium edge of the celiac and superior mesenteric branches the proximal density is comparably high, and no preferred spot orientation is observed around the brachiocephalic vessel. These results are consistent with an influence of localizing factors such as detailed hemodynamic phenomena and/or arterial wall structural and/or functional variations.