Lacunar-canalicular network in femoral cortical bone is reduced in aged women and is predominantly due to a loss of canalicular porosity

The lacunar-canalicular network (LCN) of bone contains osteocytes and their dendritic extensions, which allow for intercellular communication, and are believed to serve as the mechanosensors that coordinate the processes of bone modeling and remodeling. Imbalances in remodeling, for example, are linked to bone disease, including fragility associated with aging. We have reported that there is a reduction in scale for one component of the LCN, osteocyte lacunar volume, across the human lifespan in females. In the present study, we explore the hypothesis that canalicular porosity also declines with age. To visualize the LCN and to determine how its components are altered with aging, we examined samples from young (age: 20–23 y; n = 5) and aged (age: 70–86 y; n = 6) healthy women donors utilizing a fluorescent labelling technique in combination with confocal laser scanning microscopy. A large cross-sectional area of cortical bone spanning the endosteal to periosteal surfaces from the anterior proximal femoral shaft was examined in order to account for potential trans-cortical variation in the LCN. Overall, we found that LCN areal fraction was reduced by 40.6% in the samples from aged women. This reduction was due, in part, to a reduction in lacunar density (21.4% decline in lacunae number per given area of bone), but much more so due to a 44.6% decline in canalicular areal fraction. While the areal fraction of larger vascular canals was higher in endosteal vs. periosteal regions for both age groups, no regional differences were observed in the areal fractions of the LCN and its components for either age group. Our data indicate that the LCN is diminished in aged women, and is largely due to a decline in the canalicular areal fraction, and that, unlike vascular canal porosity, this diminished LCN is uniform across the cortex.

•

The lacunar-canalicular network (LCN) is reduced by 40.6% in aged women

•

The lacunar density (lacunae number per given area of bone) is reduced by 21.4% in aged women

•

The reduction in LCN in aged women is primarily due to a 44.6% loss of canaliculi

•

No endosteal vs. periosteal regional differences were observed in the LCN and its components in either young or aged women

•

A reduction in canaliculi with age may contribute to bone fragility in aged women

Estimation of anisotropic permeability in trabecular bone based on microCT imaging and pore-scale fluid dynamics simulations

In this paper, a comprehensive framework is proposed to estimate the anisotropic permeability matrix in trabecular bone specimens based on micro-computed tomography (microCT) imaging combined with pore-scale fluid dynamics simulations. Two essential steps in the proposed methodology are the selection of (i) a representative volume element (RVE) for calculation of trabecular bone permeability and (ii) a converged mesh for accurate calculation of pore fluid flow properties. Accurate estimates of trabecular bone porosities are obtained using a microCT image resolution of approximately 10 μm. We show that a trabecular bone RVE in the order of 2 × 2 × 2 mm³ is most suitable. Mesh convergence studies show that accurate fluid flow properties are obtained for a mesh size above 125,000 elements. Volume averaging of the pore-scale fluid flow properties allows calculation of the apparent permeability matrix of trabecular bone specimens. For the four specimens chosen, our numerical results show that the so obtained permeability coefficients are in excellent agreement with previously reported experimental data for both human and bovine trabecular bone samples. We also identified that bone samples taken from long bones generally exhibit a larger permeability in the longitudinal direction. The fact that all coefficients of the permeability matrix were different from zero indicates that bone samples are generally not harvested in the principal flow directions. The full permeability matrix was diagonalized by calculating the eigenvalues, while the eigenvectors showed how strongly the bone sample's orientations deviated from the principal flow directions. Porosity values of the four bone specimens range from 0.83 to 0.86, with a low standard deviation of ± 0.016, principal permeability values range from 0.22 to 1.45 ⋅ 10 ^-8 m², with a high standard deviation of ± 0.33. Also, the anisotropic ratio ranged from 0.27 to 0.83, with high standard deviation. These results indicate that while the four specimens are quite similar in terms of average porosity, large variability exists with respect to permeability and specimen anisotropy. The utilized computational approach compares well with semi-analytical models based on homogenization theory. This methodology can be applied in bone tissue engineering applications for generating accurate pore morphologies of bone replacement materials and to consistently select similar bone specimens in bone bioreactor studies.

The relationship between porosity and specific surface in human cortical bone is subject specific

