Fluid flow in the osteocyte mechanical environment: a fluid–structure interaction approach

The Flow of Life: Convergent Approaches to Understanding Musculoskeletal Health from Molecular- to Meso-Length Scales

In the current perspective and review article, we address the human body as a living ecosystem with collecting watersheds and draining hydrosheds; we integrate our discoveries over the past quarter of a century and pose the critical open research questions to be addressed going forward, with the aim to improve cell, tissue, organ and organismal health. First, we address the flow of fluid through the tissues of the musculoskeletal system, after which we describe the interactions of the fluid, at multiple lengths and time scales, with the molecular to macroscopic non-fluid tissue components, discussing bone and tissues in the context of "living" chromatography and/or electrophoresis columns. Thereafter, we discuss the implications of functional barrier integrity, and the effects of cytokines on active barrier function and molecular transport between organ systems, tissue compartments, and within tissues. In addition, we address the fluid and its flow and the multi-physics implications thereof for the living inhabitants of tissues, i.e., the cells. Finally, we describe the implications of the solid and fluid components and the cellular inhabitants on ecosystem health, where the tissues and organs comprise the organism form interacting ecosystems throughout life and in the context of health and disease. By taking convergent approaches to understanding musculoskeletal, human and environmental health (which themselves are interdependent), we hope to pave new paths of innovation and discovery, to improve the lives of our worlds' inhabitants, from the worlds of our bone and joints and bodies to the interacting ecosystems of our Earth to unknown worlds beyond our current understanding.

Multiscale interstitial fluid computation modeling of cortical bone to characterize the hydromechanical stimulation of lacunar-canalicular network

Bone tissue is a biological composite material with a complex hierarchical structure that could continuously adjust its internal structure to adapt to the alterations in the external load environment. The fluid flow within bone is the main route of osteocyte metabolism, and the pore pressure as well as the fluid shear stress generated by it are important mechanical stimuli perceived by osteocytes. Owing to the irregular multiscale structure of bone tissue, the fluid stimulation that lacunar-canalicular network (LCN) in different regions of the tissue underwent remained unclear. In this study, we constructed a multiscale conduction model of fluid flow stimulus signals in bone tissue based on the poroelasticity theory. We analyzed the fluid flow behaviors at the macro-scale (whole bone tissue), macro-meso scale (periosteum, interstitial bone, osteon and endosteum), and micro-scale (lacunar-osteocyte-canalicular) levels. We explored how fluid stimulation at the tissue level correlated with that at the cellular level in cortical bone and characterized the distributions of the pore pressure, fluid velocity and fluid shear stress that the osteocytes experienced across the entire tissue structure. The results showed that the initial conditions of intramedullary pressure had a significant impact on the pore pressure of Haversian systems, but had a relatively small influence on the fluid velocity. The osteocyte which were located at different positions in the bone tissue received very distinct fluid stimuli. Osteocytes in the vicinity of the Haversian Canals experienced higher fluid shear stress stimulation. When the permeability of the LCN was within the range from 10-21 m2 to 10-18 m2, the distribution of pressure, fluid velocity and fluid shear stress within the osteon near the periosteum and endosteum was significantly different from that in other parts of the bone. However, when the permeability was less than 10-22 m2, such a difference did not exist. Particularly, the flow velocity at the lacunae was markedly higher than that in the canaliculi. Meanwhile, the pore pressure and fluid shear stress were conspicuously lower than those in the canaliculi. In this study, we considered the interconnections of different biofunctional units at different scales of bone tissue, construct a more complete multiscale model of bone tissue, and propose that osteocytes at different locations receive different fluid stimuli, which provides a reference for a deeper understanding of bone mechanotransduction.

Scaffold geometries designed to promote bone ingrowth by enhancing mechanobiological stimulation and biotransportation - A multiobjective optimisation approach

In a tissue-engineered bone scaffold implant, the process of neo-tissue ingrowth and remodelling into hard lamellar bone occurs slowly; it typically requires a period of several months to a year (or more) to complete. This research considers the design optimisation of a scaffold's unit cell geometry for the purpose of accelerating the rate at which neo-tissue forms in the porous network of the scaffold (ingrowth), and hence, reduce the length of time to complete the bone ingrowth process. In this study, the basic structure of the scaffold is the Schwarz Primitive (P) surface unit cell, selected for its compelling biomechanical and permeability characteristics. The geometry of the scaffold is varied using two parameters (namely iso-value, k, and spatial period, a) within the surface equation defining the Schwarz P-surface unit cell. In total, sixteen different unit cell geometries are considered here with the porosity ranging from 50% to 82%. The design objectives for the scaffold are to (i) enhance mechanobiological stimulus conditions conducive to bone apposition and (ii) enhance permeability to improve the transport of nutrients/oxygen and metabolities to and from the sites of neo-tissue formation throughout the porous scaffold. The independent design variables (k and a) of the periodic unit cell geometry are optimised to best satisfy the design objectives of appositional mechanobiological stimulus and biotransporting permeability. First, a reconstructed sheep mandible computed tomographic (CT)-based finite element (FE) analysis model is used to relate the strain energy density and mechanobiological stimulus to the design variables. Next, a computational fluid dynamics (CFD) model of a 5 × 5 × 5 unit cell scaffold is developed to relate the distributions of pressure and fluid velocity to the design variables. Then, surrogate modelling is undertaken in which bivariate cubic polynomial response surfaces are fitted to the FE and CFD analysis output data to form mathematical functions of each objective with respect to the two design variables. Finally, a multiobjective optimisation algorithm is invoked to determine the best trade-off between the competing design objectives of mechanobiological stimulus and biofluidic permeability. The novel design of the scaffold structure is anticipated to provide a better biomechanical and biotransport environment for tissue regeneration.

A systematic study of the effect of measurement parameters on determination of osteocyte lacunar properties using laboratory X-ray micro-computed tomography

Accurate 3D characterization of osteocyte lacunae is important when investigating the role of osteocytes under various physiological and pathological conditions but remains a challenge. With the continued development of laboratory X-ray micro-computed tomography, an increasing number of studies employ these techniques beyond traditional bone morphometry to quantify osteocyte lacunae. However, there is a lack of knowledge on the effect of measurement parameters on the image quality and resolution and in turn the osteocyte lacunar quantification. Herein, we have examined the interplay between scan parameters and the resultant lacunar quantification in terms of lacunar size, shape, and density by comparison with a synchrotron benchmark dataset. We summarize our conclusions in a guide for use of μ-CT for osteocyte lacunar quantification: (1) Identification of the measurement requirements to address the research questions. (2) Collection and preparation of suitable sample(s) that fulfills these requirements. (3) Experimental considerations including determination of the required voxel size, in turn dictating the maximum FOV and by extension the maximum size of the sample(s). The experimental parameters chosen should ensure optimal image contrast, sufficient signal to noise, angular sampling etc. Usually, it is advisable to measure as well as possible within the limits of time, budget, data storage and analysis capabilities. (4) Data analysis and reporting of the results, including visual examination of the data at multiple steps in the analysis, to ensure correct feature identification and suitable reporting approaches. (5) Cross study comparisons, which may be unsuitable if the experimental conditions and analysis strategies are not comparable.

Numerical study of interstitial fluid flow behavior in osteons under dynamic loading

BACKGROUND

The porous structure in bone tissue is essential for maintaining the physiological functions and overall health of intraosseous cells. The lacunar-canalicular net (LCN), a microscopic porous structure within osteons, facilitates the transport of nutrients and signaling molecules through interstitial fluid flow. However, the transient behavior of fluid flow within these micro-pores under dynamic loading conditions remains insufficiently studied.

METHODS

The study constructs a fluid-solid coupling model including the Haversian canal, canaliculi, lacunae, and interstitial fluid, to examine interstitial fluid flow behavior within the LCN under dynamic loading with varying frequencies and amplitudes. The relationship between changes of LCN pore volume and fluid velocity, and pressure is researched.

RESULTS

