Finite element analysis of bone remodelling with piezoelectric effects using an open-source framework

Bone tissue exhibits piezoelectric properties and thus is capable of transforming mechanical stress into electrical potential. Piezoelectricity has been shown to play a vital role in bone adaptation and remodelling processes. Therefore, to better understand the interplay between mechanical and electrical stimulation during these processes, strain-adaptive bone remodelling models without and with considering the piezoelectric effect were simulated using the Python-based open-source software framework. To discretise numerical attributes, the finite element method (FEM) was used for the spatial variables and an explicit Euler scheme for the temporal derivatives. The predicted bone apparent density distributions were qualitatively and quantitatively evaluated against the radiographic scan of a human proximal femur and the bone apparent density calculated using a bone mineral density (BMD) calibration phantom, respectively. Additionally, the effect of the initial bone density on the resulting predicted density distribution was investigated globally and locally. The simulation results showed that the electrically stimulated bone surface enhanced bone deposition and these are in good agreement with previous findings from the literature. Moreover, mechanical stimuli due to daily physical activities could be supported by therapeutic electrical stimulation to reduce bone loss in case of physical impairment or osteoporosis. The bone remodelling algorithm implemented using an open-source software framework facilitates easy accessibility and reproducibility of finite element analysis made.

     

A multiscale mechanobiological model of bone remodelling predicts site-specific bone loss in the femur during osteoporosis and mechanical disuse

We propose a multiscale mechanobiological model of bone remodelling to investigate the site-specific evolution of bone volume fraction across the midshaft of a femur. The model includes hormonal regulation and biochemical coupling of bone cell populations, the influence of the microstructure on bone turnover rate, and mechanical adaptation of the tissue. Both microscopic and tissue-scale stress/strain states of the tissue are calculated from macroscopic loads by a combination of beam theory and micromechanical homogenisation. This model is applied to simulate the spatio-temporal evolution of a human midshaft femur scan subjected to two deregulating circumstances: (i) osteoporosis and (ii) mechanical disuse. Both simulated deregulations led to endocortical bone loss, cortical wall thinning and expansion of the medullary cavity, in accordance with experimental findings. Our model suggests that these observations are attributable to a large extent to the influence of the microstructure on bone turnover rate. Mechanical adaptation is found to help preserve intracortical bone matrix near the periosteum. Moreover, it leads to non-uniform cortical wall thickness due to the asymmetry of macroscopic loads introduced by the bending moment. The effect of mechanical adaptation near the endosteum can be greatly affected by whether the mechanical stimulus includes stress concentration effects or not.     

Multiscale modeling of bone tissue with surface and permeability control

Natural biological materials usually present a hierarchical arrangement with various structural levels. The biomechanical behavior of the complex hierarchical structure of bone is investigated with models that address the various levels corresponding to different scales. Models that simulate the bone remodeling process concurrently at different scales are in development. We present a multiscale model for bone tissue adaptation that considers the two top levels, whole bone and trabecular architecture. The bone density distribution is calculated at the macroscale (whole bone) level, and the trabecular structure at the microscale level takes into account its mechanical properties as well as surface density and permeability. The bone remodeling process is thus formulated as a material distribution problem at both scales. At the local level, the biologically driven information of surface density and permeability characterizes the trabecular structure. The model is tested by a three-dimensional simulation of bone tissue adaptation for the human femur. The density distribution of the model shows good agreement with the actual bone density distribution. Permeability at the microstructural level assures interconnectivity of pores, which mimics the interconnectivity of trabecular bone essential for vascularization and transport of nutrients. The importance of this multiscale model relays on the flexibility to control the morphometric parameters that characterize the trabecular structure. Therefore, the presented model can be a valuable tool to define bone quality, to assist with diagnosis of osteoporosis, and to support the development of bone substitutes.     

Low magnitude and high frequency mechanical loading prevents decreased bone formation responses of 2T3 preosteoblasts

