Modeling of biological doses and mechanical effects on bone transduction

Shear stress, hormones like parathyroid and mineral elements like calcium mediate the amplitude of stimulus signal, which affects the rate of bone remodeling. The current study investigates the theoretical effects of different metabolic doses in stimulus signal level on bone. The model was built considering the osteocyte as the sensing center mediated by coupled mechanical shear stress and some biological factors. The proposed enhanced model was developed based on previously published works dealing with different aspects of bone transduction. It describes the effects of physiological doses variations of calcium, parathyroid hormone, nitric oxide and prostaglandin E2 on the stimulus level sensed by osteocytes in response to applied shear stress generated by interstitial fluid flow. We retained the metabolic factors (parathyroid hormone, nitric oxide and prostaglandin E2) as parameters of bone cell mechanosensitivity because stimulation/inhibition of induced pathways stimulates osteogenic response in vivo. We then tested the model response in terms of stimulus signal variation versus the biological factors doses to external mechanical stimuli. Despite the limitations of the model, it is consistent and has physiological bases. Biological inputs are histologically measurable. This makes the model amenable to experimental verification.     

Osteocyte calcium signaling response to bone matrix deformation

Osteocytes embedded in calcified bone matrix have been widely believed to play important roles in mechanosensing to achieve adaptive bone remodeling in a changing mechanical environment. In vitro studies have clarified several types of mechanical stimuli such as hydrostatic pressure, fluid shear stress, and direct deformation influence osteocyte functions. However, osteocyte response to mechanical stimuli in the bone matrix has not been clearly understood. In this study, we observed the osteocyte calcium signaling response to the quantitatively applied deformation in the bone matrix. A novel experimental system was developed to apply deformation to cultured bone tissue with osteocytes on a microscope stage. As a mechanical stimulus to the osteocytes in bone matrix, in-plane shear deformation was applied using a pair of glass microneedles to bone fragments, obtained from 13-day-old embryonic chick calvariae. Deformation of bone matrix and cells was quantitatively evaluated using an image correlation method by applying for differential interference contrast images of the matrix and fluorescent images of immunolabeled osteocytes, together with imaging of the cellular calcium transient using a ratiometric method. As a result, it was confirmed that the newly developed system enables us to apply deformation to bone matrix and osteocytes successfully under the microscope without significant focal plane shift or deviation from the observation view field. The system could be a basis for further development to investigate the mechanosensing mechanism of osteocytes in bone matrix through examination of various types of rapid biochemical signaling responses and intercellular communication induced by matrix deformation.     

Decreased pericellular matrix production and selection for enhanced cell membrane repair may impair osteocyte responses to mechanical loading in the aging skeleton

Transient plasma membrane disruptions (PMD) occur in osteocytes with in vitro and in vivo loading, initiating mechanotransduction. The goal here was to determine whether osteocyte PMD formation or repair is affected by aging. Osteocytes from old (24 months) mice developed fewer PMD (?76% females, ?54% males) from fluid shear than young (3 months) mice, and old mice developed fewer osteocyte PMD (?51%) during treadmill running. This was due at least in part to decreased pericellular matrix production, as studies revealed that pericellular matrix is integral to formation of osteocyte PMD, and aged osteocytes produced less pericellular matrix (?55%). Surprisingly, osteocyte PMD repair rate was faster (+25% females, +26% males) in osteocytes from old mice, and calcium wave propagation to adjacent nonwounded osteocytes was blunted, consistent with impaired mechanotransduction downstream of PMD in osteocytes with fast PMD repair in previous studies. Inducing PMD via fluid flow in young osteocytes in the presence of oxidative stress decreased postwounding cell survival and promoted accelerated PMD repair in surviving cells, suggesting selective loss of slower‐repairing osteocytes. Therefore, as oxidative stress increases during aging, slower‐repairing osteocytes may be unable to successfully repair PMD, leading to slower‐repairing osteocyte death in favor of faster‐repairing osteocyte survival. Since PMD are an important initiator of mechanotransduction, age‐related decreases in pericellular matrix and loss of slower‐repairing osteocytes may impair the ability of bone to properly respond to mechanical loading with bone formation. These data suggest that PMD formation and repair mechanisms represent new targets for improving bone mechanosensitivity with aging.     

Model of cellular mechanotransduction via actin stress fibers

