Often, resistance to drugs is an obstacle to a successful treatment of cancer. In spite of the importance of the problem, the actual mechanisms that control the evolution of drug resistance are not fully understood. Many attempts to study drug resistance have been made in the mathematical modeling literature. Clearly, in order to understand drug resistance, it is imperative to have a good model of the underlying dynamics of cancer cells. One of the main ingredients that has been recently introduced into the rapidly growing pool of mathematical cancer models is stem cells. Surprisingly, this all-so-important subset of cells has not been fully integrated into existing mathematical models of drug resistance. In this work we incorporate the various possible ways in which a stem cell may divide into the study of drug resistance. We derive a previously undescribed estimate of the probability of developing drug resistance by the time a tumor is detected and calculate the expected number of resistant cancer stem cells at the time of tumor detection. To demonstrate the significance of this approach, we combine our previously undescribed mathematical estimates with clinical data that are taken from a recent six-year follow-up of patients receiving imatinib for the first-line treatment of chronic myelogenous leukemia. Based on our analysis we conclude that leukemia stem cells must tend to renew symmetrically as opposed to their healthy counterparts that predominantly divide asymmetrically.

Recent mathematical models have been developed to study the dynamics of chronic myelogenous leukemia (CML) under imatinib treatment. None of these models incorporates the anti-leukemia immune response. Recent experimental data show that imatinib treatment may promote the development of anti-leukemia immune responses as patients enter remission. Using these experimental data we develop a mathematical model to gain insights into the dynamics and potential impact of the resulting anti-leukemia immune response on CML. We model the immune response using a system of delay differential equations, where the delay term accounts for the duration of cell division. The mathematical model suggests that anti-leukemia T cell responses may play a critical role in maintaining CML patients in remission under imatinib therapy. Furthermore, it proposes a novel concept of an “optimal load zone” for leukemic cells in which the anti-leukemia immune response is most effective. Imatinib therapy may drive leukemic cell populations to enter and fall below this optimal load zone too rapidly to sustain the anti-leukemia T cell response. As a potential therapeutic strategy, the model shows that vaccination approaches in combination with imatinib therapy may optimally sustain the anti-leukemia T cell response to potentially eradicate residual leukemic cells for a durable cure of CML. The approach presented in this paper accounts for the role of the anti-leukemia specific immune response in the dynamics of CML. By combining experimental data and mathematical models, we demonstrate that persistence of anti-leukemia T cells even at low levels seems to prevent the leukemia from relapsing (for at least 50 months). As a consequence, we hypothesize that anti-leukemia T cell responses may help maintain remission under imatinib therapy. The mathematical model together with the new experimental data imply that there may be a feasible, low-risk, clinical approach to enhancing the effects of imatinib treatment. Recent mathematical models have been developed to study the dynamics of chronic myelogenous leukemia (CML) under imatinib treatment. None of these models incorporates the anti-leukemia immune response. Recent experimental data show that imatinib treatment may promote the development of anti-leukemia immune responses as patients enter remission. Using these experimental data, we developed a mathematical model to gain insights into the dynamics and potential impact of the resulting anti-leukemia immune response on CML. The mathematical model suggests that anti-leukemia T cell responses may play a critical role in maintaining CML patients in remission under imatinib therapy. Furthermore, it proposes a novel concept of an “optimal load zone” for leukemic cells in which the anti-leukemia immune response is most effective. Imatinib therapy may drive leukemic cell populations to enter and fall below this optimal load zone too rapidly to sustain the anti-leukemia T cell response. As a potential therapeutic strategy, the model shows that vaccination approaches in combination with imatinib therapy may optimally sustain the anti-leukemia T cell response to potentially eradicate residual leukemic cells for a durable cure of CML.

Resistance to drugs has been an ongoing obstacle to a successful treatment of many diseases. In this work we consider the problem of drug resistance in cancer, focusing on random genetic point mutations. Most previous works on mathematical models of such drug resistance have been based on stochastic methods. In contrast, our approach is based on an elementary, compartmental system of ordinary differential equations. We use our very simple approach to derive results on drug resistance that are comparable to those that were previously obtained using much more complex mathematical techniques. The simplicity of our model allows us to obtain analytic results for resistance to any number of drugs. In particular, we show that the amount of resistance generated before the start of the treatment, and present at some given time afterward, always depends on the turnover rate, no matter how many drugs are simultaneously used in the treatment.

Imatinib and other tyrosine kinase inhibitors (TKI) have improved treatment of chronic myelogenous leukemia (CML); however, most patients are not cured. Deeper mechanistic understanding may improve TKI combination therapies to better control the residual leukemic cell population. In analyzing our patients' data, we found that many patients who otherwise responded well to imatinib therapy still showed variations in their BCR–ABL transcripts. To investigate this phenomenon, we applied a mathematical model that integrates CML and an autologous immune response to the patients' data. We define an immune window or a range of leukemic loads for which the autologous immune system induces an improved response. Our modeling results suggest that, at diagnosis, a patient's leukemic load is able to partially or fully suppress the autologous immune response developed in a majority of patients, toward the CML clone(s). Imatinib therapy drives the leukemic population into the “immune window,” allowing the patient's autologous immune cells to expand and eventually mount an efficient recognition of the residual leukemic burden. This response drives the leukemic load below this immune window, allowing the leukemic population to partially recover until another weaker immune response is initiated. Thus, the autologous immune response may explain the oscillations in BCR–ABL transcripts regularly observed in patients on imatinib. Cancer Res; 75(19); 4053–62. ©2015 AACR.