A fluid-structure interaction model to characterize bone cell stimulation in parallel-plate flow chamber systems

Evaluating the influence on osteocyte mechanobiology within the lacunar-canalicular system for varying lacunar equancy and perilacunar elasticity: A multiscale fluid-structure interaction analysis

The lacunar morphology and perilacunar tissue properties of osteocytes in bone can vary under different physiological and pathological conditions. How these alterations collectively change the overall micromechanics of osteocytes in the lacunar-canalicular system (LCS) of an osteon still requires special focus. Therefore, a Haversian canal and LCS-based osteon model was established to evaluate the changes in the hydrodynamic environment around osteocytes under physiological loading using fluid-structure interaction analysis, followed by a sub-modelled finite element analysis to assess the mechanical responses of osteocytes and their components. Osteocytes were modelled with detailed configurations, including cytoplasm, nucleus, and cytoskeleton, and parametric variations in lacunar equancy (L.Eq) and perilacunar elasticity (Pl.E) were considered within the osteon model. The study aimed to conduct a comparative study among osteon models with varying L. Eq and Pl. E to check the resulting differences in osteocyte mechanobiology. The results demonstrated that the average mechanical stimulation of each subcellular component of osteocytes increased with decreases in L. Eq and Pl. E, reflecting conditions typically seen in young, healthy bone as per previous literature. However, hydrodynamic responses, such as fluid flow and fluid shear stress on osteocytes, varied proportionally with the elasticity difference between the bone matrix and the perilacunar region during Pl. E variation. Additionally, the findings revealed that a minimal percentage of energy was used to transmit mechanical responses through microtubules from the cell membrane to the nucleus, and this energy percentage increased with higher L. Eq. The outcomes of the study could help to quantify how the osteocyte microenvironment and its mechanosensitivity within cortical bone changes with L. Eq and Pl. E alterations in different bone conditions, from young to aged and healthy to diseased.

Comparing the bone regeneration potential between a trabecular bone and a porous scaffold through osteoblast migration and differentiation: A multiscale approach

Both cell migration and osteogenic differentiation are critical for successful bone regeneration. Therefore, understanding the mechanobiological aspects that govern these two processes is essential in designing effective scaffolds that promote faster bone regeneration. Studying these two factors at different locations is necessary to manage bone regeneration in various sections of a scaffold. Hence, a multiscale computational model was used to observe the mechanical responses of osteoblasts placed in different positions of the trabecular bone and gyroid scaffold. Fluid shear stresses in scaffolds at cell seeded locations (representing osteogenic differentiation) and strain energy densities in cells at cell substrate interface (representing cell migration) were observed as mechanical response parameters in this study. Comparison of these responses, as two critical factors for bone regeneration, between the trabecular bone and gyroid scaffold at different locations, is the overall goal of the study. This study reveals that the gyroid scaffold exhibits higher osteogenic differentiation and cell migration potential compared to the trabecular bone. However, the responses in the gyroid only mimic the trabecular bone in two out of nine positions. These findings can guide us in predicting the ideal cell seeded sites within a scaffold for better bone regeneration and in replicating a replaced bone condition by altering the physical parameters of a scaffold.

Computational simulation of applying mechanical vibration to mesenchymal stem cell for mechanical modulation toward bone tissue engineering

Evaluation of cell response to mechanical stimuli at in vitro conditions is known as one of the important issues for modulating cell behavior. Mechanical stimuli, including mechanical vibration and oscillatory fluid flow, act as important biophysical signals for the mechanical modulation of stem cells. In the present study, mesenchymal stem cell (MSC) consists of cytoplasm, nucleus, actin, and microtubule. Also, integrin and primary cilium were considered as mechanoreceptors. In this study, the combined effect of vibration and oscillatory fluid flow on the cell and its components were investigated using numerical modeling. The results of the FEM and FSI model showed that the cell response (stress and strain values) at the frequency of 30 H z mechanical vibration has the highest value. The achieved results on shear stress caused by the fluid flow on the cell showed that the cell experiences shear stress in the range of 0 . 1 - 10 Pa . Mechanoreceptors that bind separately to the cell surface, can be highly stimulated by hydrodynamic pressure and, therefore, can play a role in the mechanical modulation of MSCs at in vitro conditions. The results of this research can be effective in future studies to optimize the conditions of mechanical stimuli applied to the cell culture medium and to determine the mechanisms involved in mechanotransduction.

