Finite-element stress analysis of a multicomponent model of sheared and focally-adhered endothelial cells

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.

A cell-and-plasma numerical model reveals hemodynamic stress and flow adaptation in zebrafish microvessels after morphological alteration

The development of a functional cardiovascular system ensures a sustainable oxygen, nutrient and hormone delivery system for successful embryonic development and homeostasis in adulthood. While early vessels are formed by biochemical signaling and genetic programming, the onset of blood flow provides mechanical cues that participate in vascular remodeling of the embryonic vascular system. The zebrafish is a prolific animal model for studying the quantitative relationship between blood flow and vascular morphogenesis due to a combination of favorable factors including blood flow visualization in optically transparent larvae. In this study, we have developed a cell-and-plasma blood transport model using computational fluid dynamics (CFD) to understand how red blood cell (RBC) partitioning affect lumen wall shear stress (WSS) and blood pressure in zebrafish trunk blood vascular networks with altered rheology and morphology. By performing live imaging of embryos with reduced hematocrit, we discovered that cardiac output and caudal artery flow rates were maintained. These adaptation trends were recapitulated in our CFD models, which showed reduction in network WSS via viscosity reduction in the caudal artery/vein and via pressure gradient weakening in the intersegmental vessels (ISVs). Embryos with experimentally reduced lumen diameter showed reduced cardiac output and caudal artery flow rate. Factoring in this trend into our CFD models, simulations highlighted that lumen diameter reduction increased vessel WSS but this increase was mitigated by flow reduction due to the adaptive network pressure gradient weakening. Additionally, hypothetical network CFD models with different vessel lumen diameter distribution characteristics indicated the significance of axial variation in lumen diameter and cross-sectional shape for establishing physiological WSS gradients along ISVs. In summary, our work demonstrates how both experiment-driven and hypothetical CFD modeling can be employed for the study of blood flow physiology during vascular remodeling.

Disrupted Stiffness Ratio Alters Nuclear Mechanosensing

The ability of endothelial cells to sense and respond to dynamic changes in blood flow is critical for vascular homeostasis and cardiovascular health. The mechanical and geometric properties of the nuclear and cytoplasmic compartments affect mechanotransduction. We hypothesized that alterations to these parameters have resulting mechanosensory consequences. Using atomic force microscopy and mathematical modeling, we assessed how the nuclear and cytoplasmic compartment stiffnesses modulate shear stress transfer to the nucleus within aging endothelial cells. Our computational studies revealed that the critical parameter controlling shear transfer is not the individual mechanics of these compartments, but the stiffness ratio between them. Replicatively aged cells had a reduced stiffness ratio, attenuating shear transfer, while the ratio was not altered in a genetic model of accelerated aging. We provide a theoretical framework suggesting that dysregulation of the shear stress response can be uniquely imparted by relative mechanical changes in subcellular compartments.

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.

In silico stress fibre content affects peak strain in cytoplasm and nucleus but not in the membrane for uniaxial substrate stretch

Existing in silico models for single cell mechanics feature limited representations of cytoskeletal structures that contribute substantially to the mechanics of a cell. We propose a micromechanical hierarchical approach to capture the mechanical contribution of actin stress fibres. For a cell-specific fibroblast geometry with membrane, cytoplasm and nucleus, the Mori-Tanaka homogenization method was employed to describe cytoplasmic inhomogeneities and constitutive contribution of actin stress fibres. The homogenization was implemented in a finite element model of the fibroblast attached to a substrate through focal adhesions. Strain in cell membrane, cytoplasm and nucleus due to uniaxial substrate stretch was assessed for different stress fibre volume fractions and different elastic modulus of the substrate. A considerable decrease of the peak strain with increasing stress fibre content was observed in cytoplasm and nucleus but not the membrane, whereas the peak strain in cytoplasm, nucleus and membrane increased for increasing elastic modulus of the substrate. Finite element mesh of reconstructed human fibroblast and intracellular strain distribution in cell subjected to substrate stretch.

Image-derived modeling of nucleus strain amplification associated with chromatin heterogeneity