Mechanical stresses due to blood flow regulate vascular endothelial cell structure and function and play a key role in arterial physiology and pathology. In particular, the development of atherosclerosis has been shown to correlate with regions of disturbed blood flow where endothelial cells are round and have a randomly organized cytoskeleton. Thus, deciphering the relation between the mechanical environment, cell structure, and cell function is a key step toward understanding the early development of atherosclerosis. Recent experiments have demonstrated very rapid (\(\sim \)100 ms) and long-distance (\(\sim \)10 \(\upmu \)m) cellular mechanotransduction in which prestressed actin stress fibers play a critical role. Here, we develop a model of mechanical signal transmission within a cell by describing strains in a network of prestressed viscoelastic stress fibers following the application of a force to the cell surface. We find force transmission dynamics that are consistent with experimental results. We also show that the extent of stress fiber alignment and the direction of the applied force relative to this alignment are key determinants of the efficiency of mechanical signal transmission. These results are consistent with the link observed experimentally between cytoskeletal organization, mechanical stress, and cellular responsiveness to stress. Based on these results, we suggest that mechanical strain of actin stress fibers under force constitutes a key link in the mechanotransduction chain.     

Mechanosensation and fluid transport in living bone

The mechanosensory mechanisms in bone include (i) the cell system that is stimulated by external mechanical loading applied to the bone; (ii) the system that transduces that mechanical loading to a communicable signal; and (iii) the systems that transmit that signal to the effector cells for the maintenance of bone homeostasis and for strain adaptation of the bone structure. The effector cells are the osteoblasts and the osteoclasts. These systems and the mechanisms that they employ have not yet been unambiguously identified. The candidate systems will be reviewed. It will be argued that the current theoretical and experimental evidence suggests that osteocytes are the principal mechanosensory cells of bone, that they are activated by shear stress from fluid flowing through the osteocyte canaliculi, and that the electrically coupled three-dimensional network of osteocytes and lining cells is a communications system for the control of bone homeostasis and structural strain adaptation. The movement of bone fluid from the region of the bone vasculature through the canaliculi and the lacunae of the surrounding mineralized tissue accomplishes three important tasks. First, it transports nutrients to the osteocytes in the lacunae buried in the mineralized matrix. Second, it carries away the cell waste. Third, the bone fluid exerts a force on the cell process, a force that is large enough for the cell to sense. This is probably the basic mechanotransduction mechanism in bone, the way in which bone senses the mechanical load to which it is subjected. The mechanisms of bone fluid flow are described with particular emphasis on mechanotransduction. Also described is the cell to cell communication by which higher frequency signals might be transferred, a potential mechanism in bone by which the small whole tissue strain is amplified so the bone cells can respond to it. One of the conclusions is that higher frequency low amplitude strains can maintain bone as effectively as low frequency high amplitude strains. This conclusion leads to a paradigm shift in how to treat osteoporosis and how to cope with microgravity.     

Osteoblasts infill irregular pores under curvature and porosity controls: a hypothesis-testing analysis of cell behaviours

The geometric control of bone tissue growth plays a significant role in bone remodelling, age-related bone loss, and tissue engineering. However, how exactly geometry influences the behaviour of bone-forming cells remains elusive. Geometry modulates cell populations collectively through the evolving space available to the cells, but it may also modulate the individual behaviours of cells. To factor out the collective influence of geometry and gain access to the geometric regulation of individual cell behaviours, we develop a mathematical model of the infilling of cortical bone pores and use it with available experimental data on cortical infilling rates. Testing different possible modes of geometric controls of individual cell behaviours consistent with the experimental data, we find that efficient smoothing of irregular pores only occurs when cell secretory rate is controlled by porosity rather than curvature. This porosity control suggests the convergence of a large scale of intercellular signalling to single bone-forming cells, consistent with that provided by the osteocyte network in response to mechanical stimulus. After validating the mathematical model with the histological record of a real cortical pore infilling, we explore the infilling of a population of randomly generated initial pore shapes. We find that amongst all the geometric regulations considered, the collective influence of curvature on cell crowding is a dominant factor for how fast cortical bone pores infill, and we suggest that the irregularity of cement lines thereby explains some of the variability in double labelling data as well as the overall speed of osteon infilling.     

Directional dependence of osteoblastic calcium response to mechanical stimuli

