Effects of biaxial oscillatory shear stress on endothelial cell proliferation and morphology

A comprehensive MRI-based computational model of blood flow in compliant aorta using radial basis function interpolation

BACKGROUND

Properly understanding the origin and progression of the thoracic aortic aneurysm (TAA) can help prevent its growth and rupture. For a better understanding of this pathogenesis, the aortic blood flow has to be studied and interpreted in great detail. We can obtain detailed aortic blood flow information using magnetic resonance imaging (MRI) based computational fluid dynamics (CFD) with a prescribed motion of the aortic wall.

METHODS

We performed two different types of simulations-static (rigid wall) and dynamic (moving wall) for healthy control and a patient with a TAA. For the latter, we have developed a novel morphing approach based on the radial basis function (RBF) interpolation of the segmented 4D-flow MRI geometries at different time instants. Additionally, we have applied reconstructed 4D-flow MRI velocity profiles at the inlet with an automatic registration protocol.

RESULTS

The simulated RBF-based movement of the aorta matched well with the original 4D-flow MRI geometries. The wall movement was most dominant in the ascending aorta, accompanied by the highest variation of the blood flow patterns. The resulting data indicated significant differences between the dynamic and static simulations, with a relative difference for the patient of 7.47±14.18% in time-averaged wall shear stress and 15.97±43.32% in the oscillatory shear index (for the whole domain).

CONCLUSIONS

In conclusion, the RBF-based morphing approach proved to be numerically accurate and computationally efficient in capturing complex kinematics of the aorta, as validated by 4D-flow MRI. We recommend this approach for future use in MRI-based CFD simulations in broad population studies. Performing these would bring a better understanding of the onset and growth of TAA.

Analyzing the effects of helical flow in blood vessels using acoustofluidic-based dynamic flow generator

The effects of helical flow in a blood vessel are investigated in a dynamic flow generator using surface acoustic wave (SAW) in the microfluidic device. The SAW, generated by an interdigital transducer (IDT), induces acoustic streaming, resulting in a stable and consistent helical flow pattern in microscale channels. This approach allows rapid development of helical flow within the channel without directly contacting the medium. The precise design of the window enables the creation of distinct unidirectional vortices, which can be controlled by adjusting the amplitude of the SAW. Within this device, optimal operational parameters of the dynamic flow generator to preserve the integrity of endothelial cells are found, and in such settings, the actin filaments within the cells are aligned to the desired state. Our findings reveal that intracellular Ca2+ concentrations vary in response to flow conditions. Specifically, comparable maximum intensity and graphical patterns were observed between low-flow rate helical flow and high-flow rate Hagen-Poiseuille flow. These suggest that the cells respond to the helical flow through mechanosensitive ion channels. Finally, adherence of monocytes is effectively reduced under helical flow conditions in an inflammatory environment, highlighting the atheroprotective role of helical flow. STATEMENT OF SIGNIFICANCE: Helical flow in blood vessels is well known to prevent atherosclerosis. However, despite efforts to replicate helical flow in microscale channels, there is still a lack of in vitro models which can generate helical flow for analyzing its effects on the vascular system. In this study, we developed a method for generating steady and constant helical flow in microfluidic channel using acoustofluidic techniques. By utilizing this dynamic flow generator, we were able to observe the atheroprotective aspects of helical flow in vitro, including the enhancement of calcium ion flux and reduction of monocyte adhesion. This study paves the way for an in vitro model of dynamic cell culture and offers advanced investigation into helical flow in our circulatory system.

Impact of disturbed flow and arterial stiffening on mechanotransduction in endothelial cells

Disturbed flow promotes progression of atherosclerosis at particular regions of arteries where the recent studies show the arterial wall becomes stiffer. Objective of this study is to show how mechanotransduction in subcellular organelles of endothelial cells (ECs) will alter with changes in blood flow profiles applied on ECs surface and mechanical properties of arterial wall where ECs are attached to. We will examine the exposure of ECs to atherogenic flow profiles (disturbed flow) and non-atherogenic flow profiles (purely forward flow), while stiffness and viscoelasticity of arterial wall will change. A multicomponent model of endothelial cell monolayer was applied to quantify the response of subcellular organelles to the changes in their microenvironment. Our results show that arterial stiffening alters mechanotransduction in intra/inter-cellular organelles of ECs by slight increase in the transmitted stresses, particularly over central stress fibers (SFs). We also observed that degradation of glycocalyx and exposure to non-atherogenic flow profiles result in significantly higher stresses in subcellular organelles, while degradation of glycocalyx and exposure to atherogenic flow profiles result in dramatically lower stresses in the organelles. Moreover, we show that increasing the arterial wall viscoelasticity leads to slight increase in the stresses transmitted to subcellular organelles. FAs are particularly influenced with the changes in the arterial wall properties and viscoelasticity. Our study suggests that changes in viscoelasticity of arterial wall and degradation state of glycocalyx have to be considered along with arterial stiffening in designing more efficient treatment strategies for atherosclerosis. Our study provides insight into significant role of mechanotransduction in the localization of atherosclerosis by quantifying the role of ECs mechanosensors and suggests that mechanotransduction may play a key role in design of more efficient and precision therapeutics to slow down or block the progression of atherosclerosis.

Pump-Less Platform Enables Long-Term Recirculating Perfusion of 3D Printed Tubular Tissues

The direction and pattern of fluid flow affect vascular structure and function, in which vessel-lining endothelial cells exhibit variable cellular morphologies and vessel remodeling by mechanosensing. To recapitulate this microenvironment, some approaches have been reported to successfully apply unidirectional flow on endothelial cells in organ-on-a-chip systems. However, these platforms encounter drawbacks such as the dependency on pumps or confinement to closed microfluidic channels. These constraints impede their synergy with advanced biofabrication techniques like 3D bioprinting, thereby curtailing the potential to introduce greater complexity into engineered tissues. Herein, a pumpless recirculating platform (UniPlate) that enables unidirectional media recirculation through 3D printed tubular tissues, is demonstrated.The device is made of polystyrene via injection molding in combination with 3D printed sacrifical gelatin templates. Tubular blood vessels with unidirectional perfusion are firstly engineered. Then the design is expanded to incorporate duo-recirculating flow for culturing vascularized renal proximal tubules with glucose reabsorption function. In addition to media recirculation, human monocyte recirculation in engineered blood vessels is also demonstrated for over 24 h, with minimal loss of cells, cell viability, and inflammatory activation. UniPlate can be a valuable tool to more precisely control the cellular microenvironment of organ-on-a-chip systems for drug discovery.

