Bone remodelling is carried out by `bone multicellular units (BMUs) in which active osteoclasts and active osteoblasts are spatially and temporally coupled. The refilling of new bone by osteoblasts towards the back of the BMU occurs at a rate that de
pends both on the number of osteoblasts and on their secretory activity. In cortical bone, a linear phenomenological relationship between matrix apposition rate (MAR) and BMU cavity radius is found experimentally. How this relationship emerges from the combination of complex, nonlinear regulations of osteoblast number and secretory activity is unknown. Here, we extend our previous mathematical model of cell development within a single BMU to investigate how osteoblast number and osteoblast secretory activity vary along the BMUs closing cone. MARs predicted by the model are compared with data from tetracycline double labelling experiments. We find that the linear phenomenological relationship observed in these experiments between MAR and BMU cavity radius holds for most of the refilling phase simulated by our model, but not near the start and end of refilling. This suggests that at a particular bone site undergoing remodelling, bone formation starts and ends rapidly. Our model also suggests that part of the observed cross-sectional variability in tetracycline data may be due to different bone sites being refilled by BMUs at different stages of their lifetime. The different stages of a BMUs lifetime depend on whether the cell populations within the BMU are still developing or have reached a quasi-steady state while travelling through bone. We find that due to their longer lifespan, active osteoblasts reach a quasi-steady distribution more slowly than active osteoclasts. We suggest that this fact may locally enlarge the Haversian canal diameter (due to a local lack of osteoblasts compared to osteoclasts) near the BMUs point of origin.
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 struct
ured 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.
