A caution for oncologists: chemotherapy can cause chaotic dynamics

BACKGROUND AND OBJECTIVE The effect of chemotherapy in cancer models is mostly handled by using a separate equation for chemotherapeutic agent. In this study, we do not consider a separate equation for drug but rather introduce its effect in terms of a parameter m representing the fraction of tumor cells killed by chemotherapeutic drug module. The main objective of this study is to provide conditions on model parameters which when fulfilled the grave consequences of cancer can be avoided. This study also shows that chemotherapy at times can produce unexpected results. METHODS Linearization method to study the stability of model equilibria. RESULTS The results obtained in this study are governed by the trichotomy law on the number 1-a12-d1, where a12 represents the negative effect on the growth of cancer cells due to their competition with host cells for resources and d1 is rate of annihilation of cancer cells due to chemotherapy. It is seen that in case of under-dose drug module when d1<1-a12, the complete eradication of cancer is not possible. When d1=1-a12, the model suggests occurrence of chaotic dynamics. When the drug dose is properly adjusted so that d1>1-a12, the complete eradication of cancer is guaranteed. CONCLUSION The results of the model of this paper given for the post vascular stages of tumor suggest criteria to select a particular drug module (a single drug or a combination of drugs) that the chemotherapy procedure should adapt to eradicate cancer. This study injects a note of caution for oncologists that chemotherapy as cancer treatment can also cause chaotic dynamics in certain situations. This study also presents a plausible explanation to the question why sometimes a tumor grows in the body and then gets cured without any medical intervention.

[1]  A. Cerone,et al.  Mathematical modeling of drug resistance due to KRAS mutation in colorectal cancer. , 2016, Journal of theoretical biology.

[2]  Lin Shen,et al.  Mathematical Modelling and Analysis of the Tumor Treatment Regimens with Pulsed Immunotherapy and Chemotherapy , 2016, Comput. Math. Methods Medicine.

[3]  Frank Nani,et al.  A chemotherapy model for the treatment of cancer with metastasis , 2002 .

[4]  H. I. Freedman,et al.  A mathematical model of cancer treatment by immunotherapy. , 2000, Mathematical biosciences.

[5]  R. Laubenbacher,et al.  Addressing current challenges in cancer immunotherapy with mathematical and computational modeling , 2017, bioRxiv.

[6]  J. Sprott,et al.  Chaos in low-dimensional Lotka–Volterra models of competition , 2006 .

[7]  Stephen P. Banks,et al.  Chaos in a Three-Dimensional Cancer Model , 2010, Int. J. Bifurc. Chaos.

[8]  U. an der Heiden,et al.  A basic mathematical model of the immune response. , 1995, Chaos.

[9]  Xiaofeng Liao,et al.  On the global boundedness of the Lü system , 2016, Appl. Math. Comput..

[10]  Sanyi Tang,et al.  Chaos in functional response host–parasitoid ecosystem models , 2002 .

[11]  Xiaofeng Liao,et al.  Dynamical behaviors of a generalized Lorenz family , 2017 .

[12]  T. N. Mokaev,et al.  Homoclinic orbits, and self-excited and hidden attractors in a Lorenz-like system describing convective fluid motion Homoclinic orbits, and self-excited and hidden attractors , 2015 .

[13]  John A. Adam,et al.  A mathematical model of cycle-specific chemotherapy , 1995 .

[14]  John Carl Panetta,et al.  A mathematical model of periodically pulsed chemotherapy: Tumor recurrence and metastasis in a competitive environment , 1996 .

[15]  H. I. Freedman,et al.  A mathematical model for the effect of anti-angiogenic therapy in the treatment of cancer tumours by chemotherapy , 2013 .

[16]  A. Radunskayay,et al.  A Mathematical Tumor Model with Immune Resistance and Drug Therapy : An Optimal Control Approach , 1999 .

[17]  Nikolay V. Kuznetsov,et al.  Hidden attractors in Dynamical Systems. From Hidden oscillations in Hilbert-Kolmogorov, Aizerman, and Kalman Problems to Hidden Chaotic Attractor in Chua Circuits , 2013, Int. J. Bifurc. Chaos.

[18]  R. Gatenby,et al.  Models of tumor-host interaction as competing populations: implications for tumor biology and treatment. , 1995, Journal of theoretical biology.

[19]  A. Jemal,et al.  Global cancer transitions according to the Human Development Index (2008-2030): a population-based study. , 2012, The Lancet. Oncology.

[20]  Mustafa Mamat,et al.  Mathematical Model of Cancer Treatments Using Immunotherapy, Chemotherapy and Biochemotherapy , 2013 .

[21]  Sunita Gakkhar,et al.  Existence of chaos in two-prey, one-predator system , 2003 .

[22]  Giuseppe Jordão,et al.  Mathematical models in cancer therapy , 2017, Biosyst..

[23]  T. N. Mokaev,et al.  Finite-time Lyapunov dimension and hidden attractor of the Rabinovich system , 2015, 1504.04723.

[24]  Manar A. Alqudah Cancer treatment by stem cells and chemotherapy as a mathematical model with numerical simulations , 2020 .

[25]  Xiaofeng Liao,et al.  Dynamics of a New 5D Hyperchaotic System of Lorenz Type , 2018, Int. J. Bifurc. Chaos.

[26]  Kazem Nouri,et al.  Chaos in a stochastic cancer model , 2020 .

[27]  Nikolay V. Kuznetsov,et al.  Hidden attractors in electromechanical systems with and without equilibria , 2016 .

[28]  X. Liao,et al.  Some New Results for the Generalized Lorenz System , 2017 .

[29]  Gennady A. Leonov,et al.  Lyapunov's direct method in the estimation of the Hausdorff dimension of attractors , 1992 .

[30]  A. Dalgleish,et al.  The relevance of non-linear mathematics (chaos theory) to the treatment of cancer, the role of the immune response and the potential for vaccines. , 1999, QJM : monthly journal of the Association of Physicians.

[31]  G. P. Samanta,et al.  Dynamical Behaviour of a Tumor-Immune System with Chemotherapy and Optimal Control , 2013 .

[32]  H. I. Freedman,et al.  A mathematical model of vascular tumor treatment by chemotherapy , 2005, Math. Comput. Model..

[33]  E. Ahmed Fractals and chaos in cancer models , 1993 .

[34]  Abdur Raheem,et al.  A mixed therapy minimal model: Some strategies for eradication or minimization of cancer , 2020, Comput. Methods Programs Biomed..

[35]  A. Eladdadi,et al.  Modeling cancer-immune responses to therapy , 2014, Journal of Pharmacokinetics and Pharmacodynamics.

[36]  G. D. Knott,et al.  Modeling tumor regrowth and immunotherapy , 2001 .

[37]  Ricardo Gómez Aiza,et al.  Analysis of a mathematical model of periodically pulsed chemotherapy treatment , 2017 .

[38]  D. Kirschner,et al.  Modeling immunotherapy of the tumor – immune interaction , 1998, Journal of mathematical biology.

[39]  Nikolay V. Kuznetsov,et al.  Algorithms for searching for hidden oscillations in the Aizerman and Kalman problems , 2011 .