Shear-induced force transmission in a multicomponent, multicell model of the endothelium

Computational analysis of cancer cell adhesion in curved vessels affected by wall shear stress for prediction of metastatic spreading

Introduction: The dynamics of circulating tumor cells (CTCs) within blood vessels play a pivotal role in predicting metastatic spreading of cancer within the body. However, the limited understanding and method to quantitatively investigate the influence of vascular architecture on CTC dynamics hinders our ability to predict metastatic process effectively. To address this limitation, the present study was conducted to investigate the influence of blood vessel tortuosity on the behaviour of CTCs, focusing specifically on establishing methods and examining the role of shear stress in CTC-vessel wall interactions and its subsequent impact on metastasis. Methods: We computationally simulated CTC behaviour under various shear stress conditions induced by vessel tortuosity. Our computational model, based on the lattice Boltzmann method (LBM) and a coarse-grained spectrin-link membrane model, efficiently simulates blood plasma dynamics and CTC deformability. The model incorporates fluid-structure interactions and receptor-ligand interactions crucial for CTC adhesion using the immersed boundary method (IBM). Results: Our findings reveal that uniform shear stress in straight vessels leads to predictable CTC-vessel interactions, whereas in curved vessels, asymmetrical flow patterns and altered shear stress create distinct adhesion dynamics, potentially influencing CTC extravasation. Quantitative analysis shows a 25% decrease in the wall shear stress in low-shear regions and a 58.5% increase in the high-shear region. We observed high-shear regions in curved vessels to be potential sites for increased CTC adhesion and extravasation, facilitated by elevated endothelial expression of adhesion molecules. This phenomenon correlates with the increased number of adhesion bonds, which rises to approximately 40 in high-shear regions, compared to around 12 for straight vessels and approximately 5-6 in low-shear regions. The findings also indicate an optimal cellular stiffness necessary for successful CTC extravasation in curved vessels. Discussion: By the quantitative assessment of the risk of CTC extravasation as a function of vessel tortuosity, our study offers a novel tool for the prediction of metastasis risk to support the development of personalized therapeutic interventions based on individual vascular characteristics and tumor cell properties.

Modeling Physical Forces Experienced by Cancer and Stromal Cells Within Different Organ-Specific Tumor Tissue

Mechanical force exerted on cancer cells by their microenvironment have been reported to drive cells toward invasive phenotypes by altering cells' motility, proliferation, and apoptosis. These mechanical forces include compressive, tensile, hydrostatic, and shear forces. The importance of forces is then hypothesized to be an alteration of cancer cells' and their microenvironment's biophysical properties as the indicator of a tumor's malignancy state. Our objective is to investigate and quantify the correlation between a tumor's malignancy state and forces experienced by the cancer cells and components of the microenvironment. In this study, we have developed a multicomponent, three-dimensional model of tumor tissue consisting of a cancer cell surrounded by fibroblasts and extracellular matrix (ECM). Our results on three different organs including breast, kidney, and pancreas show that: A) the stresses within tumor tissue are impacted by the organ specific ECM's biophysical properties, B) more invasive cancer cells experience higher stresses, C) in pancreas which has a softer ECM (Young modulus of 1.0 kPa) and stiffer cancer cells (Young modulus of 2.4 kPa and 1.7 kPa) than breast and kidney, cancer cells experienced significantly higher stresses, D) cancer cells in contact with ECM experienced higher stresses compared to cells surrounded by fibroblasts but the area of tumor stroma experiencing high stresses has a maximum length of 40 μm when the cancer cell is surrounded by fibroblasts and 12 μm for when the cancer cell is in vicinity of ECM. This study serves as an important first step in understanding of how the stresses experienced by cancer cells, fibroblasts, and ECM are associated with malignancy states of cancer cells in different organs. The quantification of forces exerted on cancer cells by different organ-specific ECM and at different stages of malignancy will help, first to develop theranostic strategies, second to predict accurately which tumors will become highly malignant, and third to establish accurate criteria controlling the progression of cancer cells malignancy. Furthermore, our in silico model of tumor tissue can yield critical, useful information for guiding ex vivo or in vitro experiments, narrowing down variables to be investigated, understanding what factors could be impacting cancer treatments or even biomarkers to be looking for.

