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.

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.

Mechanotransduction in Endothelial Cells in Vicinity of Cancer Cells.

Introduction-Local hemodynamics impact the mechanotransduction in endothelial cells (ECs) lining the vascular network. On the other hand, cancer cells are shown to influence the local hemodynamics in their vicinity, in microvasculature. The first objective of present study is to explore how cancer cell-induced changes in local hemodynamics can impact the forces experienced by intra/inter-cellular organelles of ECs that are believed to play important roles in mechanotransduction. Moreover, extracellular matrix (ECM) stiffening has been shown to correlate with progression of most cancer types. However, it is still not well understood how ECM stiffness impacts ECs mechanosensors. The second objective of this study is to elucidate the role of ECM stiffness on mechanotransduction in ECs. Methods-A three-dimensional, multiscale, multicomponent, viscoelastic model of focally adhered ECs is developed to simulate the force transmission through ECs mechanosensors [actin cortical layer, nucleus, cytoskeleton, focal adhesions (FAs), and adherens junctions (ADJs)]. Results-Our results show that cancer cell-altered hemodynamics results in significantly high forces transmitted to subcellular organelles of ECs which are in vicinity of cancer cells. This impact is more drastic on stress fibers (SFs) both centrally located and peripheral ones. Furthermore, we demonstrate that ADJs, FAs, and SFs experience higher stresses in ECs attached to stiffer ECM. Impact of ECM stiffness is particularly significant in ECs exposed to fluid shear stresses of 2 Pa or lower. This finding reveals the role of organ-specific stiffness in promoting cancer cell transmigration even in capillaries larger than cancer cell diameter. Conclusions-ÊCancer cell-induced-changes in ECs mechanotransduction represents an important potential mechanism for cancer cell transmigration in the microvasculature particularly with stiffer ECM. The identification of ECs mechanosensors involved in early stages of EC-cancer cell interaction will help with developing more efficient therapeutic interventions to suppress cancer cell transmigration in the microvasculature.

Association of Hypertension and Organ-Specific Cancer: A Meta-Analysis.

Hypertension and cancer are two of the leading global causes of death. Hypertension, known as chronic high blood pressure, affects approximately 45% of the American population and is a growing condition in other parts of the world, particularly in Asia and Europe. On the other hand, cancer resulted in approximately 10 million deaths in 2020 worldwide. Several studies indicate a coexistence of these two conditions, specifically that hypertension, independently, is associated with an increased risk of cancer. In the present study, we conducted a meta-analysis initially to reveal the prevalence of hypertension and cancer comorbidity and then to assess which organ-specific cancers were associated with hypertension by calculating the summary relative risks (RRs) and 95% confidence intervals (CIs). Our analysis shows that hypertension plays a role in cancer initiation. Our extended analysis on how the hypertension-associated angiogenesis factors are linked to cancer demonstrated that matrix metalloproteinases 2 and 9 appear to be two key factors facilitating cancer in hypertensive patients. This work serves as an important step in the current assessment of hypertension-promoted increased risk of 19 different cancers, particularly kidney, renal cell carcinoma, breast, colorectal, endometrial, and bladder. These findings provide new insight into how to treat and prevent cancer in hypertensive patients.

Localization of Rolling and Firm-Adhesive Interactions Between Circulating Tumor Cells and the Microvasculature Wall.

INTRODUCTION: The adhesion of tumor cells to vessel wall is a critical stage in cancer metastasis. Firm adhesion of cancer cells is usually followed by their extravasation through the endothelium. Despite previous studies identifying the influential parameters in the adhesive behavior of the cancer cell to a planer substrate, less is known about the interactions between the cancer cell and microvasculature wall and whether these interactions exhibit organ specificity. The objective of our study is to characterize sizes of microvasculature where a deformable circulating cell (DCC) would firmly adhere or roll over the wall, as well as to identify parameters that facilitate such firm adherence and underlying mechanisms driving adhesive interactions. METHODS: A three-dimensional model of DCCs is applied to simulate the fluid-structure interaction between the DCC and surrounding fluid. A dynamic adhesion model, where an adhesion molecule is modeled as a spring, is employed to represent the stochastic receptor-ligand interactions using kinetic rate expressions. RESULTS: Our results reveal that both the cell deformability and low shear rate of flow promote the firm adhesion of DCC in small vessels ( < 10 µ m ). Our findings suggest that ligand-receptor bonds of PSGL-1-P-selectin may lead to firm adherence of DCC in smaller vessels and rolling-adhesion of DCC in larger ones where cell velocity drops to facilitate the activation of integrin-ICAM-1 bonds. CONCLUSIONS: Our study provides a framework to predict accurately where different DCC-types are likely to adhere firmly in microvasculature and to establish the criteria predisposing cancer cells to such firm adhesion.