Beyond the critical role of cell nuclei in gene expression and DNA replication, they also have a significant influence on cell mechanosensation and migration. Nuclear stiffness can impact force transmission and, furthermore, act as a physical barrier to translocation across tight spaces. As such, it is of wide interest to accurately characterize nucleus mechanical behavior. In this study, we present a computational investigation of the in situ deformation of a heterogeneous chondrocyte nucleus. A methodology is developed to accurately reconstruct a three-dimensional finite-element model of a cell nucleus from confocal microscopy. By incorporating the reconstructed nucleus into a chondrocyte model embedded in pericellular and extracellular matrix, we explore the relationship between spatially heterogeneous nuclear DNA content, shear stiffness, and resultant shear strain. We simulate an externally applied extracellular matrix shear deformation and compute intranuclear strain distributions, which are directly compared with corresponding experimentally measured distributions. Simulations suggest that the mechanical behavior of the nucleus is highly heterogeneous, with a nonlinear relationship between experimentally measured grayscale values and corresponding local shear moduli (μn). Three distinct phases are identified within the nucleus: a low-stiffness mRNA-rich interchromatin phase (0.17 kPa ≤ μn ≤ 0.63 kPa), an intermediate-stiffness euchromatin phase (1.48 kPa ≤ μn ≤ 2.7 kPa), and a high-stiffness heterochromatin phase (3.58 kPa ≤ μn ≤ 4.0 kPa). Our simulations also indicate that disruption of the nuclear envelope associated with lamin A/C depletion significantly increases nuclear strain in regions of low DNA concentration. We further investigate a phenotypic shift of chondrocytes to fibroblast-like cells, a signature for osteoarthritic cartilage, by increasing the contractility of the actin cytoskeleton to a level associated with fibroblasts. Peak nucleus strains increase by 35% compared to control, with the nucleus becoming more ellipsoidal. Our findings may have broad implications for current understanding of how local DNA concentrations and associated strain amplification can impact cell mechanotransduction and drive cell behavior in development, migration, and tumorigenesis.

Assessing haemorrhage-critical values of cerebral blood flow by modelling biomechanical stresses on capillaries in the immature brain

The development of intraventricular haemorrhages (IVH) in preterm newborns is triggered by a disruption of the vessels responsible for cerebral microcirculation. Analysis of the stresses exerted on vessel walls enables the identification of the critical values of cerebral blood flow (CBF) associated with the development of IVH in preterm infants. The purpose of the present study is the estimation of these critical CBF values using the biomechanical stresses obtained by the finite element modelling of immature brain capillaries. The properties of the endothelial cells and basement membranes employed were selected on the basis of published nanoindentation measurements using atomic force microscopes. The forces acting on individual capillaries were derived with a mathematical model that accounts for the peculiarities of microvascularity in the immature brain. Calculations were based on clinical measurements obtained from 254 preterm infants with the gestational age ranging from 23 to 30 weeks, with and without diagnosis of IVH. No distinction between the affected and control groups with the gestational age of 23 to 26 weeks was possible. For infants with the gestational age of 27 to 30 weeks, the CBF value of 17.03 ml/100 g/min was determined as the critical upper value, above which the likelihood of IVH increases.

Good advice for endothelial cells: Get in line, relax tension, and go with the flow

Endothelial cells (ECs) are continuously subjected to fluid wall shear stress (WSS) and cyclic strain caused by pulsatile blood flow and pressure. It is well established that these hemodynamic forces each play important roles in vascular disease, but their combined effects are not well understood. ECs remodel in response to both WSS and cyclic strain to align along the vessel axis, but in areas prone to atherogenesis, such an alignment is absent. In this perspective, experimental and clinical findings will be reviewed, which have revealed the characteristics of WSS and cyclic strain, which are associated with atherosclerosis, spanning studies on whole blood vessels to individual cells to mechanosensing molecules. Examples are described regarding the use of computational modeling to elucidate the mechanisms by which EC alignment contributes to mechanical homeostasis. Finally, the need to move toward an integrated understanding of how hemodynamic forces influence EC mechanotransduction is presented, which holds the potential to move our currently fragmented understanding to a true appreciation of the role of mechanical stimuli in atherosclerosis.

Cell focal adhesion clustering leads to decreased and homogenized basal strains

The subendothelial matrix of the artery is a complex mechanical environment where endothelial cells respond to and affect changes upon the underlying substrate. Our recent work has demonstrated that endothelial cell strain heterogeneity increases on a more heterogeneous underlying subendothelial matrix, and these cells display increased focal adhesion presence on stiffer substrate areas. However, the impact of these grouped focal adhesions on endothelial cell strains has not been explored. Here, we use finite element modeling to investigate the effects of microscale stiffness heterogeneities and focal adhesion location and stiffness on endothelial cell strains. Shear stress applied to the apical cell layer demonstrated a minimal effect on cell strain values, while substrate stretch had a greater effect on cell strain in the cell-substrate model. The addition of focal adhesions into the computational model (cell-FA-substrate model) predicted a decrease and homogenization of the cell strains. For simulations including focal adhesions, stiffer and more distributed adhesions caused increased and more heterogeneous endothelial cell strains. Overall, our data indicate that cells may group focal adhesions to minimize and homogenize their basal strains.

