Quantification of Lacunar–Canalicular Interstitial Fluid Flow Through Computational Modeling of Fluorescence Recovery After Photobleaching

Skeletal adaptation to mechanical loading has been widely hypothesized to involve the stimulation of osteocytes by interstitial fluid flow (IFF). However, direct investigation of this hypothesis has been difficult due in large part to the inability to directly measure IFF velocities within the lacunar–canalicular system. Measurements of fluorescence recovery after photobleaching (FRAP) within individual lacunae could be used to quantify lacunar–canalicular IFF when combined with mathematical modeling. In this study, we used a computational transport model to characterize the relationship between flow frequency (0.5–10 Hz), peak flow velocity (0–300 μm/s), tracer diffusion coefficient (100–300 μm2/s), and transport enhancement (i.e., (k/k0) − 1, where k and k0 are the transport rates in the presence/absence of flow) during lacunar FRAP investigations. We show that this relationship is well described by a simple power law with frequency-dependent coefficients, and is relatively insensitive to variations in lacunar geometry. Using this power law relationship, we estimated peak IFF velocities in hindlimb mice subjected to intramedullary pressurization using values of k and k0 previously obtained from ex vivo lacunar FRAP investigations. Together, our findings suggest that skeletal adaptation in hindlimb suspended mice subjected to dynamic intramedullary pressure occurred in the presence of IFF at levels associated with physiological loading.