Balloon-Mounted Stents for Treatment of Refractory Flow Diverting Device Wall Malapposition.

BACKGROUND: As indications for flow diversion (FD) have expanded, new challenges in deployment of flow diverting devices (FDDs) have been encountered. We present 4 cases with aneurysms in which deployment of FDDs were complicated by device malapposition and compromised opening in regions of parent vessel stenosis. In all 4 cases, a balloon-mounted stent was ultimately deployed within the FDD. OBJECTIVE: To describe the use of balloon-mounted stents (BMS) within FDDs for correction of flow-limiting stenosis and device malapposition. METHODS: Patients undergoing FD for treatment of aneurysms complicated by refractory flow-limiting stenosis were identified through multi-center retrospective review. Those cases requiring use of BMS were identified. Further investigation in one of the cases was performed with a simulated pulsatile blood flow model. RESULTS: After attempts to perform balloon angioplasty proved unsuccessful, BMS deployment successfully opened the stenotic parent artery and improved FDD wall apposition in all 4 cases. Simulated pulsatile blood flow modeling confirmed improvements in the distribution of velocity, wall shear stress, oscillatory shear index, and flow pattern structure after stent deployment. One case was complicated by asymptomatic in-stent thrombosis. CONCLUSION: In cases of FDD deployment complicated by flow-limiting stenosis refractory to conventional techniques, a BMS deployed within the FD can provide radial support to open both the stenotic device and parent artery. Resulting improvements in device wall apposition may portend greater long-term efficacy of FD. In-stent occlusion can occur and may reflect a thrombogenic interaction between the devices.

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.

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.

Tissue prolapse and stresses in stented coronary arteries: A computer model for multi-layer atherosclerotic plaque.

Among the many factors influencing the effectiveness of cardiovascular stents, tissue prolapse indicates the potential of a stent to cause restenosis. The deflection of the arterial wall between the struts of the stent and the tissue is known as a prolapse or draping. The prolapse is associated with injury and damage to the vessel wall due to the high stresses generated around the stent when it expands. The current study investigates the impact of stenosis severity and plaque morphology on prolapse in stented coronary arteries. A finite element method is applied for the stent, plaque, and artery set to quantify the tissue prolapse and the corresponding stresses in stenosed coronary arteries. The variable size of atherosclerotic plaques is considered. A plaque is modelled as a multi-layered medium with different thicknesses attached to the single layer of an arterial wall. The results reveal that the tissue prolapse is influenced by the degree of stenosis severity and the thickness of the plaque layers. Stresses are observed to be significantly different between the plaque layers and the arterial wall tissue. Higher stresses are concentrated in fibrosis layer of the plaque (the harder core), while lower stresses are observed in necrotic core (the softer core) and the arterial wall layer. Moreover, the morphology of the plaque regulates the magnitude and distribution of the stress. The fibrous cap between the necrotic core and the endothelium constitutes the most influential layer to alter the stresses. In addition, the thickness of the necrotic core and the stenosis severity affect the stresses. This study reveals that the morphology of atherosclerotic plaques needs to be considered a key parameter in designing coronary stents.

Effects of severity and location of stenosis on the hemodynamics in human aorta and its branches.

Pulsatile blood flow is studied in a three-dimensional model of human thoracic aorta at different stages of atherosclerotic lesion growth, taking into account the effect of atherosclerotic plaque location and peripheral symmetry. The model is reconstructed from the computed tomography images. The wall shear stress (WSS), time-averaged WSS, and the oscillatory shear index are applied to determine susceptible sites for the onset of early atherosclerosis. Then, two different degrees of stenosis severity, 50 and 80 %, are introduced to vulnerable areas of the healthy aorta geometry. The overriding issue addressed is that the WSS distribution and magnitude are strongly affected by the atherosclerotic plaque size, its symmetric features, and the location, i.e., the branch it is formed. The present study, for the first time, is capable of providing information on the high shear environment that may exist upon the rupture of plaque surface and any thrombosis due to platelet deposition. The magnitude of WSS and its distribution at the throat of 50 % stenosed aortic arch are in agreement with the previous numerical study (Huang et al. in Exp Fluids 48(3):497-508, 2010). Results show that WSS values exceed 50 Pa at the throat of 80 % stenosed left common carotid and brachiocephalic arteries.

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

