BACKGROUND: Microstructural simulations of bone remodeling are particularly relevant in the clinical management of osteoporosis. Before a model can be applied in the clinics, a validation against controlled in vivo data is crucial. Here we present a strain-adaptive feedback algorithm for the simulation of trabecular bone remodeling in response to loading and pharmaceutical treatment and report on the results of the large-scale validation against in vivo data. METHODS: The algorithm follows the mechanostat principle and incorporates mechanical feedback, based on the local strain-energy density. For the validation, simulations of bone remodeling and adaptation in 180 osteopenic mice were performed. Permutations of the conditions for early (20th week) and late (26th week) loading of 8N or 0N, and treatments with bisphosphonates, or parathyroid hormone were simulated. Static and dynamic morphometry and local remodeling sites from in vivo and in silico studies were compared. FINDINGS: For each study an individual set of model parameters was selected. Trabecular bone volume fraction was chosen as an indicator of the accuracy of the simulations. Overall errors for this parameter were 0.1-4.5%. Other morphometric indices were simulated with errors of less than 19%. Dynamic morphometry was more difficult to predict, which resulted in significant differences from the experimental data. INTERPRETATION: We validated a new algorithm for the simulation of bone remodeling in trabecular bone. The results indicate that the simulations accurately reflect the effects of treatment and loading seen in respective experimental data, and, following adaptation to human data, could be transferred into clinics.
BACKGROUND: Microstructural simulations of bone remodeling are particularly relevant in the clinical management of osteoporosis. Before a model can be applied in the clinics, a validation against controlled in vivo data is crucial. Here we present a strain-adaptive feedback algorithm for the simulation of trabecular bone remodeling in response to loading and pharmaceutical treatment and report on the results of the large-scale validation against in vivo data. METHODS: The algorithm follows the mechanostat principle and incorporates mechanical feedback, based on the local strain-energy density. For the validation, simulations of bone remodeling and adaptation in 180 osteopenicmice were performed. Permutations of the conditions for early (20th week) and late (26th week) loading of 8N or 0N, and treatments with bisphosphonates, or parathyroid hormone were simulated. Static and dynamic morphometry and local remodeling sites from in vivo and in silico studies were compared. FINDINGS: For each study an individual set of model parameters was selected. Trabecular bone volume fraction was chosen as an indicator of the accuracy of the simulations. Overall errors for this parameter were 0.1-4.5%. Other morphometric indices were simulated with errors of less than 19%. Dynamic morphometry was more difficult to predict, which resulted in significant differences from the experimental data. INTERPRETATION: We validated a new algorithm for the simulation of bone remodeling in trabecular bone. The results indicate that the simulations accurately reflect the effects of treatment and loading seen in respective experimental data, and, following adaptation to human data, could be transferred into clinics.