A mathematical multiscale model of bone remodeling, accounting for pore space-specific mechanosensation

Glucocorticoids disrupt longitudinal advance of cortical bone basic multicellular units in the rabbit distal tibia

Glucocorticoids (GCs) are the leading cause of secondary osteoporosis. The emerging perspective, derived primarily from 2D histological study of trabecular bone, is that GC-induced bone loss arises through the uncoupling of bone formation and resorption at the level of the basic multicellular unit (BMU), which carries out bone remodeling. Here we explore the impact of GCs on cortical bone remodeling in the rabbit model. Based upon the rapid reduction of bone formation and initial elevation of resorption caused by GCs, we hypothesized that the rate of advance (longitudinal erosion rate; LER) of cortical BMUs would be increased. To test this hypothesis we divided 20 female New Zealand White rabbits into four experimental groups: ovariohysterectomy (OVH), glucocorticoid (GC), OVH + GC and SHAM controls (n = 5 animals each). Ten weeks post-surgery (OVH or sham), and two weeks after the initiation of dosing (daily subcutaneous injections of 1.5 mg/kg of methylprednisolone sodium succinate in the GC-treated groups and 1 ml of saline for the others), the right tibiae were scanned in vivo using Synchrotron Radiation (SR) in-line phase contrast micro-CT at the Canadian Light Source. After an additional 2 weeks of dosing, the rabbits were euthanized and ex vivo images were collected using desktop micro-CT. The datasets were co-registered in 3D and LER was calculated as the distance traversed by BMU cutting-cones in the 14-day interval between scans. Counter to our hypothesis, LER was greatly reduced in GC-treated rabbits. Mean LER was lower in GC (4.27 μm/d; p < 0.001) and OVH + GC (4.19 μm/d; p < 0.001), while similar in OVH (40.13 μm/d; p = 0.990), compared to SHAM (40.44 μm/d). This approximately 90 % reduction in LER with GCs was also associated with an overall disruption of BMU progression, with radial expansion of the remodeling space occurring in all directions. This unexpected outcome suggests that GCs do not simply uncouple formation and resorption within cortical BMUs and highlights the value of the time-lapsed 4D approach employed.

Estrogen and estrogen receptors mediate the mechanobiology of bone disease and repair

It is well understood that the balance of bone formation and resorption is dependent on both mechanical and biochemical factors. In addition to cell-secreted cytokines and growth factors, sex hormones like estrogen are critical to maintaining bone health. Although the direct osteoprotective function of estrogen and estrogen receptors (ERs) has been reported extensively, evidence that estrogen signaling also has a role in mediating the effects of mechanical loading on maintenance of bone mass and healing of bone injuries has more recently emerged. Recent studies have underscored the role of estrogen and ERs in many pathways of bone mechanosensation and mechanotransduction. Estrogen and ERs have been shown to augment integrin-based mechanotransduction as well as canonical Wnt/b-catenin, RhoA/ROCK, and YAP/TAZ pathways. Estrogen and ERs also influence the mechanosensitivity of not only osteocytes but also osteoblasts, osteoclasts, and marrow stromal cells. The current review will highlight these roles of estrogen and ERs in cellular mechanisms underlying bone mechanobiology and discuss their implications for management of osteoporosis and bone fractures. A greater understanding of the mechanisms behind interactions between estrogen and mechanical loading may be crucial to addressing the shortcomings of current hormonal and pharmaceutical therapies. A combined therapy approach including high-impact exercise therapy may mitigate adverse side effects and allow an effective long-term solution for the prevention, treatment, and management of bone fragility in at-risk populations. Furthermore, future implications to novel local delivery mechanisms of hormonal therapy for osteoporosis treatment, as well as the effects on bone health of applications of sex hormone therapy outside of bone disease, will be discussed.

Tumor growth for remodeling process: A 2D approach

This paper aims to present a comprehensive framework for coupling tumor-bone remodeling processes in a 2-dimensional geometry. This is achieved by introducing a bio-inspired damage that represents the growing tumor, which subsequently affects the main populations involved in the remodeling process, namely, osteoclasts, osteoblasts, and bone tissue. The model is constructed using a set of differential equations based on the Komarova's and Ayati's models, modified to incorporate the bio-inspired damage that may result in tumor mass formation. Three distinct models were developed. The first two models are based on the Komarova's governing equations, with one demonstrating an osteolytic behavior and the second one an osteoblastic model. The third model is a variation of Ayati's model, where the bio-inspired damage is induced through the paracrine and autocrine parameters, exhibiting an osteolytic behavior. The obtained results are consistent with existing literature, leading us to believe that our in-silico experiments will serve as a cornerstone for paving the way towards targeted interventions and personalized treatment strategies, ultimately improving the quality of life for those affected by these conditions.

Research progress on the regulatory mechanism of integrin-mediated mechanical stress in cells involved in bone metabolism

Mechanical stress is an internal force between various parts of an object that resists external factors and effects that cause an object to deform, and mechanical stress is essential for various tissues that are constantly subjected to mechanical loads to function normally. Integrins are a class of transmembrane heterodimeric glycoprotein receptors that are important target proteins for the action of mechanical stress stimuli on cells and can convert extracellular physical and mechanical signals into intracellular bioelectrical signals, thereby regulating osteogenesis and osteolysis. Integrins play a bidirectional regulatory role in bone metabolism. In this paper, relevant literature published in recent years is reviewed and summarized. The characteristics of integrins and mechanical stress are introduced, as well as the mechanisms underlying responses of integrin to mechanical stress stimulation. The paper focuses on integrin-mediated mechanical stress in different cells involved in bone metabolism and its associated signalling mechanisms. The purpose of this review is to provide a theoretical basis for the application of integrin-mediated mechanical stress to the field of bone tissue repair and regeneration.

Cortical bone adaptation response is region specific, but not peak load dependent: insights from μ CT image analysis and mechanostat simulations of the mouse tibia loading model

