Load-adaptive bone remodeling simulations reveal osteoporotic microstructural and mechanical changes in whole human vertebrae.

Osteoporosis is a major medical burden and its impact is expected to increase in our aging society. It is associated with low bone density and microstructural deterioration. Treatments are available, but the critical factor is to define individuals at risk from osteoporotic fractures. Computational simulations investigating not only changes in net bone tissue volume, but also changes in its microstructure where osteoporotic deterioration occur might help to better predict the risk of fractures. In this study, bone remodeling simulations with a mechanical feedback loop were used to predict microstructural changes due to osteoporosis and their impact on bone fragility from 50 to 80 years of age. Starting from homeostatic bone remodeling of a group of seven, mixed sex whole vertebrae, five mechanostat models mimicking different biological alterations associated with osteoporosis were developed, leading to imbalanced bone formation and resorption with a total net loss of bone tissue. A model with reduced bone formation rate and cell sensitivity led to the best match of morphometric indices compared to literature data and was chosen to predict postmenopausal osteoporotic bone loss in the whole group. Thirty years of osteoporotic bone loss were predicted with changes in morphometric indices in agreement with experimental measurements, and only showing major deviations in trabecular number and trabecular separation. In particular, although being optimized to match to the morphometric indices alone, the predicted bone loss revealed realistic changes on the organ level and on biomechanical competence. While the osteoporotic bone was able to maintain the mechanical stability to a great extent, higher fragility towards error loads was found for the osteoporotic bones.