In adaptive bone remodeling, mechanical signals such as stress/strain caused by loading/deformation are believed to play important roles as regulators of the process in which osteoclastic resorption and osteoblastic formation are coordinated under a local mechanical environment. The mechanism by which cells sense and transduce mechanical signals to the intracellular biochemical signaling cascade is still unclear, however to address this issue, the present study investigated the characteristic response of a single osteoblastic cell, MC3T3-E1, to a well-defined mechanical stimulus and the involvement of the cytoskeletal actin fiber structure in the mechanotransduction pathway. First, by mechanically perturbing to a single cell using a microneedle, a change in the intracellular calcium ion concentration [Ca²⁺]_i was observed as a primal signaling response to a mechanical stimulus, and the threshold value of the perturbation as the mechanical stimulus was evaluated quantitatively. Second, to study directional dependence of the response to the mechanical stimulus, the effect of actin fiber orientation on the threshold value of the calcium response was investigated at various magnitudes and directions of the stimulus. It was found that the osteoblastic response to the perturbation exhibited a directional dependence. That is, the sensitivity of osteoblastic cells to a mechanical stimulus depends on the angle of the applied deformation with respect to the cytoskeletal actin fiber orientation. This finding is phenomenological evidence that cytoskeletal actin fiber structures are involved in the mechanotransduction mechanism, which may be related to cell polarization behaviors such as cellular alignment caused by mechanical stimulation.     

The effect of cytoskeletal disruption on pulsatile fluid flow-induced nitric oxide and prostaglandin E2 release in osteocytes and osteoblasts

Fluid flowing through the bone porosity might be a primary stimulus for functional adaptation of bone. Osteoblasts, and osteocytes in particular, respond to fluid flow in vitro with enhanced nitric oxide (NO) and prostaglandin E(2) (PGE(2)) release; both of these signaling molecules mediate mechanically-induced bone formation. Because the cell cytoskeleton is involved in signal transduction, we hypothesized that the pulsatile fluid flow-induced release of NO and PGE(2) in both osteoblastic and osteocytic cells involves the actin and microtubule cytoskeleton. In testing this hypothesis we found that fluid flow-induced NO response in osteoblasts was accompanied by parallel alignment of stress fibers, whereas PGE(2) response was related to fluid flow stimulation of focal adhesions formed after cytoskeletal disruption. Fluid flow-induced PGE(2) response in osteocytes was inhibited by cytoskeletal disruption, whereas in osteoblasts it was enhanced. These opposite PGE(2) responses are likely related to differences in cytoskeletal composition (osteocyte structure was more dependent on actin), but may occur via cytoskeletal modulation of shear/stretch-sensitive ion channels that are known to be dominant in osteocyte (and not osteoblast) response to mechanical loading.     

Osteocyte network; a negative regulatory system for bone mass augmented by the induction of Rankl in osteoblasts and Sost in osteocytes at unloading

Reduced mechanical stress is a major cause of osteoporosis in the elderly, and the osteocyte network, which comprises a communication system through processes and canaliculi throughout bone, is thought to be a mechanosensor and mechanotransduction system; however, the functions of osteocytes are still controversial and remain to be clarified. Unexpectedly, we found that overexpression of BCL2 in osteoblasts eventually caused osteocyte apoptosis. Osteoblast and osteoclast differentiation were unaffected by BCL2 transgene in vitro. However, the cortical bone mass increased due to enhanced osteoblast function and suppressed osteoclastogenesis at 4 months of age, when the frequency of TUNEL-positive lacunae reached 75%. In the unloaded condition, the trabecular bone mass decreased in both wild-type and BCL2 transgenic mice at 6 weeks of age, while it decreased due to impaired osteoblast function and enhanced osteoclastogenesis in wild-type mice but not in BCL2 transgenic mice at 4 months of age. Rankl and Opg were highly expressed in osteocytes, but Rankl expression in osteoblasts but not in osteocytes was increased at unloading in wild-type mice but not in BCL2 transgenic mice at 4 months of age. Sost was locally induced at unloading in wild-type mice but not in BCL2 transgenic mice, and the dissemination of Sost was severely interrupted in BCL2 transgenic mice, showing the severely impaired osteocyte network. These findings indicate that the osteocyte network is required for the upregulation of Rankl in osteoblasts and Sost in osteocytes in the unloaded condition. These findings suggest that the osteocyte network negatively regulate bone mass by inhibiting osteoblast function and activating osteoclastogenesis, and these functions are augmented in the unloaded condition at least partly through the upregulation of Rankl expression in osteoblasts and that of Sost in osteocytes, although it cannot be excluded that low BCL2 transgene expression in osteoblasts contributed to the enhanced osteoblast function.     

Functions of the osteocyte network in the regulation of bone mass

