The formation of new bone involves both the deposition of bone matrix, and the formation of a network of cells embedded within the bone matrix, called osteocytes. Osteocytes derive from bone-synthesising cells (osteoblasts) that become buried in bone
matrix during bone deposition. The generation of osteocytes is a complex process that remains incompletely understood. Whilst osteoblast burial determines the density of osteocytes, the expanding network of osteocytes regulates in turn osteoblast activity and osteoblast burial. In this paper, a spatiotemporal continuous model is proposed to investigate the osteoblast-to-osteocyte transition. The aims of the model are (i) to link dynamic properties of osteocyte generation with properties of the osteocyte network imprinted in bone, and (ii) to investigate Marottis hypothesis that osteocytes prompt the burial of osteoblasts when they become covered with sufficient bone matrix. Osteocyte density is assumed in the model to be generated at the moving bone surface by a combination of osteoblast density, matrix secretory rate, rate of entrapment, and curvature of the bone substrate, but is found to be determined solely by the ratio of the instantaneous burial rate and matrix secretory rate. Osteocyte density does not explicitly depend on osteoblast density nor curvature. Osteocyte apoptosis is also included to distinguish between the density of osteocyte lacuna and the density of live osteocytes. Experimental measurements of osteocyte lacuna densities are used to estimate the rate of burial of osteoblasts in bone matrix. These results suggest that: (i) burial rate decreases during osteonal infilling, and (ii) the control of osteoblast burial by osteocytes is likely to emanate as a collective signal from a large group of osteocytes, rather than from the osteocytes closest to the bone deposition front.
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.
Age-related bone loss and postmenopausal osteoporosis are disorders of bone remodelling, in which less bone is reformed than resorbed. Yet, this dysregulation of bone remodelling does not occur equally in all bone regions. Loss of bone is more pronou
nced near and at the endocortex, leading to cortical wall thinning and medullary cavity expansion, a process sometimes referred to as trabecularisation or cancellisation. Cortical wall thinning is of primary concern in osteoporosis due to the strong deterioration of bone mechanical properties that it is associated with. In this paper, we examine the possibility that the non-uniformity of microscopic bone surface availability could explain the non-uniformity of bone loss in osteoporosis. We use a computational model of bone remodelling in which microscopic bone surface availability influences bone turnover rate and simulate the evolution of the bone volume fraction profile across the midshaft of a long bone. We find that bone loss is accelerated near the endocortical wall where the specific surface is highest. Over time, this leads to a substantial reduction of cortical wall thickness from the endosteum. The associated expansion of the medullary cavity can be made to match experimentally observed cross-sectional data from the Melbourne Femur Collection. Finally, we calculate the redistribution of the mechanical stresses in this evolving bone structure and show that mechanical load becomes critically transferred to the periosteal cortical bone.
Bone is a biomaterial undergoing continuous renewal. The renewal process is known as bone remodelling and is operated by bone-resorbing cells (osteoclasts) and bone-forming cells (osteoblasts). Both biochemical and biomechanical regulatory mechanisms
have been identified in the interaction between osteoclasts and osteoblasts. Here we focus on an additional and poorly understood potential regulatory mechanism of bone cells, that involves the morphology of the microstructure of bone. Bone cells can only remove and replace bone at a bone surface. However, the microscopic availability of bone surface depends in turn on the ever-changing bone microstructure. The importance of this geometrical dependence is unknown and difficult to quantify experimentally. Therefore, we develop a sophisticated mathematical model of bone cell interactions that takes into account biochemical, biomechanical and geometrical regulations. We then investigate numerically the influence of bone surface availability in bone remodelling within a representative bone tissue sample. The interdependence between the bone cells activity, which modifies the bone microstructure, and changes in the microscopic bone surface availability, which in turn influences bone cell development and activity, is implemented using a remarkable experimental relationship between bone specific surface and bone porosity. Our model suggests that geometrical regulation of the activation of new remodelling events could have a significant effect on bone porosity and bone stiffness. On the other hand, geometrical regulation of late stages of osteoblast and osteoclast differentiation seems less significant. We conclude that the development of osteoporosis is probably accelerated by this geometrical regulation in cortical bone, but probably slowed down in trabecular bone.
To maintain bone mass during bone remodelling, coupling is required between bone resorption and bone formation. This coordination is achieved by a network of autocrine and paracrine signalling molecules between cells of the osteoclast lineage and cel
ls of the osteoblastic lineage. Mathematical modelling of signalling between cells of both lineages can assist in the interpretation of experimental data, clarify signalling interactions and help develop a deeper understanding of complex bone diseases. In this paper, we further develop a mathematical model of bone cell interactions by Pivonka et al. (2008) to include the proliferation of precursor osteoblasts into the model. This inclusion is important to be able to account for Wnt signalling, believed to play an important role in anabolic responses of bone. We show that an increased rate of differentiation to precursor cells or an increased rate of proliferation of precursor osteoblasts themselves both result in increased bone mass. However, modelling these different processes separately enables the new model to represent recent experimental discoveries such as the role of Wnt signalling in bone biology and the recruitment of osteoblast progenitor cells by transforming growth factor beta. Finally, we illustrate the power of the new models capabilities by applying the model to prostate cancer metastasis to bone. In the bone microenvironment, prostate cancer cells are believed to release some of the same signalling molecules used to coordinate bone remodelling (i.e. Wnt and PTHrP), enabling the cancer cells to disrupt normal signalling and coordination between bone cells. This disruption can lead to either bone gain or bone loss. We demonstrate that the new computational model developed here is capable of capturing some key observations made on the evolution of the bone mass due to metastasis of prostate cancer to the bone microenvironment
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.
