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.     

Functional adaptation of cancellous bone in human proximal femur predicted by trabecular surface remodeling simulation toward uniform stress state

Two-dimensional simulation of trabecular surface remodeling was conducted for a human proximal femur to investigate the structural change of cancellous bone toward a uniform stress state. Considering that a local mechanical stimulus plays an important role in cellular activities in bone remodeling, local stress nonuniformity was assumed to drive trabecular structural change to seek a uniform stress state. A large-scale pixel-based finite element model was used to simulate structural changes of individual trabeculae over the entire bone. As a result, the initial structure of trabeculae changed from isotropic to anisotropic due to trabecular microstructural changes caused by surface remodeling according to the mechanical environment in the proximal femur. Under a single-loading condition, it was shown that the apparent structural property evaluated by fabric ellipses corresponded to the apparent stress state in cancellous bone. As is observed in the actual bone, a distributed trabecular structure was obtained under a multiple-loading condition. Through these studies, it was concluded that trabecular surface remodeling toward a local uniform stress state at the trabecular level could naturally bring about functional adaptation phenomenon at the apparent tissue level. The proposed simulation model would be capable of providing insight into the hierarchical mechanism of trabecular surface remodeling at the microstructural level up to the apparent tissue level.     

Computer simulation of trabecular remodeling in human proximal femur using large-scale voxel FE models: Approach to understanding Wolff's law

Ever since Julius Wolff proposed the law of bone transformation in the 19th century, it has been widely known that the trabecular structure of cancellous bone adapts functionally to the loading environment. To understand the mechanism of Wolff's law, a three-dimensional (3D) computer simulation of trabecular structural changes due to surface remodeling was performed for a human proximal femur. A large-scale voxel finite element model was constructed to simulate the structural changes of individual trabeculae over the entire cancellous region. As a simple remodeling model that considers bone cellular activities regulated by the local mechanical environment, nonuniformity of local stress was assumed to drive the trabecular surface remodeling to seek a uniform stress state. Simulation results demonstrated that cell-scale (～10 μm) remodeling in response to mechanical stimulation created complex 3D trabecular structures of the entire bone-scale (～10 cm), as illustrated in the reference of Wolff. The bone remodeling reproduced the characteristic anisotropic structure in the coronal cross section and the isotropic structures in other cross sections. The principal values and axes of a structure characterized by fabric ellipsoids corresponded to those of the apparent stress of the structure. The proposed large-scale computer simulation indicates that in a complex mechanical environment of a hierarchical bone structure of over 10⁴ length scale (from ～10 μm to ～10 cm), a simple remodeling at the cellular/trabecular levels creates a highly complex and functional trabecular structure, as characterized by bone density and orientation.     

Changes in the Fabric and Compliance Tensors of Cancellous Bone due to Trabecular Surface Remodeling,Predicted by a Digital Image-based Model

Fabric and compliance tensors of a cube of cancellous bone with a complicated three-dimensional trabecular structure were obtained for trabecular surface remodeling by using a digital image-based model combined with a large-scale finite element method. Using mean intercept length and a homogenization method, the fabric and compliance tensors were determined for the trabecular structure obtained in the computer remodeling simulation. The tensorial quantities obtained indicated that anisotropic structural changes occur in cancellous bone adapting to the compressive loading condition. There were good correlations between the fabric tensor, bone volume fraction, and compliance tensor in the remodeling process. The result demonstrates that changes in the structural and mechanical properties of cancellous bone are essentially anisotropic and should be expressed by tensorial quantities.     

Changes in the fabric and compliance tensors of cancellous bone due to trabecular surface remodeling, predicted by a digital image-based model

Fabric and compliance tensors of a cube of cancellous bone with a complicated three-dimensional trabecular structure were obtained for trabecular surface remodeling by using a digital image-based model combined with a large-scale finite element method. Using mean intercept length and a homogenization method, the fabric and compliance tensors were determined for the trabecular structure obtained in the computer remodeling simulation. The tensorial quantities obtained indicated that anisotropic structural changes occur in cancellous bone adapting to the compressive loading condition. There were good correlations between the fabric tensor, bone volume fraction, and compliance tensor in the remodeling process. The result demonstrates that changes in the structural and mechanical properties of cancellous bone are essentially anisotropic and should be expressed by tensorial quantities.     