In vitro measurements of hemodynamic forces and their effects on endothelial cell mechanics at the sub-cellular level

This paper presents micro-particle tracking velocimetry measurements over cultured bovine aortic endothelial cell monolayers in microchannels. The objective was to quantify fluid forces and cell morphology at the sub-cellular scale for monolayers subjected to steady shear rates of 5, 10, and 20 dyn/cm2. The ultimate goal of this study was to develop an experimental methodology for in vitro detailed study of physiologically realistic healthy and diseased conditions. Cell topography, shear stress, and pressure distributions were calculated from sets of velocity fields made in planes parallel to the microchannel wall. For each experiment, measurements were made in 3 h intervals for 18 h. It was found that there is a three-dimensional change in cell morphology as a result of applied shear stress. That is, cells flatten and become more wedge shaped in the stream direction while conserving volume by spreading laterally, i.e., in the cross-stream direction. These changes in cell morphology are directly related to local variations in fluid loading, i.e., shear stress and pressure. This paper describes the first flow measurements over a confluent layer of endothelial cells that are spatially resolved at the sub-cellular scale with a simultaneous temporal resolution to quantify the response of cells to fluid loading.

Endothelial cells on an aged subendothelial matrix display heterogeneous strain profiles in silico

Within the artery intima, endothelial cells respond to mechanical cues and changes in subendothelial matrix stiffness. Recently, we found that the aging subendothelial matrix stiffens heterogeneously and that stiffness heterogeneities are present on the scale of one cell length. However, the impacts of these complex mechanical micro-heterogeneities on endothelial cells have not been fully understood. Here, we simulate the effects of matrices that mimic young and aged vessels on single- and multi-cell endothelial cell models and examine the resulting cell basal strain profiles. Although there are limitations to the model which prohibit the prediction of intracellular strain distributions in alive cells, this model does introduce mechanical complexities to the subendothelial matrix material. More heterogeneous basal strain distributions are present in the single- and multi-cell models on the matrix mimicking an aged artery over those exhibited on the young artery. Overall, our data indicate that increased heterogeneous strain profiles in endothelial cells are displayed in silico when there is an increased presence of microscale arterial mechanical heterogeneities in the matrix.

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 Nuclear Option: Evidence Implicating the Cell Nucleus in Mechanotransduction

Biophysical stimuli presented to cells via microenvironmental properties (e.g., alignment and stiffness) or external forces have a significant impact on cell function and behavior. Recently, the cell nucleus has been identified as a mechanosensitive organelle that contributes to the perception and response to mechanical stimuli. However, the specific mechanotransduction mechanisms that mediate these effects have not been clearly established. Here, we offer a comprehensive review of the evidence supporting (and refuting) three hypothetical nuclear mechanotransduction mechanisms: physical reorganization of chromatin, signaling at the nuclear envelope, and altered cytoskeletal structure/tension due to nuclear remodeling. Our goal is to provide a reference detailing the progress that has been made and the areas that still require investigation regarding the role of nuclear mechanotransduction in cell biology. Additionally, we will briefly discuss the role that mathematical models of cell mechanics can play in testing these hypotheses and in elucidating how biophysical stimulation of the nucleus drives changes in cell behavior. While force-induced alterations in signaling pathways involving lamina-associated polypeptides (LAPs) (e.g., emerin and histone deacetylase 3 (HDAC3)) and transcription factors (TFs) located at the nuclear envelope currently appear to be the most clearly supported mechanism of nuclear mechanotransduction, additional work is required to examine this process in detail and to more fully test alternative mechanisms. The combination of sophisticated experimental techniques and advanced mathematical models is necessary to enhance our understanding of the role of the nucleus in the mechanotransduction processes driving numerous critical cell functions.

Effects of Frequency and Acceleration Amplitude on Osteoblast Mechanical Vibration Responses: A Finite Element Study