Haemodynamic forces applied at the apical surface of vascular endothelial cells (ECs) provide the mechanical signals at intracellular organelles and through the inter-connected cellular network. The objective of this study is to quantify the intracellular and intercellular stresses in a confluent vascular EC monolayer. A novel three-dimensional, multiscale and multicomponent model of focally adhered ECs is developed to account for the role of potential mechanosensors (glycocalyx layer, actin cortical layer, nucleus, cytoskeleton, focal adhesions (FAs) and adherens junctions (ADJs)) in mechanotransmission and EC deformation. The overriding issue addressed is the stress amplification in these regions, which may play a role in subcellular localization of mechanotransmission. The model predicts that the stresses are amplified 250-600-fold over apical values at ADJs and 175-200-fold at FAs for ECs exposed to a mean shear stress of 10 dyne cm(-2). Estimates of forces per molecule in the cell attachment points to the external cellular matrix and cell-cell adhesion points are of the order of 8 pN at FAs and as high as 3 pN at ADJs, suggesting that direct force-induced mechanotransmission by single molecules is possible in both. The maximum deformation of an EC in the monolayer is calculated as 400 nm in response to a mean shear stress of 1 Pa applied over the EC surface which is in accord with measurements. The model also predicts that the magnitude of the cell-cell junction inclination angle is independent of the cytoskeleton and glycocalyx. The inclination angle of the cell-cell junction is calculated to be 6.6° in an EC monolayer, which is somewhat below the measured value (9.9°) reported previously for ECs subjected to 1.6 Pa shear stress for 30 min. The present model is able, for the first time, to cross the boundaries between different length scales in order to provide a global view of potential locations of mechanotransmission.

Finite element modelling of pulsatile blood flow in idealized model of human aortic arch: study of hypotension and hypertension.

A three-dimensional computer model of human aortic arch with three branches is reproduced to study the pulsatile blood flow with Finite Element Method. In specific, the focus is on variation of wall shear stress, which plays an important role in the localization and development of atherosclerotic plaques. Pulsatile pressure pulse is used as boundary condition to avoid flow entry development, and the aorta walls are considered rigid. The aorta model along with boundary conditions is altered to study the effect of hypotension and hypertension. The results illustrated low and fluctuating shear stress at outer and inner wall of aortic arch, proximal wall of branches, and entry region. Despite the simplification of aorta model, rigid walls and other assumptions results displayed that hypertension causes lowered local wall shear stresses. It is the sign of an increased risk of atherosclerosis. The assessment of hemodynamics shows that under the flow regimes of hypotension and hypertension, the risk of atherosclerosis localization in human aorta may increase.

The transport of LDL across the deformable arterial wall: the effect of endothelial cell turnover and intimal deformation under hypertension.

A multilayered model of the aortic wall is introduced to investigate the transport of low-density lipoprotein (LDL) under hypertension, taking into account the influences of increased endothelial cell turnover and deformation of the intima at higher pressure. Meanwhile, the thickness and properties of the endothelium, intima, internal elastic lamina (IEL), and media are affected by the transmural pressure. The LDL macromolecules enter the intima through leaky junctions over the endothelium, which are created by dying or dividing cells. Water molecules enter the intima via the paracellular pathway through breaks in tight junctions after passing the glycocalyx as well as through leaky junctions. The glycocalyx is modeled as a Brinkman porous medium to describe the fluid filtration associated with its structure. Combined Navier-Stokes and Brinkman equations are solved for the transmural flow, and the convective-diffusion equation is employed for LDL transport. The permeation of LDL over the surface of smooth muscle cells is modeled through a uniform reaction evenly distributed in the macroscopically homogeneous media layer. Simulations are performed in an axisymmetric plane centered at a leaky cell. The overriding issue addressed is that LDL fluxes across the leaky junction, the intima, fenestral pores in the IEL, and the media layer are highly affected by the transmural pressure, which affects the endothelial cell turnover rate and the compaction of intima. The present model, for the first time and with no adjustable parameters, is capable of making many realistic predictions including the proper magnitudes for the permeability of endothelium and intimal layers and the hydraulic conductivity of all layers as well as their trends with pressure. Results for the volume flux through the wall and the hydraulic conductivity of the entire arterial wall, the endothelium, and subendothelial layers at 70 and 180 mmHg are in good agreement with previous experimental studies.

The study of wall deformation and flow distribution with transmural pressure by three-dimensional model of thoracic aorta wall.