The results demonstrate that increasing strain amplitude leads to significant changes of LCN pore volume within osteons. In a complete loading cycle, with the increase of compressive strain, the pore volume in the osteon gradually shrinks, and the pressure gradient in the LCN increases, which promotes the increase of interstitial fluid velocity. When the compressive strain reaches the peak value, the flow velocity also reaches the maximum. In the subsequent unloading process, the pore volume began to recover, the pressure gradient gradually decreased, the flow rate decreased accordingly, and finally returned to the steady state level. At a loading amplitude of 1000 µε, the pore volume within LCN decreases by 1.1‰. At load amplitudes of 1500 µε, 2000 µε, and 2500 µε, the pore volume decreases by 1.6‰, 2.2‰ and 2.7‰ respectively, and the average flow velocity at the center of the superficial lacuna is 1.36 times, 1.77 times, and 2.14 times that at 1000 µε, respectively. Additionally, at a loading amplitude of 1000 µε under three different loading frequencies, the average flow velocities at the center of the superficial bone lacuna are 0.60 μm/s, 1.04 μm/s, and 1.54 μm/s, respectively. This indicates that high-frequency and high-amplitude dynamic loading can promote more vigorous fluid flow and pressure fluctuations with changes in LCN pore volume.

CONCLUSIONS

Dynamic mechanical loading can significantly enhance the interstitial fluid flow in LCN by the changes of LCN pore volume. and dynamic loading promoted fluid flow in shallow lacunae significantly higher than that in deep lacunae. The relationship between changes of LCN pore volume and interstitial fluid flow behavior has implications for drug delivery and bone tissue engineering research.

Non-Newtonian lacuno-canalicular fluid flow in bone altered by mechanical loading and magnetic field

As humans age, they experience deformity and a decrease in their bone strength, such brittleness in the bones ultimately lead to bone fracture. Magnetic field exposure combined with physical exercise may be useful in mitigating age-related bone loss by improving the canalicular fluid motion within the bone's lacuno-canalicular system (LCS). Nevertheless, an adequate amount of fluid induced shear stress is necessary for the bone mechano-transduction and solute transport in the case of brittle bone diseases. The underlying mechanisms of how magnetic fields, in combination to mechanical loading, affects the canalicular fluid motion still need to be explored. Accordingly, this study aims to develop a computer model to investigate the role of magnetic fields on loading-induced canalicular fluid flow in a curvy lacunar canalicular space with irregular osteocyte cell processes and walls. Moreover, this study considers canalicular fluid as non-Newtonian fluid, i.e., Jeffery fluid. In addition, a machine learning model was further employed for the estimation of parameters which significantly influence the canalicular fluid flow in response to loading and magnetic field. The results show that static magnetic field modulates the loading-induced canalicular fluid flow. Additionally, present study accelerates the fluid induced wall shear stress in case of osteoporosis.

Numerical Simulation of Fluid Shear Stress Distribution in Microcracks of Trabecular Bone

Bone is one of the hardest tissues in the human body, but it can undergo microcracks under long-term and periodic mechanical loads. The Newton iterative method was used to calculate the steady state, and the effects of different inlet and outlet pressures, trabecular gap width and height, and microcrack's depth and width on the fluid shear stress (FSS) were studied, and the gradient of FSS inside the microcrack was analyzed. The results show that the pressure difference and trabecular gap heigh are positively correlated with the FSS (the linear correlation coefficients R 2 were 0.9768 and 0.96542, respectively). When the trabecular gap width was 100 μm, the peak of FSS decreased by 28.57% compared with 800 and 400 μm, and the gradient of FSS inside the microcrack was 0.1-0.4 Pa/mm. This study can help people more intuitively understand the internal fluid distribution of trabecular bone and provide a reliable theoretical basis for the subsequent construction of gradient FSS devices in vitro.

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.

Multiscale simulation of the effect of low-intensity pulsed ultrasound on the mechanical properties distribution of osteocytes

Low-intensity pulsed ultrasound (LIPUS) is a potential effective means for the prevention and treatment of disuse osteoporosis. In this paper, the effect of LIPUS exposure on the mechanical properties distribution of the osteocyte system (osteocyte body contains nucleus, osteocyte process, and primary cilia) is simulated. The results demonstrate that the mechanical micro-environment of the osteocyte is significantly improved by ultrasound exposure, and the mean von Mises stress of the osteocyte system increases linearly with the excitation sound pressure amplitude. The mechanical effect of LIPUS on osteocytes is enhanced by the stress amplification mechanism of the primary cilia and osteocyte processes.

Osteocyte Sptbn1 Deficiency Alters Cell Survival and Mechanotransduction Following Formation of Plasma Membrane Disruptions (PMD) from Mechanical Loading

We and others have shown that application of high-level mechanical loading promotes the formation of transient plasma membrane disruptions (PMD) which initiate mechanotransduction. We hypothesized that increasing osteocyte cell membrane fragility, by disrupting the cytoskeleton-associated protein β2-spectrin (Sptbn1), could alter osteocytic responses and bone adaptation to loading in a PMD-related fashion. In MLO-Y4 cells, treatment with the spectrin-disrupting agent diamide or knockdown of Sptbn1 via siRNA increased the number of PMD formed by fluid shear stress. Primary osteocytes from an osteocyte-targeted DMP1-Cre Sptbn1 conditional knockout (CKO) model mimicked trends seen with diamide and siRNA treatment and suggested the creation of larger PMD, which repaired more slowly, for a given level of stimulus. Post-wounding cell survival was impaired in all three models, and calcium signaling responses from the wounded osteocyte were mildly altered in Sptbn1 CKO cultures. Although Sptbn1 CKO mice did not demonstrate an altered skeletal phenotype as compared to WT littermates under baseline conditions, they showed a blunted increase in cortical thickness when subjected to an osteogenic tibial loading protocol as well as evidence of increased osteocyte death (increased lacunar vacancy) in the loaded limb after 2 weeks of loading. The impaired post-wounding cell viability and impaired bone adaptation seen with Sptbn1 disruption support the existence of an important role for Sptbn1, and PMD formation, in osteocyte mechanotransduction and bone adaptation to mechanical loading.

Prkd1 regulates the formation and repair of plasma membrane disruptions (PMD) in osteocytes

We and others have seen that osteocytes sense high-impact osteogenic mechanical loading via transient plasma membrane disruptions (PMDs) which initiate downstream mechanotransduction. However, a PMD must be repaired for the cell to survive this wounding event. Previous work suggested that the protein Prkd1 (also known as PKCμ) may be a critical component of this PMD repair process, but the specific role of Prkd1 in osteocyte mechanobiology had not yet been tested. We treated MLO-Y4 osteocytes with Prkd1 inhibitors (Go6976, kbNB 142-70, staurosporine) and generated an osteocyte-targeted (Dmp1-Cre) Prkd1 conditional knockout (CKO) mouse. PMD repair rate was measured via laser wounding and FM1-43 dye uptake, PMD formation and post-wounding survival were assessed via fluid flow shear stress (50 dyn/cm2), and in vitro osteocyte mechanotransduction was assessed via measurement of calcium signaling. To test the role of osteocyte Prkd1 in vivo, Prkd1 CKO and their wildtype (WT) littermates were subjected to 2 weeks of unilateral axial tibial loading and loading-induced changes in cortical bone mineral density, geometry, and formation were measured. Prkd1 inhibition or genetic deletion slowed osteocyte PMD repair rate and impaired post-wounding cell survival. These effects could largely be rescued by treating osteocytes with the FDA-approved synthetic copolymer Poloxamer 188 (P188), which was previously shown to facilitate membrane resealing and improve efficiency in the repair rate of PMD in skeletal muscle myocytes. In vivo, while both WT and Prkd1 CKO mice demonstrated anabolic responses to tibial loading, the magnitude of loading-induced increases in tibial BMD, cortical thickness, and periosteal mineralizing surface were blunted in Prkd1 CKO as compared to WT mice. Prkd1 CKO mice also tended to show a smaller relative difference in the number of osteocyte PMD in loaded limbs and showed greater lacunar vacancy, suggestive of impaired post-wounding osteocyte survival. While P188 treatment rescued loading-induced increases in BMD in the Prkd1 CKO mice, it surprisingly further suppressed loading-induced increases in cortical bone thickness and cortical bone formation. Taken together, these data suggest that Prkd1 may play a pivotal role in the regulation and repair of the PMD response in osteocytes and support the idea that PMD repair processes can be pharmacologically targeted to modulate downstream responses, but suggest limited utility of PMD repair-promoting P188 in improving bone anabolic responses to loading.