Bone loss due to osteoporosis or disuse such as in paraplegia or microgravity is a significant health problem. As a treatment for osteoporosis, brief exposure of intact animals or humans to low magnitude and high frequency (LMHF) mechanical loading has been shown to normalize and prevent bone loss. However, the underlying molecular changes and the target cells by which LMHF mechanical loading alleviate bone loss are not known. Here, we hypothesized that direct application of LMHF mechanical loading to osteoblasts alters their cell responses, preventing decreased bone formation induced by disuse or microgravity conditions. To test our hypothesis, preosteoblast 2T3 cells were exposed to a disuse condition using the random positioning machine (RPM) and intervened with an LMHF mechanical load (0.1–0.4 g at 30 Hz for 10–60 min/day). Exposure of 2T3 cells to the RPM decreased bone formation responses as determined by alkaline phosphatase (ALP) activity and mineralization even in the presence of a submaximal dose of BMP4 (20 ng/ml). However, LMHF mechanical loading prevented the RPM‐induced decrease in ALP activity and mineralization. Mineralization induced by LMHF mechanical loading was enhanced by treatment with bone morphogenic protein 4 (BMP4) and blocked by the BMP antagonist noggin, suggesting a role for BMPs in this response. In addition, LMHF mechanical loading rescued the RPM‐induced decrease in gene expression of ALP, runx2, osteomodulin, parathyroid hormone receptor 1, and osteoglycin. These findings suggest that preosteoblasts may directly respond to LMHF mechanical loading to induce differentiation responses. The mechanosensitive genes identified here provide potential targets for pharmaceutical treatments that may be used in combination with low level mechanical loading to better treat osteoporosis or disuse‐induced bone loss. J. Cell. Biochem. 106: 306–316, 2009. © 2009 Wiley‐Liss, Inc.     

Altered Mechanical Environment of Bone Cells in an Animal Model of Short- and Long-Term Osteoporosis

Alterations in bone tissue composition during osteoporosis likely disrupt the mechanical environment of bone cells and may thereby initiate a mechanobiological response. It has proved challenging to characterize the mechanical environment of bone cells in vivo, and the mechanical environment of osteoporotic bone cells is not known. The objective of this research is to characterize the local mechanical environment of osteocytes and osteoblasts from healthy and osteoporotic bone in a rat model of osteoporosis. Using a custom-designed micromechanical loading device, we apply strains representative of a range of physical activity (up to 3000 με) to fluorescently stained femur samples from normal and ovariectomized rats. Confocal imaging was simultaneously performed, and digital image correlation techniques were applied to characterize cellular strains. In healthy bone tissue, osteocytes experience higher maximum strains (31,028 ± 4213 με) than osteoblasts (24,921 ± 3,832 με), whereas a larger proportion of the osteoblast experiences strains >10,000 με. Most interestingly, we show that osteoporotic bone cells experience similar or higher maximum strains than healthy bone cells after short durations of estrogen deficiency (5 weeks), and exceeded the osteogenic strain threshold (10,000 με) in a similar or significantly larger proportion of the cell (osteoblast, 12.68% vs. 13.68%; osteocyte, 15.74% vs. 5.37%). However, in long-term estrogen deficiency (34 weeks), there was no significant difference between bone cells in healthy and osteoporotic bone. These results suggest that the mechanical environment of bone cells is altered during early-stage osteoporosis, and that mechanobiological responses act to restore the mechanical environment of the bone tissue after it has been perturbed by ovariectomy.     

Peripheral blood monocyte-expressed ANXA2 gene is involved in pathogenesis of osteoporosis in humans