The two major aims of the present study were: (i) quantify localised cortical bone adaptation at the surface level using contralateral endpoint imaging data and image analysis techniques, and (ii) investigate whether cortical bone adaptation responses are universal or region specific and dependent on the respective peak load. For this purpose, we re-analyse previously published μ CT data of the mouse tibia loading model that investigated bone adaptation in response to sciatic neurectomy and various peak load magnitudes (F = 0, 2, 4, 6, 8, 10, 12 N). A beam theory-based approach was developed to simulate cortical bone adaptation in different sections of the tibia, using longitudinal strains as the adaptive stimuli. We developed four mechanostat models: universal, surface-based, strain directional-based, and combined surface and strain direction-based. Rates of bone adaptation in these mechanostat models were computed using an optimisation procedure (131,606 total simulations), performed on a single load case (F = 10 N). Subsequently, the models were validated against the remaining six peak loads. Our findings indicate that local bone adaptation responses are quasi-linear and bone region specific. The mechanostat model which accounted for differences in endosteal and periosteal regions and strain directions (i.e. tensile versus compressive) produced the lowest root mean squared error between simulated and experimental data for all loads, with a combined prediction accuracy of 76.6, 55.0 and 80.7% for periosteal, endosteal, and cortical thickness measurements (in the midshaft of the tibia). The largest root mean squared errors were observed in the transitional loads, i.e. F = 2 to 6 N, where inter-animal variability was highest. Finally, while endpoint imaging studies provide great insights into organ level bone adaptation responses, the between animal and loaded versus control limb variability make simulations of local surface-based adaptation responses challenging.

Studying the mechanical properties of the mandible and injury prediction under the effect of ossification factors

BACKGROUND AND OBJECTIVE

It is essential to know the quantitative interactions between biological tissues and external mechanical and chemical stimuli. This assists the physicians to better know the quantitative behavior of the tissue and plan more effective therapy. In the literature, the effect of the chemical and mechanical loading was investigated on the bone biological cell activities and some mechanical features, but a lack of prediction of bone injury under the chemical and mechanical factors was sensed. Therefore, the present study aims to investigate the effects of the application of major chemical factors involved in ossification, including RANKL1 (Receptor Activator of Nuclear Factor Kappa Beta Ligand), PTH2 (Parathyroid Hormone), and OPG3 (Osteoprotegerin) on the mandibular bone biological osteoblast and osteoclast activities and mechanical properties. Moreover, the study assesses the bone injury possibility under uniform mastication pressure applied on the premolar tooth in terms of the mechanostat theory undergoing the effects of the chemical factors.

METHODS

A 3D geometry of the mandible-tooth assembly was generated from the CT image dataset. The geometry was next purified, solidified, and exported to FEM4 (Finite Element Method) software to be meshed, where boundary conditions and loading were applied. Moreover, the mathematical system of differential equations to model the chemical factor effects on osteoblast and osteoclast activities as well as bone matrix volume fraction and elastic stiffness relations were applied. Next, the values of the equivalent strain were calculated to predict the injury states of the bone.

RESULTS

The results complied with the literature data. The results showed that RANKL and PTH increased the values of the equivalent strain from 450 με to 11500 με, while OPG reduced that from 450 με to 300 με.

CONCLUSIONS

Therefore, RANKL and PTH doses of this study put the bone at risk of injury compared to the baseline of no dose applied, while OPG secured the bone from injury.

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.

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.

Trabecular bone remodeling in the aging mouse: A micro-multiphysics agent-based in silico model using single-cell mechanomics

Bone remodeling is regulated by the interaction between different cells and tissues across many spatial and temporal scales. Notably, in silico models are regarded as powerful tools to further understand the signaling pathways that regulate this intricate spatial cellular interplay. To this end, we have established a 3D multiscale micro-multiphysics agent-based (micro-MPA) in silico model of trabecular bone remodeling using longitudinal in vivo data from the sixth caudal vertebra (CV6) of PolgA^{(D257A/D257A)} mice, a mouse model of premature aging. Our in silico model includes a variety of cells as single agents and receptor-ligand kinetics, mechanomics, diffusion and decay of cytokines which regulate the cells' behavior. We highlighted its capabilities by simulating trabecular bone remodeling in the CV6 of five mice over 4 weeks and we evaluated the static and dynamic morphometry of the trabecular bone microarchitecture. Based on the progression of the average trabecular bone volume fraction (BV/TV), we identified a configuration of the model parameters to simulate homeostatic trabecular bone remodeling, here named basal. Crucially, we also produced anabolic, anti-anabolic, catabolic and anti-catabolic responses with an increase or decrease by one standard deviation in the levels of osteoprotegerin (OPG), receptor activator of nuclear factor kB ligand (RANKL), and sclerostin (Scl) produced by the osteocytes. Our results showed that changes in the levels of OPG and RANKL were positively and negatively correlated with the BV/TV values after 4 weeks in comparison to basal levels, respectively. Conversely, changes in Scl levels produced small fluctuations in BV/TV in comparison to the basal state. From these results, Scl was deemed to be the main driver of equilibrium while RANKL and OPG were shown to be involved in changes in bone volume fraction with potential relevance for age-related bone features. Ultimately, this micro-MPA model provides valuable insights into how cells respond to their local mechanical environment and can help to identify critical pathways affected by degenerative conditions and ageing.

Experimental-based mechanobiological modeling of the anabolic and catabolic effects of breast cancer on bone remodeling

Bone is a biological tissue characterized by its hierarchical organization. This material has the ability to be continually renewed, which makes it highly adaptative to external loadings. Bone renewing is managed by a dynamic biological process called bone remodeling (BR), where continuous resorption of old bone and formation of new bone permits to change the bone composition and microstructure. Unfortunately, because of several factors, such as age, hormonal imbalance, and a variety of pathologies including cancer metastases, this process can be disturbed leading to various bone diseases. In this study, we have investigated the effect of breast cancer (BC) metastases causing osteolytic bone loss. BC has the ability to affect bone quantity in different ways in each of its primary and secondary stages. Based on a BR mathematical model, we modeled the BC cells' interaction with bone cells to assess their effect on bone volume fraction (BV/TV) evolution during the remodeling process. Some of the parameters used in our model have been determined experimentally using the enzyme-linked immune-sorbent assay (ELISA) and the MTT assay. Our numerical simulations show that primary BC plays a significant role in enhancing bone-forming cells' activity leading to a 6.22% increase in BV/TV over 1 year. On the other hand, secondary BC causes a noticeable decrease in BV/TV reaching 15.74% over 2 years.

Clinical Data for Parametrization of In Silico Bone Models Incorporating Cell-Cytokine Dynamics: A Systematic Review of Literature