β-glycerophosphate, not low magnitude fluid shear stress, increases osteocytogenesis in the osteoblast-to-osteocyte cell line IDG-SW3

AIM

As osteoblasts deposit a mineralized collagen network, a subpopulation of these cells differentiates into osteocytes. Biochemical and mechanical stimuli, particularly fluid shear stress (FSS), are thought to regulate this, but their relative influence remains unclear. Here, we assess both biochemical and mechanical stimuli on long-term bone formation and osteocytogenesis using the osteoblast-osteocyte cell line IDG-SW3.

METHODS

Due to the relative novelty and uncommon culture conditions of IDG-SW3 versus other osteoblast-lineage cell lines, effects of temperature and media formulation on matrix deposition and osteocytogenesis were initially characterized. Subsequently, the relative influence of biochemical (β-glycerophosphate (βGP) and ascorbic acid 2-phosphate (AA2P)) and mechanical stimulation on osteocytogenesis was compared, with intermittent application of low magnitude FSS generated by see-saw rocker.

RESULTS

βGP and AA2P supplementation were required for mineralization and osteocytogenesis, with 33°C cultures retaining a more osteoblastic phenotype and 37°C cultures undergoing significantly higher osteocytogenesis. βGP concentration positively correlated with calcium deposition, whilst AA2P stimulated alkaline phosphatase (ALP) activity and collagen deposition. We demonstrate that increasing βGP concentration also significantly enhances osteocytogenesis as quantified by the expression of green fluorescent protein linked to Dmp1. Intermittent FSS (~0.06 Pa) rocker had no effect on osteocytogenesis and matrix deposition.

CONCLUSIONS

This work demonstrates the suitability and ease with which IDG-SW3 can be utilized in osteocytogenesis studies. IDG-SW3 mineralization was only mediated through biochemical stimuli with no detectable effect of low magnitude FSS. Osteocytogenesis of IDG-SW3 primarily occurred in mineralized areas, further demonstrating the role mineralization of the bone extracellular matrix has in osteocyte differentiation.

Spatial variations in the osteocyte lacuno-canalicular network density and analysis of the connectomic parameters

Osteocyte lacuno-canalicular network (LCN) is comprised of micrometre-sized pores and submicrometric wide channels in bone. Accumulating evidence suggests multiple functions of this network in material transportation, mechanobiological signalling, mineral homeostasis and bone remodelling. Combining rhodamine staining and confocal laser scanning microscopy, the longitudinal cross-sections of six mouse tibiae were imaged, and the connectome of the network was quantified with a focus on the spatial heterogeneities of network density, connectivity and length of canaliculi. In-vivo loading and double calcein labelling on these tibiae allowed differentiating the newly formed bone from the pre-existing regions. The canalicular density of the murine cortical bone varied between 0.174 and 0.243 μm/μm3, and therefore is three times larger than the corresponding value for human femoral midshaft osteons. The spatial heterogeneity of the network was found distinctly more pronounced across the cortex than along the cortex. We found that in regions with a dense network, the LCN conserves its largely tree-like character, but increases the density by including shorter canaliculi. The current study on healthy mice should serve as a motivating starting point to study the connectome of genetically modified mice, including models of bone diseases and of reduced mechanoresponse.

Numerical calculation of streaming potential around osteocytes under human gait loading

Streaming potential is a type of stress-generated potential in bone that affects the electrical environment of osteocytes and may play a role in bone remodeling. Because the electrical environment around osteocytes has been difficult to measure experimentally until now, a numerical solid-liquid-streaming potential coupling method was proposed to analyze the streaming potential generated by bone deformation in the lacunae and canaliculus network (LCN) of the bone. Using this method, the cellular shear stress caused by liquid flow on the osteocyte surface was first calculated, and the results were consistent with those reported in the literature. Subsequently, the streaming potentials in the LCN caused by bone matrix deformation under an external gait load were calculated numerically. The results showed that the streaming potential increased slowly in the lacuna and relatively rapidly in the canaliculus and that the streaming potential increased with a decrease in the radius or an increase in the length of the canaliculus. The results also showed that relatively large gaps between the lacunae and osteocytes could induce higher streaming potentials under the same loading.

Computer model of non-Newtonian canalicular fluid flow in lacunar-canalicular system of bone tissue

Brittle bone diseases are a global healthcare problem for orthopaedic clinicians, that reduces bone strength and promotes bone fracture risk. In vivo studies reported that loading-induced fluid flow through the lacunar-canalicular channel (LCS) of bone tissue inhibit such bone loss and encourages osteogenesis i.e. new bone formation. Canalicular fluid flow converts mechanical signals into biological signals and regulates bone reconstruction by releasing signalling molecules responsible for mechanotransduction. In-silico model mostly considers canalicular fluid is Newtonian, however, physiological canalicular fluid may be non-Newtonian in nature as it contains nutrients and supplements. Accordingly, this study attempts to develop a two-dimensional in-silico model to compute loading-induced non-Newtonian canalicular fluid flow in a complex LCS of bone tissue. Moreover, canalicular fluid is considered as a Jeffery fluid, that can easily be reduced to Newtonian fluid as a special case. The results show that physiological loading modulates the canalicular fluid flow, wall shear stress (WSS) and streamline in bone LCS. Fluid velocity and WSS increases with increase in non-dimensional frequency and non-Newtonian parameter (Jeffery fluid parameters) and reduce with change in permeability. The outcomes of this study may provide new insights in the role of mechanical loading-induced non-Newtonian canalicular fluid flow dynamics in bone LCS. The key findings of this study can be used to improve the understanding of osteocyte mechanobiology involved inside the bone tissue.

Computational modeling based on confocal imaging predicts changes in osteocyte and dendrite shear stress due to canalicular loss with aging

Exercise and physical activity exert mechanical loading on the bones which induces bone formation. However, the relationship between the osteocyte lacunar-canalicular morphology and mechanical stress experienced locally by osteocytes transducing signals for bone formation is not fully understood. In this study, we used computational modeling to predict the effect of canalicular density, the number of fluid inlets, and load direction on fluid flow shear stress (FFSS) and bone strains and how these might change following the microstructural deterioration of the lacunar-canalicular network that occurs with aging. Four distinct computational models were initially generated of osteocytes with either ten or eighteen dendrites using a fluid-structure interaction method with idealized geometries. Next, a young and a simulated aged osteocyte were developed from confocal images after FITC staining of the femur of a 4-month-old C57BL/6 mouse to estimate FFSS using a computational fluid dynamics approach. The models predicted higher fluid velocities in the canaliculi versus the lacunae. Comparison of idealized models with five versus one fluid inlet indicated that with four more inlets, one-half of the dendrites experienced FFSS greater than 0.8 Pa, which has been associated with osteogenic responses. Confocal image-based models of real osteocytes indicated a six times higher ratio of canalicular to lacunar surface area in the young osteocyte model than the simulated aged model and the average FFSS in the young model (FFSS = 0.46 Pa) was three times greater than the aged model (FFSS = 0.15 Pa). Interestingly, the surface area with FFSS values above 0.8 Pa was 23 times greater in the young versus the simulated aged model. These findings may explain the impaired mechano-responsiveness of osteocytes with aging.

Fluid-solid coupling numerical simulation of entire rat caudal vertebrae under dynamic loading