Bone cells are deformed according to mechanical stimulation they receive and their mechanical characteristics. However, how osteoblasts are affected by mechanical vibration frequency and acceleration amplitude remains unclear. By developing 3D osteoblast finite element (FE) models, this study investigated the effect of cell shapes on vibration characteristics and effect of acceleration (vibration intensity) on vibrational responses of cultured osteoblasts. Firstly, the developed FE models predicted natural frequencies of osteoblasts within 6.85-48.69 Hz. Then, three different levels of acceleration of base excitation were selected (0.5, 1, and 2 g) to simulate vibrational responses, and acceleration of base excitation was found to have no influence on natural frequencies of osteoblasts. However, vibration response values of displacement, stress, and strain increased with the increase of acceleration. Finally, stress and stress distributions of osteoblast models under 0.5 g acceleration in Z-direction were investigated further. It was revealed that resonance frequencies can be a monotonic function of cell height or bottom area when cell volumes and material properties were assumed as constants. These findings will be useful in understanding how forces are transferred and influence osteoblast mechanical responses during vibrations and in providing guidance for cell culture and external vibration loading in experimental and clinical osteogenesis studies.

Cell shape-dependent shear stress on adherent cells in a micro-physiologic system as revealed by FEM

Flow-induced shear stress on adherent cells leads to biochemical signaling and mechanical responses of the cells. To determine the flow-induced shear stress on adherent cells cultured in a micro-scaled reaction chamber, we developed a suitable finite element method model. The influence of the most important parameters-cell shape, cell density, shear modulus and fluid velocity-was investigated. Notably, the cell shape strongly influences the resulting shear stress. Long and smooth cells undergo lower shear stress than more rounded cells. Changes in the curvature of the cells lead to stress peaks and single cells experience higher shear stress values than cells of a confluent monolayer. The computational results of the fluid flow simulation were validated experimentally. We also analyzed the influence of flow-induced shear stress on the metabolic activity and shape of L929, a mouse fibroblast cell line, experimentally. The results indicate that threshold stress values for continuous flow conditions cannot be transferred to quasi static flow conditions interrupted by short fluid exchange events.

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.

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.

Dynamic monitoring of mechano-sensing of cells by gold nanoslit surface plasmon resonance sensor

We demonstrated a real-time monitoring of live cells upon laminar shear stress stimulation via surface plasmon resonance (SPR) in gold nanoslit array. A large-area gold nanostructure consisted of 500-nm-period nanoslits was fabricated on a plastic film using the thermal-annealed template-stripping method. The SPR in the gold nanoslit array provides high surface sensitivity to monitor cell adhesion changes near the sensor surface. The human non-small cell lung cancer (CL1-0), human lung fibroblast (MRC-5), and human dermal fibroblast (Hs68) were cultured on the gold nanoslits and their dynamic responses to laminar shear stress were measured under different stress magnitudes from 0 to 30 dyne/cm(2). Cell adhesion was increased in CL1-0 under shear flow stimulation. No adhesion recovery was observed after stopping the flow. On the other hand, MRC-5 and Hs68 decreased adhesion and recovered from the shear stress. The degree of recovery was around 70% for MRC-5. This device provides dynamic study and early detection of cell adhesion changes under shear flow conditions.

Surface projections of titanium substrates increase antithrombotic endothelial function in response to shear stress

Despite the therapeutic benefits of both mechanical circulatory assist devices and nitinol stents with titanium (Ti) outer surfaces, problems remain with thrombosis at the blood-contacting surface. Covering these surfaces with a layer of endothelium would mimic the native lining of the cardiovascular system, potentially decreasing thrombotic complications. Since surface topography is known to affect the phenotype of a seeded cell layer and since stents and ventricular assist devices exhibit surface protrusions, we tested the hypothesis that endothelial cells (ECs) have altered function on Ti surfaces with protrusions of 1.25, 3, and 5 μm height, compared with smooth Ti surfaces. ECs and nuclei were more aligned and ECs were more elongated on all patterned surfaces. Cell area was reduced on the 3 and 5 μm features. Expression of eNOS and COX2 was not altered by patterned surfaces, but expression of KLF-2 was higher on 1.25 and 5 μm features. Nitric oxide production following exposure to flow was higher on the 5 μm features. These results show that some antithrombogenic functions of ECs are significantly enhanced for ECs cultured on surface protrusions, and no functions are diminished, informing the future design of implant surfaces for endothelialization.

TRANSDUCTION OF STRAIN TO CELLS SEEDED ONTO SCAFFOLDS EXPOSED TO UNIAXIAL STRETCHING: A THREE DIMENSIONAL FINITE ELEMENT STUDY