Vortical Structures Promote Atheroprotective Wall Shear Stress Distributions in a Carotid Artery Bifurcation Model

Carotid artery diseases, such as atherosclerosis, are a major cause of death in the United States. Wall shear stresses are known to prompt plaque formation, but there is limited understanding of the complex flow structures underlying these stresses and how they differ in a pre-disposed high-risk patient cohort. A 'healthy' and a novel 'pre-disposed' carotid artery bifurcation model was determined based on patient-averaged clinical data, where the 'pre-disposed' model represents a pathological anatomy. Computational fluid dynamic simulations were performed using a physiological flow based on healthy human subjects. A main hairpin vortical structure in the internal carotid artery sinus was observed, which locally increased instantaneous wall shear stress. In the pre-disposed geometry, this vortical structure starts at an earlier instance in the cardiac flow cycle and persists over a much shorter period, where the second half of the cardiac cycle is dominated by perturbed secondary flow structures and vortices. This coincides with weaker favorable axial pressure gradient peaks over the sinus for the 'pre-disposed' geometry. The findings reveal a strong correlation between vortical structures and wall shear stress and imply that an intact internal carotid artery sinus hairpin vortical structure has a physiologically beneficial role by increasing local wall shear stresses. The deterioration of this beneficial vortical structure is expected to play a significant role in atherosclerotic plaque formation.

Modification of the swirling well cell culture model to alter shear stress metrics

Effects of hemodynamic shear stress on endothelial cells have been extensively investigated using the "swirling well" method, in which cells are cultured in dishes or multiwell plates placed on an orbital shaker. A wave rotates around the well, producing complex patterns of shear. The method allows chronic exposure to flow with high throughput at low cost but has two disadvantages: a number of shear stress characteristics change in a broadly similar way from the center to the edge of the well, and cells at one location in the well may release mediators into the medium that affect the behavior of cells at other locations, exposed to different shears. These properties make it challenging to correlate cell properties with shear. The present study investigated simple alterations to ameliorate these issues. Flows were obtained by numerical simulation. Increasing the volume of fluid in the well-altered dimensional but not dimensionless shear metrics. Adding a central cylinder to the base of the well-forced fluid to flow in a square toroidal channel and reduced multidirectionality. Conversely, suspending a cylinder above the base of the well made the flow highly multidirectional. Increasing viscosity in the latter model increased the magnitude of dimensional but not dimensionless metrics. Finally, tilting the well changed the patterns of different wall shear stress metrics in different ways. Collectively, these methods allow similar flows over most of the cells cultured and/or allow the separation of different shear metrics. A combination of the methods overcomes the limitations of the baseline model.

Numerical and Experimental Analysis of Shear Stress Influence on Cellular Viability in Serpentine Vascular Channels

3D bioprinting has emerged as a tool for developing in vitro tissue models for studying disease progression and drug development. The objective of the current study was to evaluate the influence of flow driven shear stress on the viability of cultured cells inside the luminal wall of a serpentine network. Fluid-structure interaction was modeled using COMSOL Multiphysics for representing the elasticity of the serpentine wall. Experimental analysis of the serpentine model was performed on the basis of a desirable inlet flow boundary condition for which the most homogeneously distributed wall shear stress had been obtained from numerical study. A blend of Gelatin-methacryloyl (GelMA) and PEGDA200 PhotoInk was used as a bioink for printing the serpentine network, while facilitating cell growth within the pores of the gelatin substrate. Human umbilical vein endothelial cells were seeded into the channels of the network to simulate the blood vessels. A Live-Dead assay was performed over a period of 14 days to observe the cellular viability in the printed vascular channels. It was observed that cell viability increases when the seeded cells were exposed to the evenly distributed shear stresses at an input flow rate of 4.62 mm/min of the culture media, similar to that predicted in the numerical model with the same inlet boundary condition. It leads to recruitment of a large number of focal adhesion point nodes on cellular membrane, emphasizing the influence of such phenomena on promoting cellular morphologies.

Atherogenic potential of microgravity hemodynamics in the carotid bifurcation: a numerical investigation

Long-duration spaceflight poses multiple hazards to human health, including physiological changes associated with microgravity. The hemodynamic adaptations occurring upon entry into weightlessness have been associated with retrograde stagnant flow conditions and thromboembolic events in the venous vasculature but the impact of microgravity on cerebral arterial hemodynamics and function remains poorly understood. The objective of this study was to quantify the effects of microgravity on hemodynamics and wall shear stress (WSS) characteristics in 16 carotid bifurcation geometries reconstructed from ultrasonography images using computational fluid dynamics modeling. Microgravity resulted in a significant 21% increase in flow stasis index, a 22-23% decrease in WSS magnitude and a 16-26% increase in relative residence time in all bifurcation branches, while preserving WSS unidirectionality. In two anatomies, however, microgravity not only promoted flow stasis but also subjected the convex region of the external carotid arterial wall to a moderate increase in WSS bidirectionality, which contrasted with the population average trend. This study suggests that long-term exposure to microgravity has the potential to subject the vasculature to atheroprone hemodynamics and this effect is modulated by subject-specific anatomical features. The exploration of the biological impact of those microgravity-induced WSS aberrations is needed to better define the risk posed by long spaceflights on cardiovascular health.

NO Synthesis but Not Apoptosis, Mitosis or Inflammation Can Explain Correlations between Flow Directionality and Paracellular Permeability of Cultured Endothelium

Haemodynamic wall shear stress varies from site to site within the arterial system and is thought to cause local variation in endothelial permeability to macromolecules. Our aim was to investigate mechanisms underlying the changes in paracellular permeability caused by different patterns of shear stress in long-term culture. We used the swirling well system and a substrate-binding tracer that permits visualisation of transport at the cellular level. Permeability increased in the centre of swirled wells, where flow is highly multidirectional, and decreased towards the edge, where flow is more uniaxial, compared to static controls. Overall, there was a reduction in permeability. There were also decreases in early- and late-stage apoptosis, proliferation and mitosis, and there were significant correlations between the first three and permeability when considering variation from the centre to the edge under flow. However, data from static controls did not fit the same relation, and a cell-by-cell analysis showed that <5% of uptake under shear was associated with each of these events. Nuclear translocation of NF-κB p65 increased and then decreased with the duration of applied shear, as did permeability, but the spatial correlation between them was not significant. Application of an NO synthase inhibitor abolished the overall decrease in permeability caused by chronic shear and the difference in permeability between the centre and the edge of the well. Hence, shear and paracellular permeability appear to be linked by NO synthesis and not by apoptosis, mitosis or inflammation. The effect was mediated by an increase in transport through tricellular junctions.

