Trabecular bone remodelling simulation considering osteocytic response to fluid-induced shear stress

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

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

Bone Remodeling Process Based on Hydrostatic and Deviatoric Strain Mechano-Sensing

A macroscopic continuum model intended to provide predictions for the remodeling process occurring in bone tissue is proposed. Specifically, we consider a formulation in which two characteristic stiffnesses, namely the bulk and shear moduli, evolve independently to adapt the hydrostatic and deviatoric response of the bone tissue to environmental changes. The formulation is deliberately simplified, aiming at constituting a preliminary step toward a more comprehensive modeling approach. The evolutive process for describing the functional adaptation of the two stiffnesses is proposed based on an energetic argument. Numerical experiments reveal that it is possible to model the bone remodeling process with a different evolution for more than one material parameter, as usually done. Moreover, the results motivate further investigations into the subject.

Three-dimensional topology optimization model to simulate the external shapes of bone

Elucidation of the mechanism by which the shape of bones is formed is essential for understanding vertebrate development. Bones support the body of vertebrates by withstanding external loads, such as those imposed by gravity and muscle tension. Many studies have reported that bone formation varies in response to external loads. An increased external load induces bone synthesis, whereas a decreased external load induces bone resorption. This relationship led to the hypothesis that bone shape adapts to external load. In fact, by simulating this relationship through topology optimization, the internal trabecular structure of bones can be successfully reproduced, thereby facilitating the study of bone diseases. In contrast, there have been few attempts to simulate the external structure of bones, which determines vertebrate morphology. However, the external shape of bones may be reproduced through topology optimization because cells of the same type form both the internal and external structures of bones. Here, we constructed a three-dimensional topology optimization model to attempt the reproduction of the external shape of teleost vertebrae. In teleosts, the internal structure of the vertebral bodies is invariable, exhibiting an hourglass shape, whereas the lateral structure supporting the internal structure differs among species. Based on the anatomical observations, we applied different external loads to the hourglass-shaped part. The simulations produced a variety of three-dimensional structures, some of which exhibited several structural features similar to those of actual teleost vertebrae. In addition, by adjusting the geometric parameters, such as the width of the hourglass shape, we reproduced the variation in the teleost vertebrae shapes. These results suggest that a simulation using topology optimization can successfully reproduce the external shapes of teleost vertebrae. By applying our topology optimization model to various bones of vertebrates, we can understand how the external shape of bones adapts to external loads.

In this paper, we developed a computational method to investigate the relationship between three-dimensional bone shape and external loads imposed on bones. Many studies report that bone formation varies in response to external loads. An increased external load induces bone synthesis, whereas a decreased external load induces bone resorption. This relationship led to the hypothesis that the shape of bones adapts to external load. However, it remains unclear whether this hypothesis can explain the shape of bones. Here, we constructed a three-dimensional mathematical model that imitates the cellular activities of bone formation to attempt the reproduction of the shape of teleost vertebrae. In teleosts, the shape of the vertebrae differs among the species. We set the multiple types of external load conditions in the simulations and compared the simulation results with different teleost vertebrae. The produced structures that can resist the deformation of the surrounding tissues exhibited multiple structural features similar to the vertebrae of several teleost species. This result shows that the formation of bone shape can be explained by the adaptation to external load.

Maintaining Bone Health in the Lumbar Spine: Routine Activities Alone Are Not Enough