Trabeculae bone undergoes directional growth along the applied force under physiological loading. The growth of bone structure relies on the coordinated interplay among osteocytes, osteoblasts, and osteoclasts. Under normal circumstances, bone remodeling maintains a state of equilibrium. Excessive bone formation can lead to osteosclerosis, while excessive bone resorption can result in osteoporosis and osteonecrosis. The investigation of the structural characteristics of trabeculae and the mechanotransduction between bone cells plays a vital role in the treatment of bone-related diseases. In this study, a fluid-solid coupling model of the entire vertebral bone was established based on micro-CT images obtained from rat tail vertebrae subjected to tensile loading experiments. The flow characteristics of bone marrow and the mechanical response of osteocytes in different regions under physiological loading were investigated. The results revealed a U-shaped distribution of wall fluid shear stress (FSS) along the longitudinal axis in trabecular bone, with higher FSS regions exhibiting greater mechanical stimulation on osteocytes. These findings elucidate a positive correlation between the mechanical microenvironment among osteocytes, osteoblasts, and osteoclasts, providing potential strategies for the prevention and treatment of bone diseases.

A Novel Primary Cilium-Mediated Mechanism Through which Osteocytes Regulate Metastatic Behavior of Both Breast and Prostate Cancer Cells

Bone metastases are a common cause of suffering in breast and prostate cancer patients, however, the interaction between bone cells and cancer cells is poorly understood. Using a series of co-culture, conditioned media, human cancer spheroid, and organ-on-a-chip experiments, this study reveals that osteocytes suppress cancer cell proliferation and increase migration via tumor necrosis factor alpha (TNF-α) secretion. This action is regulated by osteocyte primary cilia and associated intraflagellar transport protein 88 (IFT88). Furthermore, it shows that cancer cells block this mechanism by secreting transforming growth factor beta (TGF-β), which disrupts osteocyte cilia and IFT88 gene expression. This bi-directional crosstalk signaling between osteocytes and cancer cells is common to both breast and prostate cancer. This study also proposes that osteocyte inhibition of cancer cell proliferation decreases as cancer cells increase, producing more TGF-β. Hence, a positive feedback loop develops accelerating metastatic tumor growth. These findings demonstrate the importance of cancer cell-osteocyte signaling in regulating breast and prostate bone metastases and support the development of therapies targeting this pathway.

Deep learning models to map osteocyte networks can successfully distinguish between young and aged bone

Osteocytes, the most abundant and mechanosensitive cells in bone tissue, play a pivotal role in bone homeostasis and mechano-responsiveness, orchestrating the intricate balance between bone formation and resorption under daily activity. Studying osteocyte connectivity and understanding their intricate arrangement within the lacunar canalicular network (LCN) is essential for unraveling bone physiology. This is particularly true as our bones age, which is associated with decreased integrity of the osteocyte network, disrupted mass transport, and lower sensitivity to the mechanical stimuli that allow the skeleton to adapt to changing demands. Much work has been carried out to investigate this relationship, often involving high resolution microscopy of discrete fragments of this network, alongside advanced computational modelling of individual cells. However, traditional methods of segmenting and measuring osteocyte connectomics are time-consuming and labour-intensive, often hindered by human subjectivity and limited throughput. In this study, we explore the application of deep learning and computer vision techniques to automate the segmentation and measurement of osteocyte connectomics, enabling more efficient and accurate analysis. We compare several state-of-the-art computer vision models (U-Nets and Vision Transformers) to successfully segment the LCN, finding that an Attention U-Net model can accurately segment and measure 81.8% of osteocytes and 42.1% of dendritic processes, when compared to manual labelling. While further development is required, we demonstrate that this degree of accuracy is already sufficient to distinguish between bones of young (2 month old) and aged (36 month old) mice, as well as capturing the degeneration induced by genetic modification of osteocytes. By harnessing the power of these advanced technologies, further developments can unravel the complexities of osteocyte networks in unprecedented detail, revolutionising our understanding of bone health and disease.

JANUS: an open-source 3D printable perfusion bioreactor and numerical model-based design strategy for tissue engineering

Bioreactors have been employed in tissue engineering to sustain longer and larger cell cultures, managing nutrient transfer and waste removal. Multiple designs have been developed, integrating sensor and stimulation technologies to improve cellular responses, such as proliferation and differentiation. The variability in bioreactor design, stimulation protocols, and cell culture conditions hampered comparison and replicability, possibly hiding biological evidence. This work proposes an open-source 3D printable design for a perfusion bioreactor and a numerical model-driven protocol development strategy for improved cell culture control. This bioreactor can simultaneously deliver capacitive-coupled electric field and fluid-induced shear stress stimulation, both stimulation systems were validated experimentally and in agreement with numerical predictions. A preliminary in vitro validation confirmed the suitability of the developed bioreactor to sustain viable cell cultures. The outputs from this strategy, physical and virtual, are openly available and can be used to improve comparison, replicability, and control in tissue engineering applications.

Osteocytes and Primary Cilia

PURPOSE OF REVIEW

The purpose of this review is to provide a background on osteocytes and the primary cilium, discussing the role it plays in osteocyte mechanosensing.

RECENT FINDINGS

Osteocytes are thought to be the primary mechanosensing cells in bone tissue, regulating bone adaptation in response to exercise, with the primary cilium suggested to be a key mechanosensing mechanism in bone. More recent work has suggested that, rather than being direct mechanosensors themselves, primary cilia in bone may instead form a key chemo-signalling nexus for processing mechanoregulated signalling pathways. Recent evidence suggests that pharmacologically induced lengthening of the primary cilium in osteocytes may potentiate greater mechanotransduction, rather than greater mechanosensing. While more research is required to delineate the specific osteocyte mechanobiological molecular mechanisms governed by the primary cilium, it is clear from the literature that the primary cilium has significant potential as a therapeutic target to treat mechanoregulated bone diseases, such as osteoporosis.

Multiscale morphological analysis of bone microarchitecture around Mg-10Gd implants

The utilization of biodegradable magnesium (Mg)-based implants for restoration of bone function following trauma represents a transformative approach in orthopaedic application. One such alloy, magnesium-10 weight percent gadolinium (Mg-10Gd), has been specifically developed to address the rapid degradation of Mg while enhancing its mechanical properties to promote bone healing. Previous studies have demonstrated that Mg-10Gd exhibits favorable osseointegration; however, it exhibits distinct ultrastructural adaptation in comparison to conventional implants like titanium (Ti). A crucial aspect that remains unexplored is the impact of Mg-10Gd degradation on the bone microarchitecture. To address this, we employed hierarchical three-dimensional imaging using synchrotron radiation in conjunction with image-based finite element modelling. By using the methods outlined, the vascular porosity, lacunar porosity and the lacunar-canaliculi network (LCN) morphology of bone around Mg-10Gd in comparison to Ti in a rat model from 4 weeks to 20 weeks post-implantation was investigated. Our investigation revealed that within our observation period, the degradation of Mg-10Gd implants was associated with significantly lower (p < 0.05) lacunar density in the surrounding bone, compared to Ti. Remarkably, the LCN morphology and the fluid flow analysis did not significantly differ for both implant types. In summary, a more pronounced lower lacunae distribution rather than their morphological changes was detected in the surrounding bone upon the degradation of Mg-10Gd implants. This implies potential disparities in bone remodelling rates when compared to Ti implants. Our findings shed light on the intricate relationship between Mg-10Gd degradation and bone microarchitecture, contributing to a deeper understanding of the implications for successful osseointegration.

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.

Biomechanical analysis of a lacunar-canalicular system under different cyclic displacement loading

The objective of this study is to use the finite element (FE) method to predict the mechanical signals (interstitial fluid velocity, strain, pore pressure, and pore fluid velocity) produced by osteocyte during physiological activities. The model predicts that the amplitude and distribution of the mechanical signals are mainly affected by the loading rate. The magnitude of mechanical signals in the lacunar-canalicular system increases as the amplitude, frequency and amount of direction of load increase. Collagen hillocks can effectively amplify strain signals at the process. The established model can be used for studying the mechanism of bone mechanotransduction at the micro-level.

