Research Papers

Graphical Hopf Bifurcation of a Filippov HTLV-1 Model With Delay in Cytotoxic T Cells Response

[+] Author and Article Information
Elham Shamsara

Management and Social Determinants of Health
Research Center,
Mashhad University of Medical Sciences,
Mashhad 9133913716, Iran;
School of Mathematics,
Institute for Research in Fundamental
Sciences (IPM),
P.O. Box 19395-5746,
Tehran, Iran
e-mail: elham.shamsara@mail.um.ac.ir

Zahra Afsharnezhad

Department of Applied Mathematics,
Ferdowsi University of Mashhad,
Mashhad 9177948974, Iran
e-mail: afsharnezhad@math.um.ac.ir

Elham Javidmanesh

Department of Applied Mathematics,
Ferdowsi University of Mashhad,
Mashhad 9177948974, Iran
e-mail: javidmanesh@um.ac.ir

Contributed by the Dynamic Systems Division of ASME for publication in the JOURNAL OF DYNAMIC SYSTEMS, MEASUREMENT,AND CONTROL. Manuscript received March 13, 2017; final manuscript received February 22, 2018; published online April 9, 2018. Assoc. Editor: Sergey Nersesov.

J. Dyn. Sys., Meas., Control 140(9), 091007 (Apr 09, 2018) Paper No: DS-17-1150; doi: 10.1115/1.4039488 History: Received March 13, 2017; Revised February 22, 2018

In this paper, we present a discontinuous cytotoxic T cells (CTLs) response for HTLV-1. Moreover, a delay parameter for the activation of CTLs is considered. In fact, a system of differential equation with discontinuous right-hand side with delay is defined for HTLV-1. For analyzing the dynamical behavior of the system, graphical Hopf bifurcation is used. In general, Hopf bifurcation theory will help to obtain the periodic solutions of a system as parameter varies. Therefore, by applying the frequency domain approach and analyzing the associated characteristic equation, the existence of Hopf bifurcation by using delay immune response as a bifurcation parameter is determined. The stability of Hopf bifurcation periodic solutions is obtained by the Nyquist criterion and the graphical Hopf bifurcation theorem. At the end, numerical simulations demonstrated our results for the system of HTLV-1.

Copyright © 2018 by ASME
Your Session has timed out. Please sign back in to continue.


