A mechanistic model of high dose irradiation damage

The main goal of our study is to develop a realistic mechanistic model of the effect of ionizing radiation on DNA in mammalian cells. We consider a population of cells structured by the number of DNA double strand breaks due to radiation. Using the system of linear differential equation, the model describes the evolution of the irradiated population of cells in time. The work is in three parts. First, we consider the effect of a single dose of radiation, while in the second part we work on the model parameter estimation using Nelder–Mead simplex algorithm which allows us to relate the clinically useful parameters of the LQ relation to aspects of cellular activity that can be manipulated experimentally. In the third part, we deal with cell killing effects of fractioned doses of radiation. Using MATLAB, we observed the cell survival fractions can be well approximated by the Linear–Quadratic relation and also show fewer cell will die if the dose is fractionated in two or more fractions