Numerical simulation of osteocyte cell in response to directional mechanical loadings and mechanotransduction analysis: Considering lacunar-canalicular interstitial fluid flow

Combining anabolic loading and raloxifene improves bone quantity and some quality measures in a mouse model of osteogenesis imperfecta

Osteogenesis imperfecta (OI) increases fracture risk due to changes in bone quantity and quality caused by mutations in collagen and its processing proteins. Current therapeutics improve bone quantity, but do not treat the underlying quality deficiencies. Male and female G610C+/- mice, a murine model of OI, were treated with a combination of raloxifene and in vivo axial tibial compressive loading starting at 10 weeks of age and continuing for 6 weeks to improve bone quantity and quality. Bone geometry and mechanical properties were measured to determine whole bone and tissue-level material properties. A colocalized Raman/nanoindentation system was used to measure chemical composition and nanomechanical properties in newly formed bone compared to old bone to determine if bone formed during the treatment regimen differed in quality compared to bone formed prior to treatment. Lastly, lacunar geometry and osteocyte apoptosis were assessed. OI mice were able to build bone in response to the loading, but this response was less robust than in control mice. Raloxifene improved some bone material properties in female but not male OI mice. Raloxifene did not alter nanomechanical properties, but loading did. Lacunar geometry was largely unchanged with raloxifene and loading. However, osteocyte apoptosis was increased with loading in raloxifene treated female mice. Overall, combination treatment with raloxifene and loading resulted in positive but subtle changes to bone quality.

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.

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.

Modeling and simulation of bone cells dynamic behavior under the late effect of breast cancer treatments

Breast Cancer (BC) treatments have been proven to interfere with the health of bones. Chemotherapy and endocrinal treatment regimens such as tamoxifen and aromatase inhibitors are frequently prescribed for women with BC. However, these drugs increase bone resorption and reduce the Bone Mineral Density (BMD), thus increasing the risk of bone fracture. In the current study, a mechanobiological bone remodeling model has been developed by coupling cellular activities, mechanical stimuli, and the effect of breast cancer treatments (chemotherapy, tamoxifen, and aromatase inhibitors). This model algorithm has been programmed and implemented on MATLAB software to simulate different treatment scenarios and their effects on bone remodeling and also predict the evolution of Bone Volume fraction (BV/TV) and the associated Bone Density Loss (BDL) over a period of time. The simulation results, achieved from different combinations of Breast Cancer treatments, allow the researchers to predict the intensity of each combination treatment on BV/TV and BMD. The combination of chemotherapy, tamoxifen, and aromatase inhibitors, followed by the combination of chemotherapy and tamoxifen remain the most harmful regimen. This is because they have a strong ability to induce the bone degradation which is represented by a decrease of 13.55% and 11.55% of the BV/TV value, respectively. These results were compared with the experimental studies and clinical observations which showed good agreement. The proposed model can be used by clinicians and physicians to choose the most appropriate combination of treatments, according to the patient's case.

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.

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.

Cumulative Effects of Strontium Ranelate and Impact Exercise on Bone Mass in Ovariectomized Rats

OBJECTIVE

To explore the effect of physical exercise (EXE), strontium ranelate (SR), or their combination on bone status in ovariectomized (OVX) rats.

DESIGN

Sixty female Wistar rats were randomized to one of five groups: sham (Sh), OVX (O), OVX+EXE (OE), OVX+SR (OSR), and OVX+EXE+SR (OESR). Animals in EXE groups were subjected to 10 drops per day (45 cm in height); rats in SR groups received 625 mg/kg/day of SR, 5 days/week for 8 weeks. Bone mineral density (BMD) and bone mineral content (BMC, dual-energy X-ray absorptiometry (DXA)), mechanical strength of the left femur (three-point bending test), and femur microarchitecture of (micro-computed tomography imaging, microCT) analyses were performed to characterize biomechanical and trabecular/cortical structure. Bone remodeling, osteocyte apoptosis, and lipid content were evaluated by ELISA and immunofluorescence tests.

RESULTS

