![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
Abstract
In this work, we consider a chaotic system that plays a vital role in the treatment of cancer by injection of a virus externally. Due to the sensitivity of this disease, most of its treatments are highly risky. Therefore, we have designed control inputs using adaptive and passive control techniques for virotherapy. Both controllers are designed to bring global stability to the cancer system with the aid of a quadratic Lyapunov function. Furthermore, we use simulations to verify our controllers. Moreover, we show that our adaptive control technique gives better results in comparison.
1. Introduction
Cancer remains a high priority of research due to its prevalence, treatment difficulties, and high mortality rate. In late 19 century, the removal of tumors by surgical operations and radiation was common but was not effective as hoped [Citation24]. The failure of these two treatments was due to the uncontrollable growth and recurring nature of cancer cells. During the last two decades, the study of cancer has accelerated not only biologists, but among mathematician as well. In Saleem et al. [Citation22] several mathematical models are used to analyse the interaction of tumor cells with healthy cells, the immune system and with viruses. These mathematical models facilitated investigations into the growth rate of cancer and its prognosis and treatment. Using a bifurcation analysis of a modified tomorrow growth model, Wei [Citation26] observed chemotherapy-induced damage may cause a microscopic tumor. In 2018, Maddali1 et al. [Citation12] studied chaotic dynamics caused by the mutual interaction of cancer cells with immune cells and healthy host cells. Eftimie et al. [Citation5] investigated oncolytic therapy with the addition of viruses in the system. A mathematical model for the growth dynamics of infected and uninfected tumors is considered by Malinzi et al. [Citation15].
Sensitivity to initial conditions is a significant characteristic of chaos [Citation11, Citation17]. Chaos plays a remarkable role in the field of biology [Citation4], engineering [Citation21], finance [Citation8], quantum mechanics [Citation20], weather forecasting [Citation18] and physical sciences [Citation13]. Regardless of any positive aspects of chaos, its occurrence in cancer patients may cause difficulties as the unpredictable growth rate of tumor cells in the body [Citation23] and predicting cardiac arrhythmia [Citation2]. These pros and cons of chaos lead to the introduction of new subject; control of chaos [Citation17]. Since the pioneering work of Ott, Grebogi and Yorke in century, the control and synchronization of chaos have become an active field for scientific research [Citation6]. A variety of productive control and synchronization [Citation7] approaches have been implemented, such as sliding [Citation3], adaptive [Citation9, Citation25], optimal [Citation16], active [Citation19], passive [Citation1], feedback linearization [Citation10] control techniques. It is important to mention that the main aim of all the control techniques appearing in literature is to overcome chaos and gain global stability. In the adaptive control technique, controllers are designed by using parameter estimation and low complexity algorithms. Virotherapy [Citation24] is the most recent and effective treatment regimen in which external viruses are injected for the removal of a tumor and to stimulate the immune system. The main hurdle in virus therapy is that all viruses have pre-existing antibodies and the challenging task is to choose the right virus. However, in 2017, Malinzi showed that cancer can be treated with a mixture of chemotherapy and virotherapy [Citation14]. However, recent research and clinical trials have shown that virotherapy has fewer side effects on patients and is less toxic towards healthy cells than chemotherapy. The target of virotherapy treatments is to lower the number of tumor cells by increasing the concentration of infected cells to the highest level. Virotherapy is extremely nonlinear because of its dependence on cancer cells disintegration and virus transmission rates. A model introduced by Soufiani [Citation24] in 2017 is described by the three differential equations
(1)
(1) In this model, tumor cells are classified into uninfected cells
, infected cells
and the virus population
, whereas r is the effective growth rate for uninfected cells, σ is the death rate of infected cells, K is the carrying capacity, p shows inter-cellular fusion, k is constant rate of virus transmission, α is the virus production rate by replication of infected cells and ω gives constant rate of virus eradication are the involved parameters in our considered model. System (Equation1
(1)
(1) ) is chaotic for parameter values; K = 2139,
, k = 0.05, r = 0.6, p = 0.02,
,
and
.
In aforementioned paragraphs, we have discussed that virotheraphy is one of the treatments for cancer. But, along with its importance it has side effects with injecting dangerous virus and consequently life of a patient can be in danger. Therefore, keeping in mind its drawback, we have introduced two control inputs for Virotherphy chaotic model to bring stability and save lives.
This paper is organized as follows. In Sections 2 and 3, we design control inputs by using adaptive and passive control techniques. We compare the results obtained from adaptive and passive control techniques in Section 4, and concluding remarks appear in Section 5.
2. Adaptive control
In this section, we use adaptive techniques to attain the global stability of the system at fixed points. The main components of this technique are the creation of an error dynamical system and the construction of a Lyapunov functions. We consider a modification of system (Equation1(1)
(1) )
(2)
(2) where
; i = 1, 2, 3 are adaptive control input rates for bringing stability to the virotherapy process in such a way that infected cell, uninfected cell and virus populations converge to the desired stable points with the passage of time. These control inputs can be designed with the help of an error dynamical system and an updating law (Equation9
(9)
(9) ). Denote an equilibria of the virotherapy system by
and introduce the error terms
(3)
(3) Now, differentiating Equation (Equation3
(3)
(3) ), we get a dynamical system for the error terms
(4)
(4) In this technique parameters involved in system are unknown and time-dependent. Therefore, arbitrary parameters
,
,
,
,
are introduced for σ, k, p, α and ω, respectively. Control terms with respect to error terms (Equation3
(3)
(3) ) are
(5)
(5) where
for i = 1, 2, 3. We get the following error dynamical system by putting Equation (Equation5
(5)
(5) ) into Equation (Equation4
(4)
(4) )
(6)
(6) where
,
,
,
and
are the error terms associated with the estimated parameters as defined by
,
,
,
,
. Taking derivative of each error term, we obtain
In order to obtain an updating law, we use a quadratic Lyapunov function
(7)
(7) Clearly,
is greater than or equal to zero and to be negative definite, we require its derivative
(8)
(8) to satisfy
. Therefore, in view of Equation (Equation8
(8)
(8) ), we can attain this by the following updating law
(9)
(9) where
,
,
,
and
are positive gain constants. There are several techniques for calculating these gain constants. However, we require only that
> 0,
for simplicity. Now, substituting Equation (Equation9
(9)
(9) ) into Equation (Equation8
(8)
(8) ) to get
(10)
(10) which clearly is a negative definite function. Hence, we get that solution of system (Equation1
(1)
(1) ) is asymptotically stable with the updating law (Equation9
(9)
(9) ).
In Figure , it can be observed that trajectories of error terms converge to zero asymptotically. In other words, when error terms approach zero, state variables start moving towards the desired point. In Figure , estimated parameters converge to the original values. To plot updated law (Equation9(9)
(9) ), initial conditions are considered as
,
,
,
,
. As these estimated parameter values approach towards their original values, chaotic trajectories start converging toward their desired point.
3. Passive feedback
In this section, chaotic virotheraphy system is controlled about equilibria using passive control function for Equation (Equation1
(1)
(1) )
(11)
(11) Here
is the control input for our considered chaotic system (Equation1
(1)
(1) ), which is designed on the basis of Feedback Linearization. In similar way as adaptive technique, the error dynamics for system (Equation11
(11)
(11) ) is defined as
(12)
(12) where error state variable
is considered as output,
=
and
=
. Error dynamics system in the form of
,
and output Y can be rewritten as
(13)
(13) According to passive control normal form, Equation (Equation13
(13)
(13) ) can be written as
(14)
(14) and
(15)
(15) Storage function for passive control method is considered as
(16)
(16) where
is Lyapunov function of
coupled with
and
(17)
(17) To ensure existence of stability for solution of system (Equation11
(11)
(11) ), we differentiate Equation (Equation16
(16)
(16) ) and get
(18)
(18) Therefore,
for all real values and hence we achieve global stability of
. Hence, zero dynamics of controlled system (Equation13
(13)
(13) ) is asymptotically stable and is a minimum phase system. According to passive control theory controller is designed to make system passive as
(19)
(19) Here β is positive real number and μ is an external signal. Using
,
into Equation (Equation19
(19)
(19) ), control function can be achieved as
For convenience, writing the system in form of error terms;
,
and
, yield
(20)
(20) Equation (Equation20
(20)
(20) ) gives the control function to control chaos that appears in virotherapy.
4. Numerical simulation results
In this section, our aim is to demonstrate graphical results of adaptive and passive control techniques obtained in sections 2 and 3. For this purpose, parameter values are considered as K = 2139, , k = 0.05, r = 0.6, p = 0.02,
,
and
along with initial conditions
,
and
.
For these initial conditions and parameter values, trajectories of system (Equation1(1)
(1) ) show unpredictable behaviour. For adaptive control, positive gain constants are selected as
for
.
Figure represents global stability of uninfected tumor cells u, infected tumor cells w and virus cells v about fixed points ,
and
when adaptive controllers are introduced into original system (Equation1
(1)
(1) ).
Figure illustrates the convergence of trajectories towards fixed points when passive feedback controller is activated at t = 0. On the other hand, constant used in control function is designed through passive technique, considered as and external signal
. Figures and verify the analytical results obtained in sections 2 and 3, respectively. Now, to analyze the control method for transient speed, control effort and convergence time, we plot each state variable using both techniques.
From Figure , it is observed that control obtained for chaotic system through adaptive technique is more efficient than passive technique. Further, it is observed that controller designed in virotheraphy becomes effective in less time using adaptive control technique, while passive control function takes more time in comparison. In Figure , control of chaos and convergence to equilibrium point is achieved faster for adaptive method instead of passive. It is also observed that, trajectories in both cases start with oscillation of large amplitude but control inputs at t = 0.1 and t = 0.9, using adaptive and passive control techniques, respectively control the systems and converge to their point. In Figures –, it it is noticed that trajectories of system (Equation1
(1)
(1) ) are stabilizing towards equilibria
is target for virotheraphy treatment, because at this point virus replication rate increases, while population of uninfected tumor cells decreases and converges to its lowest state. Third equilibrium point,
(21)
(21) represent the partial success of virotherapy. However, if virotherapy fails, then its cure state can be achieved through global stability about equilibria
. This point state can be reached when
(22)
(22) The term;
indicates infection rate or the release of virions, whereas its population in the cure state is a consequence of assuming a continuous constant dosage. Moreover, the term
represents virus elimination by death of infected cells. Hence,
reaches when the infection rate is more than virus elimination. In the present work, we have proved in Sections 2 and 3 that if system (Equation1
(1)
(1) ) is controlled about
then population of uninfected tumor cells decreases. Thus, from clinical point of view, our findings indicate that regeneration of infected cells can improve treatment.
5. Conclusion
In this paper, global stability for virotheraphy system is investigated using adaptive and passive control techniques. In both cases, we used the concept of Lyapunov stability. Passive control method coupled with feedback is used to attain convergence of virotheraphy model. Analytical work and graphical simulations confirm that both methods successfully achieve the control of virotherapy system. But, adaptive controller works faster in comparison.
Figure 5. Comparison of adaptive and passive control inputs for each state variables using equilibria; .
![Figure 5. Comparison of adaptive and passive control inputs for each state variables using equilibria; E1=(0,0,0).](/cms/asset/69463332-201c-4b94-9df6-c51d0305c15c/tjbd_a_2104391_f0005_oc.jpg)
Figure 6. Comparison of adaptive and passive control inputs for each state variables using equilibria; .
![Figure 6. Comparison of adaptive and passive control inputs for each state variables using equilibria; E2=(2139,0,0).](/cms/asset/a33acc43-f1c5-4f5a-a6c5-b946aed1616f/tjbd_a_2104391_f0006_oc.jpg)
Figure 7. Comparison of adaptive and passive control inputs for each state variables using equilibria; .
![Figure 7. Comparison of adaptive and passive control inputs for each state variables using equilibria; E3=(0.4,2.22,11.11).](/cms/asset/0464f1a9-57b3-4ace-9b7c-e480bf36631f/tjbd_a_2104391_f0007_oc.jpg)
Virotherapy shows complex behaviour due to chaotic growth of tumor cells and virus decomposition rate. In this article, we used chaos control regimens to suppress chaoticity and converge the state variables toward the target point at which tumor cells are minimized. Global stability towards steady state and
shows the complete removal of tumor and partial success of therapy, whereas convergence towards
refers to incurable stage at which population of tumor cells converges to maximum level. Moreover, stabilizing trajectories of system (Equation1
(1)
(1) ) towards equilibria
is the target state for virotheraphy treatment. At
virus replication rate increases, while population of uninfected tumor cells decreases and converges to its lowest state.
Compliance with ethical standards
Ethical approval: This paper does not contain any studies with human participants or animals performed by any of the author.
Acknowledgements
All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.
Disclosure statement
No potential conflict of interest was reported by the author(s).
Data availability statement
MATLAB is used for graphs and tedious calculations in this work; codes are available upon reasonable request from corresponding author.
References
- M. Babalola, B. Iyorzor, and I. Odesanya, Chaos control of a modified 4-d memristor chaotic oscillator via passive control technique, Niger. J. Technol. 37 (2018), pp. 480.
- S. Dai and D. Schaeffer, Chaos for cardiac arrhythmias through a one-dimensional modulation equation for alternans, Chaos 20 (2010), pp. 1–8.
- W. Ding, F. Liu, H. Chen, N. Wang, and D. Liang, Sliding mode control of fractional-order delayed memristive chaotic system with uncertainty and disturbance, Commun. Theor. Phys. 68(6) (2017), pp. 741.
- L. Dzyubak, O. Dzyubak, and J. Awrejcewicz, Controlling and stabilizing unpredictable behavior of metabolic reactions and carcinogenesis in biological systems, Nonlinear Dyn. 97(3) (2019), pp. 1853–1866.
- R. Eftimie, C.K. Macnamara, J. Dushoff, J.L. Bramson, and D.J.D Earn, Bifurcations and chaotic dynamics in a tumour-immune-virus system, Math. Model. Nat. Phenom. 11(5) (2016), pp. 65–85.
- M. Feki, Sliding mode based control and synchronization of chaotic systems in presence of parametric uncertainties, Appl. Slid. Mode Control Sci. Eng. 709 (2017), pp. 35–59.
- J.M.V. Grzybowski, E.E.N. Macau, and T. Yoneyama, Isochronal synchronization of time delay and delay-coupled chaotic systems, J. Phys. A: Math. Theor. 44(17) (2011), pp. 1–16.
- D. Guegan, Chaos in economics and finance, Annu. Rev. Control 33(1) (2009), pp. 89–93.
- R. Guo, U.E. Vincent, and B.A. Idowu, Synchronization of chaos in RCL-shunted josephson junction using a simple adaptive controller, Phys. Scr. 79(3) (2009), pp. 1–6.
- S. Hadhrami, A. Saaban, A. Ibrahim, M. Shahzad, and I. Ahmad, Linear active control algorithm to synchronize a nonlinear hiv/aids dynamical system, Asian J. Appl. Sci. Eng. 3 (2014), pp. 96–113.
- R. Klages, Chaos: A very short introduction, J. Phys. A: Math. Theor. 40(29) (2007), pp. 8604–8605.
- R.K. Maddali, D. Ahluwalia, A. Chaudhuri, and S. Sarif Hassan, Dynamics of a three dimensional chaotic cancer model, Int. J. Math. Trends Technol. 53(5) (2018), pp. 353–368.
- G.M. Mahmoud, M.A. Al-Kashif, and A.A. Farghaly, Chaotic and hyperchaotic attractors of a complex nonlinear system, J. Phys. A: Math. Theor. 41(5) (2008), pp. 1–12.
- J. Malinzi, A. Eladdadi, and P. Sibanda, Modelling the spatiotemporal dynamics of chemovirotherapy cancer treatment, J. Biol. Dyn. 11(1) (2017), pp. 244–274.
- J. Malinzi, R. Ouifki, A. Eladdadi, D.F.M Torres, and K.A. White, Enhancement of chemotherapy using oncolytic virotherapy: mathematical and optimal control analysis, Math. Biosci. Eng. 15(2) (2018), pp. 1435.
- M. Marwan and S. Ahmad, Bifurcation analysis for energy transport system and its optimal control using parameter self-tuning law, Soft. Comput. 24 (2020), pp. 17221–17231.
- M. Marwan, S. Ahmad, M. Aqeel, and M. Sabir, Control analysis of rucklidge chaotic system, J. Dyn. Syst. Meas. Control 141(4) (2019), pp. 1–7.
- S. Moon, J.M. Seo, B.-S. Han, J. Park, and J.-J. Baik, A physically extended lorenz system, Chaos: Interdiscip. J. Nonlinear Sci. 29(6) (2019), pp. 1–12.
- L. Moysis, E. Petavratzis, M. Marwan, C. Volos, H. Nistazakis, and S. Ahmad, Analysis, synchronization, and robotic application of a modified hyperjerk chaotic system, Complexity 2020 (2020), pp. 1–15.
- T. Prosen, Chaos and complexity of quantum motion, J. Phys. A: Math. Theor. 40(28) (2007), pp. 7881–7918.
- M. Sabir, M. Marwan, S. Ahmad, M. Fiaz, and F. Khan, Observer and descriptor satisfying incremental quadratic constraint for class of chaotic systems and its applications in a quadrotor chaotic system, Chaos, Solit. Fractals 137 (2020), pp. 1–13.
- M. Saleem, T. Agrawal, and A. Anees, A study of tumour growth based on stoichiometric principles: a continuous model and its discrete analogue, J. Biol. Dyn. 8(1) (2014), pp. 117–134.
- M. Shahzad, Chaos control in three dimensional cancer model by state space exact linearization based on lie algebra, Mathematics 4 (2016), pp. 33.
- B.N. Soufiani and M.U. Salamci, Chaotic behavior in virotherapy for cancer treatment, Int. Conf. Inf. Commun. Autom. Technol. (2017), pp. 1–6.
- S. Vaidyanathan, Adaptive control of a chemical chaotic reactor, Int. J. PharmTech Res. 8(3) (2015), pp. 377–382.
- H.C. Wei, A mathematical model of tumour growth with Beddington-Deangelis functional response: a case of cancer without disease, J. Biol. Dyn. 12(1) (2018), pp. 194–210.