Computational simulation of spontaneous bone straightening in growing children

Mechanoregulated trabecular bone adaptation: Progress report on in silico approaches

Adaptation is the process by which bone responds to changes in loading environment and modulates its properties and spatial organization to meet the mechanical demands. Adaptation in trabecular bone is achieved through increase in bone mass and alignment of trabecular-bone morphology along the loading direction. This transformation of internal microstructure is governed by mechanical stimuli sensed by mechanosensory cells in the bone matrix. Realisation of adaptation in the form of local bone-resorption and -formation activities as a function of mechanical stimuli is still debated. In silico modelling is a useful tool for simulation of various scenarios that cannot be investigated in vivo and particularly well suited for prediction of trabecular bone adaptation. This progress report presents the recent advances in in silico modelling of mechanoregulated adaptation at the scale of trabecular bone tissue. Four well-established bone-adaptation models are reviewed in terms of their recent improvements and validation. They consider various mechanical factors: (i) strain energy density, (ii) strain and damage, (iii) stress nonuniformity and (iv) daily stress. Contradictions of these models are discussed and their ability to describe adequately a real-life mechanoregulation process in bone is compared.

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.

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.

Biomechanical Study of the Development of Long Bones: Finite Element Structure Synthesis of the Human Second Proximal Phalanx Under Growth Conditions

Torsional loads are a possible mechanical explanation for the architecture of long bone. Finite element structure synthesis (FESS) has previously successfully been used as a deductive technique using Wolff's Law by applying expected loads to an unspecific homogeneous solid and eliminating stress free parts to verify muscle forces. The extended approach presented in this article includes further mechanobiological rules to model the development from a cartilage model to a finger bone. In contrast to former computational models, simulation of processes leading to both external growth and internal differentiation are included. Combined axial and torsional loads synthesize a complete human secondary proximal phalanx model comparable to form and internal structure to that observed in vivo. While the computational model is very sensitive to initial alterations of loads, changes after growth have a minor effect as observed in animal models. Predictions of cartilage growth and ossification during FESS showed significant similarities to ontogeny indicating the importance of mechanical factors for the morphogenesis of bone during growth. Anat Rec, 302:1389-1398, 2019. © 2018 Wiley Periodicals, Inc.

The Proprioceptive System Regulates Morphologic Restoration of Fractured Bones

Successful fracture repair requires restoration of bone morphology and mechanical integrity. Recent evidence shows that fractured bones of neonatal mice undergo spontaneous realignment, dubbed "natural reduction." Here, we show that natural reduction is regulated by the proprioceptive system and improves with age. Comparison among mice of different ages revealed, surprisingly, that 3-month-old mice exhibited more rapid and effective natural reduction than newborns. Fractured bones of null mutants for transcription factor Runx3, lacking functional proprioceptors, failed to realign properly. Blocking Runx3 expression in the peripheral nervous system, but not in limb mesenchyme, recapitulated the null phenotype, as did inactivation of muscles flanking the fracture site. Egr3 knockout mice, which lack muscle spindles but not Golgi tendon organs, displayed a less severe phenotype, suggesting that both receptor types, as well as muscle contraction, are required for this regulatory mechanism. These findings uncover a physiological role for proprioception in non-autonomous regulation of skeletal integrity.

Tooth eruption results from bone remodelling driven by bite forces sensed by soft tissue dental follicles: a finite element analysis

Intermittent tongue, lip and cheek forces influence precise tooth position, so we here examine the possibility that tissue remodelling driven by functional bite-force-induced jaw-strain accounts for tooth eruption. Notably, although a separate true 'eruptive force' is widely assumed, there is little direct evidence for such a force. We constructed a three dimensional finite element model from axial computerized tomography of an 8 year old child mandible containing 12 erupted and 8 unerupted teeth. Tissues modelled included: cortical bone, cancellous bone, soft tissue dental follicle, periodontal ligament, enamel, dentine, pulp and articular cartilage. Strain and hydrostatic stress during incisive and unilateral molar bite force were modelled, with force applied via medial and lateral pterygoid, temporalis, masseter and digastric muscles. Strain was maximal in the soft tissue follicle as opposed to surrounding bone, consistent with follicle as an effective mechanosensor. Initial numerical analysis of dental follicle soft tissue overlying crowns and beneath the roots of unerupted teeth was of volume and hydrostatic stress. To numerically evaluate biological significance of differing hydrostatic stress levels normalized for variable finite element volume, 'biological response units' in Nmm were defined and calculated by multiplication of hydrostatic stress and volume for each finite element. Graphical representations revealed similar overall responses for individual teeth regardless if incisive or right molar bite force was studied. There was general compression in the soft tissues over crowns of most unerupted teeth, and general tension in the soft tissues beneath roots. Not conforming to this pattern were the unerupted second molars, which do not erupt at this developmental stage. Data support a new hypothesis for tooth eruption, in which the follicular soft tissues detect bite-force-induced bone-strain, and direct bone remodelling at the inner surface of the surrounding bony crypt, with the effect of enabling tooth eruption into the mouth.

Ankle joint mechanics and foot proportions differ between human sprinters and non-sprinters

Recent studies of sprinters and distance runners have suggested that variations in human foot proportions and plantarflexor muscle moment arm correspond to the level of sprint performance or running economy. Less clear, however, is whether differences in muscle moment arm are mediated by altered tendon paths or by variation in the centre of ankle joint rotation. Previous measurements of these differences have relied upon assumed joint centres and measurements of bone geometry made externally, such that they would be affected by the thickness of the overlying soft tissue. Using magnetic resonance imaging, we found that trained sprinters have shorter plantarflexor moment arms (p = 0.011) and longer forefoot bones (p = 0.019) than non-sprinters. The shorter moment arms of sprinters are attributable to differences in the location of the centre of rotation (p < 0.001) rather than to differences in the path of the Achilles tendon. A simple computer model suggests that increasing the ratio of forefoot to rearfoot length permits more plantarflexor muscle work during plantarflexion that occurs at rates expected during the acceleration phase following the sprint start.

Computational modeling of bone density profiles in response to gait: a subject-specific approach

The goal of this study is to explore the potential of computational growth models to predict bone density profiles in the proximal tibia in response to gait-induced loading. From a modeling point of view, we design a finite element-based computational algorithm using the theory of open system thermodynamics. In this algorithm, the biological problem, the balance of mass, is solved locally on the integration point level, while the mechanical problem, the balance of linear momentum, is solved globally on the node point level. Specifically, the local bone mineral density is treated as an internal variable, which is allowed to change in response to mechanical loading. From an experimental point of view, we perform a subject-specific gait analysis to identify the relevant forces during walking using an inverse dynamics approach. These forces are directly applied as loads in the finite element simulation. To validate the model, we take a Dual-Energy X-ray Absorptiometry scan of the subject's right knee from which we create a geometric model of the proximal tibia. For qualitative validation, we compare the computationally predicted density profiles to the bone mineral density extracted from this scan. For quantitative validation, we adopt the region of interest method and determine the density values at fourteen discrete locations using standard and custom-designed image analysis tools. Qualitatively, our two- and three-dimensional density predictions are in excellent agreement with the experimental measurements. Quantitatively, errors are less than 3% for the two-dimensional analysis and less than 10% for the three-dimensional analysis. The proposed approach has the potential to ultimately improve the long-term success of possible treatment options for chronic diseases such as osteoarthritis on a patient-specific basis by accurately addressing the complex interactions between ambulatory loads and tissue changes.