Haemodynamic Wall Shear Stress, Endothelial Permeability and Atherosclerosis-A Triad of Controversy

A striking feature of atherosclerosis is its patchy distribution within the vascular system; certain arteries and certain locations within each artery are preferentially affected. Identifying the local risk factors underlying this phenomenon may lead to new therapeutic strategies. The large variation in lesion prevalence in areas of curvature and branching has motivated a search for haemodynamic triggers, particular those related to wall shear stress (WSS). The fact that lesions are rich in blood-derived lipids has motivated studies of local endothelial permeability. However, the location of lesions, the underlying haemodynamic triggers, the role of permeability, the routes by which lipids cross the endothelium, and the mechanisms by which WSS affects permeability have all been areas of controversy. This review presents evidence for and against the current consensus that lesions are triggered by low and/or oscillatory WSS and that this type of shear profile leads to elevated entry of low density lipoprotein (LDL) into the wall via widened intercellular junctions; it also evaluates more recent evidence that lesion location changes with age, that multidirectional shear stress plays a key role, that LDL dominantly crosses the endothelium by transcytosis, and that the link between flow and permeability results from hitherto unrecognised shear-sensitive mediators.

Simple add-on devices to enhance the efficacy of conventional surface immunoassays implemented on standard labware

We propose microfluidic add-ons easily placed on standard assay labware such as microwells and slides to enhance the kinetics of immunoassays. The devices are compatible with mass production, well-established assay protocols and automated platforms.

Comparison of arteriovenous fistulas constructed with main or internal branch of the cephalic vein: a retrospective analysis of 32 cases

OBJECTIVE

To evaluate the effect on the maturation of arteriovenous fistulas (AVFs) when using the internal branch of the cephalic vein compared with the main branch of the cephalic vein.

METHODS

The study enrolled patients with end-stage renal disease and divided them into an internal branch group (AVF constructed using the internal branch of the cephalic vein) or a main branch group (AVF constructed using the main branch of the cephalic vein). The surgical outcomes including complications were observed in these patients after 12 weeks.

RESULTS

Thirty-two patients with end-stage renal disease were included in the study. There were 16 patients in each group. The demographic and clinical characteristics were not significantly different between the two groups. The diameter of the arteries and veins were not significantly different between the two groups before the operation. In the internal branch group, significantly more (n = 7) patients failed to mature or required surgical intervention compared with the main branch group (n = 1).

CONCLUSION

For veins of the same diameter, these findings suggest that constructing AVFs using the main branch of the cephalic vein instead of the internal branch was more suitable for patients with end-stage renal disease requiring haemodialysis.

Improvement and validation of a computational model of flow in the swirling well cell culture model

Effects of fluid dynamics on cells are often studied by growing the cells on the base of cylindrical wells or dishes that are swirled on the horizontal platform of an orbital shaker. The swirling culture medium applies a shear stress to the cells that varies in magnitude and directionality from the center to the edge of the vessel. Computational fluid dynamics methods are used to simulate the flow and hence calculate shear stresses at the base of the well. The shear characteristics at each radial location are then compared with cell behavior at the same position. Previous simulations have generally ignored effects of surface tension and wetting, and results have only occasionally been experimentally validated. We investigated whether such idealized simulations are sufficiently accurate, examining a commonly-used swirling well configuration. The breaking wave predicted by earlier simulations was not seen, and the edge-to-center difference in shear magnitude (but not directionality) almost disappeared, when surface tension and wetting were included. Optical measurements of fluid height and velocity agreed well only with the computational model that incorporated surface tension and wetting. These results demonstrate the importance of including accurate fluid properties in computational models of the swirling well method.

Vascular dysfunction and pathology: focus on mechanical forces

The role of mechanical forces is emerging as a new player in the pathophysiologic programming of the cardiovascular system. The ability of the cells to 'sense' mechanical forces does not relate only to perception of movement or flow, as intended traditionally, but also to the biophysical properties of the extracellular matrix, the geometry of the tissues, and the force distribution inside them. This is also supported by the finding that cells can actively translate mechanical cues into discrete gene expression and epigenetic programming. In the present review, we will contextualize these new concepts in the vascular pathologic programming.

Endothelial cells do not align with the mean wall shear stress vector

The alignment of arterial endothelial cells (ECs) with the mean wall shear stress (WSS) vector is the prototypical example of their responsiveness to flow. However, evidence for this behaviour rests on experiments where many WSS metrics had the same orientation or where they were incompletely characterized. In the present study, we tested the phenomenon more rigorously. Aortic ECs were cultured in cylindrical wells on the platform of an orbital shaker. In this system, orientation would differ depending on the WSS metric to which the cells aligned. Variation in flow features and hence in possible orientations was further enhanced by altering the viscosity of the medium. Orientation of endothelial nuclei was compared with WSS characteristics obtained by computational fluid dynamics. At low mean WSS magnitudes, ECs aligned with the modal WSS vector, while at high mean WSS magnitudes they aligned so as to minimize the shear acting across their long axis (transverse WSS). Their failure to align with the mean WSS vector implies that other aspects of endothelial behaviour attributed to this metric require re-examination. The evolution of a mechanism for minimizing transverse WSS is consistent with it having detrimental effects on the cells and with its putative role in atherogenesis.

Long-term biofouling formation mediated by extracellular proteins in Nannochloropsis gaditana microalga cultures at different medium N/P ratios