Public health organisations typically recommend a minimum amount of moderate intensity activities such as walking or cycling for two and a half hours a week, combined with some more demanding physical activity on at least 2 days a week to maintain a healthy musculoskeletal condition. For populations at risk of bone loss in the lumbar spine, these guidelines are particularly relevant. However, an understanding of how these different activities are influential in maintaining vertebral bone health is lacking. A predictive structural finite element modelling approach using a strain-driven algorithm was developed to study mechanical stimulus and bone adaptation in the lumbar spine under various physiological loading conditions. These loading conditions were obtained with a previously developed full-body musculoskeletal model for a range of daily living activities representative of a healthy lifestyle. Activities of interest for the simulations include moderate intensity activities involving limited spine movements in all directions such as, walking, stair ascent and descent, sitting down and standing up, and more demanding activities with large spine movements during reaching and lifting tasks. For a combination of moderate and more demanding activities, the finite element model predicted a trabecular and cortical bone architecture representative of a healthy vertebra. When more demanding activities were removed from the simulations, areas at risk of bone degradation were observed at all lumbar levels in the anterior part of the vertebral body, the transverse processes and the spinous process. Moderate intensity activities alone were found to be insufficient in providing a mechanical stimulus to prevent bone degradation. More demanding physical activities are essential to maintain bone health in the lumbar spine.

Remodelling of trabecular bone in human distal tibia: A model based on an in-vivo HR-pQCT study

An abnormal remodelling process of bones can lead to various bone disorders, such as osteoporosis, making them prone to fracture. Simulations of load-induced remodelling of trabecular bone were used to investigate its response to mechanical signal. However, the role of mechanostat in trabecular-bone remodelling has not yet been investigated in simulations underpinned by a longitudinal in-vivo study in humans. In this work, a finite-element model based on a 6-month longitudinal in-vivo HR-pQCT study was developed and validated to investigate the effect of mechanical stimuli on bone remodelling. The simulated changes in microstructural parameters and density of trabecular bone were compared with respective experimental results. A maximum principal strain (MPS) and a maximum principal strain gradient (∇MPS) were used as mechanical signals to drive a five-stage mechanostat remodelling model, including additional over-strain and damage stages. It was found that the density distribution varied with the studied mechanical signals, along with decreasing with time levels of bone volume fraction BV/TV, trabecular thickness Tb.Th and bone surface area Tb.BS as well as increased trabecular separation Tb.Sp. Among these parameters, BV/TV and Tb.Th together with the bone-remodelling parameters from the MPS model demonstrated a significant correlation with the experimental data. The developed model provides a good foundation for further development and investigation of the relationships between mechanical loading and human-bone microarchitecture.

Mathematical modeling of bone in-growth into undegradable porous periodic scaffolds under mechanical stimulus

Undegradable scaffolds, as a key element in bone tissue engineering, prevail in the present clinical applications, and the bone in-growth into such scaffolds under mechanical stimulus is an important issue to evaluate the bone-repair effect. This work aims to develop a mathematical framework to investigate the effect of mechanical stimulus on the bone in-growth into undegradable scaffolds. First, the osteoclast and osteoblast activities were coupled by their autocrine and paracrine effects. Second, the mechanical stimulus was empirically incorporated into the coupling cell activities on the basis of experimental observations. Third, the effect of mechanical stimulus including intensity and duration on the bone in-growth process was numerically studied; moreover, the homeostasis of scaffold–bone system under the mechanical stimulus was also treated. The results showed that the numbers of osteoblasts and osteoclasts in the scaffold–bone system tended to constants representing the system homeostasis. Both the mechanical intensity and duration optimized the final bone formation. The numerical results of the bone formation were comparable to the experimental results in rats. The findings from this modeling study could be used to explain many physiological phenomena and clinical observations. The developed model integrates both cell and tissue scales, which can be used as a platform to investigate bone remodeling under mechanical stimulus.

Microscale poroelastic metamodel for efficient mesoscale bone remodelling simulations

Bone functional tissue adaptation is a multiaspect physiological process driven by interrelated mechanical and biological stimuli which requires the combined activity of osteoclasts and osteoblasts. In previous work, the authors developed a phenomenological mesoscale structural modelling approach capable of predicting internal structure of the femur based on daily activity loading, which relied on the iterative update of the cross-sectional areas of truss and shell elements representative of trabecular and cortical bones, respectively. The objective of this study was to introduce trabecular reorientation in the phenomenological model at limited computational cost. To this aim, a metamodel derived from poroelastic microscale continuum simulations was used to predict the functional adaptation of a simplified proximal structural femur model. Clear smooth trabecular tracts are predicted to form in the regions corresponding to the main trabecular groups identified in literature, at minimal computational cost.