In silico simulations aim to provide fast, inexpensive, and ethical alternatives to years of costly experimentation on animals and humans for studying bone remodeling, its deregulation during osteoporosis and the effect of therapeutics. Within the varied spectrum of in silico modeling techniques, bone cell population dynamics and agent-based multiphysics simulations have recently emerged as useful tools to simulate the effect of specific signaling pathways. In these models, parameters for cell and cytokine behavior are set based on experimental values found in literature; however, their use is currently limited by the lack of clinical in vivo data on cell numbers and their behavior as well as cytokine concentrations, diffusion, decay and reaction rates. Further, the settings used for these parameters vary across research groups, prohibiting effective cross-comparisons. This review summarizes and evaluates the clinical trial literature that can serve as input or validation for in silico models of bone remodeling incorporating cells and cytokine dynamics in post-menopausal women in treatment, and control scenarios. The GRADE system was used to determine the level of confidence in the reported data, and areas lacking in reported measures such as binding site occupancy, reaction rates and cell proliferation, differentiation and apoptosis rates were highlighted as targets for further research. We propose a consensus for the range of values that can be used for the cell and cytokine settings related to the RANKL-RANK-OPG, TGF-β and sclerostin pathways and a Levels of Evidence-based method to estimate parameters missing from clinical trial literature.

Effects of Mechanical Stress Stimulation on Function and Expression Mechanism of Osteoblasts

Osteoclasts and osteoblasts play a major role in bone tissue homeostasis. The homeostasis and integrity of bone tissue are maintained by ensuring a balance between osteoclastic and osteogenic activities. The remodeling of bone tissue is a continuous ongoing process. Osteoclasts mainly play a role in bone resorption, whereas osteoblasts are mainly involved in bone remodeling processes, such as bone cell formation, mineralization, and secretion. These cell types balance and restrict each other to maintain bone tissue metabolism. Bone tissue is very sensitive to mechanical stress stimulation. Unloading and loading of mechanical stress are closely related to the differentiation and formation of osteoclasts and bone resorption function as well as the differentiation and formation of osteoblasts and bone formation function. Consequently, mechanical stress exerts an important influence on the bone microenvironment and bone metabolism. This review focuses on the effects of different forms of mechanical stress stimulation (including gravity, continuously compressive pressure, tensile strain, and fluid shear stress) on osteoclast and osteoblast function and expression mechanism. This article highlights the involvement of osteoclasts and osteoblasts in activating different mechanical transduction pathways and reports changings in their differentiation, formation, and functional mechanism induced by the application of different types of mechanical stress to bone tissue. This review could provide new ideas for further microscopic studies of bone health, disease, and tissue damage reconstruction.

Fiber Rearrangement and Matrix Compression in Soft Tissues: Multiscale Hypoelasticity and Application to Tendon

It is widely accepted that the nonlinear macroscopic mechanical behavior of soft tissue is governed by fiber straightening and re-orientation. Here, we provide a quantitative assessment of this phenomenon, by means of a continuum micromechanics approach. Given the negligibly small bending stiffness of crimped fibers, the latter are represented through a number of hypoelastic straight fiber phases with different orientations, being embedded into a hypoelastic matrix phase. The corresponding representative volume element (RVE) hosting these phases is subjected to “macroscopic” strain rates, which are downscaled to fiber and matrix strain rates on the one hand, and to fiber spins on the other hand. This gives quantitative access to the fiber decrimping (or straightening) phenomenon under non-affine conditions, i.e. in the case where the fiber orientations cannot be simply linked to the macroscopic strain state. In the case of tendinous tissue, such an RVE relates to the fascicle material with 50 μm characteristic length, made up of crimped collagen bundles and a gel-type matrix in-between. The fascicles themselves act as parallel fibers in a similar matrix at the scale of a tissue-related RVE with 500 μm characteristic length. As evidenced by a sensitivity analysis and confirmed by various mechanical tests, it is the initial crimping angle which drives both the degree of straightening and the shape of the macroscopic stress-strain curve, while the final linear portion of this curve depends almost exclusively on the collagen bundle elasticity. Our model also reveals the mechanical cooperation of the tissue’s key microstructural components: while the fibers carry tensile forces, the matrices undergo hydrostatic pressure.

Multiscale modeling of bone tissue mechanobiology

Mechanical environment has a crucial role in our organism at the different levels, ranging from cells to tissues and our own organs. This regulatory role is especially relevant for bones, given their importance as load-transmitting elements that allow the movement of our body as well as the protection of vital organs from load impacts. Therefore bone, as living tissue, is continuously adapting its properties, shape and repairing itself, being the mechanical loads one of the main regulatory stimuli that modulate this adaptive behavior. Here we review some key results of bone mechanobiology from computational models, describing the effect that changes associated to the mechanical environment induce in bone response, implant design and scaffold-driven bone regeneration.

Parameter reduction, sensitivity studies, and correlation analyses applied to a mechanobiologically regulated bone cell population model of the bone metabolism

When striving for reconstructing and predicting bone remodeling processes by means of mathematical models, cell population models have become a popular option. From a conceptual point of view, these models are able to take into account an arbitrary amount of regulatory mechanisms driving the development of bone cells and their activities. However, in most cases, the models include a large number of parameters; and most of those parameters cannot be measured, which certainly compromises the credibility of cell population models. Here, new insights are presented as to the potential improvement of this unsatisfactory situation. In particular, a previously published bone remodeling model was considered, and based on combination and merging of the original parameters, the total number of parameters could be reduced from 28 to 18, without impairing the model's versatility and significance. Furthermore, a comprehensive number of one- and two-variable sensitivity studies were performed, pointing out which parameters (alone and in combination with other parameters) influence the model predictions significantly - for that purpose, the mean squared relative error (MSRE) between simulations based on the original parameters and based on varied parameters was considered as failure measure. It has turned out that the model is significantly more sensitive to parameters which can be considered as phenomenological (such as differentiation, proliferation, and apoptosis rates) than to parameters which are directly related to specific processes (such as dissociation rate constants, and maximum concentrations of the involved factors). Using common correlation measures (such as Pearson, Spearman, and partial ranked correlation coefficients), correlation studies revealed that the correlations between most parameters and the MSRE are weak, while a few parameters exhibited moderate correlations. In conclusion, the results shown in this paper provide valuable insights concerning the design of new experiments allowing for measurement of the parameters which are most influential in the context of bone remodeling simulation.

A General Mechano-Pharmaco-Biological Model for Bone Remodeling Including Cortisol Variation