Fluid flow shear stress and tissue remodeling-an orthodontic perspective: evidence synthesis and differential gene expression network analysis

Introduction: This study aimed to identify and analyze in vitro studies investigating the biological effect of fluid-flow shear stress (FSS) on cells found in the periodontal ligament and bone tissue. Method: We followed the PRISMA guideline for systematic reviews. A PubMed search strategy was developed, studies were selected according to predefined eligibility criteria, and the risk of bias was assessed. Relevant data related to cell source, applied FSS, and locus-specific expression were extracted. Based on this evidence synthesis and, as an original part of this work, analysis of differential gene expression using over-representation and network-analysis was performed. Five relevant publicly available gene expression datasets were analyzed using gene set enrichment analysis (GSEA). Result: A total of 6,974 articles were identified. Titles and abstracts were screened, and 218 articles were selected for full-text assessment. Finally, 120 articles were included in this study. Sample size determination and statistical analysis related to methodological quality and the ethical statement item in reporting quality were most frequently identified as high risk of bias. The analyzed studies mostly used custom-made fluid-flow apparatuses (61.7%). FSS was most frequently applied for 0.5 h, 1 h, or 2 h, whereas FSS magnitudes ranged from 6 to 20 dyn/cm2 depending on cell type and flow profile. Fluid-flow frequencies of 1 Hz in human cells and 1 and 5 Hz in mouse cells were mostly applied. FSS upregulated genes/metabolites responsible for tissue formation (AKT1, alkaline phosphatase, BGLAP, BMP2, Ca2+, COL1A1, CTNNB1, GJA1, MAPK1/MAPK3, PDPN, RUNX2, SPP1, TNFRSF11B, VEGFA, WNT3A) and inflammation (nitric oxide, PGE-2, PGI-2, PTGS1, PTGS2). Protein-protein interaction networks were constructed and analyzed using over-representation analysis and GSEA to identify shared signaling pathways. Conclusion: To our knowledge, this is the first review giving a comprehensive overview and discussion of methodological technical details regarding fluid flow application in 2D cell culture in vitro experimental conditions. Therefore, it is not only providing valuable information about cellular molecular events and their quantitative and qualitative analysis, but also confirming the reproducibility of previously published results.

Numerical simulations of fluid flow in trabecular-lacunar cavities under cyclic loading

BACKGROUND

Under external loading, the fluid shear stress (FSS) in the porous structures of bones, such as trabecular or lacunar-canalicular cavity, can influence the biological response of bone cells. However, few studies have considered both cavities. The present study investigated the characteristics of fluid flow at different scales in cancellous bone in rat femurs, as well as the effects of osteoporosis and loading frequency.

METHODS

Sprague Dawley rats (3 months old) were divided into normal and osteoporotic groups. A multiscale 3D fluid-solid coupling finite element model considering trabecular system and lacunar-canalicular system was established. Cyclic displacement loadings with frequencies of 1, 2, and 4 Hz were applied.

FINDINGS

Results showed that the wall FSS around the adhesion complexes of osteocyte on the canaliculi was higher than that on the osteocyte body. Under the same loading conditions, the wall FSS of the osteoporotic group was smaller than that of the normal group. The fluid velocity and FSS in trabecular pores exhibited a linear relationship with loading frequency. Similarly, the FSS around osteocytes also showed the loading frequency-dependent phenomenon.

INTERPRETATION

The high cadence in movement can effectively increase the FSS level on osteocytes for osteoporotic bone, i.e., expand the space within the bone with physiological load. This study might help in understanding the process of bone remodeling under cyclic loading and provide the fundamental data for the development of strategies for osteoporosis treatment.

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.

Computation of physiological loading induced interstitial fluid motion in muscle standardized femur: Healthy vs. osteoporotic bone

BACKGROUND AND OBJECTIVES

Physiological loading-induced mechanical environments regulate bone modeling and remodeling. Thus, loading-induced normal strain is typically considered a stimulus to osteogenesis. However, several studies noticed new bone formation near the sites of minimal normal strain, e.g., the neutral axis of bending in long bones, which raises a question on how bone mass is maintained near these sites. Secondary mechanical components such as shear strain and interstitial fluid flow also stimulate bone cells and regulate bone mass. However, the osteogenic potential of these components is not well established. Accordingly, the present study estimates the distribution of physiological muscle loading-induced mechanical environments such as normal strain, shear strain, pore pressure, and interstitial fluid flow in long bones.

METHODS

A poroelastic finite element muscle standardized femur (MuscleSF) model is developed to compute the distribution of the mechanical environment as a function of bone porosities associated with osteoporotic and disuse bone loss.

RESULTS

The results indicate the presence of higher shear strain and interstitial fluid motion near the minimal strain sites, i.e., the neutral axis of bending of femoral cross-sections. This suggests that secondary stimuli may maintain the bone mass at these locations. Pore pressure and interstitial fluid motion reduce with the increased porosity associated with bone disorders, possibly resulting in diminished skeletal mechano-sensitivity to exogenous loading.

CONCLUSIONS

These outcomes present a better understanding of mechanical environment-mediated regulation of site-specific bone mass, which can be beneficial in developing prophylactic exercise to prevent bone loss in osteoporosis and muscle disuse.

Computational modeling for osteogenic potential assessment of physical exercises based on loading-induced mechanobiological environments in cortical bone remodeling

Osteoporosis and disuse can cause bone loss which reduces the weight-bearing strength of long bones. Physical exercise or mechanical loading prevents bone loss as it promotes bone modeling through osteogenesis, i.e., new bone formation. Several studies have observed distinct bone remodeling responses to physical exercises; nevertheless, the underlying mechanism behind such responses is not well established. Loading-induced pore-pressure and fluid motion act as mechanobiological stimuli to bone cells namely osteocytes which further initiate osteoactivities. The shape of loading waveforms also affects the poromechanical environment of bone. Accordingly, the present study hypothesizes that loading waveforms associated with physiological exercises may expose the bone to different mechanobiological stimuli resulting in distinct bone remodeling. A poromechanical finite element model is developed to compute pore-pressure and interstitial fluid velocity in femoral cortical bone tissue (healthy and osteoporotic) subjected to loading waveforms of three physiological exercises namely walking, running, and jumping. The model also computes the mechanobiological stimulus as a function of fluid velocity. The outcomes indicate that pore-pressure and fluid velocity decrease significantly in osteoporotic bone tissue in comparison with healthy tissue. Jumping and running both improve pore-pressure and fluid velocity in healthy and osteoporotic tissues, whereas running significantly enhances mechanobiological stimulus in both the tissues which indicates a possible explanation for distinct bone remodeling to different physical exercises. The present work also suggests that running may be recommended as a potential biomechanical therapeutic to prevent bone loss. Overall, the present work contributes to the area of orthopedic research to develop effective designs of prophylactic exercises to improve bone health.

A review on design of scaffold for osteoinduction: Toward the unification of independent design variables

Biophysical stimulus quantifies the osteoinductivity of the scaffold concerning the mechanoregulatory mathematical models of scaffold-assisted cellular differentiation. Consider a set of independent structural variables ($) that comprises bulk porosity levels ([Formula: see text]) and a set of morphological features of the micro-structure ([Formula: see text]) associated with scaffolds, i.e., [Formula: see text]. The literature suggests that biophysical stimulus ([Formula: see text]) is a function of independent structural variables ($). Limited understanding of the functional correlation between biophysical stimulus and structural features results in the lack of the desired osteoinductivity in a scaffold. Consequently, it limits their broad applicability to assist bone tissue regeneration for treating critical-sized bone fractures. The literature indicates the existence of multi-dimensional independent design variable space as a probable reason for the general lack of osteoinductivity in scaffolds. For instance, known morphological features are the size, shape, orientation, continuity, and connectivity of the porous regions in the scaffold. It implies that the number of independent variables ([Formula: see text]) is more than two, i.e., [Formula: see text], which interact and influence the magnitude of [Formula: see text] in a unified manner. The efficiency of standard engineering design procedures to analyze the correlation between dependent variable ([Formula: see text]) and independent variables ($) in 3D mutually orthogonal Cartesian coordinate system diminishes proportionally with the increase in the number of independent variables ([Formula: see text]) (Deb in Optimization for engineering design-algorithms and examples, PHI Learning Private Limited, New Delhi, 2012). Therefore, there is an immediate need to devise a framework that has the potential to quantify the micro-structural's morphological features in a unified manner to increase the prospects of scaffold-assisted bone tissue regeneration.