Low bone mineral density (BMD) is a risk factor of osteoporosis and has strong genetic determination. Genes influencing BMD and fundamental mechanisms leading to osteoporosis have yet to be fully determined. Peripheral blood monocytes (PBM) are potential osteoclast precursors, which could access to bone resorption surfaces and differentiate into osteoclasts to resorb bone. Herein, we attempted to identify osteoporosis susceptibility gene(s) and characterize their function(s), through an initial proteomics discovery study on PBM in vivo, and multiscale validation studies in vivo and in vitro. Utilizing the quantitative proteomics methodology LC-nano-ESI-MS(E), we discovered that a novel protein, i.e. ANXA2, was up-regulated twofold in PBM in vivo in Caucasians with extremely low BMD (cases) versus those with extremely high BMD (controls) (n = 28, p < 0.05). ANXA2 gene up-regulation in low BMD subjects was replicated at the mRNA level in PBM in vivo in a second and independent case-control sample (n = 80, p < 0.05). At the DNA level, we found that SNPs in the ANXA2 gene were associated with BMD variation in a 3(rd) and independent case-control sample (n = 44, p < 0.05), as well as in a random population sample (n = 997, p < 0.05). The above integrative evidence strongly supports the concept that ANXA2 is involved in the pathogenesis of osteoporosis in humans. Through a follow-up cellular functional study, we found that ANXA2 protein significantly promoted monocyte migration across an endothelial barrier in vitro (p < 0.001). Thus, elevated ANXA2 protein expression level, as detected in low BMD subjects, probably stimulates more PBM migration through the blood vessel walls to bone resorption surfaces in vivo, where they differentiate into higher number of osteoclasts and resorb bone at higher rates, thereby decreasing BMD. In conclusion, this study identified a novel osteoporosis susceptibility gene ANXA2, and suggested a novel pathophysiological mechanism, mediated by ANXA2, for osteoporosis in humans.     

Multiscale characterization and micromechanical modeling of crop stem materials

An essential prerequisite for the efficient biomechanical tailoring of crops is to accurately relate mechanical behavior to compositional and morphological properties across different length scales. In this article, we develop a multiscale approach to predict macroscale stiffness and strength properties of crop stem materials from their hierarchical microstructure. We first discuss the experimental multiscale characterization based on microimaging (micro-CT, light microscopy, transmission electron microscopy) and chemical analysis, with a particular focus on oat stems. We then derive in detail a general micromechanics-based model of macroscale stiffness and strength. We specify our model for oats and validate it against a series of bending experiments that we conducted with oat stem samples. In the context of biomechanical tailoring, we demonstrate that our model can predict the effects of genetic modifications of microscale composition and morphology on macroscale mechanical properties of thale cress that is available in the literature.

     

New Treatment Modalities in Osteoporosis

ObjectiveTo describe recently discovered agents for the management of osteoporosis.MethodsA literature review (PubMed search) was conducted to identify agents at various stages of development for osteoporosis treatment. Agents under study or review for approval were included.ResultsIn menopause, bone remodeling is increased, and agents that suppress bone resorption can stabilize bone mass. In contrast, agents that target the osteoblast can increase bone formation and bone mass. Novel antiresorptive agents can target the formation or the activity of osteoclasts. They include denosumab, an antibody to receptor activated nuclear factor kB; new selective estrogen receptor modulators, such as bazedoxifene; and cathepsin K inhibitors, such as odanacatib. Src kinase inhibitors are in the early phases of development. Parathyroid hormone is the only approved anabolic agent for the treatment of osteoporosis. Novel anabolic therapies for osteoporosis may include the use of factors with anabolic properties for bone or the neutralization of growth factor antagonists. Recent discoveries have demonstrated that the Wnt/β-catenin signaling pathway has a central role in osteoblastic cell differentiation. Antibodies to Wnt antagonists, such as sclerostin, are under development as new therapeutic approaches for osteoporosis. Anabolic therapies have the potential to enhance bone mass, but their long-term safety must be proven.ConclusionsNew developments in the treatment of osteoporosis include novel antiresorptive and anabolic agents. Their success will depend on their long-term effectiveness and safety profile. (Endocr Pract. 2010;16:855-863)     

氧化损伤和微重力致骨质疏松症的研究进展

伴随着人口老龄化日益严重,骨质疏松症作为"悄无声息的流行病"逐渐引起人们的注意。氧化损伤和力学刺激是造成骨质疏松的两个主要原因。一方面氧化损伤可通过刺激FoxOs信号通路抑制成骨细胞分化,造成骨质疏松,另一方面机体在长期缺乏负荷力刺激时也会发生废用性骨丢失,二者之间存在着紧密的联系。Nrf2作为细胞应对氧化损伤的主要防御机制,可调控多种抗氧化蛋白酶转录,在氧化损伤所造成的骨质疏松中扮演着重要角色。本文综述了氧化损伤和微重力造成骨质疏松的机制以及Nrf2对抗氧化损伤的调节和对修复骨质发育的影响。     