When preparing tissue engineering and regenerative medicine constructs, a commonly encountered problem is the failure of seeded cells to infiltrate the scaffold. In an increasing number of cases, constructs are being mechanically preconditioned with the expectation that preconditioning will enhance the construct's maturation and effectiveness by pre-exposing seeded cells to stimuli the tissue of interest experiences in vivo. However, whether or not mechanostimulation of a scaffold actually results in transmission of stimuli to the seeded cells remains poorly understood. The purpose of this research was to develop a model that quantifies how strain is transmitted to cells layered on a scaffold's surface compared to cells embedded within a scaffold. Three-dimensional finite element models representative of these conditions were created. When 10% strain was applied to the construct, embedded cells received the full imposed strain. However, cells growing on top of the scaffold received 5% strain within the first layer of cells, and the strain transmitted to cells in subsequent layers decreased exponentially with increasing distance from the scaffold's surface. When experimentally testing the model, strain-induced biological responses were muted in conditions where cell to scaffold contact was reduced. This research illustrates the importance of achieving cellular penetration and cell-to-scaffold contacts when mechanically conditioning tissue engineering constructs.

Coordinated Mechanosensitivity of Membrane Rafts and Focal Adhesions

Endothelial cells sense mechanical forces of blood flow through mechanisms that involve focal adhesions (FAs). The mechanosensitive pathways that originate from FA-associated integrin activation may involve membrane rafts, small cholesterol- and sphigolipid-rich domains that are either immobilized, by virtue of their attachment to the cytoskeleton, or highly mobile in the plane of the plasma membrane. In this study, we fluorescently labeled non-mobile and mobile populations of GM1, a ganglioside associated with lipid rafts, and transfected cells with the red fluorescent protein-(RFP-) talin, an indicator of integrin activation at FAs, in order to determine the kinetics and sequential order of raft and talin mechanosensitivity. Cells were imaged under confocal microscopy during mechanical manipulation of a FA induced by a fibronectin (FN)-functionalized nanoelectrode with feedback control of position. First, FA deformation led to long range deformation of immobile rafts followed by active recoil of a subpopulation of displaced rafts. Second, initial adhesion between the FN-probe and the cell induced rapid accumulation of GM1 at the probe site with a time constant of 1.7 s. Talin accumulated approximately 20 s later with a time constant of 0.6 s. Third, a 1 μm deformation of the FA lead to immediate (0.3 s) increase in GM1 fluorescence and a later (6 s) increase in talin. Fourth, long term deformation of FAs led to continual GM1 accumulation at the probe site that was reversed upon removal of the deformation. These results demonstrate that rafts are directly mechanosensitive and that raft mobility may enable the earliest events related to FA mechanosensing and reinforcement upon force application.

A model of localised Rac1 activation in endothelial cells due to fluid flow

Endothelial cells respond to fluid flow by elongating in the direction of flow. Cytoskeletal changes and activation of signalling molecules have been extensively studied in this response, including: activation of receptors by mechano-transduction, actin filament alignment in the direction of flow, changes to cell-substratum adhesions, actin-driven lamellipodium extension, and localised activation of Rho GTPases. To study this process we model the force over a single cell and couple this to a model of the Rho GTPases, Rac and Rho, via a Kelvin-body model of mechano-transduction. It is demonstrated that a mechano-transducer can respond to the normal component of the force is likely to be a necessary component of the signalling network in order to establish polarity. Furthermore, the rate-limiting step of Rac1 activation is predicted to be conversion of Rac-GDP to Rac-GTP, rather than activation of upstream components. Modelling illustrates that the aligned endothelial cell morphology could attenuate the signalling network.

Theoretical models for coronary vascular biomechanics: progress & challenges

A key aim of the cardiac Physiome Project is to develop theoretical models to simulate the functional behaviour of the heart under physiological and pathophysiological conditions. Heart function is critically dependent on the delivery of an adequate blood supply to the myocardium via the coronary vasculature. Key to this critical function of the coronary vasculature is system dynamics that emerge via the interactions of the numerous constituent components at a range of spatial and temporal scales. Here, we focus on several components for which theoretical approaches can be applied, including vascular structure and mechanics, blood flow and mass transport, flow regulation, angiogenesis and vascular remodelling, and vascular cellular mechanics. For each component, we summarise the current state of the art in model development, and discuss areas requiring further research. We highlight the major challenges associated with integrating the component models to develop a computational tool that can ultimately be used to simulate the responses of the coronary vascular system to changing demands and to diseases and therapies.