Comparison of the linear finite element prediction of deformation and strain of human cancellous bone to 3D digital volume correlation measurements

The mechanical properties of cancellous bone and the biological response of the tissue to mechanical loading are related to deformation and strain in the trabeculae during function. Due to the small size of trabeculae, their motion is difficult to measure. To avoid the need to measure trabecular motions during loading the finite element method has been used to estimate trabecular level mechanical deformation. This analytical approach has been empirically successful in that the analytical models are solvable and their results correlate with the macroscopically measured stiffness and strength of bones. The present work is a direct comparison of finite element predictions to measurements of the deformation and strain at near trabecular level. Using the method of digital volume correlation, we measured the deformation and calculated the strain at a resolution approaching the trabecular level for cancellous bone specimens loaded in uniaxial compression. Smoothed results from linearly elastic finite element models of the same mechanical tests were correlated to the empirical three-dimensional (3D) deformation in the direction of loading with a coefficient of determination as high as 97% and a slope of the prediction near one. However, real deformations in the directions perpendicular to the loading direction were not as well predicted by the analytical models. Our results show, that the finite element modeling of the internal deformation and strain in cancellous bone can be accurate in one direction but that this does not ensure accuracy for all deformations and strains.     

The discrete nature of trabecular bone microarchitecture affects implant stability

Small endosseous implants, such as screws, are important components of modern orthopedics and dentistry. Hence they have to reliably fulfill a variety of requirements, which makes the development of such implants challenging. Finite element analysis is a widely used computational tool used to analyze and optimize implant stability in bone. For these purposes, bone is generally modeled as a continuum material. However, bone failure and bone adaptation processes are occurring at the discrete level of individual trabeculae; hence the assessment of stresses and strains at this level is relevant. Therefore, the aim of the present study was to investigate how peri-implant strain distribution and load transfer between implant and bone are affected by the continuum assumption. We performed a computational study in which cancellous screws were inserted in continuum and discrete models of trabecular bone; axial loading was simulated. We found strong differences in bone-implant stiffness between the discrete and continuum bone model. They depended on bone density and applied boundary conditions. Furthermore, load transfer from the screw to the surrounding bone differed strongly between the continuum and discrete models, especially for low-density bone. Based on our findings we conclude that continuum bone models are of limited use for finite element analysis of peri-implant mechanical loading in trabecular bone when a precise quantification of peri-implant stresses and strains is required. Therefore, for the assessment and improvement of trabecular bone implants, finite element models which accurately represent trabecular microarchitecture should be used.     

Age-related changes in microarchitecture and mineralization of cancellous bone in the porcine mandibular condyle

The mandibular condyle is considered a good model for developing cancellous bone because of its rapid growth and high rate of remodeling. The aim of the present study was to analyze the simultaneous changes in microarchitecture and mineralization of cancellous bone during development in a three-dimensional fashion. Eight mandibular condyles of pigs aged 8 weeks prepartum to 108 weeks postpartum were scanned using microCT with an isotropic spatial resolution of 10 microm. The number of trabeculae decreased during development, whereas both the trabecular thickness and the distance between the trabeculae increased. The bone surface to volume ratio decreased during development, possibly limiting the amount of (re)modeling. Both the mean degree of mineralization and intratrabecular differences in mineralization between the surfaces and cores of trabecular elements increased during development. The trabecular surfaces were more highly mineralized in the older condyles compared to the younger ones. Together with the observed decrease in the relative size of trabecular surface, this finding suggests a decrease in (re)modeling activity during development. In accordance with the general growth and development of the pig, it was concluded that most developmental changes in cancellous bone occur until the age of 40 weeks postpartum.     

Trabecular bone's mechanical properties are affected by its non-uniform mineral distribution