Basal endothelial glycocalyx's response to shear stress: a review of structure, function, and clinical implications

The endothelial glycocalyx encompasses the entire endothelial cell, transducing extracellular signals and regulating vascular permeability and barrier functions. The apical glycocalyx, which forms the lumen of the vessel, and the basal glycocalyx, at the smooth muscle cell interface, are often investigated separately as they are exposed to vastly different stimuli. The apical glycocalyx directly senses fluid shear forces transmitting them intracellularly through connection to the cytoskeleton of the endothelial cell. The basal glycocalyx has demonstrated sensitivity to shear due to blood flow transmitted through the cytoskeleton, promoting alternate signaling processes. In this review, we discuss current literature on the basal glycocalyx's response to shear stress in the context of mechanotransduction and remodeling. The possible implications of basal glycocalyx degradation in pathologies are also explored. Finally, this review seeks to highlight how addressing the gaps discussed would improve our wholistic understanding of the endothelial glycocalyx and its role in maintaining vascular homeostasis.

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.

Impact of Wound Dressing on Mechanotransduction within Tissues of Chronic Wounds

Chronic wounds are significant public health problems impacting the health-related quality of individuals' lives (due to disability, decreased productivity, and loss of independence) and an immense economic burden to healthcare systems around the world. In this study, our main objective is to investigate how mechanotransduction can impact the healing process in chronic wounds. We have developed new three-dimensional models of wound tissue to study the distribution of forces within these tissues exerted by wound dressings with different characteristics. The roles of mechanical forces on wound healing have gained significant clinical attention; the application of mechanical forces is expected to influence the physiology of tissue surrounding a wound. We aim to investigate whether the force transmission within wound tissue is impacted by the dressing characteristics and whether this impact may differ with wound tissue's properties. Our results show that wound dressings with lower stiffnesses promote force transmission within a wound tissue. This impact is even more significant on stiffer wound tissues. Furthermore, we show that size of wound dressing alters forces that transmit within the wound tissue where dressings with 9 cm length show higher stresses. The wound tissue stiffening has been associated with healing of a wound. Our results demonstrate that wounds with stiffer tissue experience higher stresses. Taken all together, our findings suggest that low stiffness of wound dressing and its size may be introduced as a criterion to explain parameters predisposing a chronic wound to heal. This study's findings on the role of dressings and tissue characteristics demonstrate that precision dressings are required for wound management and understanding how a dressing impacts mechanotransduction in wound tissue will lead to design of novel dressings promoting healing in chronic wounds.

Biomechanics of cells and subcellular components: A comprehensive review of computational models and applications

Cells are a fundamental structural, functional and biological unit for all living organisms. Up till now, considerable efforts have been made to study the responses of single cells and subcellular components to an external load, and understand the biophysics underlying cell rheology, mechanotransduction and cell functions using experimental and in silico approaches. In the last decade, computational simulation has become increasingly attractive due to its critical role in interpreting experimental data, analysing complex cellular/subcellular structures, facilitating diagnostic designs and therapeutic techniques, and developing biomimetic materials. Despite the significant progress, developing comprehensive and accurate models of living cells remains a grand challenge in the 21st century. To understand current state of the art, this review summarises and classifies the vast array of computational biomechanical models for cells. The article covers the cellular components at multi-spatial levels, that is, protein polymers, subcellular components, whole cells and the systems with scale beyond a cell. In addition to the comprehensive review of the topic, this article also provides new insights into the future prospects of developing integrated, active and high-fidelity cell models that are multiscale, multi-physics and multi-disciplinary in nature. This review will be beneficial for the researchers in modelling the biomechanics of subcellular components, cells and multiple cell systems and understanding the cell functions and biological processes from the perspective of cell mechanics.

Hydrodynamic instability and flow reduction in polymer brush coated channels