A characteristic relationship for bone between bone volume fraction (BV/TV) and specific surface (BS/TV) has previously been proposed based on 2D histological measurements. This relationship has been suggested to be bone intrinsic, i.e., to not depend on bone type, bone site and health state. In these studies, only limited data comes from cortical bone. The aim of this paper was to investigate the relationship between BV/TV and BS/TV in human cortical bone using high-resolution micro-CT imaging and the correlations with subject-specific biometric data such as height, weight, age and sex. Images from femoral cortical bone samples of the Melbourne Femur Collection were obtained using synchrotron radiation micro-CT (SPring8, Japan). Sixteen bone samples from thirteen individuals were analysed in order to find bone volume fraction values ranging from 0.20 to 1. Finally, morphological models of the tissue microstructure were developed to help explain the relationship between BV/TV and BS/TV. Our experimental findings indicate that the BV/TV vs BS/TV relationship is subject specific rather than intrinsic. Sex and pore density were statistically correlated with the individual curves. However no correlation was found with body height, weight or age. Experimental cortical data points deviate from interpolating curves previously proposed in the literature. However, these curves are largely based on data points from trabecular bone samples. This finding challenges the universality of the curve: highly porous cortical bone is significantly different to trabecular bone of the same porosity. Finally, our morphological models suggest that changes in BV/TV within the same sample can be explained by an increase in pore area rather than in pore density. This is consistent with the proposed mechanisms of age-related endocortical bone loss. In addition, these morphological models highlight that the relationship between BV/TV and BS/TV is not linear at high BV/TV as suggested in the literature but is closer to a square root function.

Mathematical modeling of postmenopausal osteoporosis and its treatment by the anti-catabolic drug denosumab

Denosumab, a fully human monoclonal antibody, has been approved for the treatment of postmenopausal osteoporosis. The therapeutic effect of denosumab rests on its ability to inhibit osteoclast differentiation. Here, we present a computational approach on the basis of coupling a pharmacokinetics model of denosumab with a pharmacodynamics model for quantifying the effect of denosumab on bone remodeling. The pharmacodynamics model comprises an integrated systems biology-continuum micromechanics approach, including a bone cell population model, considering the governing biochemical factors of bone remodeling (including the action of denosumab), and a multiscale micromechanics-based bone mechanics model, for implementing the mechanobiology of bone remodeling in our model. Numerical studies of postmenopausal osteoporosis show that denosumab suppresses osteoclast differentiation, thus strongly curtailing bone resorption. Simulation results also suggest that denosumab may trigger a short-term bone volume gain, which is, however, followed by constant or decreasing bone volume. This evolution is accompanied by a dramatic decrease of the bone turnover rate by more than one order of magnitude. The latter proposes dominant occurrence of secondary mineralization (which is not anymore impeded through cellular activity), leading to higher mineral concentration per bone volume. This explains the overall higher bone mineral density observed in denosumab-related clinical studies.

Spatio-temporal structure of cell distribution in cortical bone multicellular units: a mathematical model

Bone remodelling maintains the functionality of skeletal tissue by locally coordinating bone-resorbing cells (osteoclasts) and bone-forming cells (osteoblasts) in the form of Bone Multicellular Units (BMUs). Understanding the emergence of such structured units out of the complex network of biochemical interactions between bone cells is essential to extend our fundamental knowledge of normal bone physiology and its disorders. To this end, we propose a spatio-temporal continuum model that integrates some of the most important interaction pathways currently known to exist between cells of the osteoblastic and osteoclastic lineage. This mathematical model allows us to test the significance and completeness of these pathways based on their ability to reproduce the spatio-temporal dynamics of individual BMUs. We show that under suitable conditions, the experimentally observed structured cell distribution of cortical BMUs is retrieved. The proposed model admits travelling-wave-like solutions for the cell densities with tightly organised profiles, corresponding to the progression of a single remodelling BMU. The shapes of these spatial profiles within the travelling structure can be linked to the intrinsic parameters of the model such as differentiation and apoptosis rates for bone cells. In addition to the cell distribution, the spatial distribution of regulatory factors can also be calculated. This provides new insights on how different regulatory factors exert their action on bone cells leading to cellular spatial and temporal segregation, and functional coordination.