Biofouling represents an important limitation in photobioreactor cultures. The biofouling propensity of different materials (polystyrene, borosilicate glass, polymethyl methacrylate, and polyethylene terephthalate glycol-modified) and coatings (two spray-applied and nanoparticle-based superhydrophobic coatings and a hydrogel-based fouling release coating) was evaluated by means of a short-term protein test, using bovine serum albumin (BSA) as a model protein, and by the long-term culture of the marine microalga Nannochloropsis gaditana under practical conditions. The results from both methods were similar, confirming that the BSA test predicts microalgal biofouling on surfaces exposed to microalgae cultures whose cells secrete macromolecules, such as proteins, with a high capacity for forming a conditioning film before cell adhesion. The hydrogel-based coating showed significantly reduced BSA and N. gaditana adhesion, whereas the other surfaces failed to control biofouling. Microalgal biofouling was associated with an increased concentration of sticky extracellular proteins at low N/P ratios (below 15).

Shake It or Shrink It: Mass Transport and Kinetics in Surface Bioassays Using Agitation and Microfluidics

Surface assays, such as ELISA and immunofluorescence, are nothing short of ubiquitous in biotechnology and medical diagnostics today. The development and optimization of these assays generally focuses on three aspects: immobilization chemistry, ligand-receptor interaction, and concentrations of ligands, buffers, and sample. A fourth aspect, the transport of the analyte to the surface, is more rarely delved into during assay design and analysis. Improving transport is generally limited to the agitation of reagents, a mode of flow generation inherently difficult to control, often resulting in inconsistent reaction kinetics. However, with assay optimization reaching theoretical limits, the role of transport becomes decisive. This perspective develops an intuitive and practical understanding of transport in conventional agitation systems and in microfluidics, the latter underpinning many new life science technologies. We give rules of thumb to guide the user on system behavior, such as advection regimes and shear stress, and derive estimates for relevant quantities that delimit assay parameters. Illustrative cases with examples of experimental results are used to clarify the role of fundamental concepts such as boundary and depletion layers, mass diffusivity, or surface tension.

Shear stress rosettes capture the complex flow physics in diseased arteries

Wall shear stress (WSS) is an important parameter in arterial mechanobiology. Various flow metrics, such as time averaged WSS (TAWSS), oscillatory shear index (OSI), and transWSS, have been used to characterize and relate possible WSS variations in arterial diseases like aneurysms and atherosclerosis. We use a graphical representation of WSS using shear rosettes to map temporal changes in the flow dynamics during a cardiac cycle at any spatial location on the vessel surface. The presence of secondary flows and flow reversals can be interpreted directly from the shape of the shear rosette. The mean WSS is given by the rosette centroid, the OSI by the splay around the rosette origin, and the transWSS by its width. We define a new metric, anisotropy ratio (AR), based on the ratio of the length to width of the shear rosette, to capture flow bi-directionality. We characterized the flow physics in controls and patient specific geometries of the ascending aorta (AA) and internal carotid artery (ICA) that have fundamentally different flow dynamics due to differences in the Reynolds and Womersley numbers. The differences in the flow dynamics are well reflected in the shapes of the WSS rosettes and the corresponding flow metrics.

Development of a two-stage limb ischemia model to better simulate human peripheral artery disease

Peripheral arterial disease (PAD) develops due to the narrowing or blockage of arteries supplying blood to the lower limbs. Surgical and endovascular interventions are the main treatments for advanced PAD but alternative and adjunctive medical therapies are needed. Currently the main preclinical experimental model employed in PAD research is based on induction of acute hind limb ischemia (HLI) by a 1-stage procedure. Since there are concerns regarding the ability to translate findings from this animal model to patients, we aimed to develop a novel clinically relevant animal model of PAD. HLI was induced in male Apolipoprotein E (ApoE-/-) deficient mice by a 2-stage procedure of initial gradual femoral artery occlusion by ameroid constrictors for 14 days and subsequent excision of the femoral artery. This 2-stage HLI model was compared to the classical 1-stage HLI model and sham controls. Ischemia severity was assessed using Laser Doppler Perfusion Imaging (LDPI). Ambulatory ability was assessed using an open field test, a treadmill test and using established scoring scales. Molecular markers of angiogenesis and shear stress were assessed within gastrocnemius muscle tissue samples using quantitative polymerase chain reaction. HLI was more severe in mice receiving the 2-stage compared to the 1-stage ischemia induction procedure as assessed by LDPI (p = 0.014), and reflected in a higher ischemic score (p = 0.004) and lower average distance travelled on a treadmill test (p = 0.045). Mice undergoing the 2-stage HLI also had lower expression of angiogenesis markers (vascular endothelial growth factor, p = 0.004; vascular endothelial growth factor- receptor 2, p = 0.008) and shear stress response mechano-transducer transient receptor potential vanilloid 4 (p = 0.041) within gastrocnemius muscle samples, compared to animals having the 1-stage HLI procedure. Mice subjected to the 2-stage HLI receiving an exercise program showed significantly greater improvement in their ambulatory ability on a treadmill test than a sedentary control group. This study describes a novel model of HLI which leads to more severe and sustained ischemia than the conventionally used model. Exercise therapy, which has established efficacy in PAD patients, was also effective in this new model. This new model maybe useful in the evaluation of potential novel PAD therapies.

Design and Computational Validation of a Novel Bioreactor for Conditioning Vascular Tissue to Time-Varying Multidirectional Fluid Shear Stress

PURPOSE

The cardiovascular endothelium experiences pulsatile and multidirectional fluid wall shear stress (WSS). While the effects of non-physiologic WSS magnitude and pulsatility on cardiovascular function have been studied extensively, the impact of directional abnormalities remains unknown due to the challenge to replicate this characteristic in vitro. To address this gap, this study aimed at designing a bioreactor capable of subjecting cardiovascular tissue to time-varying WSS magnitude and directionality.

METHODS

The device consisted of a modified cone-and-plate bioreactor. The cone rotation generates a fluid flow subjecting tissue to desired WSS magnitude, while WSS directionality is achieved by altering the alignment of the tissue relative to the flow at each instant of time. Computational fluid dynamics was used to verify the device ability to replicate the native WSS of the proximal aorta. Cone and tissue mount velocities were determined using an iterative optimization procedure.

RESULTS

Using conditions derived from cone-and-plate theory, the initial simulations yielded root-mean-square errors of 22.8 and 8.4% in WSS magnitude and angle, respectively, between the predicted and the target signals over one cycle, relative to the time-averaged target values. The conditions obtained after two optimization iterations reduced those errors to 3.5 and 0.5%, respectively, and generated 0.2% and 0.01% difference in time-averaged WSS magnitude and angle, respectively, relative to the target waveforms.