Evaluating the cell migration potential of TiO2 nanorods incorporated in a Ti6Al4V scaffold: A multiscale approach

Improvement of cell migration by the nano-topographical modification of implant surface can directly or indirectly accelerate wound healing and osseointegration between bone and implant. Therefore, modification of the implant surface was done with TiO2 nanorod (NR) arrays to develop a more osseointegration-friendly implant in this study. Modulating the migration of a cell, adhered to a scaffold, by the variations of NR diameter, density and tip diameter in vitro is the primary objective of the study. The fluid structure interaction method was used, followed by the submodelling technique in this multiscale analysis. After completing a simulation over a global model, fluid structure interaction data was applied to the sub-scaffold finite element model to predict the mechanical response over cells at the cell-substrate interface. Special focus was given to strain energy density at the cell interface as a response parameter due to its direct correlation with the migration of an adherent cell. The results showed a huge rise in strain energy density after the addition of NRs on the scaffold surface. It also highlighted that variation in NR density plays a more effective role than the variation in NR diameter to control cell migration over a substrate. However, the effect of NR diameter becomes insignificant when the NR tip was considered. The findings of this study could be used to determine the best nanostructure parameters for better osseointegration.

Finite element study of stem cells under fluid flow for mechanoregulation toward osteochondral cells

Investigating the effects of mechanical stimuli on stem cells under in vitro and in vivo conditions is a very important issue to reach better control on cellular responses like growth, proliferation, and differentiation. In this regard, studying the effects of scaffold geometry, steady, and transient fluid flow, as well as influence of different locations of the cells lodged on the scaffold on effective mechanical stimulations of the stem cells are of the main goals of this study. For this purpose, collagen-based scaffolds and implicit surfaces of the pore architecture was used. In this study, computational fluid dynamics and fluid-structure interaction method was used for the computational simulation. The results showed that the scaffold microstructure and the pore architecture had an essential effect on accessibility of the fluid to different portions of the scaffold. This leads to the optimization of shear stress and hydrodynamic pressure in different surfaces of the scaffold for better transportation of oxygen and growth factors as well as for optimized mechanoregulative responses of cell-scaffold interactions. Furthermore, the results indicated that the HP scaffold provides more optimizer surfaces to culture stem cells rather than Gyroid and IWP scaffolds. The results of exerting oscillatory fluid flow into the HP scaffold showed that the whole surface of the HP scaffold expose to the shear stress between 0.1 and 40 mPa and hydrodynamics factors on the scaffold was uniform. The results of this study could be used as an aid for experimentalists to choose optimist fluid flow conditions and suitable situation for cell culture.

Design and research of bone repair scaffold based on two-way fluid-structure interaction

BACKGROUND AND OBJECTIVE

Porous bone repair scaffolds are an important method of repairing bone defects. Fluid flow in the scaffold plays a vital role in tissue differentiation and permeability and fluid shear stress (FSS) are two important factors. The differentiation of bone tissue depends on the osteogenic differentiation of cells, FSS affects cell proliferation and differentiation, and permeability affects the transportation of nutrients and metabolic waste. Therefore, it is necessary to better understand and analyze the FSS on the cell surface and the permeability of the scaffold to obtain better osteogenic performance.

METHODS

In this study, computational fluid dynamics (CFD) was used to analyze fluid flow in the scaffold. Three structures and nine scaffold unit cell models were designed and the cell models were loaded onto the scaffold surface. Considering cell deformability, the two-way fluid-structure interaction (FSI) method was used to evaluate the FSS on the cell surface.