Informing phenomenological structural bone remodelling with a mechanistic poroelastic model

Studies suggest that fluid motion in the extracellular space may be involved in the cellular mechanosensitivity at play in the bone tissue adaptation process. Previously, the authors developed a mesoscale predictive structural model of the femur using truss elements to represent trabecular bone, relying on a phenomenological strain-based bone adaptation algorithm. In order to introduce a response to bending and shear, the authors considered the use of beam elements, requiring a new formulation of the bone adaptation drivers. The primary goal of the study presented here was to isolate phenomenological drivers based on the results of a mechanistic approach to be used with a beam element representation of trabecular bone in mesoscale structural modelling. A single-beam model and a microscale poroelastic model of a single trabecula were developed. A mechanistic iterative adaptation algorithm was implemented based on fluid motion velocity through the bone matrix pores to predict the remodelled geometries of the poroelastic trabecula under 42 different loading scenarios. Regression analyses were used to correlate the changes in poroelastic trabecula thickness and orientation to the initial strain outputs of the beam model. Linear (R(2) > 0.998) and third-order polynomial (R(2) > 0.98) relationships were found between change in cross section and axial strain at the central axis, and between beam reorientation and ratio of bending strain to axial strain, respectively. Implementing these relationships into the phenomenological predictive algorithm for the mesoscale structural femur has the potential to produce a model combining biofidelic structure and mechanical behaviour with computational efficiency.

Quantifying the osteocyte network in the human skeleton

Osteocytes form an extensive cellular network throughout the hard tissue matrix of the skeleton, which is known to regulate skeletal structure. However due to limitations in imaging techniques, the magnitude and complexity of this network remain undefined. We have used data from recent papers obtained by new imaging techniques, in order to estimate absolute and relative quantities of the human osteocyte network and form a more complete understanding of the extent and nature of this network. We estimate that the total number of osteocytes within the average adult human skeleton is ~42 billion and that the total number of osteocyte dendritic projections from these cells is ~3.7 trillion. Based on prior measurements of canalicular density and a mathematical model of osteocyte dendritic process branching, we calculate that these cells form a total of 23 trillion connections with each other and with bone surface cells. We estimate the total length of all osteocytic processes connected end-to-end to be 175,000 km. Furthermore, we calculate that the total surface area of the lacuno-canalicular system is 215 m(2). However, the residing osteocytes leave only enough space for 24 mL of extracellular fluid. Calculations based on measurements in lactation-induced murine osteocytic osteolysis indicate a potential total loss of ~16,000 mm(3) (16 mL) of bone by this process in the human skeleton. Finally, based on the average speed of remodelling in the adult, we calculate that 9.1 million osteocytes are replenished throughout the skeleton on a daily basis, indicating the dynamic nature of the osteocyte network. We conclude that the osteocyte network is a highly complex communication network, and is much more vast than commonly appreciated. It is at the same order of magnitude as current estimates of the size of the neural network in the brain, even though the formation of the branched network differs between neurons and osteocytes. Furthermore, continual replenishment of large numbers of osteocytes in the process of remodelling allows therapeutic changes to the continually renewed osteoblast population to be rapidly incorporated into the skeleton.

Biomimetic approaches in bone tissue engineering: Integrating biological and physicomechanical strategies

The development of responsive biomaterials capable of demonstrating modulated function in response to dynamic physiological and mechanical changes in vivo remains an important challenge in bone tissue engineering. To achieve long-term repair and good clinical outcomes, biologically responsive approaches that focus on repair and reconstitution of tissue structure and function through drug release, receptor recognition, environmental responsiveness and tuned biodegradability are required. Traditional orthopedic materials lack biomimicry, and mismatches in tissue morphology, or chemical and mechanical properties ultimately accelerate device failure. Multiple stimuli have been proposed as principal contributors or mediators of cell activity and bone tissue formation, including physical (substrate topography, stiffness, shear stress and electrical forces) and biochemical factors (growth factors, genes or proteins). However, optimal solutions to bone regeneration remain elusive. This review will focus on biological and physicomechanical considerations currently being explored in bone tissue engineering.