Changes in interstitial fluid flow, mass transport and the bone cell response in microgravity and normogravity

In recent years, our scientific interest in spaceflight has grown exponentially and resulted in a thriving area of research, with hundreds of astronauts spending months of their time in space. A recent shift toward pursuing territories farther afield, aiming at near-Earth asteroids, the Moon, and Mars combined with the anticipated availability of commercial flights to space in the near future, warrants continued understanding of the human physiological processes and response mechanisms when in this extreme environment. Acute skeletal loss, more severe than any bone loss seen on Earth, has significant implications for deep space exploration, and it remains elusive as to why there is such a magnitude of difference between bone loss on Earth and loss in microgravity. The removal of gravity eliminates a critical primary mechano-stimulus, and when combined with exposure to both galactic and solar cosmic radiation, healthy human tissue function can be negatively affected. An additional effect found in microgravity, and one with limited insight, involves changes in dynamic fluid flow. Fluids provide the most fundamental way to transport chemical and biochemical elements within our bodies and apply an essential mechano-stimulus to cells. Furthermore, the cell cytoplasm is not a simple liquid, and fluid transport phenomena together with viscoelastic deformation of the cytoskeleton play key roles in cell function. In microgravity, flow behavior changes drastically, and the impact on cells within the porous system of bone and the influence of an expanding level of adiposity are not well understood. This review explores the role of interstitial fluid motion and solute transport in porous bone under two different conditions: normogravity and microgravity.

Effects of Osteocyte Shape on Fluid Flow and Fluid Shear Stress of the Loaded Bone

This study was conducted to better understand the specific behavior of the intraosseous fluid flow. We calculated the number and distribution of bone canaliculi around the osteocytes based on the varying shapes of osteocytes. We then used these calculated parameters and other bone microstructure data to estimate the anisotropy permeability of the lacunar-canalicular network. Poroelastic finite element models of the osteon were established, and the influence of the osteocyte shape on the fluid flow properties of osteons under an axial displacement load was analyzed. Two types of boundary conditions (BC) that might occur in physiological environments were considered on the cement line of the osteon. BC1 allows free fluid passage from the outer elastic restraint boundary, and BC2 is impermeable and allows no free fluid passage from outer displacement constrained boundary. They both have the same inner boundary conditions that allow fluid to pass through. Changes in the osteocyte shape altered the maximum value of pressure gradient (PG), pore pressure (PP), fluid velocity (FV), and fluid shear stress (FSS) relative to the reference model (spherical osteocytes). The maximum PG, PP, FV, and FSS in BC2 were nearly 100% larger than those in BC1, respectively. It is found that the BC1 was closer to the real physiological environment. The fluid flow along different directions in the elongated osteocyte model was more evident than that in other models, which may have been due to the large difference in permeability along different directions. Changes in osteocyte shape significantly affect the degrees of anisotropy of fluid flow and porous media of the osteon. The model presented in this study can accurately quantify fluid flow in the lacunar-canalicular network.

Interactive effects of various loading parameters on the fluid dynamics within the lacunar-canalicular system for a single osteocyte

The osteocyte lacunar-canalicular system (LCS) serves as a mechanotransductive core where external loading applied to the skeleton is transduced into mechanical signals (e.g., fluid shear) that can be sensed by mechanosensors (osteocytes). The fluid velocity and shear stress within the LCS are affected by various loading parameters. However, the interactive effect of distinct loading parameters on the velocity and shear stress in the LCS remains unclear. To address this issue, we developed a multiscale modeling approach, combining a poroelastic finite element (FE) model with a single osteocytic LCS unit model to calculate the flow velocity and shear stress within the LCS. Next, a sensitivity analysis was performed to investigate individual and interactive effects of strain magnitude, strain rate, number of cycles, and intervening short rests between loading cycles on the velocity and shear stress around the osteocyte. Lastly, we developed a relatively simple regression model to predict those outcomes. Our results demonstrated that the strain magnitude or rate alone were the main factors affecting the velocity and shear stress; however, the combination of these two was not directly additive, and addition of a short rest between cycles could enhance the combination of these two related factors. These results show highly interactive effects of distinct loading parameters on fluid velocity and shear stress in the LCS. Specifically, our results suggest that an enhanced fluid dynamics environment in the LCS can be achieved with a brief number of load cycles combined with short rest insertion and high strain magnitude and rate.

Bone permeability and mechanotransduction: Some current insights into the function of the lacunar-canalicular network

Lacunar-canalicular (LC) permeability involves the passage of fluids, nutrients, oxygen, ions, and signalling molecules through bone tissue, facilitating the maintenance of bone vitality and function and responses to various physiological conditions and diseases. LC permeability and fluid flow-shear stress/drag force play important roles in mechanotransduction in bone tissue by inducing mechanical stimuli in osteocytes, modulating cellular functions, and determining bone adaptation. Alterations in LC structure may therefore influence the fluid flow pattern through the LC network, thereby affecting the ability of osteocytes to sense and translate mechanical signals and possibly contributing to bone remodelling. Several bone-health conditions are associated with changes in LC structure and function and may affect mechanotransduction and responses, although the mechanisms underlying these associations are still not fully understood. In this review, recent studies of LC networks, their formation and transfer mechanical stimuli, and changes in structure, functional permeability, and mechanotransduction that result from age, pathology, and mechanical loading are discussed. Additionally, applications of vibration and low-intensity pulsed ultrasound in bone healthcare and regeneration fields are also presented.

Finite Element Models of Osteocytes and Their Load-Induced Activation

PURPOSE OF REVIEW

Osteocytes are the conductors of bone adaptation and remodelling. Buried inside the calcified matrix, they sense mechanical cues and signal osteoclasts in case of low activity, and osteoblasts when stresses are high. How do osteocytes detect mechanical stress? What physical signal do they perceive? Finite element analysis is a useful tool to address these questions as it allows calculating stresses, strains and fluid flow where they cannot be measured. The purpose of this review is to evaluate the capabilities and challenges of finite element models of bone, in particular the osteocytes and load-induced activation mechanisms.

RECENT FINDINGS

High-resolution imaging and increased computational power allow ever more detailed modelling of osteocytes, either in isolation or embedded within the mineralised matrix. Over the years, homogeneous models of bone and osteocytes got replaced by heterogeneous and microstructural models, including, e.g. the lacuno-canalicular network and the cytoskeleton. The lacuno-canalicular network induces strain amplifications and the osteocyte protrusions seem to be stimulated much more than the cell body, both by strain and fluid flow. More realistic cell geometries, like minute constrictions of the canaliculi, increase this effect. Microstructural osteocyte models describe the transduction of external stimuli to the nucleus. Supracellular multiscale models (e.g. of a tunnelling osteon) allow to study differential loading of osteocytes and to distinguish between strain and fluid flow as the pivotal stimulatory cue. In the future, the finite element models may be enhanced by including chemical transport and intercellular communication between osteocytes, osteoclasts and osteoblasts.

Simulation study on the effect of resistance exercise on the hydrodynamic microenvironment of osteocytes in microgravity