CONCLUSIONS

A bioreactor capable of generating simultaneously desired time-varying WSS magnitude and directionality was designed and validated computationally. The ability to subject tissue to in vivo-like WSS will provide new insights into cardiovascular mechanobiology and disease.

Understanding mechanobiology in cultured endothelium: A review of the orbital shaker method

A striking feature of atherosclerosis is its highly non-uniform distribution within the arterial tree. This has been attributed to variation in the haemodynamic wall shear stress (WSS) experienced by endothelial cells, but the WSS characteristics that are important and the mechanisms by which they lead to disease remain subjects of intensive investigation despite decades of research. In vivo evidence suggests that multidirectional WSS is highly atherogenic. This possibility is increasingly being studied by culturing endothelial cells in wells that are swirled on an orbital shaker. The method is simple and cost effective, has high throughput and permits chronic exposure, but interpretation of the results can be difficult because the fluid mechanics are complex; hitherto, their description has largely been restricted to the engineering literature. Here we review the findings of such studies, which indicate that putatively atherogenic flow characteristics occur at the centre of the well whilst atheroprotective ones occur towards the edge, and we describe simple mathematical methods for choosing experimental variables that avoid resonance, wave breaking and uncovering of the cells. We additionally summarise a large number of studies showing that endothelium cultured at the centre of the well expresses more pro-inflammatory and fewer homeostatic genes, has higher permeability, proliferation, apoptosis and senescence, and shows more endothelial-to-mesenchymal transition than endothelium at the edge. This simple method, when correctly interpreted, has the potential to greatly increase our understanding of the homeostatic and pathogenic mechanobiology of endothelial cells and may help identify new therapeutic targets in vascular disease.

Hydrodynamic shear stress promotes epithelial-mesenchymal transition by downregulating ERK and GSK3β activities

BACKGROUND

Epithelial-mesenchymal transition (EMT) occurs in the tumor microenvironment and presents an important mechanism of tumor cell intravasation, stemness acquisition, and metastasis. During metastasis, tumor cells enter the circulation to gain access to distant tissues, but how this fluid microenvironment influences cancer cell biology is poorly understood.

METHODS AND RESULTS

Here, we present both in vivo and in vitro evidence that EMT-like transition also occurs in circulating tumor cells (CTCs) as a result of hydrodynamic shear stress (+SS), which promotes conversion of CD24middle/CD44high/CD133middle/CXCR4low/ALDH1low primary patient epithelial tumor cells into specific high sphere-forming CD24low/CD44low/CD133high/CXCR4high/ALDH1high cancer stem-like cells (CSLCs) or tumor-initiating cells (TICs) with elevated tumor progression and metastasis capacity in vitro and in vivo. We demonstrate that conversion of CSLCs/TICs from epithelial tumor cells via +SS is dependent on reactive oxygen species (ROS)/nitric oxide (NO) generation, and suppression of extracellular signal-related kinase (ERK)/glycogen synthase kinase (GSK)3β, a mechanism similar to that operating in embryonic stem cells to prevent their differentiation while promoting self-renewal.

CONCLUSION

Fluid shear stress experienced during systemic circulation of human breast tumor cells can lead to specific acquisition of mesenchymal stem cell (MSC)-like potential that promotes EMT, mesenchymal-epithelial transition, and metastasis to distant organs. Our data revealed that biomechanical forces appeared to be important microenvironmental factors that not only drive hematopoietic development but also lead to acquisition of CSLCs/TIC potential in cancer metastasis. Our data highlight that +SS is a critical factor that promotes the conversion of CTCs into distinct TICs in blood circulation by endowing plasticity to these cells and by maintaining their self-renewal signaling pathways.

Vascular Endothelial Cell Behavior in Complex Mechanical Microenvironments

The vascular mechanical microenvironment consists of a mixture of spatially and temporally changing mechanical forces. This exposes vascular endothelial cells to both hemodynamic forces (fluid flow, cyclic stretching, lateral pressure) and vessel forces (basement membrane mechanical and topographical properties). The vascular mechanical microenvironment is "complex" because these forces are dynamic and interrelated. Endothelial cells sense these forces through mechanosensory structures and transduce them into functional responses via mechanotransduction pathways, culminating in behavior directly affecting vascular health. Recent in vitro studies have shown that endothelial cells respond in nuanced and unique ways to combinations of hemodynamic and vessel forces as compared to any single mechanical force. Understanding the interactive effects of the complex mechanical microenvironment on vascular endothelial behavior offers the opportunity to design future biomaterials and biomedical devices from the bottom-up by engineering for the cellular response. This review describes and defines (1) the blood vessel structure, (2) the complex mechanical microenvironment of the vascular endothelium, (3) the process in which vascular endothelial cells sense mechanical forces, and (4) the effect of mechanical forces on vascular endothelial cells with specific attention to recent works investigating the influence of combinations of mechanical forces. We conclude this review by providing our perspective on how the field can move forward to elucidate the effects of the complex mechanical microenvironment on vascular endothelial cell behavior.

Mechanotransmission in endothelial cells subjected to oscillatory and multi-directional shear flow

Local haemodynamics are linked to the non-uniform distribution of atherosclerosic lesions in arteries. Low and oscillatory (reversing in the axial flow direction) wall shear stress (WSS) induce inflammatory responses in endothelial cells (ECs) mediating disease localization. The objective of this study is to investigate computationally how the flow direction (reflected in WSS variation on the EC surface over time) influences the forces experienced by structural components of ECs that are believed to play important roles in mechanotransduction. A three-dimensional, multi-scale, multi-component, viscoelastic model of focally adhered ECs is developed, in which oscillatory WSS (reversing or non-reversing) parallel to the principal flow direction, or multi-directional oscillatory WSS with reversing axial and transverse components are applied over the EC surface. The computational model includes the glycocalyx layer, actin cortical layer, nucleus, cytoskeleton, focal adhesions (FAs), stress fibres and adherens junctions (ADJs). We show the distinct effects of atherogenic flow profiles (reversing unidirectional flow and reversing multi-directional flow) on subcellular structures relative to non-atherogenic flow (non-reversing flow). Reversing flow lowers stresses and strains due to viscoelastic effects, and multi-directional flow alters stress on the ADJs perpendicular to the axial flow direction. The simulations predict forces on integrins, ADJ filaments and other substructures in the range that activate mechanotransduction.

