Mathematical model of radiation effects on thrombopoiesis in rhesus macaques and humans

A mathematical model that describes the effects of acute radiation exposure on thrombopoiesis in primates and humans is presented. Thrombopoiesis is a complex multistage dynamic process with potential differences between species. Due to known differences in cellular radiosensitivities, nadir times, and cytopenia durations, direct extrapolation from rhesus to human platelet dynamics is unrealistic. Developing mathematical models of thrombopoiesis for both humans and primates allows for the comparison of the system's response across species. Thus, data obtained in primate experiments can be extrapolated to predictions in humans. Parameter values for rhesus macaques and humans were obtained either from direct experimental measurements or through optimization procedures using dynamic data on platelet counts following radiation exposure. Model simulations accurately predict trends observed in platelet dynamics: at low radiation doses platelet counts decline after a time lag, and nadir depth is dose dependent. The models were validated using data that was not used during the parameterization process. In particular, additional experimental data was used for rhesus, and accident and platelet donor data was used for humans. The model aims to simulate the average response in rhesus and humans following irradiation. Variation in platelet dynamics due to individual variability can be modeled using Monte Carlo simulations in which parameter values are sampled from distributions. This model provides insight into the time course of the physiological effects of radiation exposure, information which could be valuable for disaster planning and survivability analysis and help in drug development of radiation medical countermeasures.