A polymer brush is a passive medium. At equilibrium the knowledge of its chemical composition and thickness is enough for a full system characterization. However, when the brush is exposed to fluid flow it reveals a much more intriguing nature, in which filamentous protrusions and the way they interact among themselves and with the surrounding fluid are of outmost importance. Here we investigate such a rich behavior via numerical simulations. We focus on the brush hydrodynamic response at low Reynolds numbers, observing a significant fluid flow reduction inside a polymer-brush coated channel. We find that the reduction of the flow inside the channel is significantly larger than what would happen if the brush effect consisted only in reducing the effective channel width. This amplified reduction is understood as being due to the morphological instability of the brush-liquid interface which is shown to have an elastic origin: the mechanical stress acting on the brush due to the imposed flow is partially released by the interface modulation. In turn, this modulation dissipates more energy than a flat interface in the surrounding fluid, causing a reduction of flow velocity. Our results and interpretations provide an explanation for recent experimental measurements.

Understanding the Role of Endothelial Glycocalyx in Mechanotransduction via Computational Simulation: A Mini Review

Endothelial glycocalyx (EG) is a forest-like structure, covering the lumen side of blood vessel walls. EG is exposed to the mechanical forces of blood flow, mainly shear, and closely associated with vascular regulation, health, diseases, and therapies. One hallmark function of the EG is mechanotransduction, which means the EG senses the mechanical signals from the blood flow and then transmits the signals into the cells. Using numerical modelling methods or in silico experiments to investigate EG-related topics has gained increasing momentum in recent years, thanks to tremendous progress in supercomputing. Numerical modelling and simulation allows certain very specific or even extreme conditions to be fulfilled, which provides new insights and complements experimental observations. This mini review examines the application of numerical methods in EG-related studies, focusing on how computer simulation contributes to the understanding of EG as a mechanotransducer. The numerical methods covered in this review include macroscopic (i.e., continuum-based), mesoscopic [e.g., lattice Boltzmann method (LBM) and dissipative particle dynamics (DPD)] and microscopic [e.g., molecular dynamics (MD) and Monte Carlo (MC) methods]. Accounting for the emerging trends in artificial intelligence and the advent of exascale computing, the future of numerical simulation for EG-related problems is also contemplated.

Early Atherosclerotic Changes in Coronary Arteries are Associated with Endothelium Shear Stress Contraction/Expansion Variability

Although unphysiological wall shear stress (WSS) has become the consensus hemodynamic mechanism for coronary atherosclerosis, the complex biomechanical stimulus affecting atherosclerosis evolution is still undetermined. This has motivated the interest on the contraction/expansion action exerted by WSS on the endothelium, obtained through the WSS topological skeleton analysis. This study tests the ability of this WSS feature, alone or combined with WSS magnitude, to predict coronary wall thickness (WT) longitudinal changes. Nine coronary arteries of hypercholesterolemic minipigs underwent imaging with local WT measurement at three time points: baseline (T1), after 5.6 ± 0.9 (T2), and 7.6 ± 2.5 (T3) months. Individualized computational hemodynamic simulations were performed at T1 and T2. The variability of the WSS contraction/expansion action along the cardiac cycle was quantified using the WSS topological shear variation index (TSVI). Alone or combined, high TSVI and low WSS significantly co-localized with high WT at the same time points and were significant predictors of thickening at later time points. TSVI and WSS magnitude values in a physiological range appeared to play an atheroprotective role. Both the variability of the WSS contraction/expansion action and WSS magnitude, accounting for different hemodynamic effects on the endothelium, (1) are linked to WT changes and (2) concur to identify WSS features leading to coronary atherosclerosis.

Hemodynamic and morphological characteristics of a growing cerebral aneurysm