A new biosensor for osteoporosis detection

Abstract

Osteoporosis is a disease that is characterized by deterioration of bone tissue and increased risk of fracture as it leads to a decrease in bone mineral density, which is an important public health problem. Today, bone mineral density is measured by radiological techniques. Alternative techniques are needed because of the disadvantages such as excessive radiation intake, the cost of radiological techniques, and the necessity for specialist personnel for the devices. The quantitative determination of biochemical markers that play a role in bone mineralization may be a good alternative for the osteoporosis diagnosis and especially in the follow-up of treatment.

In this study, a specific and sensitive immunological biosensor, which quantitatively determines the osteocalcin molecule, has been developed to be used in the early osteoporosis diagnosis and to evaluate the response to the drug treatment. Anti-osteocalcin antibody was immobilized onto gold electrode surface via covalent immobilization method by using 6-mercaptohexanol, 1,4-butanedioldiglycidyl ether, ethanolamine, and glutaraldehyde. Immobilization steps and biosensor characterization were specified by cyclic voltammetry and electrochemical impedance spectroscopy. The detection time and range of Ocn biosensor were determined as 45?min and 10–60?pg µL^?1 Ocn concentration, respectively. The Ocn biosensor was successfully applied in artificial serum samples spiked with Ocn.     

Top down and bottom up engineering of bone

The goal of this retrospective article is to place the body of my lab's multiscale mechanobiology work in context of top-down and bottom-up engineering of bone. We have used biosystems engineering, computational modeling and novel experimental approaches to understand bone physiology, in health and disease, and across time (in utero, postnatal growth, maturity, aging and death, as well as evolution) and length scales (a single bone like a femur, m; a sample of bone tissue, mm-cm; a cell and its local environment, μm; down to the length scale of the cell's own skeleton, the cytoskeleton, nm). First we introduce the concept of flow in bone and the three calibers of porosity through which fluid flows. Then we describe, in the context of organ-tissue, tissue-cell and cell-molecule length scales, both multiscale computational models and experimental methods to predict flow in bone and to understand the flow of fluid as a means to deliver chemical and mechanical cues in bone. Addressing a number of studies in the context of multiple length and time scales, the importance of appropriate boundary conditions, site specific material parameters, permeability measures and even micro-nanoanatomically correct geometries are discussed in context of model predictions and their value for understanding multiscale mechanobiology of bone. Insights from these multiscale computational modeling and experimental methods are providing us with a means to predict, engineer and manufacture bone tissue in the laboratory and in the human body.     

Benzofuran pyran hybrid prevents glucocorticoid induced osteoporosis in mice via modulation of canonical Wnt/β-catenin signaling

Glucocorticoid induced osteoporosis (GIOP) is the second most leading cause of osteoporosis. We have identified a compound, a benzofuran pyran hybrid compound 4e that has osteogenic potential and we wanted to assess its efficacy in GIOP in male mice. We assessed the effect of dexamethasone and compound 4e on primary osteoblasts using various cell based and immunofluorescence assays. For in vivo studies we administered methylprednisolone and compound 4e as a prophylactic measure in male Balb/c mice for 28 days and then evaluated the effect on bone microarchitecture by microCT, bone formation by histology along with clinically relevant bone markers. Compound 4e preserved osteoblast differentiation as evident by higher ALP positive cells and mineralization in compound treated groups. Compound 4e also increased the expression of osteogenic genes. This compound guarded β-catenin expression both in vitro and in vivo as confirmed by western blot and immunofluorescence assays. This led to the preservation of bone microarchitecture and cortical thickness at 2.5 mg kg^?1 and 5 mg kg^?1 doses. Further compound 4e enhanced bone formation rate and regulated osteocyte death. The osteogenic potential of compound 4e was reflected by an increased level of serum marker osteocalcin and decreased levels of SOST and CTX-I. Overall, Compound 4e is able to overcome the catabolic effect of dexamethasone on bone by targeting the canonical WNT/β-catenin signaling as evidenced by both in vitro and in vivo studies.

     