Osteocytes establish an extensive intracellular and extracellular communication system via gap-junction-coupled cell processes and canaliculi throughout bone and the communication system is extended to osteoblasts on the bone surface. The osteocyte network is an ideal mechanosensory system and suitable for mechanotransduction. However, the overall function of the osteocyte network remains to be clarified, since bone resorption is enhanced by osteocyte apoptosis, which is followed by a process of secondary necrosis attributable to the lack of scavengers. The enhanced bone resorption is caused by the release of intracellular content, including immunostimulatory molecules that activate osteoclastogenesis through the canaliculi. Therefore, a mouse model is required in which the osteocyte network is disrupted but in which no bone resorption is induced, in order to evaluate the overall functions of the osteocyte network. One such model is the BCL2 transgenic mouse, in which the osteocyte network, including both intracellular and extracellular networks, is disrupted. Another model is the osteocyte-specific Gja1 knockout mouse, in which intercellular communication through gap junctions is impaired but the canalicular system is intact. Combining the findings from these mouse models with previous histological observations showing the inverse linkage between osteocyte density and bone formation, we conclude that the osteocyte network enhances bone resorption and inhibits bone formation under physiological conditions. Further, studies with BCL2 transgenic mice show that these osteocyte functions are augmented in the unloaded condition. In this condition, Rankl upregulation in osteoblasts and Sost upregulation in osteocytes are, at least in part, responsible for enhanced bone resorption and suppressed bone formation, respectively.     

Bone remodelling algorithms incorporating both strain and microdamage stimuli

Biomechanical theories to predict bone remodelling have used either mechanical strain or microdamage as the stimulus driving cellular responses. Even though experimental data have implicated both stimuli in bone cell regulation, a mechano-regulatory system incorporating both stimuli has not yet been proposed. In this paper, we test the hypothesis that bone remodelling may be regulated by signals due to both strain and microdamage. Four mechano-regulation algorithms are studied where the stimulus is: strain, damage, combined strain/damage, and either strain or damage with damage-adaptive remodelling prioritised when damage is above a critical level. Each algorithm is implemented with both bone lining cell (surface) sensors and osteocyte cell (internal) sensors. Each algorithm is applied to prediction of a bone multicellular unit (BMU) remodelling on the surface of a bone trabecula. It is predicted that a regulatory system capable of responding to changes in either strain or microdamage but which prioritises removal of damaged bone when damage is above a critical level, is the only one that provides a plausible prediction of BMU behaviour. A mechanism for this may be that, below a certain damage threshold, osteocyte processes can sense changes in strain and fluid flow but above the threshold damage interferes with the signalling mechanism, or causes osteocyte apoptosis so that a remodelling response occurs to remove the dead osteocytes.     

Mechanisms for osteogenic differentiation of human mesenchymal stem cells induced by fluid shear stress

Mechanical stimuli can improve bone function by promoting the proliferation and differentiation of bone cells and osteoblasts. As precursors of osteoblasts, human mesenchymal stem cells (hMSCs) are sensitive to mechanical stimuli. In recent years, fluid shear stress (FSS) has been widely used as a method of mechanical stimulation in bone tissue engineering to induce the osteogenic differentiation of hMSCs. However, the mechanism of this differentiation is not completely clear. Several signaling pathways are involved in the mechanotransduction of hMSCs responding to FSS, such as MAPK, NO/cGMP/PKG and Ca²⁺ signaling pathway. Here, we briefly review how hMSCs respond to fluid flow stimuli and focus on the signal molecules involved in this mechanotransduction.     

Structural differences in the osteocyte network between the calvaria and long bone revealed by three-dimensional fluorescence morphometry, possibly reflecting distinct mechano-adaptations and sensitivities

The structural features of osteocytes and their cellular process network are thought to allow for mechanotransduction from the bone tissue to these cells. This study applied three-dimensional fluorescence microscopy to fixed and decalcified bone specimens to quantitatively compare the osteocytes and their networks between mouse parietal bone and tibia that are physiologically enforced by distinct mechanical loads. The subsequent morphometric analysis by the surface rendering of osteocyte cell bodies revealed the tibia to have relatively enriched cytoplasm in the osteocyte cell body in comparison to the parietal bone. Furthermore, quantitative tracing of the cellular processes in silico demonstrated that the numbers of the cellular processes and their bifurcation points per osteocyte in the tibia were significantly higher than those in the parietal bone. Though the total length of the processes per osteocyte in the tibia was two times longer, its total surface area and total volume were smaller than those in the parietal bone, due to its thinner diameter. These architectural differences in the osteocytes and their networks are thus implicated in the adaptation to physiologically different loading, and may also induce distinct mechanosensitivities.     

The stress reaction and its molecular events: splicing variants

