Holistic View on Cell Survival and DNA Damage: How Model-Based Data Analysis Supports Exploration of Dynamics in Biological Systems

In this work, a method is established to calibrate a model that describes the basic dynamics of DNA damage and repair. The model can be used to extend planning for radiotherapy and hyperthermia in order to include the biological effects. In contrast to “syntactic” models (e.g., describing molecular kinetics), the model used here describes radiobiological semantics, resulting in a more powerful model but also in a far more challenging calibration. Model calibration is attempted from clonogenic assay data (doses of 0–6 Gy) and from time-resolved comet assay data obtained within 6 h after irradiation with 6 Gy. It is demonstrated that either of those two sources of information alone is insufficient for successful model calibration, and that both sources of information combined in a holistic approach are necessary to find viable model parameters. Approximate Bayesian computation (ABC) with simulated annealing is used for parameter search, revealing two aspects that are beneficial to resolving the calibration problem: (1) assessing posterior parameter distributions instead of point-estimates and (2) combining calibration runs from different assays by joining posterior distributions instead of running a single calibration run with a combined, computationally very expensive objective function.