Strain energy density gradients in bone marrow predict osteoblast and osteoclast activity: a finite element study

Huiskes et al. hypothesized that mechanical strains sensed by osteocytes residing in trabecular bone dictate the magnitude of load-induced bone formation. More recently, the mechanical environment in bone marrow has also been implicated in bone׳s response to mechanical stimulation. In this study, we hypothesize that trabecular load-induced bone formation can be predicted by mechanical signals derived from an integrative µFE model, incorporating a description of both the bone and marrow phase. Using the mouse tail loading model in combination with in vivo micro-computed tomography (µCT) we tracked load induced changes in the sixth caudal vertebrae of C57BL/6 mice to quantify the amount of newly mineralized and eroded bone volumes. To identify the mechanical signals responsible for adaptation, local morphometric changes were compared to micro-finite element (µFE) models of vertebrae prior to loading. The mechanical parameters calculated were strain energy density (SED) on trabeculae at bone forming and resorbing surfaces, SED in the marrow at the boundary between bone forming and resorbing surfaces, along with SED in the trabecular bone and marrow volumes. The gradients of each parameter were also calculated. Simple regression analysis showed mean SED gradients in the trabecular bone matrix to significantly correlate with newly mineralized and eroded bone volumes R(2)=0.57 and 0.41, respectively, p<0.001). Nevertheless, SED gradients in the marrow were shown to be the best predictor of osteoblastic and osteoclastic activity (R(2)=0.83 and 0.60, respectively, p<0.001). These data suggest that the mechanical environment of the bone marrow plays a significant role in determining osteoblast and osteoclast activity.

Connecting mechanics and bone cell activities in the bone remodeling process: an integrated finite element modeling

Bone adaptation occurs as a response to external loadings and involves bone resorption by osteoclasts followed by the formation of new bone by osteoblasts. It is directly triggered by the transduction phase by osteocytes embedded within the bone matrix. The bone remodeling process is governed by the interactions between osteoblasts and osteoclasts through the expression of several autocrine and paracrine factors that control bone cell populations and their relative rate of differentiation and proliferation. A review of the literature shows that despite the progress in bone remodeling simulation using the finite element (FE) method, there is still a lack of predictive models that explicitly consider the interaction between osteoblasts and osteoclasts combined with the mechanical response of bone. The current study attempts to develop an FE model to describe the bone remodeling process, taking into consideration the activities of osteoclasts and osteoblasts. The mechanical behavior of bone is described by taking into account the bone material fatigue damage accumulation and mineralization. A coupled strain-damage stimulus function is proposed, which controls the level of autocrine and paracrine factors. The cellular behavior is based on Komarova et al.'s (2003) dynamic law, which describes the autocrine and paracrine interactions between osteoblasts and osteoclasts and computes cell population dynamics and changes in bone mass at a discrete site of bone remodeling. Therefore, when an external mechanical stress is applied, bone formation and resorption is governed by cells dynamic rather than adaptive elasticity approaches. The proposed FE model has been implemented in the FE code Abaqus (UMAT routine). An example of human proximal femur is investigated using the model developed. The model was able to predict final human proximal femur adaptation similar to the patterns observed in a human proximal femur. The results obtained reveal complex spatio-temporal bone adaptation. The proposed FEM model gives insight into how bone cells adapt their architecture to the mechanical and biological environment.

Interstitial fluid flow in canaliculi as a mechanical stimulus for cancellous bone remodeling: in silico validation