Osteoporosis occurs in astronauts after long-term space flight owing to the lack of gravity. The mechanical microenvironment of osteocytes in load-bearing bone are changed during resistance exercise, which prevents massive bone loss in the human body. A cylindrical fluid-structure coupling finite element model for osteons with a two-stage pore structure (i.e., Haversian canal, lacunar-canalicular system) was established with the software COMSOL. In the Earth's gravity field and in microgravity, considering the effects of pulsating pressure of arterioles, a comparative study was performed on the changes in hydrodynamic microenvironment of osteocytes during human body high-intensity exercise at different frequencies (defined as causing bone to produce 3000 με) and the body is at rest. Positive and negative liquid pressure (with respect to one atmosphere pressure) alternately acted on osteocytes during human exercising, but only positive pressure acted on osteocytes during human resting. The variation range of liquid pressure acted on osteocytes during human exercising was significantly higher than that during resting. The liquid flow velocity around osteocytes during body exercise was about four orders of magnitude higher than that during resting. In microgravity, moderate physical exercise can obviously improve the hydrodynamic microenvironment of osteocytes in load-bearing bone, which could compensate for the lack of mechanical stimulation to osteocytes caused by the lack of gravity, thereby promoting the normal physiological function of osteocytes. To a certain extent, these results revealed the biomechanical mechanism by which exercise has an effect in fighting osteoporosis in astronauts.

Computational framework for analyzing flow-induced strain on osteocyte as modulated by microenvironment

Osteocytes buried in bone matrix are major mechanosensory cells that regulate bone remodeling in response to interstitial fluid flow in a lacuno-canalicular porosity. To gain an understanding of the mechanism of osteocyte mechanosensing, it is important to be able to evaluate the local strain on the osteocyte process membrane induced by the interstitial fluid flow. The microenvironment of the osteocytes, including the pericellular matrix (PCM) and canalicular ultrastructure, is a key modulator of the flow-induced strain on the osteocyte process membrane because it produces heterogeneous flow patterns in the pericellular space. To investigate the effect of changes in the microenvironment of osteocytes on the flow-induced strain, we developed a novel computational framework for analyzing the fluid-structure interaction. Computer simulations based on the proposed framework enabled evaluation of the spatial distribution of flow-induced strain on the osteocyte process membrane according to changes in the PCM density and canalicular curvature. The simulation results reveal that a decrease in PCM density and an increase in canalicular curvature, each of which is associated with aging and bone disease, have the notable effect of enhancing local flow-induced strain on the osteocyte process membrane. We believe that the proposed computational framework is a promising framework for investigating cell-specific mechanical stimuli and that it has the potential to accelerate the mechanobiological study of osteocytes by providing a deeper understanding of their mechanical environment in living bone tissue.

Fluid-Solid Coupling Simulation of Wall Fluid Shear Stress on Cells under Gradient Fluid Flow

Fluid shear stress (FSS) plays a crucial role for cell migration within bone cavities filled with interstitial fluid. Whether the local wall FSS distribution on cell surface depends on the global gradient FSS of flow field should be clarified to explain our previous experimental observation. In this study, finite element models of discretely distributed or hexagonal closely packed cells adherent on the bottom plate in a modified plate flow chamber with different global FSS gradient were constructed. Fluid-solid coupling simulation of wall fluid shear stress on cells was performed, and two types of data analysis methods were used. The results showed that the profile of local FSS distribution on cell surface coincides with the angle of cell migration determined in the previous study, suggesting that RAW264.7 osteoclast precursors may sense the global FSS gradient and migrate toward the low-FSS region under a high gradient. For hexagonal closely packed cells, this profile on the surface of central cells decreased along with the increase of cell spacing, which may be caused by the higher local FSS difference along the direction of FSS gradient in the regions close to the bottom plate. This study may explain the phenomenon of the targeted migration of osteoclast precursors under gradient FSS field and further provide insights into the mechanism of mechanical stimulation-induced bone remodeling.

MULTI-SCALE MECHANICAL BEHAVIOR ANALYSIS ON FLUID–SOLID COUPLING FOR OSTEONS IN VARIOUS GRAVITATIONAL FIELDS

Background: Compact bone mainly consists of cylindrical osteon structures. In microgravity, the change in the mechanical microenvironment of osteocytes might be the root cause of astronauts’ bone loss during space flights.

Methods: A multi-scale three-dimensional (3D) fluid–solid coupling finite element model of osteons with a two-stage pore structure was developed using COMSOL software based on the natural structure of osteocytes. Gradients in gravitational fields of [Formula: see text]1, 0, 1, 2.5, and 3.7[Formula: see text]g were used to investigate the changes in the mechanical microenvironment on osteocyte structure. The difference in arteriole pulsating pressure and static compression stress caused by each gravity gradient was investigated.

Results: The mechanical response of osteocytes increased with the value of g, compared with the Earth’s gravitational field. For instance, the fluid pressure of osteocytes and the von Mises stress of bone matrix near lacunae decreased by 31.3% and 99.9%, respectively, in microgravity. Under static loading, only about 16.7% of osteocytes in microgravity and 58.3% of osteocytes in the Earth’s gravitational field could reach the fluid shear stress threshold of biological reactions in cell culture experiments. Compared with the Earth’s gravitational field, the pressure gradient inside osteocytes severely decreased in microgravity.

Conclusion: The mechanical microenvironment of osteocytes in microgravity might cause significant changes in the mechanical microenvironment of osteocytes, which may lead to disuse osteoporosis in astronauts.

Osteocytes and Estrogen Deficiency

PURPOSE OF REVIEW

Postmenopausal osteoporosis reduces circulating estrogen levels, which leads to osteoclast resorption, bone loss, and fracture. This review addresses emerging evidence that osteoporosis is not simply a disease of bone loss but that mechanosensitive osteocytes that regulate both osteoclasts and osteoblasts are also impacted by estrogen deficiency.

RECENT FINDINGS

At the onset of estrogen deficiency, the osteocyte mechanical environment is altered, which coincides with temporal changes in bone tissue composition. The osteocyte microenvironment is also altered, apoptosis is more prevalent, and hypermineralization occurs. The mechanobiological responses of osteocytes are impaired under estrogen deficiency, which exacerbates osteocyte paracrine regulation of osteoclasts. Recent research reveals changes in osteocytes during estrogen deficiency that may play a critical role in the etiology of the disease. A paradigm change for osteoporosis therapy requires an advanced understanding of such changes to establish the efficacy of osteocyte-targeted therapies to inhibit resorption and secondary mineralization.

Control of osteoblast arrangement by osteocyte mechanoresponse through prostaglandin E2 signaling under oscillatory fluid flow stimuli

Anisotropic collagen/apatite microstructure is a prominent determinant of bone tissue functionalization; in particular, bone matrix modulates its anisotropic microstructure depending on the surrounding mechanical condition. Although mechanotransduction in bones is governed by osteocyte function, the precise mechanisms linking mechanical stimuli and anisotropic formation of collagen/apatite microstructure are poorly understood. Here we developed a novel anisotropic mechano-coculture system which enables the understanding of the biological mechanisms regulating the oriented bone matrix formation, which is constructed by aligned osteoblasts. The developed model provides bone-mimetic coculture platform that enables simultaneous control of mechanical condition and osteoblast-osteocyte communication with an anisotropic culture scaffold. The engineered coculture device helps in understanding the relationship between osteocyte mechanoresponses and osteoblast arrangement, which is a significant contributor to anisotropic organization of bone tissue. Our study showed that osteocyte responses to oscillatory flow stimuli regulated osteoblast arrangement through soluble molecular interactions. Importantly, we found that prostaglandin E2 is a novel determinant for oriented collagen/apatite organization of bone matrix, through controlling osteoblast arrangement.

High-resolution image-based simulation reveals membrane strain concentration on osteocyte processes caused by tethering elements