The bone remodeling process takes place at the surface of trabeculae and results in a non-uniform mineral distribution. This will affect the mechanical properties of cancellous bone, because the properties of bone tissue depend on its mineral content. We investigated how large this effect is by simulating several non-uniform mineral distributions in 3D finite element models of human trabecular bone and calculating the apparent stiffness of these models. In the ‘linear model’ we assumed a linear relation between mineral content and Young's modulus of the tissue. In the ‘exponential model’ we included an empirical exponential relation in the model. When the linear model was used the mineral distribution slightly changed the apparent stiffness, the difference varied between an 8% decrease and a 4% increase compared to the uniform model with the same BMD. The exponential model resulted in up to 20% increased apparent stiffness in the main load-bearing direction. A thin less mineralized surface layer (28 μm) and highly mineralized interstitial bone (mimicking mineralization resulting from anti-resorptive treatment) resulted in the highest stiffness. This could explain large reductions in fracture risk resulting from small increases in BMD. The non-uniform mineral distribution could also explain why bone tissue stiffness determined using nano-indentation is usually higher than finite element (FE)-determined stiffness. We conclude that the non-uniform mineral distribution in trabeculae does affect the mechanical properties of cancellous bone and that the tissue stiffness determined using FE-modeling could be improved by including detailed information about mineral distribution in trabeculae in the models.     

Effect of trabecular curvature on the stiffness of trabecular bone

Simplified structural models of trabecular bone have been used to model various forms of trabecular variability. The structural effects of variability of direction, length and thickness of the trabeculae have been studied using 'lattice-type' finite element models. However, many of the trabeculae are not perfectly straight, and have a small degree of curvature. The objective of this study is to quantify the influence of small curvatures of the trabeculae on the effective modulus of trabecular bone, in the principal material direction. An analytical analysis of the effect of curvature on a single trabecula is performed, utilizing the concept of cellular-solid models. Closed-form expressions are derived for the effect of curvature on the flexibility in the principal material direction. For comparison, expressions are derived for the flexibility of a straight oblique element, representing angular variability. A quantitative comparison is presented, which is dependent on the thickness of the trabeculae. It was found that small curvatures have a large effect on the stiffness of the trabecular structure. This effect is largest for thin trabeculae, and decreases for thick trabeculae. The stiffness of the trabecular structure can be reduced by a factor of up to four for thin trabeculae and up to two for thick trabeculae, even for small curvatures. The flexibility of curved elements is found to be larger than the flexibility of oblique elements with similar eccentricities. Thus it seems that curvature might play a role in determining the effective modulus of trabecular bone.     

Analysis of microstructural and mechanical alterations of trabecular bone in a simulated three-dimensional remodeling process

Bone remodeling is a complex dynamic process, which modulates both bone mass and bone microstructure. In addition to bone mass, bone microstructure is an important contributor to bone quality in osteoporosis and fragility fractures. However, the quantitative knowledge of evolution of three-dimensional (3D) trabecular microstructure in adaptation to the external forces is currently limited. In this study, a new 3D simulation method of remodeling of human trabecular bone was developed to quantitatively study the dynamic evolution of bone mass and trabecular microstructure in response to different external loading conditions. The morphological features of trabecular plate and rod, such as thickness and number density in different orientations were monitored during the remodeling process using a novel imaging analysis technique, namely Individual Trabecula Segmentation (ITS). We showed that the volume fraction and microstructures of trabecular bone including, trabecular type and orientation, were determined by the applied mechanical load. Particularly, the morphological parameters of trabecular plates were more sensitive to the applied load, indicating that they played the major role in the mechanical properties of the trabecular bone. Reducing the applied load caused severe microstructural deteriorations of trabecular bone, such as trabecular plate perforation, rod breakage, and a conversion from plates to rods.     

Finite Element Analysis of Simulations of Cancellous Bone Resorption

A stochastic simulation of the resorption of cancellous bone has been developed and integrated with a finite element model to predict the resultant change in structural properties of bone as bone density decreases. The resorption represents the net imbalance of osteoclast and osteoblast activity that occurs in osteoporosis. A simple lattice structure of trabecular bone is considered, with an examination of the lattice geometry and discretization indicating that just five trabeculae need to be modelled. The results from the analysis show how the mechanical properties of the cancellous bone degrade with osteoporosis and demonstrate how the method can be used to predict the relationships between stiffness and density or porosity.     