The Effect of Arterial Curvature on Blood Flow in Arterio-Venous Fistulae: Realistic Geometries and Pulsatile Flow

Arterio-Venous Fistulae (AVF) are regarded as the "gold standard" method of vascular access for patients with End-Stage Renal Disease (ESRD) who require haemodialysis. However, up to 60% of AVF do not mature, and hence fail, as a result of Intimal Hyperplasia (IH). Unphysiological flow and oxygen transport patterns, associated with the unnatural and often complex geometries of AVF, are believed to be implicated in the development of IH. Previous studies have investigated the effect of arterial curvature on blood flow in AVF using idealized planar AVF configurations and non-pulsatile inflow conditions. The present study takes an important step forwards by extending this work to more realistic non-planar brachiocephalic AVF configurations with pulsatile inflow conditions. Results show that forming an AVF by connecting a vein onto the outer curvature of an arterial bend does not, necessarily, suppress unsteady flow in the artery. This finding is converse to results from a previous more idealized study. However, results also show that forming an AVF by connecting a vein onto the inner curvature of an arterial bend can suppress exposure to regions of low wall shear stress and hypoxia in the artery. This finding is in agreement with results from a previous more idealized study. Finally, results show that forming an AVF by connecting a vein onto the inner curvature of an arterial bend can significantly reduce exposure to high WSS in the vein. The results are important, as they demonstrate that in realistic scenarios arterial curvature can be leveraged to reduce exposure to excessively low/high levels of WSS and regions of hypoxia in AVF. This may in turn reduce rates of IH and hence AVF failure.

Characterizations and Correlations of Wall Shear Stress in Aneurysmal Flow

Wall shear stress (WSS) is one of the most studied hemodynamic parameters, used in correlating blood flow to various diseases. The pulsatile nature of blood flow, along with the complex geometries of diseased arteries, produces complicated temporal and spatial WSS patterns. Moreover, WSS is a vector, which further complicates its quantification and interpretation. The goal of this study is to investigate WSS magnitude, angle, and vector changes in space and time in complex blood flow. Abdominal aortic aneurysm (AAA) was chosen as a setting to explore WSS quantification. Patient-specific computational fluid dynamics (CFD) simulations were performed in six AAAs. New WSS parameters are introduced, and the pointwise correlation among these, and more traditional WSS parameters, was explored. WSS magnitude had positive correlation with spatial/temporal gradients of WSS magnitude. This motivated the definition of relative WSS gradients. WSS vectorial gradients were highly correlated with magnitude gradients. A mix WSS spatial gradient and a mix WSS temporal gradient are proposed to equally account for variations in the WSS angle and magnitude in single measures. The important role that WSS plays in regulating near wall transport, and the high correlation among some of the WSS parameters motivates further attention in revisiting the traditional approaches used in WSS characterizations.

Impact of Bi-Axial Shear on Atherogenic Gene Expression by Endothelial Cells

This study demonstrated the effects of the directionality of oscillatory wall shear stress (WSS) on proliferation and proatherogenic gene expression (I-CAM, E-Selectin, and IL-6) in the presence of inflammatory mediators leukotriene B4 (LTB4) and bacterial lipopolysaccharide (LPS) from endothelial cells grown in an orbiting culture dish. Computational fluid dynamics (CFD) was applied to quantify the flow in the dish, while an analytical solution representing an extension of Stokes second problem was used for validation. Results indicated that WSS magnitude was relatively constant near the center of the dish and oscillated significantly (0-0.9 Pa) near the side walls. Experiments showed that LTB4 dominated the shear effects on cell proliferation and area. Addition of LPS didn't change proliferation, but significantly affected cell area. The expression of I-CAM1, E-Selectin and IL-6 were altered by directional oscillatory shear index (DOSI, a measure of the biaxiality of oscillatory shear), but not shear magnitude. The significance of DOSI was further reinforced by the strength of its interactions with other atherogenic factors. Hence, directionality of shear appears to be an important factor in regulating gene expression and provides a potential explanation of the propensity for increased vascular lesions in regions in the arteries with oscillating biaxial flow.

An orbital shear platform for real-time, in vitro endothelium characterization

Electrical impedance techniques have been used to characterize endothelium morphology, permeability, and motility in vitro. However, these impedance platforms have been limited to either static endothelium studies and/or induced laminar fluid flow at a constant, single shear stress value. In this work, we present a microfabricated impedance sensor for real-time, in vitro characterization of human umbilical vein endothelial cells (HUVECs) undergoing oscillatory hydrodynamic shear. Oscillatory shear was applied with an orbital shaker and the electrical impedance was measured by a microfabricated impedance chip with discrete electrodes positioned at radial locations of 0, 2.5, 5.0, 7.5, 10.0, and 12.5 mm from the center of the chip. Depending on their radial position within the circular orbital platform, HUVECs were exposed to shear values ranging between 0.6 and 6.71 dyne/cm(2) (according to numerical simulations) for 22 h. Impedance spectra were fit to an equivalent circuit model and the trans-endothelial resistance and monolayer's capacitance were extracted. Results demonstrated that, compared to measurements acquired before the onset of shear, cells at the center of the platform that experienced low steady shear stress (∼2.2 dyne/cm(2) ) had an average change in trans-endothelial resistance of 6.99 ± 4.06% and 1.78 ± 2.40% change in cell capacitance after 22 hours of shear exposure; cells near the periphery of the well (r = 12.5 mm) experienced transient shears (2.5-6.7 dyne/cm(2) ) and exhibited a greater change in trans-endothelial resistance (24.2 ± 10.8%) and cell capacitance (4.57 ± 5.39%). This study, demonstrates that the orbital shear platform provides a simple system that can capture and quantify the real-time cellular morphology as a result of induced shear stress. The orbital shear platform presented in this work, compared to traditional laminar platforms, subjects cells to more physiologically relevant oscillatory shear as well as exposes the sample to several shear values simultaneously. Biotechnol. Bioeng. 2016;113: 1336-1344. © 2015 Wiley Periodicals, Inc.

Device-Based In Vitro Techniques for Mechanical Stimulation of Vascular Cells: A Review