RESULTS

The simulation results showed that as the pore size of the scaffold increases, its permeability increases and the FSS decreases. The FSS received on the cell surface was much larger than scaffold surface. Moreover the FSS on the cell surface was distributed in steps.

CONCLUSIONS

The results showed the permeability of all models matches that of human bone tissue. Based on the cell surface FSS as the criterion, it was found that the spherical-560 scaffold exhibited the best osteogenic performance. This provided a strategy to design a better bone repair scaffold from biological aspects.

Stem Cell Mechanobiology and the Role of Biomaterials in Governing Mechanotransduction and Matrix Production for Tissue Regeneration

Mechanobiology has underpinned many scientific advances in understanding how biophysical and biomechanical cues regulate cell behavior by identifying mechanosensitive proteins and specific signaling pathways within the cell that govern the production of proteins necessary for cell-based tissue regeneration. It is now evident that biophysical and biomechanical stimuli are as crucial for regulating stem cell behavior as biochemical stimuli. Despite this, the influence of the biophysical and biomechanical environment presented by biomaterials is less widely accounted for in stem cell-based tissue regeneration studies. This Review focuses on key studies in the field of stem cell mechanobiology, which have uncovered how matrix properties of biomaterial substrates and 3D scaffolds regulate stem cell migration, self-renewal, proliferation and differentiation, and activation of specific biological responses. First, we provide a primer of stem cell biology and mechanobiology in isolation. This is followed by a critical review of key experimental and computational studies, which have unveiled critical information regarding the importance of the biophysical and biomechanical cues for stem cell biology. This review aims to provide an informed understanding of the intrinsic role that physical and mechanical stimulation play in regulating stem cell behavior so that researchers may design strategies that recapitulate the critical cues and develop effective regenerative medicine approaches.

Cancer cells optimize elasticity for efficient migration

Cancer progression is associated with alternations in the cytoskeletal architecture of cells and, consequently, their mechanical properties such as stiffness. Changing the mechanics of cells enables cancer cells to migrate and invade to distant organ sites. This process, metastasis, is the main reason for cancer-related mortality. Cell migration is an essential step towards increasing the invasive potential of cells. Although many studies have shown that the migratory speed and the invasion of cells can be inversely correlated to the stiffness of cells, some other investigations indicate opposing results. In the current work, based on the strain energy stored in cells due to the contractile forces, we defined an energy-dependent term, migratory index, to approximate how changes in the mechanical properties of cells influence cell migration required for cancer progression. Cell migration involves both cell deformation and force transmission within cells. The effects of these two parameters can be represented equally by the migratory index. Our mechanical modelling and computational study show that cells depending on their shape, size and other physical parameters have a maximum migratory index taking place at a specific range of cell bulk elasticity, indicating the most favourable conditions for invasive mobility. This approximate model could be used to explain why the stiffness of cells varies during cancer progression. We believe that the stiffness of cancer or malignant cells depending on the stiffness of their normal or non-malignant counterparts is either decreased or increased to reach the critical condition in which the mobility potential of cells is approximated to be maximum.

Versatile and High-throughput Force Measurement Platform for Dorsal Cell Mechanics

We present a high-throughput microfluidics technique facilitating in situ measurements of cell mechanics parameters at the dorsal side of the cell, including molecular binding strengths, local traction forces, and viscoelastic properties. By adjusting the flow rate, the force magnitude exerted on the cell can be modulated ranging from ~14 pN to 2 nN to perturb various force-dependent processees in cells. Time-lapse images were acquired to record events due to such perturbation. The values of various mechanical parameters are subsequently obtained by single particle tracking. Up to 50 events can be measured simultaneously in a single experiment. Integrating the microfluidic techniques with the analytic framework established in computational fluid dynamics, our method is physiologically relevant, reliable, economic and efficient.

Efficient Low Shear Flow-based Trapping of Biological Entities