The process of bone remodeling requires a strict coordination of bone resorption and formation in time and space in order to maintain consistent bone quality and quantity. Bone-resorbing osteoclasts and bone-forming osteoblasts are the two major players in the remodeling process. Their coordination is achieved by generating the appropriate number of osteoblasts since osteoblastic-lineage cells govern the bone mass variation and regulate a corresponding number of osteoclasts. Furthermore, diverse hormones, cytokines and growth factors that strongly link osteoblasts to osteoclasts coordinated these two cell populations. The understanding of this complex remodeling process and predicting its evolution is crucial to manage bone strength under physiologic and pathologic conditions. Several mathematical models have been suggested to clarify this remodeling process, from the earliest purely phenomenological to the latest biomechanical and mechanobiological models. In this current article, a general mathematical model is proposed to fill the gaps identified in former bone remodeling models. The proposed model is the result of combining existing bone remodeling models to present an updated model, which also incorporates several important parameters affecting bone remodeling under various physiologic and pathologic conditions. Furthermore, the proposed model can be extended to include additional parameters in the future. These parameters are divided into four groups according to their origin, whether endogenous or exogenous, and the cell population they affect, whether osteoclasts or osteoblasts. The model also enables easy coupling of biological models to pharmacological and/or mechanical models in the future.

Integration of mechanics and biology in computer simulation of bone remodeling

Bone remodeling is a complex physiological process that spans across multiple spatial and temporal scales and is regulated by both mechanical and hormonal cues. An imbalance between bone resorption and bone formation in the process of bone remodeling may lead to various bone pathologies. One powerful and non-invasive approach to gain new insights into mechano-adaptive bone remodeling is computer modeling and simulation. Recent findings in bone physiology and advances in computer modeling have provided a unique opportunity to study the integration of mechanics and biology in bone remodeling. Our objective in this review is to critically appraise recent advances and developments and discuss future research opportunities in computational bone remodeling approaches that enable integration of mechanics and cellular and molecular pathways. Based on the critical appraisal of the relevant recent published literature, we conclude that multiscale in silico integration of personalized bone mechanics and mechanobiology combined with data science and analytics techniques offer the potential to deepen our knowledge of bone remodeling and provide ample opportunities for future research.

Load adaptation through bone remodeling: a mechanobiological model coupled with the finite element method

This work proposes a novel tissue-scale mechanobiological model of bone remodeling to study bone's adaptation to distinct loading conditions. The devised algorithm describes the mechanosensitivity of bone and its impact on bone cells' functioning through distinct signaling factors. In this study, remodeling is mechanically ruled by variations of the strain energy density (SED) of bone, which is determined by performing a linear elastostatic analysis combined with the finite element method. Depending on the SED levels and on a set of biological signaling factors ([Formula: see text] parameters), osteoclasts and osteoblasts can be mechanically triggered. To reproduce this phenomenon, this work proposes a new set of [Formula: see text] parameters. The combined response of osteoclasts and osteoblasts will then affect bone's apparent density, which is correlated with other mechanical properties of bone, through a phenomenological law. Thus, this novel model proposes a constant interplay between the mechanical and biological components of the process. The spatiotemporal simulation used to validate this new approach is a benchmark example composed by two distinct phases: (1) pre-orientation and (2) load adaptation. On both of them, bone is able to adapt its morphology according to the loading condition, achieving the required trabecular distribution to withstand the applied loads. Moreover, the equilibrium morphology reflects the orientation of the load. These preliminary results support the new approach proposed in this study.

Subject-specific multiscale modeling of aortic valve biomechanics

A Finite Element workflow for the multiscale analysis of the aortic valve biomechanics was developed and applied to three physiological anatomies with the aim of describing the aortic valve interstitial cells biomechanical milieu in physiological conditions, capturing the effect of subject-specific and leaflet-specific anatomical features from the organ down to the cell scale. A mixed approach was used to transfer organ-scale information down to the cell-scale. Displacement data from the organ model were used to impose kinematic boundary conditions to the tissue model, while stress data from the latter were used to impose loading boundary conditions to the cell level. Peak of radial leaflet strains was correlated with leaflet extent variability at the organ scale, while circumferential leaflet strains varied over a narrow range of values regardless of leaflet extent. The dependency of leaflet biomechanics on the leaflet-specific anatomy observed at the organ length-scale is reflected, and to some extent emphasized, into the results obtained at the lower length-scales. At the tissue length-scale, the peak diastolic circumferential and radial stresses computed in the fibrosa correlated with the leaflet surface area. At the cell length-scale, the difference between the strains in two main directions, and between the respective relationships with the specific leaflet anatomy, was even more evident; cell strains in the radial direction varied over a relatively wide range (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.36-0.87$$\end{document}0.36-0.87) with a strong correlation with the organ length-scale radial strain (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R^{2}= 0.95$$\end{document}R2=0.95); conversely, circumferential cell strains spanned a very narrow range (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.75-0.88$$\end{document}0.75-0.88) showing no correlation with the circumferential strain at the organ level (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R^{2}= 0.02$$\end{document}R2=0.02). Within the proposed simulation framework, being able to account for the actual anatomical features of the aortic valve leaflets allowed to gain insight into their effect on the structural mechanics of the leaflets at all length-scales, down to the cell scale.

A mathematical biomechanical model for bone remodeling integrated with a radial point interpolating meshless method

Bone remodeling is a highly complex process, in which bone cells interact and regulate bone's apparent density as a response to several external and internal stimuli. In this work, this process is numerically described using a novel 2D biomechanical model. Some of the new features in this model are (i) the mathematical parameters used to determine bone's apparent density and cellular density; (ii) an automatic boundary recognition step to spatially control bone remodeling and (iii) an approach to mimic the mechanical transduction to osteoclasts and osteoblasts. Moreover, this model is combined with a meshless approach - the Radial Point Interpolation Method (RPIM). The use of RPIM is an asset for this application, especially in the construction of the boundary maps. This work studies bone's adaptation to a certain loading regime through bone resorption. The signaling pathways of bone cells are dependent on the level of strain energy density (SED) in bone. So, when SED changes, bone cells' functioning is affected, causing also changes on bone's apparent density. With this model, bone is able to achieve an equilibrium state, optimizing its structure to withstand the applied loads. Results suggest that this model has the potential to provide high quality solutions while being a simpler alternative to more complex bone remodeling models in the literature.

Toward a Mathematical Modeling of Diseases' Impact on Bone Remodeling: Technical Review

