KNOWLEDGEBASE - ARTICLE #2039

Fitting the linear-quadratic model of cell death by radiation using GraphPad Prism

What is the linear-quadratic model?

Mammalian cell death after exposure to radiation after a dose D of radiation often follows the linear-quadratic model (1-3).

Fraction of cells surviving =  e-(A*D + B*D^2)

The linear-quadratic model was derived by Chadwick and Leenhouts (2). They proposed that the linear component (A*X) represents cell death  due to a single lethal hit to the DNA, and that the quadratic (B*X2) component represents cell death that only happens with two hits.  However, it is now clear that this mechanism is not correct, and the biological interpretation of the LQ parameters, A and B, is unclear (3). Even though we don't know the biological basis for the model, it does a reasonable job of describing radiation induced cell death, except perhaps at very high radiation doses. Bodgi et. al. review the linear-quadratic equation and some alternatives (3) and propose a biological model that corresponds to the linear-quadratic model (4). 

Fitting the linear-quadratic model in GraphPad Prism

The equation below can be entered into Prism, where Y is the fraction survival and X is the dose:

Y = exp(-1*(A*X + B*X^2))

The graph of this model looks like this:

The graph on the left uses a linear Y axis. The graph on the right uses a logarithmic Y axis. This was done in GraphPad Prism by making the change in the Format Axis dialog. The nonlinear regression "saw" the original data where the Y values were fractions (not logarithm of fractions; that would require a modified model). 

Nonlinear regression requires that the user set the initial values of each parameter. For the fits above, I made the initial value of A equal -1 times the initial slope, and the initial value of B equal -0.1 times the initial slope. These rules for initial values seem to work fine with a variety of data. 

This attached file creates the graphs above, plus similar ones where Y is percentage of cells that survived, or the number of cells surviving. It also includes a set of analyses for when Y is fraction, percentage or number of dead cells. 

References

  1. Brenner, D. J. (2008). The linear-quadratic model is an appropriate methodology for determining isoeffective doses at large doses per fraction, 18/4: 234–9. Elsevier.
  2. Chadwick, K.H., Leenhouts, H.P., 1973. A molecular theory of cell survival. Phys. Med. Biol. 13, 78–87.
  3. Bodgi, L., Canet, A., Pujo-Menjouet, L., Lesne, A., Victor, J.-M., & Foray, N. (2016). Mathematical models of radiation action on living cells: From the target theory to the modern approaches. A historical and critical review. Journal of Theoretical Biology, 394: 93–101. DOI: 10.1016/j.jtbi.2016.01.018
  4. Bodgi, L., & Foray, N. (2016). The nucleo-shuttling of the ATM protein as a basis for a novel theory of radiation response: resolution of the linear-quadratic model. International Journal of Radiation Biology, 92/3: 117–31. DOI: 10.3109/09553002.2016.1135260

 

Explore the Knowledgebase

Analyze, graph and present your scientific work easily with GraphPad Prism. No coding required.