Osteocytes are vital for regulating bone remodeling by sensing the flow-induced mechanical stimuli applied to their cell processes. In this mechanosensing mechanism, tethering elements (TEs) connecting the osteocyte process with the canalicular wall potentially amplify the strain on the osteocyte processes. The ultrastructure of the osteocyte processes and canaliculi can be visualized at a nanometer scale using high-resolution imaging via ultra-high voltage electron microscopy (UHVEM). Moreover, the irregular shapes of the osteocyte processes and the canaliculi, including the TEs in the canalicular space, should considerably influence the mechanical stimuli applied to the osteocytes. This study aims to characterize the roles of the ultrastructure of osteocyte processes and canaliculi in the mechanism of osteocyte mechanosensing. Thus, we constructed a high-resolution image-based model of an osteocyte process and a canaliculus using UHVEM tomography and investigated the distribution and magnitude of flow-induced local strain on the osteocyte process by performing fluid–structure interaction simulation. The analysis results reveal that local strain concentration in the osteocyte process was induced by a small number of TEs with high tension, which were inclined depending on the irregular shapes of osteocyte processes and canaliculi. Therefore, this study could provide meaningful insights into the effect of ultrastructure of osteocyte processes and canaliculi on the osteocyte mechanosensing mechanism.

Physiological Loading-Induced Interstitial Fluid Dynamics in Osteon of Osteogenesis Imperfecta Bone

Osteogenesis imperfecta (OI), also known as "brittle bone disease," is a genetic bone disorder. OI bones experience frequent fractures. Surgical procedures are usually followed by clinicians in the management of OI. It has been observed physical activity is equally beneficial in reducing OI bone fractures in both children and adults as mechanical stimulation improves bone mass and strength. Loading-induced mechanical strain and interstitial fluid flow stimulate bone remodeling activities. Several studies have characterized strain environment in OI bones, whereas very few studies attempted to characterize the interstitial fluid flow. OI significantly affects bone micro-architecture. Thus, this study anticipates that canalicular fluid flow reduces in OI bone in comparison to the healthy bone in response to physiological loading due to altered poromechanical properties. This work attempts to understand the canalicular fluid distribution in single osteon models of OI and healthy bone. A poromechanical model of osteon is developed to compute pore-pressure and interstitial fluid flow as a function of gait loading pattern reported for OI and healthy subjects. Fluid distribution patterns are compared at different time-points of the stance phase of the gait cycle. It is observed that fluid flow significantly reduces in OI bone. Additionally, flow is more static than dynamic in OI osteon in comparison to healthy subjects. This work attempts to identify the plausible explanation behind the diminished mechanotransduction capability of OI bone. This work may further be extended for designing better biomechanical therapies to enhance the fluid flow in order to improve osteogenic activities in OI bone.

Plasma membrane disruption (PMD) formation and repair in mechanosensitive tissues

Mammalian cells employ an array of biological mechanisms to detect and respond to mechanical loading in their environment. One such mechanism is the formation of plasma membrane disruptions (PMD), which foster a molecular flux across cell membranes that promotes tissue adaptation. Repair of PMD through an orchestrated activity of molecular machinery is critical for cell survival, and the rate of PMD repair can affect downstream cellular signaling. PMD have been observed to influence the mechanical behavior of skin, alveolar, and gut epithelial cells, aortic endothelial cells, corneal keratocytes and epithelial cells, cardiac and skeletal muscle myocytes, neurons, and most recently, bone cells including osteoblasts, periodontal ligament cells, and osteocytes. PMD are therefore positioned to affect the physiological behavior of a wide range of vertebrate organ systems including skeletal and cardiac muscle, skin, eyes, the gastrointestinal tract, the vasculature, the respiratory system, and the skeleton. The purpose of this review is to describe the processes of PMD formation and repair across these mechanosensitive tissues, with a particular emphasis on comparing and contrasting repair mechanisms and downstream signaling to better understand the role of PMD in skeletal mechanobiology. The implications of PMD-related mechanisms for disease and potential therapeutic applications are also explored.

Disrupted osteocyte connectivity and pericellular fluid flow in bone with aging and defective TGF-β signaling

Skeletal fragility in the elderly does not simply result from a loss of bone mass. However, the mechanisms underlying the concurrent decline in bone mass, quality, and mechanosensitivity with age remain unclear. The important role of osteocytes in these processes and the age-related degeneration of the intricate lacunocanalicular network (LCN) in which osteocytes reside point to a primary role for osteocytes in bone aging. Since LCN complexity severely limits experimental dissection of these mechanisms in vivo, we used two in silico approaches to test the hypothesis that LCN degeneration, due to aging or an osteocyte-intrinsic defect in transforming growth factor beta (TGF-β) signaling (TβRIIocy-/-), is sufficient to compromise essential osteocyte responsibilities of mass transport and exposure to mechanical stimuli. Using reconstructed confocal images of bone with fluorescently labeled osteocytes, we found that osteocytes from aged and TβRIIocy-/- mice had 33 to 45% fewer, and more tortuous, canaliculi. Connectomic network analysis revealed that diminished canalicular density is sufficient to impair diffusion even with intact osteocyte numbers and overall LCN architecture. Computational fluid dynamics predicts that the corresponding drop in shear stress experienced by aged or TβRIIocy-/- osteocytes is highly sensitive to canalicular surface area but not tortuosity. Simulated expansion of the osteocyte pericellular space to mimic osteocyte perilacunar/canalicular remodeling restored predicted shear stress for aged osteocytes to young levels. Overall, these models show how loss of LCN volume through LCN pruning may lead to impaired fluid dynamics and osteocyte exposure to mechanostimulation. Furthermore, osteocytes emerge as targets of age-related therapeutic efforts to restore bone health and function.

Closing the osteon: Do osteocytes sense strain rate rather than fluid flow?

Osteons are cylindrical structures of bone created by matrix resorbing osteoclasts, followed by osteoblasts that deposit new bone. Osteons align with the principal loading direction and it is thought that the osteoclasts are directed by osteocytes, the mechanosensitive cells that reside inside the bone matrix. These osteocytes are presumably controlled by interstitial fluid flow, induced by the physiological loading of bones. Here I consider the stimulation of osteocytes while the osteon is closed by osteoblasts. In a conceptual finite element model, bone is considered a poro-elastic material and subjected to locomotion-induced loading conditions. It appears that the magnitude of flow is constant along the closing cone, while shear strain rate in the bone matrix diminishes linearly with the deposition of bone. This suggests that shear strain rate, rather than fluid flow, is the physical cue that controls osteocytes and bone deposition in newly formed osteons.

Mechanical Stimulation Modulates Osteocyte Regulation of Cancer Cell Phenotype

Breast and prostate cancers preferentially metastasise to bone tissue, with metastatic lesions forming in the skeletons of most patients. On arriving in bone tissue, disseminated tumour cells enter a mechanical microenvironment that is substantially different to that of the primary tumour and is largely regulated by bone cells. Osteocytes, the most ubiquitous bone cell type, orchestrate healthy bone remodelling in response to physical exercise. However, the effects of mechanical loading of osteocytes on cancer cell behaviour is still poorly understood. The aim of this study was to characterise the effects of osteocyte mechanical stimulation on the behaviour of breast and prostate cancer cells. To replicate an osteocyte-controlled environment, this study treated breast (MDA-MB-231 and MCF-7) and prostate (PC-3 and LNCaP) cancer cell lines with conditioned media from MLO-Y4 osteocyte-like cells exposed to mechanical stimulation in the form of fluid shear stress. We found that osteocyte paracrine signalling acted to inhibit metastatic breast and prostate tumour growth, characterised by reduced proliferation and invasion and increased migration. In breast cancer cells, these effects were largely reversed by mechanical stimulation of osteocytes. In contrast, conditioned media from mechanically stimulated osteocytes had no effect on prostate cancer cells. To further investigate these interactions, we developed a microfluidic organ-chip model using the Emulate platform. This new organ-chip model enabled analysis of cancer cell migration, proliferation and invasion in the presence of mechanical stimulation of osteocytes by fluid shear stress, resulting in increased invasion of breast and prostate cancer cells. These findings demonstrate the importance of osteocytes and mechanical loading in regulating cancer cell behaviour and the need to incorporate these factors into predictive in vitro models of bone metastasis.

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.