A stochastic model of DNA fragments rejoining

When cells are exposed to ionizing radiation, DNA damages in the form of single strand breaks (SSBs), double strand breaks (DSBs), base damage or their combinations are frequent events. It is known that the complexity and severity of DNA damage depends on the quality of radiation, and the microscopic dose deposited in small segments of DNA, which is often related to the linear transfer energy (LET) of the radiation. Experimental studies have suggested that under the same dose, high LET radiation induces more small DNA fragments than low-LET radiation, which affects Ku efficiently binding with DNA end and might be a main reason for high-LET radiation induced RBE [1] since DNA DSB is a major cause for radiation-induced cell death. In this work, we proposed a mathematical model of DNA fragments rejoining according to non-homologous end joining (NHEJ) mechanism. By conducting Gillespie's stochastic simulation, we found several factors that impact the efficiency of DNA fragments rejoining. Our results demonstrated that aberrant DNA damage repair can result predominantly from the occurrence of a spatial distribution of DSBs leading to short DNA fragments. Because of the low efficiency that short DNA fragments recruit repair protein and release the protein residue after fragments rejoining, Ku-dependent NHEJ is significantly interfered with short fragments. Overall, our work suggests that inhibiting the Ku-dependent NHEJ may significantly contribute to the increased efficiency for cell death and mutation observed for high LET radiation.