A wide variety of bone diseases have hitherto been discovered, such as osteoporosis, Paget's disease, osteopetrosis, and metastatic bone disease, which are not well defined in terms of changes in biochemical and mechanobiological regulatory factors. Some of these diseases are secondary to other pathologies, including cancer, or to some clinical treatments. To better understand bone behavior and prevent its deterioration, bone biomechanics have been the subject of mathematical modeling that exponentially increased over the last years. These models are becoming increasingly complex. The current paper provides a timely and critical analysis of previously developed bone remodeling mathematical models, particularly those addressing bone diseases. Besides, mechanistic pharmacokinetic/pharmacodynamic (PK/PD) models, which englobe bone disease and its treatment's effect on bone health. Therefore, the review starts by presenting bone remodeling cycle and mathematical models describing this process, followed by introducing some bone diseases and discussing models of pathological mechanisms affecting bone, and concludes with exhibiting the available bone treatment procedures considered in the PK/PD models.

Modeling Progressive Damage Accumulation in Bone Remodeling Explains the Thermodynamic Basis of Bone Resorption by Overloading

Computational modeling of skeletal tissue seeks to predict the structural adaptation of bone in response to mechanical loading. The theory of continuum damage-repair, a mathematical description of structural adaptation based on principles of damage mechanics, continues to be developed and utilized for the prediction of long-term peri-implant outcomes. Despite its technical soundness, CDR does not account for the accumulation of mechanical damage and irreversible deformation. In this work, a nonlinear mathematical model of independent damage accumulation and plastic deformation is developed in terms of the CDR formulation. The proposed model incorporates empirical correlations from uniaxial experiments. Supporting elements of the model are derived, including damage and yielding criteria, corresponding consistency conditions, and the basic, necessary forms for integration during loading. Positivity of mechanical dissipation due to damage is proved, while strain-based, associative plastic flow and linear hardening describe post-yield behavior. Calibration of model parameters to the empirical correlations from which the model was derived is then accomplished. Results of numerical experiments on a point-wise specimen show that damage and plasticity inhibit bone formation by dissipation of energy available to biological processes, leading to material failure that would otherwise be predicted to experience a net gain of bone.

Simulation of bone remodeling around a femoral prosthesis using a model that accounts for biological and mechanical interactions

The present study focuses on a model for three-dimensional bone remodeling of the human femur that considers cellular dynamics to determine the volume fraction of new bone. The model considers the interaction among responsive osteoblasts, active osteoblasts, and osteoclasts, as well as signaling molecules and parathyroid hormone (PTH). The stimulus of the model has a systemic origin due to the PTH effect, and a local origin due to the action of cytokines, growth factors, and mechanical stimuli near the site of the bone cells. The present work considers that the mechanical stimulus that activates cellular activity is obtained from stresses acting on the bone tissue and the number of daily loading cycles. In addition to simulating the bone modeling process in an intact femur, the numerical model is used to simulate bone adaptation in relation to the stress shielding phenomenon after the implantation of a femoral prosthesis. The results showed that the simulations provide a distribution of bone density that is similar to a radiograph and, in addition, allows the visualization of osteoblast and osteoclast dynamics in bone adaptation response after prosthesis implantation.

"Variances" and "in-variances" in hierarchical porosity and composition, across femoral tissues from cow, horse, ostrich, emu, pig, rabbit, and frog

It is very well known that bone is a hierarchically organized material produced by bone cells residing in the fluid environments filling (larger) vascular pores and (smaller) lacunar pores. The extracellular space consists of hydroxyapatite crystals, collagen type I molecules, and water with non-collageneous organics. It is less known to which extent the associated quantities (mineral, organic, and water concentrations; vascular, lacunar, and extracellular porosities) vary across species, organs, and ages. We here investigate the aforementioned quantities across femoral shaft tissues from cow, horse, emu, frog, ostrich, pig, and rabbit; by means of light microscopy and dehydration-demineralization tests; thereby revealing interesting invariances: The extracellular volume fractions of organic matter turn out to be similar across all tested non-amphibian tissues; as do the extracellular volume fractions of hydroxyapatite across all tested mammals. Hence, the chemical composition of the femoral extracellular bone matrix is remarkably "invariant" across differently aged mammals; while the water content shows significant variations, as does the partitions of water between the different pore spaces. The latter exhibit strikingly varying morphologies as well. This finding adds to the ample "universal patterns" in the sense of evolutionary developmental biology; and it provides interesting design requirements for the development of novel biomimetic tissue engineering solutions.

A mechano-chemo-biological model for bone remodeling with a new mechano-chemo-transduction approach

Bone remodeling is a fundamental biological process that develops in bone tissue along its whole lifetime. It refers to a continuous bone transformation with new bone formation and old bone resorption that changes the internal microstructure and composition of the tissue. The main objectives of bone remodeling are: repair of the internal microcracks; adaptation of the macroscopic stiffness and strength to the actual changing mechanical demands; and control of the calcium homeostasis. Understanding this process and predicting its evolution is critical to reduce the effects of long-term disuse as happens during periods of reduced mobility. It is also important in the design of bone implants to avoid long-term stress shielding. Many mathematical models have been proposed from the earliest purely phenomenological to the latest that include biological knowledge. However, there still exists a lack of connection between the mechanical driving force and the biochemical and cell processes it triggers. Here, and following previous works that model independently the mechanobiological and biochemical processes in bone remodeling, we present a more complete model, useful for both cortical and trabecular bone, that uses a new mechanotransduction approach based on the effect of strains onto the bonding-unbonding rate of RANK/RANKL/OPG receptor-ligand reactions. We compare the results of this model with previous ones, showing a good agreement in similar conditions. We also apply it to realistic situations such as a femoral bone after implantation of a hip prosthesis, getting similar results to the clinical ones in the final bone density distribution. Finally, we extend this approach to the anisotropic case, getting not only the mean density, but also the directional homogenization of the microstructure. This biochemical approach permits, not only to predict the bone evolution under changes in the mechanical loads, but also, to consider anabolic and catabolic drugs to control bone density, such as those used in osteoporosis.

Strain shielding for cemented hip implants

BACKGROUND

Long-term survival of hip implants is of increasing relevance due to the rising life expectancy. The biomechanical effect of strain shielding as a result of implant insertion may lead to bone resorption, thus increasing risk for implant loosening and periprosthetic fractures. Patient-specific quantification of strain shielding could assist orthopedic surgeons in choosing the biomechanically most appropriate prosthesis.

METHODS

Validated quantitative CT-based finite element models of five femurs in intact and implanted states were considered to propose a systematic algorithm for strain shielding quantification. Three different strain measures were investigated and the most appropriate measure for strain shielding quantification is recommended. It is used to demonstrate a practical femur-specific implant selection among three common designs.