The sensitivity of shear stress over smooth muscle cells (SMCs) to the deformability of media layer due to pressure is investigated in thoracic aorta wall using three-dimensional simulations. A biphasic, anisotropic model assuming the radius, thickness, and hydraulic conductivity of vessel wall as functions of transmural pressure is employed in numerical simulations. The leakage of interstitial fluid from intima to media layer is only possible through fenestral pores on the internal elastic lamina (IEL). The media layer is assumed a heterogeneous medium containing SMCs embedded in a porous extracellular matrix of elastin, proteoglycan, and collagen fibers. The applicable pressures for the deformation of media layer are varied from 0 to 180 mmHg. The SMCs are cylindrical objects of circular cross section at zero pressure. The cross sectional shape of SMCs changes from circle to ellipse as the media is compressed. The local shear stress over the nearest SMC to the IEL profoundly depends on pressure, SMCs configurations, and the corresponding distance to the IEL. The consideration of various SMC configurations, namely the staggered and square arrays, mimics various physiological conditions that can happen in positioning of an SMC. The results of our simulations show that even the second nearest SMCs to the IEL can significantly change their functions due to high shear stress levels. This is in contrast to earlier studies suggesting the highest vulnerability to shear stress for the innermost layer of SMCs at the intimal-medial interface.

Distribution of shear stress over smooth muscle cells in deformable arterial wall.

A biphasic, anisotropic model of the deformable aortic wall in combination with computational fluid dynamics is used to investigate the variation of shear stress over smooth muscle cells (SMCs) with transmural pressure. The media layer is modeled as a porous medium consisting of SMCs and a homogeneous porous medium of interstitial fluid and elastin, collagen and proteoglycans fibers. Interstitial fluid enters the media through fenestral pores, which are distributed over the internal elastic lamina (IEL). The IEL is considered as an impermeable barrier to fluid flow except at fenestral pores. The thickness and the radius of aortic wall vary with transmural pressure ranging from 10 to 180 mm Hg. It is assumed that SMCs are cylinders with a circular cross section at 0 mm Hg. As the transmural pressure increases, SMCs elongate with simultaneous change of cross sectional shape into ellipse according to the strain field in the media. Results demonstrate that the variation of shear stress within the media layer is significantly dependent on the configuration and cross sectional shape of SMCs. In the staggered array of SMCs, the shear stress over the first SMC nearest to the IEL is about 2.2 times lower than that of the square array. The shear stress even over the second nearest SMC to the IEL is considerably higher (about 15%) in the staggered array. In addition to configuration and cross sectional shape of SMCs, the variation of structural properties of the media layer with pressure and the sensitivity of the local shear stress to the minimum distance between SMCs and the IEL (reducing with transmural pressure) between SMCs and the IEL are studied. At 180 mm Hg, the ratio of the local shear stress of the nearest SMC to that of the second nearest SMC is 4.8 in the square array, whereas it reduces to about 1.8 in the staggered array. The importance of the fluid shear stress is associated with its role in the biomolecular state of smooth muscle cells bearing the shear stress.

Effect of the shape and configuration of smooth muscle cells on the diffusion of ATP through the arterial wall.

In this study, the shape and the configuration of smooth muscle cells (SMCs) within the arterial wall are altered to investigate their influence on molecular transport across the media layer of the thoracic aorta wall. In a 2D geometry of the media layer containing SMCs, the finite-element method has been employed to simulate the diffusion of solutes through the media layer. The media is modeled as a heterogeneous system composed of SMCs having elliptic or circular cross sections embedded in a homogeneous porous medium made of proteoglycan and collagen fibers with an interstitial fluid filling the void. The arrangement of SMCs is in either ordered or disordered fashion for different volume fractions of SMCs. The interstitial fluid enters the media through fenestral pores, which are assumed to be distributed uniformly over the internal elastic lamina (IEL). Results revealed that in an ordered arrangement of SMCs, the concentration of adenosine 5'-triphosphate (ATP) over the surface of SMCs with an elliptic cross section is 5-8% more than those of circular SMCs in volume fractions of 0.4-0.7. The ATP concentration at the SMC surface decreases with volume fraction in the ordered configuration of SMCs. In a disordered configuration, the local ATP concentration at the SMC surface and in the bulk are strongly dependent on the distance between neighboring SMCs, as well as the minimum distance between SMCs and fenestral pores. Moreover, the SMCs in farther distances from the IEL are as important as those just beneath the IEL in disordered configurations. The results of this study lead us to better understanding of the role of SMCs in controlling the diffusion of important species such as ATP within the arterial wall.