Fast and accurate specimen-specific simulation of trabecular bone elastic modulus using novel beam-shell finite element models

Elastic modulus and strength of trabecular bone are negatively affected by osteoporosis and other metabolic bone diseases. Micro-computed tomography-based beam models have been presented as a fast and accurate way to determine bone competence. However, these models are not accurate for trabecular bone specimens with a high number of plate-like trabeculae. Therefore, the aim of this study was to improve this promising methodology by representing plate-like trabeculae in a way that better reflects their mechanical behavior. Using an optimized skeletonization and meshing algorithm, voxel-based models of trabecular bone samples were simplified into a complex structure of rods and plates. Rod-like and plate-like trabeculae were modeled as beam and shell elements, respectively, using local histomorphometric characteristics. To validate our model, apparent elastic modulus was determined from simulated uniaxial elastic compression of 257 cubic samples of trabecular bone (4mm×4mm×4mm; 30μm voxel size; BIOMED I project) in three orthogonal directions using the beam-shell models and using large-scale voxel models that served as the gold standard. Excellent agreement (R(2)=0.97) was found between the two, with an average CPU-time reduction factor of 49 for the beam-shell models. In contrast to earlier skeleton-based beam models, the novel beam-shell models predicted elastic modulus values equally well for structures from different skeletal sites. It allows performing detailed parametric analyses that cover the entire spectrum of trabecular bone microstructures.     

Mechanical properties of cancellous bone in the human mandibular condyle are anisotropic

The objective of the present study was (1) to test the hypothesis that the elastic and failure properties of the cancellous bone of the mandibular condyle depend on the loading direction, and (2) to relate these properties to bone density parameters. Uniaxial compression tests were performed on cylindrical specimens (n=47) obtained from the condyles of 24 embalmed cadavers. Two loading directions were examined, i.e., a direction coinciding with the predominant orientation of the plate-like trabeculae (axial loading) and a direction perpendicular to the plate-like trabeculae (transverse loading). Archimedes' principle was applied to determine bone density parameters. The cancellous bone was in axial loading 3.4 times stiffer and 2.8 times stronger upon failure than in transverse loading. High coefficients of correlation were found among the various mechanical properties and between them and the apparent density and volume fraction. The anisotropic mechanical properties can possibly be considered as a mechanical adaptation to the loading of the condyle in vivo.     

Altered tissue properties induce changes in cancellous bone architecture in aging and diseases

The mechanical properties of cancellous bone depend on its architecture and the tissue properties of the mineralized matrix. The architecture is continuously adapted to external loads. In this paper, it was assumed that changes in tissue properties leading to changes in tissue deformation can induce adaptation of the architecture. We asked whether changes in cancellous bone architecture with aging and in e.g. early osteoarthrosis can be explained from changes in tissue properties.This was investigated using computer models in which the cancellous architecture was adapted to external loads. Bone tissue with deformations below a certain threshold was resorbed, deformations above another threshold induced formation. Deformations between these two boundaries, in the 'lazy zone', did not induce bone adaptation. The effects of changes in bone tissue stiffness on bone mass, global stiffness and architecture were investigated. The bone gain (40-60%) resulting from a 50% decrease in tissue stiffness (simulating diseased tissue) was much larger than the bone loss (2-30%) resulting from a 50% increase in tissue stiffness (simulating highly mineralized, old tissue). The adaptation induced by a decrease in tissue stiffness resulted in an almost constant stiffness in the main load bearing direction, but the transversal stiffness decreased. An increased tissue stiffness resulted in a higher stiffness in the main direction and overcompensation in the transversal directions: the global stiffness could become even smaller than the stiffness of the original model. Concluding, we showed that changes in trabecular bone in aging and diseases can be partly explained from changes in tissue properties.     

Hip bone trabecular architecture shows uniquely distinctive locomotor behaviour in South African australopithecines