In OVX rats, whole-body BMD, trabecular parameters, and osteocalcin (OCN) levels decreased, while weight, lean/fat mass, osteocyte apoptosis, and lipid content all increased. EXE after ovariectomy improved BMD and BMC, trabecular parameters, cross-sectional area (CSA), moment of inertia, and OCN levels while decreasing osteocyte apoptosis and lipid content. SR treatment increased BMD and BMC, trabecular parameters, CSA, stiffness, OCN, and alkaline phosphatase (ALP) levels. Furthermore, fat mass, N-telopeptide (NTX) level, osteocyte apoptosis, and lipid content significantly decreased. The combination of both EXE and SR improved bone parameters compared with EXE or SR alone.

CONCLUSION

EXE and SR had positive and synergistic effects on bone formation and resorption.

A modular programmed biphasic dual-delivery system on 3D ceramic scaffolds for osteogenesis in vitro and in vivo

Single-factor delivery is the most common characteristic of bone tissue engineering techniques. However, bone regeneration is a complex process requiring multiple factors and specialized release mechanisms. Therefore, the development of a dual-delivery system allowing for programmed release kinetics would be highly desirable. Improvement of the molarity and versatility of the delivery system has rarely been studied. Herein, we report the development of a novel, modular programmed biphasic dual-release system (SCB), carrying a BMP2 and an engineered collagen I-derived recognition motif (Stath-DGEA), with a self-remodification feature on hydroxyapatite (HA)-based materials. The SCB system was loaded onto an additive manufactured (AM) scaffold in order to evaluate its bifactor osteogenic potential and its biphasic release behavior. Further, the biomechanical properties of the scaffold were studied by using the fluid-structure interaction (FSI) method. Section fluorescent labeling revealed that the HA scaffold has a relatively higher density and efficiency. Additionally, the results of the release and inhibition experiment suggested that the SCB system could facilitate the sustained release of therapeutic levels of two factors during the initial stage of implantation, thereby exhibiting a rapid high-dose release pattern at a specific time point during the second stage. The FSI prediction model indicated that the scaffold provides an excellent biomimetic mechanical and fluid dynamic microenvironment to promote osteogenesis. Our results indicated that incorporation of BMP2 with Stath-DGEA in the biphasic SCB system could have a synergetic effect in promoting the adhesion, proliferation, and differentiation of bone marrow mesenchymal stem cells (BMSCs) in vitro, under staged stimulations. Further, in vivo studies in both ectopic and orthotopic rat models showed that the SCB system loaded onto an AM scaffold could enhance osteointegration and osteoinduction throughout the osteogenic process. Thus, the novel synthetic SCB system described herein used on an AM scaffold provides a biomimetic extracellular environment that enhances bone regeneration and is a promising multifunctional, dual-release platform.

Multiscale finite element modeling of mechanical strains and fluid flow in osteocyte lacunocanalicular system

Osteocytes form over 90% of the bone cells and are postulated to be mechanosensors responsible for regulating the function of osteoclasts and osteoblasts in bone modeling and remodeling. Physical activity results in mechanical loading on the bones. Osteocytes are thought to be the main mechanosensory cells in bone. Upon load osteocytes secrete key factors initiating downstream signaling pathways that regulate skeletal metabolism including the Wnt/β-catenin signaling pathway. Osteocytes have dendritic structures and are housed in the lacunae and canaliculi within the bone matrix. Mechanical loading is known to have two primary effects, namely a mechanical strain (membrane disruption by stretching) on the lacunae/cells, and fluid flow, in the form of fluid flow shear stress (FFSS), in the space between the cell membranes and the lacuna-canalicular walls. In response, osteocytes get activated via a process called mechanotransduction in which mechanical signals are transduced to biological responses. The study of mechanotransduction is a complex subject involving principles of engineering mechanics as well as biological signaling pathway studies. Several length scales are involved as the mechanical loading on macro sized bones are converted to strain and FFSS responses at the micro-cellular level. Experimental measurements of strain and FFSS at the cellular level are very difficult and correlating them to specific biological activity makes this a very challenging task. One of the methods commonly adopted is a multi-scale approach that combines biological and mechanical experimentation with in silico numerical modeling of the engineering aspects of the problem. Finite element analysis along with fluid-structure interaction methodologies are used to compute the mechanical strain and FFSS. These types of analyses often involve a multi-length scale approach where models of both the macro bone structure and micro structure at the cellular length scale are used. Imaging modalities play a crucial role in the development of the models and present their own challenges. This paper reviews the efforts of various research groups in addressing this problem and presents the work in our research group. A clear understanding of how mechanical stimuli affect the lacunae and perilacunar tissue strains and shear stresses on the cellular membranes may ultimately lead to a better understanding of the process of osteocyte activation.