The growth of cerebral aneurysms is linked to local hemodynamic conditions, but the driving mechanisms of the growth are poorly understood. The goal of this study was to examine the association between intraaneurysmal hemodynamic features and areas of aneurysm growth, to present the key hemodynamic parameters essential for an accurate prediction of the growth, and to gain a deeper understanding of the underlying mechanisms. Patient-specific images of a growing cerebral aneurysm in 3 different growth stages acquired over a period of 40 months were segmented and reconstructed. A unique aspect of this patient-specific case study was that while one side of the aneurysm stayed stable, the other side continued to grow. This unique case enabled the authors to examine their aims in the same patient with parent and daughter arteries under the same inlet flow conditions. Pulsatile flow in the aneurysm models was simulated using computational fluid dynamics and was validated with in vitro experiments using particle image velocimetry measurements. The authors' detailed analysis of intrasaccular hemodynamics linked the growing regions of aneurysms to flow instabilities and complex vortex structures. Extremely low velocities were observed at or around the center of the unstable vortex structure, which matched well with the growing regions of the studied cerebral aneurysm. Furthermore, the authors observed that the aneurysm wall regions with a growth greater than 0.5 mm coincided with wall regions of lower (< 0.5 Pa) time-averaged wall shear stress (TAWSS), lower instantaneous (< 0.5 Pa) wall shear stress (WSS), and high (> 0.1) oscillatory shear index (OSI). To determine which set of parameters can best identify growing and nongrowing aneurysms, the authors performed statistical analysis for consecutive stages of the growing CA. The results demonstrated that the combination of TAWSS and the distance from the center of the vortical structure has the highest sensitivity and positive predictive value, and relatively high specificity and negative predictive value. These findings suggest that an unstable, recirculating flow structure within the aneurysm sac created in the region adjacent to the aneurysm wall with low TAWSS may be introduced as an accurate criterion to explain the hemodynamic conditions predisposing the aneurysm to growth. The authors' findings are based on one patient's data set, but the study lays out the justification for future large-scale verification. The authors' findings can assist clinicians in differentiating stable and growing aneurysms during preinterventional planning.

Fluid Shear Stress Sensing by the Endothelial Layer

Blood flow produces mechanical frictional forces, parallel to the blood flow exerted on the endothelial wall of the vessel, the so-called wall shear stress (WSS). WSS sensing is associated with several vascular pathologies, but it is first a physiological phenomenon. Endothelial cell sensitivity to WSS is involved in several developmental and physiological vascular processes such as angiogenesis and vascular morphogenesis, vascular remodeling, and vascular tone. Local conditions of blood flow determine the characteristics of WSS, i.e., intensity, direction, pulsatility, sensed by the endothelial cells that, through their effect of the vascular network, impact WSS. All these processes generate a local-global retroactive loop that determines the ability of the vascular system to ensure the perfusion of the tissues. In order to account for the physiological role of WSS, the so-called shear stress set point theory has been proposed, according to which WSS sensing acts locally on vessel remodeling so that WSS is maintained close to a set point value, with local and distant effects of vascular blood flow. The aim of this article is (1) to review the existing literature on WSS sensing involvement on the behavior of endothelial cells and its short-term (vasoreactivity) and long-term (vascular morphogenesis and remodeling) effects on vascular functioning in physiological condition; (2) to present the various hypotheses about WSS sensors and analyze the conceptual background of these representations, in particular the concept of tensional prestress or biotensegrity; and (3) to analyze the relevance, explanatory value, and limitations of the WSS set point theory, that should be viewed as dynamical, and not algorithmic, processes, acting in a self-organized way. We conclude that this dynamic set point theory and the biotensegrity concept provide a relevant explanatory framework to analyze the physiological mechanisms of WSS sensing and their possible shift toward pathological situations.

Mechanobiology of dynamic enzyme systems

This Perspective paper advances a hypothesis of mechanosensation by endothelial cells in which the cell is a dynamic crowded system, driven by continuous enzyme activity, that can be shifted from one non-equilibrium state to another by external force. The nature of the shift will depend on the direction, rate of change, and magnitude of the force. Whether force induces a pathophysiological or physiological change in cell biology will be determined by whether the dynamics of a cellular system can accommodate the dynamics and magnitude of the force application. The complex interplay of non-static cytoskeletal structures governs internal cellular rheology, dynamic spatial reorganization, and chemical kinetics of proteins such as integrins, and a flaccid membrane that is dynamically supported; each may constitute the necessary dynamic properties able to sense external fluid shear stress and reorganize in two and three dimensions. The resulting reorganization of enzyme systems in the cell membrane and cytoplasm may drive the cell to a new physiological state. This review focuses on endothelial cell mechanotransduction of shear stress, but may lead to new avenues of investigation of mechanobiology in general requiring new tools for interrogation of mechanobiological systems, tools that will enable the synthesis of large amounts of spatial and temporal data at the molecular, cellular, and system levels.

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.

Role of deformable cancer cells on wall shear stress-associated-VEGF secretion by endothelium in microvasculature