Cancellous bone retains structural and behavioural properties which are time and strain-rate dependent. As the orientation of the trabeculae (trajectories) follows the direction of the principal strains imposed by daily loadings, habitual postural and locomotor behaviours are responsible for a variety of trabecular architectures and site-specific textural arrangements of the pelvic cancellous network. With respect to the great ape condition, the human trabecular pattern is characterized by a distinctive ilioischial bundle, an undivided sacropubic bundle, and a full diagonal crossing (approximately 100 degrees) over the acetabulum between the ilioischial and the sacropubic bundles. Advanced digital image processing (DIP) of hip bone radiographs has revealed that adolescent and adult South African australopithecines retained an incompletely developed human-like trabecular pattern associated with gait-related features that are unique among the extant primates.     

Stress-concentrating effect of resorption lacunae in trabecular bone

Analyses of the distributions of stress and strain within individual bone trabeculae have not yet been reported. In this study, four trabeculae were imaged and finite elements models were generated in an attempt to quantify the variability of stress/strain in real trabeculae. In three of these trabeculae, cavities were identified with depths comparable to values reported for resorption lacunae ( approximately 50 microm)-although we cannot be certain, it is most probable that they are indeed resorption lacunae. A tensile load was applied to each trabeculum to simulate physiological loading and to ensure that bending was minimized. The force carried by each trabecula was calculated from this value using the average cross sectional area of each trabecula. The analyses predict that very high stresses (>100 MPa) existed within bone trabecular tissue. Stress and strain distributions were highly heterogeneous in all cases, more so in trabeculae with the presumptive resorption lacunae where at least 30% of the tissue had a strain greater than 4000 micoepsilon in all cases. Stresses were elevated at the pit of the lacunae, and peak stress concentrations were located in the longitudinal direction ahead of the lacunae. Given these high strains, we suggest that microdamage is inevitable around resorption lacunae in trabecular bone, and may cause the bone multicellular unit to proceed to resorb a packet of bone in the trabeculum rather than just resorb whatever localized area was initially targeted.     

Variability of tissue mineral density can determine physiological creep of human vertebral cancellous bone

Creep is a time-dependent viscoelastic deformation observed under a constant prolonged load. It has been indicated that progressive vertebral deformation due to creep may increase the risk of vertebral fracture in the long-term. The objective of this study was to examine the relationships of creep with trabecular architecture and tissue mineral density (TMD) parameters in human vertebral cancellous bone at a physiological static strain level. Architecture and TMD parameters of cancellous bone were analyzed using microcomputerized tomography (micro-CT) in specimens cored out of human vertebrae. Then, creep and residual strains of the specimens were measured after a two-hour physiological compressive constant static loading and unloading cycle. Creep developed (3877 ± 2158 με) resulting in substantial levels of non-recoverable post-creep residual strain (1797 ± 1391 με). A strong positive linear correlation was found between creep and residual strain (r = 0.94, p < 0.001). The current results showed that smaller thickness, larger surface area, greater connectivity of trabeculae, less mean tissue mineral density (TMD, represented by gray levels) and higher variability of TMD are associated with increasing logarithmic creep rate. The TMD variability (GL(COV)) was the strongest correlate of creep rate (r = 0.79, p < 0.001). This result suggests that TMD variability may be a useful parameter for estimating the long-term deformation of a whole vertebral body. The results further suggest that the changes in TMD variability resulting from bone remodeling are of importance and may provide an insight into the understanding of the mechanisms underlying progressive failure of vertebral bodies and development of a clinical fracture.     

Adaptive diffuse domain approach for calculating mechanically induced deformation of trabecular bone

Remodelling of trabecular bone is essentially affected by the mechanical load of the trabeculae. Mathematical modelling and simulation of the remodelling process have to include time-consuming calculations of the displacement field within the complex trabecular structure under loading. We present an adaptive diffuse domain approach for calculating the elastic bone deformation based on micro computer tomogram data of real trabecular bone structures and compared it with a conventional voxel-based finite element method. In addition to allowing for higher computational efficiency, the adaptive approach is characterised by a very smooth representation of the bone surface, which suggests that this approach would be suitable as a basis for future simulations of bone resorption and formation processes within the trabecular structure.