Cancellous bone has a dynamic 3-dimensional architecture of trabeculae, the arrangement of which is continually reorganized via bone remodeling to adapt to the mechanical environment. Osteocytes are currently believed to be the major mechanosensory cells and to regulate osteoclastic bone resorption and osteoblastic bone formation in response to mechanical stimuli. We previously developed a mathematical model of trabecular bone remodeling incorporating the possible mechanisms of cellular mechanosensing and intercellular communication in which we assumed that interstitial fluid flow activates the osteocytes to regulate bone remodeling. While the proposed model has been validated by the simulation of remodeling of a single trabecula, it remains unclear whether it can successfully represent in silico the functional adaptation of cancellous bone with its multiple trabeculae. In the present study, we demonstrated the response of cancellous bone morphology to uniaxial or bending loads using a combination of our remodeling model with the voxel finite element method. In this simulation, cancellous bone with randomly arranged trabeculae remodeled to form a well-organized architecture oriented parallel to the direction of loading, in agreement with the previous simulation results and experimental findings. These results suggested that our mathematical model for trabecular bone remodeling enables us to predict the reorganization of cancellous bone architecture from cellular activities. Furthermore, our remodeling model can represent the phenomenological law of bone transformation toward a locally uniform state of stress or strain at the trabecular level.

Design Optimization of Scaffold Microstructures Using Wall Shear Stress Criterion Towards Regulated Flow-Induced Erosion

Tissue scaffolds aim to provide a cell-friendly biomechanical environment for facilitating cell growth. Existing studies have shown significant demands for generating a certain level of wall shear stress (WSS) on scaffold microstructural surfaces for promoting cellular response and attachment efficacy. Recently, its role in shear-induced erosion of polymer scaffold has also drawn increasing attention. This paper proposes a bi-directional evolutionary structural optimization (BESO) approach for design of scaffold microstructure in terms of the WSS uniformity criterion, by downgrading highly-stressed solid elements into fluidic elements and/or upgrading lowly-stressed fluidic elements into solid elements. In addition to this, a computational model is presented to simulate shear-induced erosion process. The effective stiffness and permeability of initial and optimized scaffold microstructures are characterized by the finite element based homogenization technique to quantify the variations of mechanical properties of scaffold during erosion. The illustrative examples show that a uniform WSS is achieved within the optimized scaffold microstructures, and their architectural and biomechanical features are maintained for a longer lifetime during shear-induced erosion process. This study provides a mathematical means to the design optimization of cellular biomaterials in terms of the WSS criterion towards controllable shear-induced erosion.

Experimental and finite element analysis of the mouse caudal vertebrae loading model: prediction of cortical and trabecular bone adaptation

In this study, we attempt to predict cortical and trabecular bone adaptation in the mouse caudal vertebrae loading model using knowledge of bone's local mechanical environment at the onset of loading. In a previous study, we demonstrated appreciable 25.9 and 11% increases in both trabecular and cortical bone volume density, respectively, when subjecting the fifth caudal vertebrae (C5) of C57BL/6 (B6) mice to an acute loading regime (amplitude of 8N, 3000 cycles, 10 Hz, 3 times a week for 4 weeks). We have also established a validated finite element (FE) model of the C5 vertebra using micro-computed tomography (micro-CT), which characterizes, in 3D, the micro-mechanical strains present in both cortical and trabecular compartments due to the applied loads. To investigate the relationship between load-induced bone adaptation and mechanical strains in-vivo and in-silico data sets were compared. Using data from the previous cross-sectional study, we divided cortical and trabecular compartments into 15 subregions and determined, for each region, a bone formation parameter ΔBV/BS (a cross-sectional measure of the bone volume added to cortical and trabecular surfaces following the described loading regime). Linear regression was then used to correlate mean regional values of ΔBV/BS with mean values of mechanical strains derived from the FE models which were similarly regionalized. The mechanical parameters investigated were strain energy density (SED), the orthogonal strains (e (x), e (y), e (z)) and the three shear strains (e (xy), e (yz), e (zx)). For cortical regions, regression analysis showed SED to correlate extremely well with ΔBV/BS (R (2) = 0.82) and e (z) (R (2) = 0.89). Furthermore, SED was found to predict expansion of the cortical shell correlating significantly with the regional percentage increases in cortical tissue volume (R (2) = 0.92), cortical marrow volume (R (2) = 0.91) and cortical thickness (R (2) = 0.56). For trabecular regions, FE parameters were found not to correlate with load-induced trabecular bone morphology. These results indicate that load-induced cortical morphology can be predicted from population data, whereas the prediction of trabecular morphology requires subject-specific micro- architecture.