Moderate-Intensity Rotating Magnetic Fields Do Not Affect Bone Quality and Bone Remodeling in Hindlimb Suspended Rats

Abundant evidence has substantiated the positive effects of pulsed electromagnetic fields (PEMF) and static magnetic fields (SMF) on inhibiting osteopenia and promoting fracture healing. However, the osteogenic potential of rotating magnetic fields (RMF), another common electromagnetic application modality, remains poorly characterized thus far, although numerous commercial RMF treatment devices have been available on the market. Herein the impacts of RMF on osteoporotic bone microarchitecture, bone strength and bone metabolism were systematically investigated in hindlimb-unloaded (HU) rats. Thirty two 3-month-old male Sprague-Dawley rats were randomly assigned to the Control (n = 10), HU (n = 10) and HU with RMF exposure (HU+RMF, n = 12) groups. Rats in the HU+RMF group were subjected to daily 2-hour exposure to moderate-intensity RMF (ranging from 0.60 T to 0.38 T) at 7 Hz for 4 weeks. HU caused significant decreases in body mass and soleus muscle mass of rats, which were not obviously altered by RMF. Three-point bending test showed that the mechanical properties of femurs in HU rats, including maximum load, stiffness, energy absorption and elastic modulus were not markedly affected by RMF. µCT analysis demonstrated that 4-week RMF did not significantly prevent HU-induced deterioration of femoral trabecular and cortical bone microarchitecture. Serum biochemical analysis showed that RMF did not significantly change HU-induced decrease in serum bone formation markers and increase in bone resorption markers. Bone histomorphometric analysis further confirmed that RMF showed no impacts on bone remodeling in HU rats, as evidenced by unchanged mineral apposition rate, bone formation rate, osteoblast numbers and osteoclast numbers in cancellous bone. Together, our findings reveal that RMF do not significantly affect bone microstructure, bone mechanical strength and bone remodeling in HU-induced disuse osteoporotic rats. Our study indicates potentially obvious waveform-dependent effects of electromagnetic fields-stimulated osteogenesis, suggesting that RMF, at least in the present form, might not be an optimal modality for inhibiting disuse osteopenia/osteoporosis.     

Effects of electromagnetic fields on osteoporosis: A systematic literature review

ABSTRACT

Electromagnetic fields (EMFs) as a safe, effective and noninvasive treatment have been researched and used for many years in orthopedics, and the common use clinically is to promote fracture healing. The effects of EMFs on osteoporosis have not been well concerned. The balance between osteoblast and osteoclast activity as well as the balance between osteogenic differentiation and adipogenic differentiation of bone marrow mesenchymal stem cells plays an important role in the process of osteoporosis. A number of recent reports suggest that EMFs have a positive impact on the balances. In this review, we discuss the recent advances of EMFs in the treatment of osteoporosis from basic research to clinical study and introduce the possible mechanism. In addition, we presented future perspectives of application of EMFs for osteoporosis.     

Predicting changes in mechanical properties of trabecular bone by adaptive remodeling

Because changes in the mechanical properties of bone are closely related to trabecular bone remodeling, methods that consider the temporal morphological changes induced by adaptive remodeling of trabecular bone are needed to estimate long-term fracture risk and bone quality in osteoporosis. We simulated bone remodeling using simplified and pig trabecular bone models and estimated the morphology of healthy and osteoporotic cases. We then displayed the fracture risk of the remodeled models based on a cumulative histogram from high stress. The histogram showed more elements had higher stresses in the osteoporosis model, indicating that the osteoporosis model had a greater risk.     

The Co-effect of Cordyceps sinensis and Strontium on Osteoporosis in Ovariectomized Osteopenic Rats