The growth of cells and tissues is regulated by stress. When body is injured, it manifests a large spectrum of metabolic, endocrine, and immune alterations, which is named stress reaction. Among them, the production of growth factors may play a critical role. For osteoblasts and myoblasts, IGF-I has been shown to be involved in the process of cells in response to overloads. There are two splicing forms, one is IGF-Ea, the other is the IGF-IEb in the rodents and corresponds to IGF-IEc in humans. The latter is markedly up-regulated in response to overloads. Therefore, it has been named mechanogrowth factor. The link between the mechanical stimulus and the gene expression represents a new and important area in cell science. Understanding the process of splicing in IGF-I helps one to investigate the mechanotransduction of cells in response to mechanical stimulation at molecular level.     

Label‐free quantitative proteome analysis of skeletal tissues under mechanical load

Skeletal tissue has the capability to adapt its mass and structure in response to mechanical stress. However, the molecular mechanism of bone and cartilage to respond to mechanical stress are not fully understood. A label‐free quantitative proteome approach was used for the first time to obtain a global perspective of the response of skeletal tissue to mechanical stress. Label‐free quantitative analysis of 1D‐PAGE‐LC/MS/MS based proteomics was applied to identify differentially expressed proteins. Differential expression analysis in the experimental groups and control group showed significant changes for 248 proteins including proteins related to proliferation, differentiation, regulation of signal transduction and energy metabolic pathways. Fluorescence labeling by incorporation of alizarin/calcein in newly formed bone minerals qualitatively demonstrated new bone formation. Skeletal tissues under mechanical load evoked marked new bone formation in comparison with the control group. Bone material apposition was evident. Our data suggest that 39 proteins were assigned a role in anabolic process. Comparisons of anabolic versus catabolic features of the proteomes show that 42 proteins were related to catabolic. In addition, some proteins were related to regulation of signal transduction and energy pathways, such as tropomyosin 4, fibronectin 1, and laminin, might be new molecular targets that are responsive to mechanical force. Differentially expressed proteins identified in this model may offer a useful starting point for elucidating novel aspects of the effects of mechanical force on skeletal tissue. J. Cell. Biochem. 108: 600–611, 2009. © 2009 Wiley‐Liss, Inc.     

Comparison of strain measurement in the mouse forearm using subject-specific finite element models,strain gaging,and digital image correlation

Mechanical loading in bone leads to the activation of bone-forming pathways that are most likely associated with a minimum strain threshold being experienced by the osteocyte. To investigate the correlation between cellular response and mechanical stimuli, researchers must develop accurate ways to measure/compute strain both externally on the bone surface and internally at the osteocyte level. This study investigates the use of finite element (FE) models to compute bone surface strains on the mouse forearm. Strains from three FE models were compared to data collected experimentally through strain gaging and digital image correlation (DIC). Each FE model was assigned subject-specific bone properties and consisted of one-dimensional springs representing the interosseous membrane. After three-point bending was performed on the ulnae and radii, moment of inertia was determined from microCT analysis of the bone region between the supports and then used along with standard beam analyses to calculate the Young’s modulus. Non-contact strain measurements from DIC were determined to be more suitable for validating numerical results than experimental data obtained through conventional strain gaging. When comparing strain responses in the three ulnae, we observed a 3–14% difference between numerical and DIC strains while the strain gage values were 37–56% lower than numerical values. This study demonstrates a computational approach for capturing bone surface strains in the mouse forearm. Ultimately, strains from these macroscale models can be used as inputs for microscale and nanoscale FE models designed to analyze strains directly in the osteocyte lacunae.     

A theoretical framework for strain-related trabecular bone maintenance and adaptation

It is assumed that density and morphology of trabecular bone is partially controlled by mechanical forces. How these effects are expressed in the local metabolic functions of osteoclast resorption and osteoblast formation is not known. In order to investigate possible mechano-biological pathways for these mechanisms we have proposed a mathematical theory (Nature 405 (2000) 704). This theory is based on hypothetical osteocyte stimulation of osteoblast bone formation, as an effect of elevated strain in the bone matrix, and a role for microcracks and disuse in promoting osteoclast resorption. Applied in a 2-D Finite Element Analysis model, the theory explained the formation of trabecular patterns. In this article we present a 3-D FEA model based on the same theory and investigated its potential morphological predictability of metabolic reactions to mechanical loads. The computations simulated the development of trabecular morphological details during growth, relative to measurements in growing pigs, reasonably realistic. They confirmed that the proposed mechanisms also inherently lead to optimal stress transfer. Alternative loading directions produced new trabecular orientations. Reduction of load reduced trabecular thickness, connectivity and mass in the simulation, as is seen in disuse osteoporosis. Simulating the effects of estrogen deficiency through increased osteoclast resorption frequencies produced osteoporotic morphologies as well, as seen in post-menopausal osteoporosis. We conclude that the theory provides a suitable computational framework to investigate hypothetical relationships between bone loading and metabolic expressions.