The development of predictive mathematical models can contribute to a deeper understanding of the specific stages of bone mechanobiology and the process by which bone adapts to mechanical forces. The objective of this work was to predict, with spatial accuracy, cortical bone adaptation to mechanical load, in order to better understand the mechanical cues that might be driving adaptation. The axial tibial loading model was used to trigger cortical bone adaptation in C57BL/6 mice and provide relevant biological and biomechanical information. A method for mapping cortical thickness in the mouse tibia diaphysis was developed, allowing for a thorough spatial description of where bone adaptation occurs. Poroelastic finite-element (FE) models were used to determine the structural response of the tibia upon axial loading and interstitial fluid velocity as the mechanical stimulus. FE models were coupled with mechanobiological governing equations, which accounted for non-static loads and assumed that bone responds instantly to local mechanical cues in an on–off manner. The presented formulation was able to simulate the areas of adaptation and accurately reproduce the distributions of cortical thickening observed in the experimental data with a statistically significant positive correlation (Kendall'sτrank coefficientτ= 0.51,p< 0.001). This work demonstrates that computational models can spatially predict cortical bone mechanoadaptation to a time variant stimulus. Such models could be used in the design of more efficient loading protocols and drug therapies that target the relevant physiological mechanisms.
Bone is a dynamic tissue and adapts its architecture in response to biological and mechanical factors. Here we investigate how cortical bone formation is spatially controlled by the local mechanical environment in the murine tibia axial loading model (C57BL/6). We obtained 3D locations of new bone formation by performing ‘slice and view’ 3D fluorochrome mapping of the entire bone and compared these sites with the regions of high fluid velocity or strain energy density estimated using a finite element model, validated with ex-vivo bone surface strain map acquired ex-vivo using digital image correlation. For the comparison, 2D maps of the average bone formation and peak mechanical stimulus on the tibial endosteal and periosteal surface across the entire cortical surface were created. Results showed that bone formed on the periosteal and endosteal surface in regions of high fluid flow. Peak strain energy density predicted only the formation of bone periosteally. Understanding how the mechanical stimuli spatially relates with regions of cortical bone formation in response to loading will eventually guide loading regime therapies to maintain or restore bone mass in specific sites in skeletal pathologies.
Bones' functionally adaptive responses to mechanical loading can usefully be studied in the tibia by the application of loads between the knee and ankle in normal and genetically modified mice. Such loading also deforms the fibula. Our present study was designed to ascertain whether the fibula adapts to loading in a similar way to the tibia and could thus provide an additional bone in which to study functional adaptation. The right tibiae/fibulae in C57BL/6 mice were subjected to a single period of axial loading (40 cycles at 10 Hz with 10-second intervals between each cycle; approximately 7 min/day, 3 alternate days/week, 2 weeks). The left tibiae/fibulae were used as non-loaded, internal controls. Both left and right fibulae and tibiae were analyzed by micro-computed tomography at the levels of the mid-shaft of the fibula and 25% from its proximal and distal ends. We also investigated the effects of intermittent parathyroid hormone (iPTH) on the (re)modelling response to 2-weeks of loading and the effect of 2-consecutive days of loading on osteocytes' sclerostin expression. These in vivo experiments confirmed that the fibula showed similar loading-related (re)modelling responses to those previously documented in the tibia and similar synergistic increases in osteogenesis between loading and iPTH. The numbers of sclerostin-positive osteocytes at the proximal and middle fibulae were markedly decreased by loading. Collectively, these data suggest that the mouse fibula, as well as the tibia and ulna, is a useful bone in which to assess bone cells' early responses to mechanical loading and the adaptive (re)modelling that this engenders.