The co-effect of Cordyceps sinensi (CS; caterpillar fungus) and strontium on ovariectomized osteopenic rats was studied in this paper. After the rats were treated orally with CS, strontium (SR), and CS rich in strontium (CSS), respectively, the urine calcium, plasma calcium, plasma phosphorus, bone mineral content, mechanical testing, and the mass of uterus, thymus, and body were examined. Both CSS and SR have a positive effect on mechanical strength and mineral content of ovariectomized osteopenic rats. However, femoral neck strength in the CSS-treated group was higher than those in the SR-treated groups. CSS and SR significantly decreased urinary calcium excretion and plasma total calcium and inorganic phosphate concentrations. On the contrary, CS and CSS significantly increased weights of atrophic uteri and weights of body and also decreased the thymus mass in animals, whereas SR did not exhibit any such effects. Our experiments have demonstrated that CSS possess a preferable effect against the decrease of bone strength and bone mineral mass caused by osteoporosis. It was caused by the co-effect of CS and strontium. The mechanism of it includes decreases bone resorption, increases bone formation, increases in body weight, and enhances 17β-estradiol-producing as well as enhancing the immune functions in animals. The data provide an important proof of concept that CSS might be a new potential therapy for the management of postmenopausal osteoporosis in humans.     

Multiscale methodology for bone remodelling simulation using coupled finite element and neural network computation

The aim of this paper is to develop a multiscale hierarchical hybrid model based on finite element analysis and neural network computation to link mesoscopic scale (trabecular network level) and macroscopic (whole bone level) to simulate the process of bone remodelling. As whole bone simulation, including the 3D reconstruction of trabecular level bone, is time consuming, finite element calculation is only performed at the macroscopic level, whilst trained neural networks are employed as numerical substitutes for the finite element code needed for the mesoscale prediction. The bone mechanical properties are updated at the macroscopic scale depending on the morphological and mechanical adaptation at the mesoscopic scale computed by the trained neural network. The digital image-based modelling technique using μ-CT and voxel finite element analysis is used to capture volume elements representativeof 2 mm³ at the mesoscale level of the femoral head. The input data for the artificial neural network are a set of bone material parameters, boundary conditions and the applied stress. The output data are the updated bone properties and some trabecular bone factors. The current approach is the first model, to our knowledge, that incorporates both finite element analysis and neural network computation to rapidly simulate multilevel bone adaptation.     

Regulation of Osteoblast Differentiation and Iron Content in MC3T3-E1 Cells by Static Magnetic Field with Different Intensities

Many studies have indicated that static magnetic fields (SMFs) have positive effects on bone tissue, including bone formation and bone healing process. Evaluating the effects of SMFs on bone cell (especially osteoblast) function and exploring the mechanism, which is critical for understanding the possible risks or benefits from SMFs to the balance of bone remodeling. Iron and magnetic fields have the natural relationship, and iron is an essential element for normal bone metabolism. Iron overload or deficiency can cause severe bone disorders including osteoporosis. However, there are few reports regarding the role of iron in the regulation of bone formation under SMFs. In this study, hypomagnetic field (HyMF) of 500 nT, moderate SMF (MMF) of 0.2 T, and high SMF (HiMF) of 16 T were used to investigate how osteoblast (MC3T3-E1) responses to SMFs and iron metabolism of osteoblast under SMFs. The results showed that SMFs did not pose severe toxic effects on osteoblast growth. During cell proliferation, iron content of osteoblast MC3T3-E1 cells was decreased in HyMF, but was increased in MMF and HiMF after exposure for 48 h. Compared to untreated control (i.e., geomagnetic field, GMF), HyMF and MMF exerted deleterious effects on osteoblast differentiation by simultaneously retarding alkaline phosphatase (ALP) activity, mineralization and calcium deposition. However, when exposed to HiMF of 16 T, the differentiation potential showed the opposite tendency with enhanced mineralization. Iron level was increased in HyMF, constant in MMF and decreased in HiMF during cell differentiation. In addition, the mRNA expression of transferrin receptor 1 (TFR1) was promoted by HyMF but was inhibited by HiMF. At the same time, HiMF of 16 T and MMF of 0.2 T increased the expression of ferroportin 1 (FPN1). In conclusion, these results indicated that osteoblast differentiation can be regulated by altering the strength of the SMF, and iron is possibly involved in this process.