Mechanistic mathematical modelling of mercaptopurine effects on cell cycle of human acute lymphoblastic leukaemia cells

The antimetabolite mercaptopurine (MP) is widely used to treat childhood acute lymphoblastic leukaemia (ALL). To study the dynamics of MP on the cell cycle, we incubated human T-cell leukaemia cell lines (Molt-4 sensitive and resistant subline and P12 resistant) with 10 μM MP and measured total cell count, cell cycle distribution, percent viable, percent apoptotic, and percent dead cells serially over 72 h. We developed a mathematical model of the cell cycle dynamics after treatment with MP and used it to show that the Molt-4 sensitive controls had a significantly higher rate of cells entering apoptosis (2.7-fold, P<0.00001) relative to the resistant cell lines. Additionally, when treated with MP, the sensitive cell line showed a significant increase in the rate at which cells enter apoptosis compared to its controls (2.4-fold, P<0.00001). Of note, the resistant cell lines had a higher rate of antimetabolite incorporation into the DNA of viable cells (>1.4-fold, P<0.01). Lastly, in contrast to the other cell lines, the Molt-4 resistant subline continued to cycle, though at a rate slower relative to its control, rather than proceed to apoptosis. This led to a larger S-phase block in the Molt-4 resistant cell line, but not a higher rate of cell death. Gene expression of apoptosis, cell cycle, and repair genes were consistent with mechanistic dynamics described by the model. In summary, the mathematical model provides a quantitative assessment to compare the cell cycle effects of MP in cells with varying degrees of MP resistance.