FINDINGS

Strain shielding measures demonstrated similar trends in all Gruen zones except zone 1, where the volumetric strain measure differed from von-Mises and maximum principal strains. The volumetric strain measure is in better agreement with clinical bone resorption records. It is also consistent with the biological mechanism of bone remodeling so it is recommended for strain shielding quantification. Applying the strain shielding algorithm on three different implants for a specific femur suggests that the collared design is preferable. Such quantitative biomechanical input is valuable for practical patient specific implant selection.

INTERPRETATION

Volumetric strain should be considered for strain shielding examination. The presented methodology may potentially enable patient-specific pre-operative strain shielding evaluation so to minimize strain shielding. It should be further used in a longitudinal study so to correlate between strain shielding predictions and clinical bone resorption.

Porous Scaffold Design for Additive Manufacturing in Orthopedics: A Review

With the increasing application of orthopedic scaffolds, a dramatically increasing number of requirements for scaffolds are precise. The porous structure has been a fundamental design in the bone tissue engineering or orthopedic clinics because of its low Young's modulus, high compressive strength, and abundant cell accommodation space. The porous structure manufactured by additive manufacturing (AM) technology has controllable pore size, pore shape, and porosity. The single unit can be designed and arrayed with AM, which brings controllable pore characteristics and mechanical properties. This paper presents the current status of porous designs in AM technology. The porous structures are stated from the cellular structure and the whole structure. In the aspect of the cellular structure, non-parametric design and parametric design are discussed here according to whether the algorithm generates the structure or not. The non-parametric design comprises the diamond, the body-centered cubic, and the polyhedral structure, etc. The Voronoi, the Triply Periodic Minimal Surface, and other parametric designs are mainly discussed in parametric design. In the discussion of cellular structures, we emphasize the design, and the resulting biomechanical and biological effects caused by designs. In the aspect of the whole structure, the recent experimental researches are reviewed on uniform design, layered gradient design, and layered gradient design based on topological optimization, etc. These parts are summarized because of the development of technology and the demand for mechanics or bone growth. Finally, the challenges faced by the porous designs and prospects of porous structure in orthopedics are proposed in this paper.

Mechanical regulation of bone formation and resorption around implants in a mouse model of osteopenic bone

Although mechanical stimulation is considered a promising approach to accelerate implant integration, our understanding of load-driven bone formation and resorption around implants is still limited. This lack of knowledge may delay the development of effective loading protocols to prevent implant loosening, especially in osteoporosis. In healthy bone, formation and resorption are mechanoregulated processes. In the intricate context of peri-implant bone regeneration, it is not clear whether bone (re)modelling can still be load-driven. Here, we investigated the mechanical control of peri-implant bone (re)modelling with a well-controlled mechanobiological experiment. We applied cyclic mechanical loading after implant insertion in tail vertebrae of oestrogen depleted mice and we monitored peri-implant bone response by in vivo micro-CT. Experimental data were combined with micro-finite element simulations to estimate local tissue strains in (re)modelling locations. We demonstrated that a substantial increase in bone mass around the implant could be obtained by loading the entire bone. This augmentation could be attributed to a large reduction in bone resorption rather than to an increase in bone formation. We also showed that following implantation, mechanical regulation of bone (re)modelling was transiently lost. Our findings should help to clarify the role of mechanical stimulation on the maintenance of peri-implant bone mass.

Endothelial cells produce angiocrine factors to regulate bone and cartilage via versatile mechanisms

Blood vessels are conduits distributed throughout the body, supporting tissue growth and homeostasis by the transport of cells, oxygen and nutrients. Endothelial cells (ECs) form the linings of the blood vessels, and together with pericytes, are essential for organ development and tissue homeostasis through producing paracrine signalling molecules, called angiocrine factors. In the skeletal system, ECs - derived angiocrine factors, combined with bone cells-released angiogenic factors, orchestrate intercellular crosstalk of the bone microenvironment, and the coupling of angiogenesis-to-osteogenesis. Whilst the involvement of angiogenic factors and the blood vessels of the skeleton is relatively well established, the impact of ECs -derived angiocrine factors on bone and cartilage homeostasis is gradually emerging. In this review, we survey ECs - derived angiocrine factors, which are released by endothelial cells of the local microenvironment and by distal organs, and act specifically as regulators of skeletal growth and homeostasis. These may potentially include angiocrine factors with osteogenic property, such as Hedgehog, Notch, WNT, bone morphogenetic protein (BMP), fibroblast growth factor (FGF), insulin-like growth factor (IGF), and platelet-derived growth factor (PDGF). Understanding the versatile mechanisms by which ECs-derived angiocrine factors orchestrate bone and cartilage homeostasis, and pathogenesis, is an important step towards the development of therapeutic potential for skeletal diseases.

Assessment of romosozumab efficacy in the treatment of postmenopausal osteoporosis: Results from a mechanistic PK-PD mechanostat model of bone remodeling

This paper introduces a theoretical framework for the study of the efficacy of romosozumab, a humanized monoclonal antibody targeting sclerostin for the treatment of osteoporosis. We developed a comprehensive mechanistic pharmacokinetic-pharmacodynamic (PK-PD) model of the effect of drug treatment on bone remodeling in postmenopausal osteoporosis (PMO). We utilized a one-compartment PK model to represent subcutaneous injections of romosozumab and subsequent absorption into serum. The PD model is based on a recently-developed bone cell population model describing the bone remodeling process at the tissue scale. The latter accounts for mechanical feedback by incorporating nitric oxide (NO) and sclerostin (Scl) as biochemical feedback molecules. Utilizing a competitive binding model, where Wnt and Scl compete for binding to LRP5/6, allows to regulate anabolic bone remodeling responses. Here, we extended this model with respect to romosozumab binding to sclerostin. For the currently approved monthly injections of 210 mg, the model predicted a 6.59%, 10.38% and 15.25% increase in BMD at the lumbar spine after 6, 12 and 24 months, respectively. These results are in good agreement with the data reported in the literature. Our model is also able to distinguish the bone-site specific drug effects. For instance, at the femoral neck, our model predicts a BMD increase of 3.85% after 12 months of 210 mg injections, which is consistent with literature observations. Finally, our simulations indicate rapid bone loss after treatment discontinuation, indicating that some additional interventions such as use of bisphosphonates are required to maintain bone mass.

In silico experiments of bone remodeling explore metabolic diseases and their drug treatment