Endothelial surface layer (glycocalyx) is the major physiological regulator of tumor cell adhesion to endothelium. Cancer cells express vascular endothelial growth factor (VEGF) which increases the permeability of a microvessel wall by degrading glycocalyx. Endothelial cells lining large arteries have also been reported, in vitro and in vivo, to mediate VEGF expression significantly under exposure to high wall shear stress (WSS) > 0.6 Pa. The objective of the present study is to explore whether local hemodynamic conditions in the vicinity of a migrating deformable cancer cell can influence the function of endothelial cells to express VEGF within the microvasculature. A three-dimensional model of deformable cancer cells (DCCs) migrating within a capillary is developed by applying a massively parallel hemodynamics application to simulate the fluid-structure interaction between the DCC and fluid surrounding the DCC. We study how dynamic interactions between the DCC and its local microenvironment affect WSS exposed on endothelium, under physiological conditions of capillaries with different diameters and flow conditions. Moreover, we quantify the area of endothelium affected by the DCC. Our results show that the DCC alters local hemodynamics in its vicinity up to an area as large as 40 times the cancer cell lateral surface. In this area, endothelium experiences high WSS values in the range of 0.6–12 Pa. Endothelial cells exposed to this range of WSS have been reported to express VEGF. Furthermore, we demonstrate that stiffer cancer cells expose higher WSS on the endothelium. A strong impact of cell stiffness on its local microenvironment is observed in capillaries with diameters <16 μm. WSS-induced-VEGF by endothelium represents an important potential mechanism for cancer cell adhesion and metastasis in the microvasculature. This work serves as an important first step in understanding the mechanisms driving VEGF-expression by endothelium and identifying the underlying mechanisms of glycocalyx degradation by endothelium in microvasculature. The identification of angiogenesis factors involved in early stages of cancer cell-endothelium interactions and understanding their regulation will help, first to develop anti-angiogenic strategies applied to diagnostic studies and therapeutic interventions, second to predict accurately where different cancer cell types most likely adhere in microvasculature, and third to establish accurate criteria predisposing the cancer metastasis.

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.

Fluid shear, intercellular stress, and endothelial cell alignment

Endothelial cell alignment along the direction of laminar fluid flow is widely understood to be a defining morphological feature of vascular homeostasis. While the role of associated signaling and structural events have been well studied, associated intercellular stresses under laminar fluid shear have remained ill-defined and the role of these stresses in the alignment process has remained obscure. To fill this gap, we report here the tractions as well as the complete in-plane intercellular stress fields measured within the human umbilical vein endothelial cell (HUVEC) monolayer subjected to a steady laminar fluid shear of 1 Pa. Tractions, intercellular stresses, as well as their time course, heterogeneity, and anisotropy, were measured using monolayer traction microscopy and monolayer stress microscopy. Prior to application of laminar fluid flow, intercellular stresses were largely tensile but fluctuated dramatically in space and in time (317 ± 122 Pa). Within 12 h of the onset of laminar fluid flow, the intercellular stresses decreased substantially but continued to fluctuate dramatically (142 ± 84 Pa). Moreover, tractions and intercellular stresses aligned strongly and promptly (within 1 h) along the direction of fluid flow, whereas the endothelial cell body aligned less strongly and substantially more slowly (12 h). Taken together, these results reveal that steady laminar fluid flow induces prompt reduction in magnitude and alignment of tractions and intercellular stress tensor components followed by the retarded elongation and alignment of the endothelial cell body. Appreciably smaller intercellular stresses supported by cell-cell junctions logically favor smaller incidence of gap formation and thus improved barrier integrity.

Computational models of the primary cilium and endothelial mechanotransmission

In endothelial cells (ECs), the mechanotransduction of fluid shear stress is partially dependent on the transmission of force from the fluid into the cell (mechanotransmission). The role of the primary cilium in EC mechanotransmission is not yet known. To motivate a framework towards quantifying cilia contribution to EC mechanotransmission, we have reviewed mechanical models of both (1) the primary cilium (three-dimensional and lower-dimensional) and (2) whole ECs (finite element, non-finite element, and tensegrity). Both the primary cilia and whole EC models typically incorporate fluid-induced wall shear stress and spatial geometry based on experimentally acquired images of cells. This paper presents future modelling directions as well as the major goals towards integrating primary cilium models into a multi-component EC mechanical model. Finally, we outline how an integrated cilium-EC model can be used to better understand mechanotransduction in the endothelium.