The most common cause of death in the developed world is cardiovascular disease. For decades, this has provided a powerful motivation to study the effects of mechanical forces on vascular cells in a controlled setting, since these cells have been implicated in the development of disease. Early efforts in the 1970 s included the first use of a parallel-plate flow system to apply shear stress to endothelial cells (ECs) and the development of uniaxial substrate stretching techniques (Krueger et al., 1971, “An in Vitro Study of Flow Response by Cells,” J. Biomech., 4(1), pp. 31–36 and Meikle et al., 1979, “Rabbit Cranial Sutures in Vitro: A New Experimental Model for Studying the Response of Fibrous Joints to Mechanical Stress,” Calcif. Tissue Int., 28(2), pp. 13–144). Since then, a multitude of in vitro devices have been designed and developed for mechanical stimulation of vascular cells and tissues in an effort to better understand their response to in vivo physiologic mechanical conditions. This article reviews the functional attributes of mechanical bioreactors developed in the 21st century, including their major advantages and disadvantages. Each of these systems has been categorized in terms of their primary loading modality: fluid shear stress (FSS), substrate distention, combined distention and fluid shear, or other applied forces. The goal of this article is to provide researchers with a survey of useful methodologies that can be adapted to studies in this area, and to clarify future possibilities for improved research methods.

The effect of oscillatory mechanical stimulation on osteoblast attachment and proliferation

The aim of this paper is to investigate the effect of the magnitude and duration of oscillatory mechanical stimulation on osteoblast attachment and proliferation as well as the time gap between seeding and applying the stimulation. Cells were exposed to three levels of speed at two different conditions. For the first group, mechanical shear stress was applied after 20 min of cell seeding. For the second group there was no time gap between cell seeding and applying mechanical stimulation. The total area subjected to shear stress was divided into three parts and for each part a comparative study was conducted at defined time points. Our results showed that both shear stress magnitude and the time gap between cell seeding and applying shear stress, are important in further cell proliferation and attachment. The effect of shear stress was not significant at lower speeds for both groups at earlier time points. However, a higher percentage of area was covered by cells at later time points under shear stress. In addition, the time gap can also improve osteoblast attachment. For the best rate of cell attachment and proliferation, the magnitude of shear stress and time gap should be optimized. The results of this paper can be utilized to improve cell attachment and proliferation in bioreactors.

The effect of in-plane arterial curvature on blood flow and oxygen transport in arterio-venous fistulae

Arterio-Venous Fistulae (AVF) are the preferred method of vascular access for patients with end stage renal disease who need hemodialysis. In this study, simulations of blood flow and oxygen transport were undertaken in various idealized AVF configurations. The objective of the study was to understand how arterial curvature affects blood flow and oxygen transport patterns within AVF, with a focus on how curvature alters metrics known to correlate with vascular pathology such as Intimal Hyperplasia (IH). If one subscribes to the hypothesis that unsteady flow causes IH within AVF, then the results suggest that in order to avoid IH, AVF should be formed via a vein graft onto the outer-curvature of a curved artery. However, if one subscribes to the hypothesis that low wall shear stress and/or low lumen-to-wall oxygen flux (leading to wall hypoxia) cause IH within AVF, then the results suggest that in order to avoid IH, AVF should be formed via a vein graft onto a straight artery, or the inner-curvature of a curved artery. We note that the recommendations are incompatible-highlighting the importance of ascertaining the exact mechanisms underlying development of IH in AVF. Nonetheless, the results clearly illustrate the important role played by arterial curvature in determining AVF hemodynamics, which to our knowledge has been overlooked in all previous studies.

Disturbed flow promotes endothelial senescence via a p53-dependent pathway

OBJECTIVE

Although atherosclerosis is associated with systemic risk factors such as age, high cholesterol, and obesity, plaque formation occurs predominately at branches and bends that are exposed to disturbed patterns of blood flow. The molecular mechanisms that link disturbed flow-generated mechanical forces with arterial injury are uncertain. To illuminate them, we investigated the effects of flow on endothelial cell (EC) senescence.

APPROACH AND RESULTS

LDLR(-/-) (low-density lipoprotein receptor(-/-)) mice were exposed to a high-fat diet for 2 to 12 weeks (or to a normal chow diet as a control) before the assessment of cellular senescence in aortic ECs. En face staining revealed that senescence-associated β-galactosidase activity and p53 expression were elevated in ECs at sites of disturbed flow in response to a high-fat diet. By contrast, ECs exposed to undisturbed flow did not express senescence-associated β-galactosidase or p53. Studies of aortae from healthy pigs (aged 6 months) also revealed enhanced senescence-associated β-galactosidase staining at sites of disturbed flow. These data suggest that senescent ECs accumulate at disturbed flow sites during atherogenesis. We used in vitro flow systems to examine whether a causal relationship exists between flow and EC senescence. Exposure of cultured ECs to flow (using either an orbital shaker or a syringe-pump flow bioreactor) revealed that disturbed flow promoted EC senescence compared with static conditions, whereas undisturbed flow reduced senescence. Gene silencing studies demonstrated that disturbed flow induced EC senescence via a p53-p21 signaling pathway. Disturbed flow-induced senescent ECs exhibited reduced migration compared with nonsenescent ECs in a scratch wound closure assay, and thus may be defective for arterial repair. However, pharmacological activation of sirtuin 1 (using resveratrol or SRT1720) protected ECs from disturbed flow-induced senescence.

CONCLUSIONS

Disturbed flow promotes endothelial senescence via a p53-p21-dependent pathway which can be inhibited by activation of sirtuin 1. These observations support the principle that pharmacological activation of sirtuin 1 may promote cardiovascular health by suppressing EC senescence at atheroprone sites.

Synergistic effects of orbital shear stress on in vitro growth and osteogenic differentiation of human alveolar bone-derived mesenchymal stem cells

