A mathematical model of cortical bone remodeling at cellular level under mechanical stimulus

A bone cell population dynamics model for cortical bone remodeling under mechanical stimulus is developed in this paper. The external experiments extracted from the literature which have not been used in the creation of the model are used to test the validity of the model. Not only can the model compare reasonably well with these experimentalresults such as the increase percentage of finalvalues of bone mineral content (BMC) and bone fracture energy(BFE) among different loading schemes (which provesthe validity of the model), but also predict the realtime developmentpattern of BMC and BFE, as well as the dynamics ofosteoblasts (OBA), osteoclasts (OCA), nitric oxide (NO) and prostaglandin E2 (PGE2) for each loading scheme, which can hardly be monitored through experiment. In conclusion, themodel is the first of its kind that is able to provide an insight into the quantitative mechanism of bone remodeling at cellular level by which bone cells are activated by mechanical stimulus in order to start resorption/formation of bone mass. More importantly, this model has laid a solid foundation based on which future work such as systemic control theory analysis of bone remodeling under mechanical stimulus  can be investigated. The to-be identified control mechanismwill help to develop effective drugs and combined nonpharmacological therapies to combat bone loss pathologies. Also this deeper understanding of how mechanical forces quantitatively interact with skeletal tissue is essential for the generation of bone tissue for tissue replacement purposes in tissue engineering.