The effect of noisy flow on endothelial cell mechanotransduction: a computational study

Flow in the arterial system is mostly laminar, but turbulence occurs in vivo under both normal and pathological conditions. Turbulent and laminar flow elicit significantly different responses in endothelial cells (ECs), but the mechanisms allowing ECs to distinguish between these different flow regimes remain unknown. The authors present a computational model that describes the effect of turbulence on mechanical force transmission within ECs. Because turbulent flow is inherently "noisy" with random fluctuations in pressure and velocity, our model focuses on the effect of signal noise (a stochastically changing force) on the deformation of intracellular transduction sites including the nucleus, cell-cell adhesion proteins (CCAPs), and focal adhesion sites (FAS). The authors represent these components of the mechanical signaling pathway as linear viscoelastic structures (Kelvin bodies) connected to the cell surface via cytoskeletal elements. The authors demonstrate that FAS are more sensitive to signal noise than the nucleus or CCAP. The relative sensitivity of these various structures to noise is affected by the nature of the cytoskeletal connections within the cell. Finally, changes in the compliance of the nucleus dramatically affect nuclear sensitivity to noise, suggesting that pathologies that alter nuclear mechanical properties will be associated with abnormal EC responsiveness to turbulent flow.

PAM2 (piston assisted microsyringe): a new rapid prototyping technique for biofabrication of cell incorporated scaffolds

Rapid prototyping techniques are widely used to fabricate well-defined three-dimensional structures of tissue homologs. The piston-assisted microsyringe (PAM2) is a rapid prototyping technology specifically developed for low-shear stress extrusion of viscous hydrogel solutions containing cells. In this article the working parameters of the system were established to guarantee the realization of spatially controlled hydrogel scaffolds. Moreover the shear stresses acting on the cell membrane during extrusion was investigated through a computational fluid-dynamic analysis. The computational models show that the shear stress on the cells is of the order of 100 Pa during the extrusion process. HepG2 cells encapsulated in alginate were then extruded into spatially organized hepatic lobule-like architectures and their viability and function were evaluated. The results show that the metabolic fingerprint of the cells is preserved with respect to controls and the cells are uniformly distributed through the gel scaffold.

Endothelial Cell Membrane Sensitivity to Shear Stress is Lipid Domain Dependent

Blood flow-associated shear stress causes physiological and pathophysiological biochemical processes in endothelial cells that may be initiated by alterations in plasma membrane lipid domains characterized as liquid-ordered (l(o)), such as rafts or caveolae, or liquid-disordered (l(d)). To test for domain-dependent shear sensitivity, we used time-correlated single photon counting instrumentation to assess the photophysics and dynamics of the domain-selective lipid analogues DiI-C(12) and DiI-C(18) in endothelial cells subjected to physiological fluid shear stress. Under static conditions, DiI-C(12) fluorescence lifetime was less than DiI-C(18) lifetime and the diffusion coefficient of DiI-C(12) was greater than the DiI-C(18) diffusion coefficient, confirming that DiI-C(12) probes l(d), a more fluid membrane environment, and DiI-C(18) probes the l(o) phase. Domains probed by DiI-C(12) exhibited an early (10 s) and transient decrease of fluorescence lifetime after the onset of shear while domains probed by DiI-C(18) exhibited a delayed decrease of fluorescence lifetime that was sustained for the 2 min the cells were subjected to flow. The diffusion coefficient of DiI-C(18) increased after shear imposition, while that of DiI-C(12) remained constant. Determination of the number of molecules (N) in the control volume suggested that DiI-C(12)-labeled domains increased in N immediately after step-shear, while N for DiI-C(18)-stained membrane transiently decreased. These results demonstrate that membrane microdomains are differentially sensitive to fluid shear stress.

Engineering high-density endothelial cell monolayers on soft substrates