Cellular behavior is dependent on a variety of physical cues required for normal tissue function. In order to mimic native tissue environments, human alveolar bone-derived mesenchymal stem cells (hABMSCs) were exposed to orbital shear stress (OSS) in a low-speed orbital shaker. The synergistic effects of OSS on proliferation and differentiation of hABMSCs were investigated. In particular, we induced the osteoblastic differentiation of hABMSCs cultured in the absence of OM by exposing hABMSCs to OSS (0.86-1.51 dyne/cm(2)). Activation of Cx43 was associated with exposure of hABMSCs to OSS. The viability of cells stimulated for 10, 30, 60, 120, and 180 min/day increased by approximately 10% compared with that of control. The OSS groups with stimulation of 10, 30, and 60 min/day had more intense mineralized nodules compared with the control group. In quantification of vascular endothelial growth factor (VEGF) and bone morphogenetic protein-2 (BMP-2) protein, VEGF protein levels under stimulation for 10, 60, and 180 min/day and BMP-2 levels under stimulation for 60, 120, and 180 min/day were significantly different compared with those of the control. In conclusion, the results indicated that exposing hABMSCs to OSS enhanced their differentiation and maturation.

Flow-mediated stem cell labeling with superparamagnetic iron oxide nanoparticle clusters

This study presents a strategy to enhance the uptake of superparamagnetic iron oxide nanoparticle (SPIO) clusters by manipulating the cellular mechanical environment. Specifically, stem cells exposed to an orbital flow ingested almost a 2-fold greater amount of SPIO clusters than those cultured statically. Improvements in magnetic resonance (MR) contrast were subsequently achieved for labeled cells in collagen gels and a mouse model. Overall, this strategy will serve to improve the efficiency of cell tracking and therapies.

Computation in the rabbit aorta of a new metric - the transverse wall shear stress - to quantify the multidirectional character of disturbed blood flow

Spatial variation of the haemodynamic stresses acting on the arterial wall is commonly assumed to explain the focal development of atherosclerosis. Disturbed flow in particular is thought to play a key role. However, widely-used metrics developed to quantify its extent are unable to distinguish between uniaxial and multidirectional flows. We analysed pulsatile flow fields obtained in idealised and anatomically-realistic arterial geometries using computational fluid dynamics techniques, and in particular investigated the multidirectionality of the flow fields, capturing this aspect of near-wall blood flow with a new metric - the transverse wall shear stress (transWSS) - calculated as the time-average of wall shear stress components perpendicular to the mean flow direction. In the idealised branching geometry, multidirectional flow was observed downstream of the branch ostium, a region of flow stagnation, and to the sides of the ostium. The distribution of the transWSS was different from the pattern of traditional haemodynamic metrics and more dependent on the velocity waveform imposed at the branch outlet. In rabbit aortas, transWSS patterns were again different from patterns of traditional metrics. The near-branch pattern varied between intercostal ostia, as is the case for lesion distribution; for some branches there were striking resemblances to the age-dependent patterns of disease seen in rabbit and human aortas. The new metric may lead to improved understanding of atherogenesis.

A material-independent cell-environment niche based on microreciprocating motion for cell growth enhancement

In tissue engineering, cell-cell, cell-scaffold and cell-environment communication balances regulate how cell populations participate in tissue generation, maintenance and repair. These communication balances are called niches. In this study, an easily implemented and material-independent cell-environment niche based on microreciprocating motions is developed to enhance cell growth. A micropositioning piezoelectric lead zirconate titanate stage is used to provide precise microreciprocating shear stress motions. Various shear stresses were applied to bovine endothelial cells (BECs) that were cultured on the artificially synthesized materials to obtain the suitable shear stress for growth enhancement. It was found that the suitable shear stress for apparent enhancement of BEC growth ranges from 1.8 to 2.2 N m(-2). Biopolymers were further used to verify the feasibility of the proposed approach using the optimized shear stress obtained from the culture on artificially synthesized polymers. The results further confirmed that the growth of BECs was enhanced as expected under the calculated reciprocating frequencies based on the suitable shear stress. It is hoped that the proposed microshear-stress-based niche could be a more cost- and time-effective solution for the enhancement of cell growth in tissue engineering applications.

Does low and oscillatory wall shear stress correlate spatially with early atherosclerosis? A systematic review

Low and oscillatory wall shear stress is widely assumed to play a key role in the initiation and development of atherosclerosis. Indeed, some studies have relied on the low shear theory when developing diagnostic and treatment strategies for cardiovascular disease. We wished to ascertain if this consensus is justified by published data. We performed a systematic review of papers that compare the localization of atherosclerotic lesions with the distribution of haemodynamic indicators calculated using computational fluid dynamics. The review showed that although many articles claim their results conform to the theory, it has been interpreted in different ways: a range of metrics has been used to characterize the distribution of disease, and they have been compared with a range of haemodynamic factors. Several studies, including all of those making systematic point-by-point comparisons of shear and disease, failed to find the expected relation. The various pre- and post-processing techniques used by different groups have reduced the range of shears over which correlations were sought, and in some cases are mutually incompatible. Finally, only a subset of the known patterns of disease has been investigated. The evidence for the low/oscillatory shear theory is less robust than commonly assumed. Longitudinal studies starting from the healthy state, or the collection of average flow metrics derived from large numbers of healthy vessels, both in conjunction with point-by-point comparisons using appropriate statistical techniques, will be necessary to improve our understanding of the relation between blood flow and atherogenesis.

Endothelial cell sensing of flow direction

OBJECTIVE

Atherosclerosis-prone regions of arteries are characterized by complex flow patterns where the magnitude of shear stress is low and direction rapidly changes, termed disturbed flow. How endothelial cells sense flow direction and how it impacts inflammatory effects of disturbed flow are unknown. We therefore aimed to understand how endothelial cells respond to changes in flow direction.

APPROACH AND RESULTS

Using a recently developed flow system capable of changing flow direction to any angle, we show that responses of aligned endothelial cells are determined by flow direction relative to their morphological and cytoskeletal axis. Activation of the atheroprotective endothelial nitric oxide synthase pathway is maximal at 180° and undetectable at 90°, whereas activation of proinflammatory nuclear factor-κB is maximal at 90° and undetectable at 180°. Similar effects were observed in randomly oriented cells in naive monolayers subjected to onset of shear. Cells aligned on micropatterned substrates subjected to oscillatory flow were also examined. In this system, parallel flow preferentially activated endothelial nitric oxide synthase and production of nitric oxide, whereas perpendicular flow preferentially activated reactive oxygen production and nuclear factor-κB.

CONCLUSIONS

These data show that the angle between flow and the cell axis defined by their shape and cytoskeleton determines endothelial cell responses. The data also strongly suggest that the inability of cells to align in low and oscillatory flow is a key determinant of the resultant inflammatory activation.