Mathematical modelling of bone adaptation of the metacarpal subchondral bone in racehorses

In Thoroughbred racehorses, fractures of the distal limb are commonly catastrophic. Most of these fractures occur due to the accumulation of fatigue damage from repetitive loading, as evidenced by microdamage at the predilection sites for fracture. Adaptation of the bone in response to training loads is important for fatigue resistance. In order to better understand the mechanism of subchondral bone adaptation to its loading environment, we utilised a square root function defining the relationship between bone volume fraction (fBM) and specific surface (Sv) of the subchondral bone of the lateral condyles of the third metacarpal bone (MCIII) of the racehorse, and using this equation, developed a mathematical model of subchondral bone that adapts to loading conditions observed in vivo. The model is expressed as an ordinary differential equation incorporating a formation rate that is dependent on strain energy density. The loading conditions applied to a selected subchondral region, i.e. volume of interest, were estimated based on joint contact forces sustained by racehorses in training. For each of the initial conditions of fBM we found no difference between subsequent homoeostatic fBM at any given loading condition, but the time to reach equilibrium differed by initial fBM and loading condition. We found that the observed values for fBM from the mathematical model output were a good approximation to the existing data for racehorses in training or at rest. This model provides the basis for understanding the effect of changes to training strategies that may reduce the risk of racehorse injury. © 2018, Springer-Verlag GmbH Germany, part of Springer Nature