Welsh, J. S. , 2011, “ Contagious Cancer,” Oncologist, 16(1), pp. 1–4. [CrossRef] [PubMed]
Koyanagi, Y. , Itoyama, Y. , Nakamura, N. , Takamatsu, K. , Kira, J. I. , Iwamasa, T. , Goto, I. , and Yamamoto, N. , 1993, “ In Vivo Infection of Human T-Cell Leukemia Virus Type I in Non-T Cells,” Virology, 196(1), pp. 25–33. [CrossRef] [PubMed]
Bai, Z. , and Zhou, Y. , 2012, “ Dynamics of a Viral Infection Model With Delayed CTL Response and Immune Circadian Rhythm,” Chaos, Solitons Fractals, 45(9–10), pp. 1133–1139. [CrossRef]
Asquith, B. , and Bangham, C. R. , 2007, “ Quantifying HTLV-I Dynamics,” Immunol. Cell Biol., 85(4), pp. 280–286. [CrossRef] [PubMed]
Bangham, C. R. , 2000, “ The Immune Response to HTLV-I,” Curr. Opinion Immunol., 12(4), pp. 397–402. [CrossRef]
Lang, J. , and Li, M. Y. , 2012, “ Stable and Transient Periodic Oscillations in a Mathematical Model for CTL Response to HTLV-I Infection,” J. Math. Biol., 65(1), pp. 181–199. [CrossRef] [PubMed]
Bangham, C. R. , 2003, “ Human T-Lymphotropic Virus Type 1 (HTLV-1): Persistence and Immune Control,” Int. J. Hematol., 78(4), pp. 297–303. [CrossRef] [PubMed]
Bangham, C. R. , 2003, “ The Immune Control and Cell-to-Cell Spread of Human T-Lymphotropic Virus Type 1,” J. Gen. Virol., 84(Pt. 12), pp. 3177–3189. [CrossRef] [PubMed]
Fan, R. , Dong, Y. , Huang, G. , and Takeuchi, Y. , 2014, “ Apoptosis in Virus Infection Dynamics Models,” J. Biol. Dyn., 8(1), pp. 20–41. [CrossRef] [PubMed]
Shamsara, E. , Afsharnejad, Z. , and Mostolizadeh, R. , 2017, “ Hopf Bifurcation for a Discontinuous Htlv-1 Model,” FILOMAT, 31(20), pp. 6247–6267.
Nowak, M. , and May, R. M. , 2000, Virus Dynamics: Mathematical Principles of Immunology and Virology, Oxford University Press, Oxford, UK.
Gómez-Acevedo, H. , Li, M. Y. , and Jacobson, S. , 2010, “ Multistability in a Model for CTL Response to HTLV-I Infection and Its Implications to HAM/TSP Development and Prevention,” Bull. Math. Biol., 72(3), pp. 681–696. [CrossRef] [PubMed]
Li, M. Y. , and Shu, H. , 2012, “ Global Dynamics of a Mathematical Model for HTLV-I Infection of CD4+ T Cells With Delayed CTL Response,” Nonlinear Anal.: Real World Appl., 13(3), pp. 1080–1092. [CrossRef]
Shamsara, E. , Mostolizadeh, R. , and Afsharnezhad, Z. , 2016, “ Transcritical Bifurcation of an Immunosuppressive Infection Model,” Iranian J. Numer. Anal. Optim., 6(2), pp. 1–16.
Wodarz, D. , 2007, Killer Cell Dynamics, Vol. 32, Springer, Cham, Switzerland. [CrossRef]
Balasubramaniam, P. , Tamilalagan, P. , and Prakash, M. , 2015, “ Bifurcation Analysis of HIV Infection Model With Antibody and Cytotoxic T-Lymphocyte Immune Responses and Beddington–DeAngelis Functional Response,” Math. Methods Appl. Sci., 38(7), pp. 1330–1341. [CrossRef]
Egami, C. , 2009, “ Bifurcation Analysis of the Nowak–Bangham Model in CTL Dynamics,” Math. Biosci., 221(1), pp. 33–42. [CrossRef] [PubMed]
Hethcote, H. W. , 2000, “ The Mathematics of Infectious Diseases,” SIAM Rev., 42(4), pp. 599–653. [CrossRef]
Z. Hu , X. Z., Liu , H. X., Wang , and Ma, W. , 2010, “ Analysis of the Dynamics of a Delayed HIV Pathogenesis Model,” J. Comput. Appl. Math., 234(2), pp. 461–476. [CrossRef]
Li, M. Y. , and Shu, H. , 2011, “ Multiple Stable Periodic Oscillations in a Mathematical Model of CTL Response to HTLV-I Infection,” Bull. Math. Biol., 73(8), pp. 1774–1793. [CrossRef] [PubMed]
Muroya, Y. , Enatsu, Y. , and Li, H. , 2013, “ Global Stability of a Delayed HTLV-I Infection Model With a Class of Nonlinear Incidence Rates and CTLs Immune Response,” Appl. Math. Comput., 219(21), pp. 10559–10573.
Shamsara, E. , Afsharnezhad, Z. , and Efatti, S. , 2018, “Optimal Drug Control of a Four Dimensional HIV Infection Model,” Acta Biotheoritica, in press.
Shu, H. , Wang, L. , and Watmough, J. , 2013, “ Global Stability of a Nonlinear Viral Infection Model With Infinitely Distributed Intracellular Delays and CTL Immune Responses,” SIAM J. Appl. Math., 73(3), pp. 1280–1302. [CrossRef]
Tian, X. , and Xu, R. , 2014, “ Global Stability and Hopf Bifurcation of an HIV-1 Infection Model With Saturation Incidence and Delayed CTL Immune Response,” Appl. Math. Comput., 237, pp. 146–154.
Xie, Q. , Huang, D. , Zhang, S. , and Cao, J. , 2010, “ Analysis of a Viral Infection Model With Delayed Immune Response,” Appl. Math. Modell., 34(9), pp. 2388–2395. [CrossRef]
Afsharnezhad, Z. , and Amaleh, M. K. , 2011, “ Continuation of the Periodic Orbits for the Differential Equation With Discontinuous Right Hand Side,” J. Dyn. Differ. Equations, 23(1), pp. 71–92. [CrossRef]
Akhmet, M. , 2010, Principles of Discontinuous Dynamical Systems, Springer Science & Business Media, New York. [CrossRef] [PubMed] [PubMed]
Gouzé, J. L. , and Sari, T. , 2002, “ A Class of Piecewise Linear Differential Equations Arising in Biological Models,” Dyn. Syst., 17(4), pp. 299–316. [CrossRef]
Elaiw, A. M. , 2010, “ Global Properties of a Class of HIV Models,” Nonlinear Anal.: Real World Appl., 11(4), pp. 2253–2263. [CrossRef]
Shamsara, E. , Shamsara, J. , and Afsharnezhad, Z. , 2016, “ Optimal Control Therapy and Vaccination for Special HIV-1 Model With Delay,” Theory Biosci., 135(4), pp. 217–230. [CrossRef] [PubMed]
Rihan, F. A. , and Rihan, B. F. , 2015, “ Numerical Modelling of Biological Systems With Memory Using Delay Differential Equations,” Appl. Math., 9(3), pp. 1645–1658.
Di Bernardo, M. , Budd, C. J. , Champneys, A. R. , Kowalczyk, P. , Nordmark, A. B. , Tost, G. O. , and Piiroinen, P. T. , 2008, “ Bifurcations in Nonsmooth Dynamical Systems,” SIAM Rev., 50(4), pp. 629–701. [CrossRef]
Filippov, A. F. , 1988, Differential Equations With Discontinuous Right Hand Sides: Control Systems, Springer, Dordrecht, The Netherlands.
Kunze, M. , and Küpper, T. , 2001, Non-Smooth Dynamical Systems: An Overview, Springer, Berlin. [CrossRef]
Leine, R. I. , and Nijmeijer, H. , 2013, Dynamics and Bifurcations of Non-Smooth Mechanical Systems, Vol. 18, Springer Science & Business Media, Berlin.
Kuznetsov, Y. A. , 2013, Elements of Applied Bifurcation Theory, Vol. 112, Springer Science & Business Media, New York.
Moiola, J. L. , and Chen, G. , 1996, Hopf Bifurcation Analysis, World Scientific, Singapore. [CrossRef]
Allwright, D. J. , 1977, “Harmonic Balance and the Hopf Bifurcation,” Math. Proc. Cambridge Philos. Soc., 82(3), pp. 453–467.
Mees, A. , and Chua, L. , 1979, “ The Hopf Bifurcation Theorem and Its Applications to Nonlinear Oscillations in Circuits and Systems,” IEEE Trans. Circuits Syst., 26(4), pp. 235–254. [CrossRef]
Mees, A. I. , 1981, Dynamics of Feedback Systems, Wiley, New York.
Xu, C. , Tang, X. , and Liao, M. , 2010, “ Frequency Domain Analysis for Bifurcation in a Simplified Tri-Neuron BAM Network Model With Two Delays,” Neural Networks, 23(7), pp. 872–880. [CrossRef] [PubMed]
Moiola, J. L. , and Chen, G. , 1993, “ Frequency Domain Approach to Computation and Analysis of Bifurcations and Limit Cycles: A Tutorial,” Int. J. Bifurcation Chaos, 3(4), pp. 843–867. [CrossRef]
Yu, W. , and Cao, J. , 2007, “ Stability and Hopf Bifurcation on a Two-Neuron System With Time Delay in the Frequency Domain,” Int. J. Bifurcation Chaos, 17(4), pp. 1355–1366. [CrossRef]
Irving, R. S. , 2003, Integers, Polynomials, and Rings: A Course in Algebra, Springer Science & Business Media, New York.


Grahic Jump Location
Fig. 5

Behavior and phase portrait at τ=0.6667582188>τ0=0.5667582188

Grahic Jump Location
Fig. 6

The typical solutions converge to a periodic solution when τ>τ0 with the initial value (12,0.6,1)

Grahic Jump Location
Fig. 7

Behavior and phase portrait at τ=0.0667582188<τ0=0.5667582188 with the initial value (12,0.6,1)

Grahic Jump Location
Fig. 8

The typical solutions are asymptotically stable when τ<τ0 with the initial value (12, 2, 1)

Grahic Jump Location
Fig. 1

Behavior and phase portrait at τ=0.1392164589<τ0=0.3182800220

Grahic Jump Location
Fig. 2

The typical solutions are asymptotically stable when τ<τ0 with the initial value (23,0.6,1.5). The unit is cell/mm3 on the vertical axes and days for time t.

Grahic Jump Location
Fig. 3

Behavior and phase portrait at τ=0.5392164589>τ0=0.3182800220

Grahic Jump Location
Fig. 4

The periodic solutions when τ=0.5392164589>τ0=0.3182800220 with the initial value (23,0.6,1.5). The unit is cell/mm3 on the vertical axes and days for time t.



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In