Bone structure and function are maintained by well-regulated bone metabolism and remodeling. Although the underlying molecular and cellular mechanisms are now being understood, physiological and pathological states of bone are still difficult to predict due to the complexity of intercellular signaling. We have now developed a novel in silico experimental platform, V-Bone, to integratively explore bone remodeling by linking complex microscopic molecular/cellular interactions to macroscopic tissue/organ adaptations. Mechano-biochemical couplings modeled in V-Bone relate bone adaptation to mechanical loading and reproduce metabolic bone diseases such as osteoporosis and osteopetrosis. V-Bone also enables in silico perturbation on a specific signaling molecule to observe bone metabolic dynamics over time. We also demonstrate that this platform provides a powerful way to predict in silico therapeutic effects of drugs against metabolic bone diseases. We anticipate that these in silico experiments will substantially accelerate research into bone metabolism and remodeling.

Study of the combined effects of PTH treatment and mechanical loading in postmenopausal osteoporosis using a new mechanistic PK-PD model

One of only a few approved and available anabolic treatments for severe osteoporosis is daily injections of PTH (1-34). This drug has a specific dual action which can act either anabolically or catabolically depending on the type of administration, i.e. intermittent or continuous, respectively. In this paper, we present a mechanistic pharmacokinetic-pharmacodynamic model of the action of PTH in postmenopausal osteoporosis. This model accounts for anabolic and catabolic activities in bone remodelling under intermittent and continuous administration of PTH. The model predicts evolution of common bone biomarkers and bone volume fraction (BV/TV) over time. We compared the relative changes in BV/TV resulting from a daily injection of 20 [Formula: see text]g of PTH with experimental data from the literature. Simulation results indicate a site-specific bone gain of 8.66[Formula: see text] (9.4 ± 1.13[Formula: see text]) at the lumbar spine and 3.14[Formula: see text] (2.82 ± 0.72[Formula: see text]) at the femoral neck. Bone gain depends nonlinearly on the administered dose, being, respectively, 0.68[Formula: see text], 3.4[Formula: see text] and 6.16[Formula: see text] for a 10, 20 and 40 [Formula: see text]g PTH dose at the FN over 2 years. Simulations were performed also taking into account a bone mechanical disuse to reproduce elderly frail subjects. The results show that mechanical disuse ablates the effects of PTH and leads to a 1.08% reduction of bone gain at the FN over a 2-year treatment period for the 20 [Formula: see text]g of PTH. The developed model can simulate a range of pathological conditions and treatments in bones including different PTH doses, different mechanical loading environments and combinations. Consequently, the model can be used for testing and generating hypotheses related to synergistic action between PTH treatment and physical activity.

Patient-Specific Bone Multiscale Modelling, Fracture Simulation and Risk Analysis-A Survey

This paper provides a starting point for researchers and practitioners from biology, medicine, physics and engineering who can benefit from an up-to-date literature survey on patient-specific bone fracture modelling, simulation and risk analysis. This survey hints at a framework for devising realistic patient-specific bone fracture simulations. This paper has 18 sections: Section 1 presents the main interested parties; Section 2 explains the organzation of the text; Section 3 motivates further work on patient-specific bone fracture simulation; Section 4 motivates this survey; Section 5 concerns the collection of bibliographical references; Section 6 motivates the physico-mathematical approach to bone fracture; Section 7 presents the modelling of bone as a continuum; Section 8 categorizes the surveyed literature into a continuum mechanics framework; Section 9 concerns the computational modelling of bone geometry; Section 10 concerns the estimation of bone mechanical properties; Section 11 concerns the selection of boundary conditions representative of bone trauma; Section 12 concerns bone fracture simulation; Section 13 presents the multiscale structure of bone; Section 14 concerns the multiscale mathematical modelling of bone; Section 15 concerns the experimental validation of bone fracture simulations; Section 16 concerns bone fracture risk assessment. Lastly, glossaries for symbols, acronyms, and physico-mathematical terms are provided.

Study on the biomechanical responses of the loaded bone in macroscale and mesoscale by multiscale poroelastic FE analysis

BACKGROUND

Bone is a hierarchically structured composite material, and different hierarchical levels exhibit diverse material properties and functions. The stress and strain distribution and fluid flow in bone play an important role in the realization of mechanotransduction and bone remodeling.

METHODS

To investigate the mechanotransduction and fluid behaviors in loaded bone, a multiscale method was developed. Based on poroelastic theory, we established the theoretical and FE model of a segment bone to provide basis for researching more complex bone model. The COMSOL Multiphysics software was used to establish different scales of bone models, and the properties of mechanical and fluid behaviors in each scale were investigated.

RESULTS

FE results correlated very well with analytical in macroscopic scale, and the results for the mesoscopic models were about less than 2% different compared to that in the macro-mesoscale models, verifying the correctness of the modeling. In macro-mesoscale, results demonstrated that variations in fluid pressure (FP), fluid velocity (FV), von Mises stress (VMS), and maximum principal strain (MPS) in the position of endosteum, periosteum, osteon, and interstitial bone and these variations can be considerable (up to 10, 8, 4 and 3.5 times difference in maximum FP, FV, VMS, and MPS between the highest and the lowest regions, respectively). With the changing of Young's modulus (E) in each osteon lamella, the strain and stress concentration occurred in different positions and given rise to microscale spatial variations in the fluid pressure field. The heterogeneous distribution of lacunar-canalicular permeability (klcp) in each osteon lamella had various influence on the FP and FV, but had little effect on VMS and MPS.

CONCLUSION

Based on the idealized model presented in this article, the presence of endosteum and periosteum has an important influence on the fluid flow in bone. With the hypothetical parameter values in osteon lamellae, the bone material parameters have effect on the propagation of stress and fluid flow in bone. The model can also incorporate alternative material parameters obtained from different individuals. The suggested method is expected to provide dependable biological information for better understanding the bone mechanotransduction and signal transduction.

Computational model of the dual action of PTH - Application to a rat model of osteoporosis