Capturing cells or biological entities is an important and challenging step toward in-vitro studies of cells under a precisely controlled microscale environment. In this work, we have developed a compact and efficient microdevice for on-chip trapping of micro-sized particles. This hydrodynamics-based trapping system allows the isolation of polystyrene micro-particles with a shorter time while inducing a less hydrodynamic deformation and stress on the particles or cells both after and before trapping. A numerical simulation was carried out to design a hydrodynamic trapping mechanism and optimize the geometric and fluidic parameters affecting the trapping efficiency of the microfluidic network. By using the finite element analysis, the velocity field, pressure field, and hydrodynamic force on the micro particles were studied. Finally, a PDMS microfluidic device was fabricated to test the device's ability to trap polystyrene microspheres. Computational fluid analysis and experimental testing showed a high trapping efficiency that is more than 90%. This microdevice can be used for single cell studies including their biological, physical and chemical characterization.

Estrogen deficiency impairs integrin αvβ3-mediated mechanosensation by osteocytes and alters osteoclastogenic paracrine signalling

The integrin α_vβ₃ has been shown to play an important role in osteocyte mechanotransduction. It has been reported that there are fewer β₃ integrin-containing cells in osteoporotic bone cells. Osteocytes cultured in vitro under estrogen deficient conditions demonstrate altered mechanotransduction. However, it is unknown whether the altered mechanotransduction in estrogen deficient osteocytes is directly associated with defective α_vβ₃ expression or signalling. The objective of this study is to investigate the role of estrogen deficiency for regulating MLO-Y4 cell morphology, α_vβ₃ expression, focal adhesion formation and mechanotransduction by osteocytes. Here, we report that estrogen withdrawal leads to a smaller focal adhesion area and reduced α_vβ₃ localisation at focal adhesion sites, resulting in an increased Rankl/Opg ratio and defective Cox-2 responses to oscillatory fluid flow. Interestingly, α_vβ₃ antagonism had a similar effect on focal adhesion assembly, Rankl/Opg ratio, and Cox-2 responses to oscillatory fluid flow. Taken together, our results provide the first evidence for a relationship between estrogen withdrawal and defective α_vβ₃-mediated signalling. Specifically, this study implicates estrogen withdrawal as a putative mechanism responsible for altered α_vβ₃ expression and resultant changes in downstream signalling in osteocytes during post-menopausal osteoporosis, which might provide an important, but previously unidentified, contribution to the bone loss cascade.

Mechanical Stimuli in the Local In Vivo Environment in Bone: Computational Approaches Linking Organ-Scale Loads to Cellular Signals

PURPOSE OF REVIEW

Connecting organ-scale loads to cellular signals in their local in vivo environment is a current challenge in the field of bone (re)modelling. Understanding this critical missing link would greatly improve our ability to anticipate mechanotransduction during different modes of stimuli and the resultant cellular responses. This review characterises computational approaches that could enable coupling links across the multiple scales of bone.

RECENT FINDINGS

Current approaches using strain and fluid shear stress concepts have begun to link organ-scale loads to cellular signals; however, these approaches fail to capture localised micro-structural heterogeneities. Furthermore, models that incorporate downstream communication from osteocytes to osteoclasts, bone-lining cells and osteoblasts, will help improve the understanding of (re)modelling activities. Incorporating this potentially key information in the local in vivo environment will aid in developing multiscale models of mechanotransduction that can predict or help describe resultant biological events related to bone (re)modelling. Progress towards multiscale determination of the cell mechanical environment from organ-scale loads remains elusive. Construction of organ-, tissue- and cell-scale computational models that include localised environmental variation, strain amplification and intercellular communication mechanisms will ultimately help couple the hierarchal levels of bone.

Cellular fluid shear stress on implant surfaces-establishment of a novel experimental set up

Background

Mechanostimuli of different cells can affect a wide array of cellular and inter-cellular biological processes responsible for dental implant healing. The purpose of this in vitro study was to establish a new test model to create a reproducible flow-induced fluid shear stress (FSS) of osteoblast cells on implant surfaces.

