Modeling the effects of inflammation in bone fracture healing

A new mathematical model is presented to study the early inflammatory effects in bone healing. It consists of a system of nonlinear ordinary differential equations that represents the interactions among macrophages, mesenchymal stem cells, and osteoblasts. A qualitative analysis of the model is performed to determine the equilibria and their corresponding stability properties. A set of numerical simulations is performed to support the theoretical results. The model is also used to numerically monitor the evolution of a broken bone for different types of fractures and to explore possible treatments to accelerate bone healing by administrating anti-inflammatory drugs.