This paper presents a pharmacokinetic/pharmacodynamic (PK/PD) model of the action of PTH(1-34) on bone modelling and remodelling, developed for quantitatively investigating the dose- and administration pattern-dependency of the bone tissue response to this drug. Firstly, a PK model of PTH(1-34) was developed, accounting for administration via subcutaneous injections. Subsequently, the PK model was coupled to a (mechanistic) bone cell population model of bone modelling and remodelling, taking into account the effects of PTH(1-34) on the differentiation of lining cells into active osteoblasts, on the apoptosis of active osteoblasts, and on proliferation of osteoblast precursors, as well as on the key regulatory pathways of bone cell activities. Numerical simulations show that the coupled PK/PD model is able to distinguish between continuous and intermittent administration patterns of PTH(1-34), in terms of yielding both catabolic bone responses (if drug administration is carried out continuously) and anabolic bone responses (if drug administration is carried out intermittently). The model also features a non-linear relation between bone gain and drug dose (as known from experiments); doubling the dose from 80 μg/kg/day to 160 μg/kg/day induced a 1.3-fold increase of the bone volume-to-total volume ratio. Furthermore, the model presented in this paper confirmed that bone modelling represents an essential mechanism of the anabolic response of bone to PTH(1-34) administration in rat models, and that the large amount of bone formation observed in such models cannot be explained via remodelling alone.

A review of recent developments in mathematical modeling of bone remodeling

In this article, we summarize the developments in the mathematical modeling of the mechanics of bone and related biological phenomena. We will devote special attention to the results of the last 10–15 years, although we will cover some relevant classical work to better frame the more recent researches. We will propose a division of the literature based on the main aim of the model (mechanical/biomathematical) and the type of biological phenomena considered (stimulus, growth, cell population dynamics). Finally, we will suggest some possible directions for future investigations.

A surrogate mechanostatistical microstructural model to inform whole hip cortical bone remodelling

Multiscale models of bone remodelling which incorporate clinically important features at a wide range of spatial scales are computationally intractable. In addition, there is a large gap in existing bone remodelling literature which deals with Haversian-level cortical bone remodelling behaviour and its effect on whole organ bone strength. In this study, we present an integration between continuum mechanics and surrogate modelling as an efficient circumvention to the problem of intractability, linked with an intricate 3D cortical bone remodelling algorithm which captures realistic microanatomical features at the Haversian level. The surrogate model predicted clear differentiation in long-term bone strength changes in clinically significant cases with extremely high computational efficiency.

Mathematical and computational models for bone tissue engineering in bioreactor systems

Research into cellular engineered bone grafts offers a promising solution to problems associated with the currently used auto- and allografts. Bioreactor systems can facilitate the development of functional cellular bone grafts by augmenting mass transport through media convection and shear flow-induced mechanical stimulation. Developing successful and reproducible protocols for growing bone tissue in vitro is dependent on tuning the bioreactor operating conditions to the specific cell type and graft design. This process, largely reliant on a trial-and-error approach, is challenging, time-consuming and expensive. Modelling can streamline the process by providing further insight into the effect of the bioreactor environment on the cell culture, and by identifying a beneficial range of operational settings to stimulate tissue production. Models can explore the impact of changing flow speeds, scaffold properties, and nutrient and growth factor concentrations. Aiming to act as an introductory reference for bone tissue engineers looking to direct their experimental work, this article presents a comprehensive framework of mathematical models on various aspects of bioreactor bone cultures and overviews modelling case studies from literature.

Bone remodeling: A tissue-level process emerging from cell-level molecular algorithms

The human skeleton undergoes constant remodeling throughout the lifetime. Processes occurring on microscopic and molecular scales degrade bone and replace it with new, fully functional tissue. Multiple bone remodeling events occur simultaneously, continuously and independently throughout the body, so that the entire skeleton is completely renewed about every ten years.Bone remodeling is performed by groups of cells called Bone Multicellular Units (BMU). BMUs consist of different cell types, some specialized in the resorption of old bone, others encharged with producing new bone to replace the former. These processes are tightly regulated so that the amount of new bone produced is in perfect equilibrium with that of old bone removed, thus maintaining bone microscopic structure.To date, many regulatory molecules involved in bone remodeling have been identified, but the precise mechanism of BMU operation remains to be fully elucidated. Given the complexity of the signaling pathways already known, one may question whether such complexity is an inherent requirement of the process or whether some subset of the multiple constituents could fulfill the essential role, leaving functional redundancy to serve an alternative safety role. We propose in this work a minimal model of BMU function that involves a limited number of signals able to account for fully functional BMU operation. Our main assumptions were i) at any given time, any cell within a BMU can select only one among a limited choice of decisions, i.e. divide, die, migrate or differentiate, ii) this decision is irreversibly determined by depletion of an appropriate internal inhibitor and iii) the dynamics of any such inhibitor are coupled to that of specific external mediators, such as hormones, cytokines, growth factors. It was thus shown that efficient BMU operation manifests as an emergent process, which results from the individual and collective decisions taken by cells within the BMU unit in the absence of any external planning.

X-ray physics-based CT-to-composition conversion applied to a tissue engineering scaffold, enabling multiscale simulation of its elastic behavior

Nowadays, the assessment of the mechanical competence of tissue engineering scaffolds based on computer simulations is a well-accepted technology. Typically, such simulations are performed by means of the Finite Element (FE) method, with the underlying structural model being created based on micro-computed tomography (microCT). Here, this analysis modality is applied to a new, ternary composite, consisting of PHBV, i.e. poly(3-hydroxybutyrate-co-3-hydroxyvalerate), PLGA, i.e. poly(lactic-co-glycolide), as well as of TCP, i.e. tricalcium phosphate hydrate. The studied scaffold structure is made up by fibers of this new composite material, manufactured by means of the rapid prototyping method. The data collected from microCT is utilized for adequately defining the mechanical properties of the FE model. In particular, the three-dimensional field of grey values is interpreted in terms of the underlying field of attenuation coefficients, taking into account the photon energy employed in microCT imaging, eventually allowing for calculation of the three-dimensionally distributed, voxel-specific composition of the studied material. For the sake of keeping the FE simulations as efficient as possible, groups of voxels are combined into one finite element; the grey value of the latter is obtained by volume averaging. Employing a two-step micromechanical homogenization scheme, the experimentally accessible stiffness of the three constituents (PHBV, PLGA, and TCP) is then, finite element by finite element, upscaled to the composition-dependent stiffness of the composite material. The plausibility and adequacy of the FE model is demonstrated by simulating the effects of uniaxial compression on the scaffold structure, in terms of resulting stress and strain fields, highlighting the importance of the fiber junctions (as they are the mechanically most stressed regions), and that neglecting the material heterogeneity would lead to a potentially significant underestimation of stresses and strains. Finally, a comparison is made of the employed analysis modality of microCT data with a previously pursued, simplified analysis strategy, highlighting the conceptual superiority of the former, and pointing out the application limits of the latter.