Methods

As FSS effects on osteoblasts are detectable at 10 dyn/cm², a custom-made flow chamber was created. Computer-aided verification of circulation processes was performed. In order to verify FSS effects, cells were analysed via light and fluorescence microscopy.

Results

Utilising computer-aided simulations, the underside of the upper plate was considered to have optimal conditions for cell culturing. At this site, a flow-induced orientation of osteoblast cell clusters and an altered cell morphology with cellular elongation and alteration of actin fibres in the fluid flow direction was detected.

Conclusions

FSS simulation using this novel flow chamber might mimic the peri-implant situation in the phase of loaded implant healing. With this FSS flow chamber, osteoblast cells’ sensitivity to FSS was verified in the form of morphological changes and cell re-clustering towards the direction of the flow. Different shear forces can be created simultaneously in a single experiment.

Oscillatory fluid flow induces the osteogenic lineage commitment of mesenchymal stem cells: The effect of shear stress magnitude, frequency, and duration

A potent regulator of bone anabolism is physical loading. However, it is currently unclear whether physical stimuli such as fluid shear within the marrow cavity is sufficient to directly drive the osteogenic lineage commitment of resident mesenchymal stem cells (MSC). Therefore, the objective of the study is to employ a systematic analysis of oscillatory fluid flow (OFF) parameters predicted to occur in vivo on early MSC osteogenic responses and late stage lineage commitment. MSCs were exposed to OFF of 1Pa, 2Pa and 5Pa magnitudes at frequencies of 0.5Hz, 1Hz and 2Hz for 1h, 2h and 4h of stimulation. Our findings demonstrate that OFF elicits a positive osteogenic response in MSCs in a shear stress magnitude, frequency, and duration dependent manner that is gene specific. Based on the mRNA expression of osteogenic markers Cox2, Runx2 and Opn after short-term fluid flow stimulation, we identified that a regime of 2Pa shear magnitude and 2Hz frequency induces the most robust and reliable upregulation in osteogenic gene expression. Furthermore, long-term mechanical stimulation utilising this regime, elicits a significant increase in collagen and mineral deposition when compared to static control demonstrating that mechanical stimuli predicted within the marrow is sufficient to directly drive osteogenesis.

In silico bone mechanobiology: modeling a multifaceted biological system

Mechanobiology, the study of the influence of mechanical loads on biological processes through signaling to cells, is fundamental to the inherent ability of bone tissue to adapt its structure in response to mechanical stimulation. The immense contribution of computational modeling to the nascent field of bone mechanobiology is indisputable, having aided in the interpretation of experimental findings and identified new avenues of inquiry. Indeed, advances in computational modeling have spurred the development of this field, shedding new light on problems ranging from the mechanical response to loading by individual cells to tissue differentiation during events such as fracture healing. To date, in silico bone mechanobiology has generally taken a reductive approach in attempting to answer discrete biological research questions, with research in the field broadly separated into two streams: (1) mechanoregulation algorithms for predicting mechanobiological changes to bone tissue and (2) models investigating cell mechanobiology. Future models will likely take advantage of advances in computational power and techniques, allowing multiscale and multiphysics modeling to tie the many separate but related biological responses to loading together as part of a larger systems biology approach to shed further light on bone mechanobiology. Finally, although the ever‐increasing complexity of computational mechanobiology models will inevitably move the field toward patient‐specific models in the clinic, the determination of the context in which they can be used safely for clinical purpose will still require an extensive combination of computational and experimental techniques applied to in vitro and in vivo applications. WIREs Syst Biol Med 2016, 8:485–505. doi: 10.1002/wsbm.1356

For further resources related to this article, please visit the WIREs website.

Cell morphology and focal adhesion location alters internal cell stress

