Computational modeling based on confocal imaging predicts changes in osteocyte and dendrite shear stress due to canalicular loss with aging.

Exercise and physical activity exert mechanical loading on the bones which induces bone formation. However, the relationship between the osteocyte lacunar-canalicular morphology and mechanical stress experienced locally by osteocytes transducing signals for bone formation is not fully understood. In this study, we used computational modeling to predict the effect of canalicular density, the number of fluid inlets, and load direction on fluid flow shear stress (FFSS) and bone strains and how these might change following the microstructural deterioration of the lacunar-canalicular network that occurs with aging. Four distinct computational models were initially generated of osteocytes with either ten or eighteen dendrites using a fluid-structure interaction method with idealized geometries. Next, a young and a simulated aged osteocyte were developed from confocal images after FITC staining of the femur of a 4-month-old C57BL/6 mouse to estimate FFSS using a computational fluid dynamics approach. The models predicted higher fluid velocities in the canaliculi versus the lacunae. Comparison of idealized models with five versus one fluid inlet indicated that with four more inlets, one-half of the dendrites experienced FFSS greater than 0.8 Pa, which has been associated with osteogenic responses. Confocal image-based models of real osteocytes indicated a six times higher ratio of canalicular to lacunar surface area in the young osteocyte model than the simulated aged model and the average FFSS in the young model (FFSS = 0.46 Pa) was three times greater than the aged model (FFSS = 0.15 Pa). Interestingly, the surface area with FFSS values above 0.8 Pa was 23 times greater in the young versus the simulated aged model. These findings may explain the impaired mechano-responsiveness of osteocytes with aging.