An insight into tumor dormancy equilibrium via the analysis of its domain of attraction

Abstract The trajectories of the dynamic system which regulates the competition between the populations of malignant cells and immune cells may tend to an asymptotically stable equilibrium in which the sizes of these populations do not vary, which is called tumor dormancy. Especially for lower steady-state sizes of the population of malignant cells, this equilibrium represents a desirable clinical condition since the tumor growth is blocked. In this context, it is of mandatory importance to analyze the robustness of this clinical favorable state of health in the face of perturbations. To this end, the paper presents an optimization technique to determine whether an assigned rectangular region, which surrounds an asymptotically stable equilibrium point of a quadratic systems, is included into the domain of attraction of the equilibrium itself. The biological relevance of the application of this technique to the analysis of tumor growth dynamics is shown on the basis of a recent quadratic model of the tumor–immune system competition dynamics. Indeed the application of the proposed methodology allows to ensure that a given safety region, determined on the basis of clinical considerations, belongs to the domain of attraction of the tumor blocked equilibrium; therefore for the set of perturbed initial conditions which belong to such region, the convergence to the healthy steady state is guaranteed. The proposed methodology can also provide an optimal strategy for cancer treatment.

[1]  Graziano Chesi,et al.  Computing optimal quadratic lyapunov functions for polynomial nonlinear systems via lmis , 2002 .

[2]  J. Leith,et al.  Interlocking triads of growth control in tumors. , 1995, Bulletin of mathematical biology.

[3]  A. Radunskaya,et al.  Mixed Immunotherapy and Chemotherapy of Tumors: Modeling, Applications and Biological Interpretations , 2022 .

[4]  J S Spratt,et al.  Decelerating growth and human breast cancer , 1993, Cancer.

[5]  Claudio De Persis,et al.  Proceedings of the 15th IFAC World Congress , 2002 .

[6]  Alberto Tesi,et al.  On the stability domain estimation via a quadratic Lyapunov function: convexity and optimality properties for polynomial systems , 1996, IEEE Trans. Autom. Control..

[7]  J. Leith,et al.  Autocrine and paracrine growth factors in tumor growth: a mathematical model. , 1991, Bulletin of mathematical biology.

[8]  R. K. Miller,et al.  Stability analysis of complex dynamical systems , 1982 .

[9]  J. Leith,et al.  Growth factors and growth control of heterogeneous cell populations. , 1993, Bulletin of mathematical biology.

[10]  A. Garulli,et al.  LMI‐based computation of optimal quadratic Lyapunov functions for odd polynomial systems , 2005 .

[11]  J. Panetta,et al.  A mathematical model of periodically pulsed chemotherapy: tumor recurrence and metastasis in a competitive environment. , 1996, Bulletin of mathematical biology.

[12]  A. Vicino,et al.  On the estimation of asymptotic stability regions: State of the art and new proposals , 1985 .

[13]  Ami Radunskaya,et al.  A mathematical tumor model with immune resistance and drug therapy: an optimal control approach , 2001 .

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

[15]  J. Blattman,et al.  Cancer Immunotherapy: A Treatment for the Masses , 2004, Science.

[16]  E. Davison,et al.  A computational method for determining quadratic lyapunov functions for non-linear systems , 1971 .

[17]  Graziano Chesi Estimating the domain of attraction: a light LMI technique for a class of polynomial systems , 2003, 42nd IEEE International Conference on Decision and Control (IEEE Cat. No.03CH37475).

[18]  A. Vicino,et al.  Some results on the asymptotic stability of second-order nonlinear systems , 1984 .

[19]  Francesco Amato,et al.  Robust Control of Linear Systems Subject to Uncertain Time-Varying Parameters , 2006 .

[20]  J. Leith,et al.  Positive feedback and angiogenesis in tumor growth control. , 1997, Bulletin of mathematical biology.

[21]  J. Leith,et al.  Dormancy, regression, and recurrence: towards a unifying theory of tumor growth control. , 1994, Journal of theoretical biology.

[22]  Luigi Preziosi,et al.  From population dynamics to modelling the competition between tumors and immune system , 1996 .

[23]  Davison,et al.  A computational method for determining quadratic Lyapunov Functions for nonlinear systems , 1970 .

[24]  A. Perelson,et al.  Nonlinear dynamics of immunogenic tumors: parameter estimation and global bifurcation analysis. , 1994, Bulletin of mathematical biology.

[25]  J. Thorp,et al.  Stability regions of nonlinear dynamical systems: a constructive methodology , 1989 .

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

[27]  Felix F. Wu,et al.  Stability regions of nonlinear autonomous dynamical systems , 1988 .

[28]  Z. Bajzer,et al.  Conceptual frameworks for mathematical modeling of tumor growth dynamics , 1996 .

[29]  L. D. Pillis,et al.  A Validated Mathematical Model of Cell-Mediated Immune Response to Tumor Growth , 2005 .

[30]  A theoretical explanation of "concomitant resistance". , 1995, Bulletin of mathematical biology.

[31]  E. Yaz Linear Matrix Inequalities In System And Control Theory , 1998, Proceedings of the IEEE.

[32]  G. Spiga,et al.  Analysis of a kinetic cellular model for tumor-immune system interaction , 1999 .

[33]  B. Tibken,et al.  Estimation of the domain of attraction for polynomial systems , 2005, The Fourth International Workshop on Multidimensional Systems, 2005. NDS 2005..

[34]  M. J. Chapman,et al.  A mathematical model for the in vitro kinetics of the anti-cancer agent topotecan. , 2004, Mathematical biosciences.

[35]  Robert A Gatenby,et al.  Application of quantitative models from population biology and evolutionary game theory to tumor therapeutic strategies. , 2003, Molecular cancer therapeutics.

[36]  Stephen P. Boyd,et al.  Linear Matrix Inequalities in Systems and Control Theory , 1994 .

[37]  Ram Rup Sarkar,et al.  Cancer self remission and tumor stability-- a stochastic approach. , 2005, Mathematical biosciences.