This study demonstrates that a confluent monolayer of endothelial cells (ECs) can be tissue engineered on a soft substrate with a cell density and morphology that approximates in vivo conditions. We achieved formation of a confluent EC monolayer on polydimethylsiloxane (PDMS) elastomer by microcontact printing of fibronectin (FN) in a square lattice array of 3microm diameter circular islands at a 6microm pitch. Uniform coatings of FN or serum proteins on PDMS or on tissue-culture-treated polystyrene failed to support the equivalent EC density and/or confluence. The ECs on the FN micropatterned PDMS achieved a density of 1,536+/-247cellsmm(-2), close to the 3,215+/-336cellsmm(-2) observed in vivo from porcine pulmonary artery and significantly higher (2- to 5-fold) than EC density on other materials. The probable mechanism for enhanced EC adhesion, growth and density is increased focal adhesion (FA) formation between the ECs and the substrate. After 14days culture, the micropatterned FN surface increased the average number of FAs per cell to 35+/-10, compared to 7+/-6 for ECs on PDMS uniformly coated with FN. Thus, microscale patterning of FN into FA-sized, circular islands on PDMS elastomer promotes the formation of EC monolayers with in vivo-like cell density and morphology.

Image-based finite element modeling of alveolar epithelial cell injury during airway reopening

The acute respiratory distress syndrome (ARDS) is characterized by fluid accumulation in small pulmonary airways. The reopening of these fluid-filled airways involves the propagation of an air-liquid interface that exerts injurious hydrodynamic stresses on the epithelial cells (EpC) lining the airway walls. Previous experimental studies have demonstrated that these hydrodynamic stresses may cause rupture of the plasma membrane (i.e., cell necrosis) and have postulated that cell morphology plays a role in cell death. However, direct experimental measurement of stress and strain within the cell is intractable, and limited data are available on the mechanical response (i.e., deformation) of the epithelium during airway reopening. The goal of this study is to use image-based finite element models of cell deformation during airway reopening to investigate how cell morphology and mechanics influence the risk of cell injury/necrosis. Confocal microscopy images of EpC in subconfluent and confluent monolayers were used to generate morphologically accurate three-dimensional finite element models. Hydrodynamic stresses on the cells were calculated from boundary element solutions of bubble propagation in a fluid-filled parallel-plate flow channel. Results indicate that for equivalent cell mechanical properties and hydrodynamic load conditions, subconfluent cells develop higher membrane strains than confluent cells. Strain magnitudes were also found to decrease with increasing stiffness of the cell and membrane/cortex region but were most sensitive to changes in the cell's interior stiffness. These models may be useful in identifying pharmacological treatments that mitigate cell injury during airway reopening by altering specific biomechanical properties of the EpC.

Effects of particulates and lipids on the hydraulic conductivity of Matrigel

The hydraulic conductivity of a connective tissue is determined both by the fine ultrastructure of the extracellular matrix and the effects of larger particles in the interstitial space. In this study, we explored this relationship by examining the effects of 30- or 90-nm-diameter latex nanospheres or low-density lipoproteins (LDL) on the hydraulic conductivity of Matrigel, a basement membrane matrix. The hydraulic conductivity of Matrigel with latex nanospheres or LDL particles added at 4.8% weight fraction was measured and compared with the hydraulic conductivity of Matrigel alone. The LDL-derived lipids in the gel were visualized by transmission electron microscopy and were seen to have aggregated into particles up to 500 nm in size. The addition of these materials to the medium markedly decreased its hydraulic conductivity, with the LDL-derived lipids having a much larger effect than did the latex nanospheres. Debye-Brinkman theory was used to predict the effect of addition of particles to the hydraulic conductivity of the medium. The theoretical predictions matched well with the results from adding latex nanospheres to the medium. However, LDL decreased hydraulic conductivity much more than was predicted by the theory. The validation of the theoretical model for rigid particles embedded in extracellular matrix suggests that it could be used to make predictions about the influence of particulates (e.g., collagen, elastin, cells) on the hydraulic conductivity of the fine filamentous matrix (the proteoglycans) in connective tissues. In addition, the larger-than-predicted effects of lipidlike particles on hydraulic conductivity may magnify the pathology associated with lipid accumulation, such as in Bruch's membrane of the retina during macular degeneration and the blood vessel wall in atherosclerosis.

Bioengineering and Bioinformatics Summer Institutes: meeting modern challenges in undergraduate summer research

Summer undergraduate research programs in science and engineering facilitate research progress for faculty and provide a close-ended research experience for students, which can prepare them for careers in industry, medicine, and academia. However, ensuring these outcomes is a challenge when the students arrive ill-prepared for substantive research or if projects are ill-defined or impractical for a typical 10-wk summer. We describe how the new Bioengineering and Bioinformatics Summer Institutes (BBSI), developed in response to a call for proposals by the National Institutes of Health (NIH) and the National Science Foundation (NSF), provide an impetus for the enhancement of traditional undergraduate research experiences with intense didactic training in particular skills and technologies. Such didactic components provide highly focused and qualified students for summer research with the goal of ensuring increased student satisfaction with research and mentor satisfaction with student productivity. As an example, we focus on our experiences with the Penn State Biomaterials and Bionanotechnology Summer Institute (PSU-BBSI), which trains undergraduates in core technologies in surface characterization, computational modeling, cell biology, and fabrication to prepare them for student-centered research projects in the role of materials in guiding cell biology.