Trabecular bone remodelling simulated by a stochastic exchange of discrete bone packets from the surface

Human bone is constantly renewed through life via the process of bone remodelling, in which individual packets of bone are removed by osteoclasts and replaced by osteoblasts. Remodelling is mechanically controlled, where osteocytes embedded within the bone matrix are thought to act as mechanical sensors. In this computational work, a stochastic model for bone remodelling is used in which the renewal of bone material occurs by exchange of discrete bone packets. We tested different hypotheses of how the mechanical stimulus for bone remodelling is integrated by osteocytes and sent to actor cells on the bone's surface. A collective (summed) signal from multiple osteocytes as opposed to an individual (maximal) signal from a single osteocyte was found to lead to lower inner porosity and surface roughness of the simulated bone structure. This observation can be interpreted in that collective osteocyte signalling provides an effective surface tension to the remodelling process. Furthermore, the material heterogeneity due to remodelling was studied on a network of trabeculae. As the model is discrete, the age of individual bone packets can be monitored with time. The simulation results were compared with experimental data coming from quantitative back scattered electron imaging by transforming the information about the age of the bone packet into a mineral content. Discrepancies with experiments indicate that osteoclasts preferentially resorb low mineralized, i.e. young, bone at the bone's surface.

Effects of loading frequency on the functional adaptation of trabeculae predicted by bone remodeling simulation

The process of bone remodeling is regulated by metabolic activities of many bone cells. While osteoclasts and osteoblasts are responsible for bone resorption and formation, respectively, activities of these cells are believed to be controlled by a mechanosensory system of osteocytes embedded in the extracellular bone matrix. Several experimental and theoretical studies have suggested that the strain-derived interstitial fluid flow in lacuno-canalicular porosity serves as the prime mover for bone remodeling. Previously, we constructed a mathematical model for trabecular bone remodeling that interconnects the microscopic cellular activities with the macroscopic morphological changes in trabeculae through the mechanical hierarchy. This model assumes that fluid-induced shear stress acting on osteocyte processes is a driving force for bone remodeling. The validity of this model has been demonstrated with a remodeling simulation using a two-dimensional trabecular model. In this study, to investigate the effects of loading frequency, which is thought to be a significant mechanical factor in bone remodeling, we simulated morphological changes of a three-dimensional single trabecula under cyclic uniaxial loading with various frequencies. The results of the simulation show the trabecula reoriented to the loading direction with the progress of bone remodeling. Furthermore, as the imposed loading frequency increased, the diameter of the trabecula in the equilibrium state was enlarged by remodeling. These results indicate that our simulation model can successfully evaluate the relationship between loading frequency and trabecular bone remodeling.

In silico design of treatment strategies in wound healing and bone fracture healing

Wound and bone fracture healing are natural repair processes initiated by trauma. Over the last decade, many mathematical models have been established to investigate the healing processes in silico, in addition to ongoing experimental work. In recent days, the focus of the mathematical models has shifted from simulation of the healing process towards simulation of the impaired healing process and the in silico design of treatment strategies. This review describes the most important causes of failure of the wound and bone fracture healing processes and the experimental models and methods used to investigate and treat these impaired healing cases. Furthermore, the mathematical models that are described address these impaired healing cases and investigate various therapeutic scenarios in silico. Examples are provided to illustrate the potential of these in silico experiments. Finally, limitations of the models and the need for and ability of these models to capture patient specificity and variability are discussed.