Osteocyte lacunar strain determination using multiscale finite element analysis

Osteocytes are thought to be the primary mechanosensory cells within bone, regulating both osteoclasts and osteoblasts to control load induced changes in bone resorption and formation. Osteocytes initiate intracellular responses including activating the Wnt/β-catenin signaling pathway after experiencing mechanical forces. In response to changing mechanical loads (strain) the osteocytes signal to cells on the bone surface. However, this process of osteocyte activation appears heterogeneous since it occurs in sub-populations of osteocytes, even within regions predicted to be experiencing similar global strain magnitudes determined based on traditional finite element modeling approaches. Several studies have investigated the strain responses of osteocyte lacunae using finite element (FE) models, but many were limited by the use of idealized geometries (e.g., ellipsoids) and analysis of a single osteocyte. Finite element models by other groups included more details, such as canaliculi, but all were done on models consisting of a single osteocyte. We hypothesized that variation in size and orientation of the osteocyte lacunae within bone would give rise to micro heterogeneity in the strain fields that could better explain the observed patterns of osteocyte activation following load. The osteocytes in our microscale and nanoscale models have an idealized oval shape and some are based on confocal scans. However, all the FE models in this preliminary study consist of multiple osteocytes. The number of osteocytes in the 3D confocal scan models ranged from five to seventeen. In this study, a multi-scale computational approach was used to first create an osteocyte FE model at the microscale level to examine both the theoretical lacunar and perilacunar strain responses based on two parameters: 1) lacunar orientation and 2) lacunar size. A parametric analysis was performed by steadily increasing the perilacunar modulus (5, 10, 15, and 20 GPa). Secondly, a nanoscale FE model was built using known osteocyte dimensions to determine the predicted strains in the perilacunar matrix, fluid space, and cell body regions. Finally, 3-D lacunar models were created using confocal image stacks from mouse femurs to determine the theoretical strain in the lacunae represented by realistic geometries. Overall, lacunar strains decreased by 14% in the cell body, 15% in the fluid space region and 25% in the perilacunar space as the perilacunar modulus increased, indicating a stress shielding effect. Lacunar strains were lower for the osteocytes aligned along the loading axis compared to those aligned perpendicular to axis. Increases in lacuna size also led to increased lacunar strains. These finite element model findings suggest that orientation and lacunar size may contribute to the heterogeneous initial pattern of osteocyte strain response observed in bone following in vivo applied mechanical loads. A better understanding of how mechanical stimuli directly affect the lacunae and perilacunar tissue strains may ultimately lead to a better understanding of the process of osteocyte activation in response to mechanical loading.

•

A multi-scale computational approach used to first create multiple osteocyte FE model at the microscale level

•

3-D Lacuna model created using confocal image stacks from a mouse femur to determine the theoretical strain in the lacunae.

•

Lacunar strains decreased as the perilacunar modulus increased.

•

Lacunar strains were lower for the osteocytes aligned along the loading axis compared to those aligned perpendicular to axis.

•

Increases in lacuna size also led to increased lacunar strains

Control of Bone Matrix Properties by Osteocytes

Osteocytes make up 90-95% of the cellular content of bone and form a rich dendritic network with a vastly greater surface area than either osteoblasts or osteoclasts. Osteocytes are well positioned to play a role in bone homeostasis by interacting directly with the matrix; however, the ability for these cells to modify bone matrix remains incompletely understood. With techniques for examining the nano- and microstructure of bone matrix components including hydroxyapatite and type I collagen becoming more widespread, there is great potential to uncover novel roles for the osteocyte in maintaining bone quality. In this review, we begin with an overview of osteocyte biology and the lacunar-canalicular system. Next, we describe recent findings from in vitro models of osteocytes, focusing on the transitions in cellular phenotype as they mature. Finally, we describe historical and current research on matrix alteration by osteocytes in vivo, focusing on the exciting potential for osteocytes to directly form, degrade, and modify the mineral and collagen in their surrounding matrix.