Extracellular mechanical cues have been shown to have a profound effect on osteogenic cell behaviour. However, it is not known precisely how these cues alter intracellular mechanics to initiate changes in cell behaviour. In this study, a combination of in vitro culture of MC3T3-E1 cells and finite-element modelling was used to investigate the effects of passive differences in substrate stiffness on intracellular mechanics. Cells on collagen-based substrates were classified based on the presence of cell processes and the dimensions of various cellular features were quantified. Focal adhesion (FA) density was quantified from immunohistochemical staining, while cell and substrate stiffnesses were measured using a live-cell atomic force microscope. Computational models of cell morphologies were developed using an applied contraction of the cell body to simulate active cell contraction. The results showed that FA density is directly related to cell morphology, while the effect of substrate stiffness on internal cell tension was modulated by both cell morphology and FA density, as investigated by varying the number of adhesion sites present in each morphological model. We propose that the cells desire to achieve a homeostatic stress state may play a role in osteogenic cell differentiation in response to extracellular mechanical cues.

The role of integrin α(V)β(3) in osteocyte mechanotransduction

Recent in vivo studies have proposed that integrin αvβ3 attachments between osteocyte cell processes and the extracellular matrix may facilitate mechanosensation in bone. However the role of these attachments in osteocyte biochemical response to mechanical stimulus has yet to be investigated. With this in mind, the objective of this study was to determine the effect of blocking integrin αvβ3 function on the biochemical response of osteocytes to mechanical stimulus. Antagonists specific to integrin subunit β3 were used to block integrin αvβ3 on MLO-Y4 mouse osteocytes. After treatment, cells were subjected to laminar oscillatory fluid flow stimulus (1 Pa, 1 Hz) for one hour. Fluorescent staining was performed to visualise cell morphology. Prostaglandin E2 (PGE2) release was assayed using an enzyme immunoassay and qRT-PCR was used to analyse the relative expression of cyclooxygenase-2 (COX-2), receptor activator of NF-κB ligand (RANKL) and osteoprotegerin (OPG). Our results show that blocking integrin αvβ3 disrupts osteocyte morphology, causing a reduction in spread area and process retraction. Integrin αvβ3 blocking also disrupted COX-2 expression and PGE2 release in response to fluid shear stress. Taken together, the results of this study indicate that integrin αvβ3 is essential for the maintenance of osteocyte cell processes and also for mechanosensation and mechanotransduction by osteocytes. A better understanding of this process may lead to the development of novel treatments for bone pathologies where mechanosensitivity is thought to be compromised.

Multiscale fluid-structure interaction modelling to determine the mechanical stimulation of bone cells in a tissue engineered scaffold

Recent studies have shown that mechanical stimulation, by means of flow perfusion and mechanical compression (or stretching), enhances osteogenic differentiation of mesenchymal stem cells and bone cells within biomaterial scaffolds in vitro. However, the precise mechanisms by which such stimulation enhances bone regeneration is not yet fully understood. Previous computational studies have sought to characterise the mechanical stimulation on cells within biomaterial scaffolds using either computational fluid dynamics or finite element (FE) approaches. However, the physical environment within a scaffold under perfusion is extremely complex and requires a multiscale and multiphysics approach to study the mechanical stimulation of cells. In this study, we seek to determine the mechanical stimulation of osteoblasts seeded in a biomaterial scaffold under flow perfusion and mechanical compression using multiscale modelling by two-way fluid-structure interaction and FE approaches. The mechanical stimulation, in terms of wall shear stress (WSS) and strain in osteoblasts, is quantified at different locations within the scaffold for cells of different attachment morphologies (attached, bridged). The results show that 75.4 % of scaffold surface has a WSS of 0.1-10 mPa, which indicates the likelihood of bone cell differentiation at these locations. For attached and bridged osteoblasts, the maximum strains are 397 and 177,200 με, respectively. Additionally, the results from mechanical compression show that attached cells are more stimulated (maximum strain = 22,600 με) than bridged cells (maximum strain = 10.000 με)Such information is important for understanding the biological response of osteoblasts under in vitro stimulation. Finally, a combination of perfusion and compression of a tissue engineering scaffold is suggested for osteogenic differentiation.