Chapter 5: Imaging in depth: controversies and opportunities

One enduring challenge of biological imaging is achieving depth of penetration-into cells, tissues, and animals. How deeply can we probe and with what resolution and efficacy? These are critical issues as microscopists seek to push ever deeper, while resolving structural details and observing specific molecular events. In this guide to depth-appropriate modalities, standard optical platforms such as confocal and two-photon microscopes are considered along with complementary imaging modalities that range in depth of penetration. After an introduction to basic techniques, the trade-offs and limitations that distinguish competing technologies are considered, with emphasis on the visualization of subcellular structures and dynamic events. Not surprisingly, there are differences of opinion regarding imaging technologies, as highlighted in a section on point-scanning and Nipkow-disk style confocal microscopes. Confocal microscopy is then contrasted with deconvolution and multi-photon imaging modalities. It is also important to consider the detectors used by current instruments (such as PMTs and CCD cameras). Ultimately specimen properties, in conjunction with instrumentation, determine the depth at which subcellular operations and larger-scale biological processes can be visualized. Relative advantages are mentioned in the context of experiment planning and instrument-purchase decisions. Given the rate at which new optical techniques are being invented, this report should be viewed as a snapshot of current capabilities, with the goal of providing a framework for thinking about new developments.

Macrorheology and adaptive microrheology of endothelial cells subjected to fluid shear stress

Vascular endothelial cells (ECs) respond to temporal and spatial characteristics of hemodynamic forces by alterations in their adhesiveness to leukocytes, secretion of vasodilators, and permeability to blood-borne constituents. These physiological and pathophysiological changes are tied to adaptation of cell mechanics and mechanotransduction, the process by which cells convert forces to intracellular biochemical signals. The exact time scales of these mechanical adaptations, however, remain unknown. We used particle-tracking microrheology to study adaptive changes in intracellular mechanics in response to a step change in fluid shear stress, which simulates both rapid temporal and steady features of hemodynamic forces. Results indicate that ECs become significantly more compliant as early as 30 s after a step change in shear stress from 0 to 10 dyn/cm(2) followed by recovery of viscoelastic parameters within 4 min of shearing, even though shear stress was maintained. After ECs were sheared for 5 min, return of shear stress to 0 dyn/cm(2) in a stepwise manner did not result in any further rheological adaptation. Average vesicle displacements were used to determine time-dependent cell deformation and macrorheological parameters by fitting creep function to a linear viscoelastic liquid model. Characteristic time and magnitude for shear-induced deformation were 3 s and 50 nm, respectively. We conclude that ECs rapidly adapt their mechanical properties in response to shear stress, and we provide the first macrorheological parameters for time-dependent deformations of ECs to a physiological forcing function. Such studies provide insight into pathologies such as atherosclerosis, which may find their origins in EC mechanics.

Integrated multimodal microscopy, time-resolved fluorescence, and optical-trap rheometry: toward single molecule mechanobiology

Cells respond to forces through coordinated biochemical signaling cascades that originate from changes in single-molecule structure and dynamics and proceed to large-scale changes in cellular morphology and protein expression. To enable experiments that determine the molecular basis of mechanotransduction over these large time and length scales, we construct a confocal molecular dynamics microscope (CMDM). This system integrates total-internal-reflection fluorescence (TIRF), epifluorescence, differential interference contrast (DIC), and 3-D deconvolution imaging modalities with time-correlated single-photon counting (TCSPC) instrumentation and an optical trap. Some of the structures hypothesized to be involved in mechanotransduction are the glycocalyx, plasma membrane, actin cytoskeleton, focal adhesions, and cell-cell junctions. Through analysis of fluorescence fluctuations, single-molecule spectroscopic measurements [e.g., fluorescence correlation spectroscopy (FCS) and time-resolved fluorescence] can be correlated with these subcellular structures in adherent endothelial cells subjected to well-defined forces. We describe the construction of our multimodal microscope in detail and the calibrations necessary to define molecular dynamics in cell and model membranes. Finally, we discuss the potential applications of the system and its implications for the field of mechanotransduction.