Next Article in Journal
A Novel Solid Nanocrystals Self-Stabilized Pickering Emulsion Prepared by Spray-Drying with Hydroxypropyl-β-cyclodextrin as Carriers
Next Article in Special Issue
Bioconvection Due to Gyrotactic Microorganisms in Couple Stress Hybrid Nanofluid Laminar Mixed Convection Incompressible Flow with Magnetic Nanoparticles and Chemical Reaction as Carrier for Targeted Drug Delivery through Porous Stretching Sheet
Previous Article in Journal
Concentration of Bioactive Phenolic Compounds in Olive Mill Wastewater by Direct Contact Membrane Distillation
Previous Article in Special Issue
The Finite Pore Volume GAB Adsorption Isotherm Model as a Simple Tool to Estimate a Diameter of Cylindrical Nanopores
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fractional Dynamics of HIV with Source Term for the Supply of New CD4+ T-Cells Depending on the Viral Load via Caputo–Fabrizio Derivative

1
Department of Mathematics, University of Lakki Marwat, Lakki Marwat 28420, Pakistan
2
Center of Excellence in Theoretical and Computational Science (TaCS-CoE), Faculty of Science, King Mongkut’s, University of Technology Thonburi (KMUTT), 126 Pracha Uthit Rd., Bang Mod, Thung Khru, Bangkok 10140, Thailand
3
School of Mathematics and Statistics, Xi’an Jiaotong University, Xi’an 710049, China
4
KMUTT Fixed Point Research Laboratory, KMUTT-Fixed Point Theory and Applications Research Group, SCL 802 Fixed Point Laboratory, Department of Mathematics, Faculty of Science, King Mongkut’s, University of Technology Thonburi (KMUTT), 126 Pracha Uthit Rd., Bang Mod, Thung Khru, Bangkok 10140, Thailand
5
Department of Medical Research, China Medical University Hospital, China Medical University, Taichung 40402, Taiwan
6
Department of Mathematics, College of Science & Arts, King Abdulaziz University, P.O. Box 344, Rabigh 21911, Saudi Arabia
*
Author to whom correspondence should be addressed.
Molecules 2021, 26(6), 1806; https://doi.org/10.3390/molecules26061806
Submission received: 21 February 2021 / Revised: 5 March 2021 / Accepted: 5 March 2021 / Published: 23 March 2021

Abstract

:
Human immunodeficiency virus (HIV) is a life life-threatening and serious infection caused by a virus that attacks CD 4 + T-cells, which fight against infections and make a person susceptible to other diseases. It is a global public health problem with no cure; therefore, it is highly important to study and understand the intricate phenomena of HIV. In this article, we focus on the numerical study of the path-tracking damped oscillatory behavior of a model for the HIV infection of CD 4 + T-cells. We formulate fractional dynamics of HIV with a source term for the supply of new CD4+ T-cells depending on the viral load via the Caputo–Fabrizio derivative. In the formulation of fractional HIV dynamics, we replaced the constant source term for the supply of new CD 4 + T-cells from the thymus with a variable source term depending on the concentration of the viral load, and introduced a term that describes the incidence of the HIV infection of CD 4 + T-cells. We present a novel numerical scheme for fractional view analysis of the proposed model to highlight the solution pathway of HIV. We inspect the periodic and chaotic behavior of HIV for the given values of input factors using numerical simulations.

1. Introduction

Mathematical biology has a wide range of applications in genetics, environmental sciences, population dynamics, medical sciences, etc. The field may be expressed as biomathematics depending on whether dealing with mathematical aspects or biological aspects [1]. Several biological phenomena and processes can be formulating in the language of mathematics in the framework of ordinary differential equations, fractional differential equations, impulsive differential equations, stochastic differential equations, delay differential equations, etc. [2,3,4,5,6]. Different assumptions, laws, axioms governing these processes are used in the formulation of these mathematical models to show the intricate dynamics of biological phenomena. Presently, the most critical and dangerous global health issue is human immunodeficiency virus (HIV), and numerous mathematical models have been developed of the human immune system to highlight the overall picture of the infection that describes the infection of HIV and its interaction with the immune system. It is reported that HIV is the causative agent for obtaining acquired immune deficiency syndrome (AIDS), which damages the body’s ability to fight other infections. HIV is an incurable deadly disease, and millions of people have lost their lives from it. However, people are living long and healthy lives with HIV due to effective HIV treatment, care, diagnosis, and prevention. When HIV enters the body of a healthy individual, it quickly reproduces and damages the CD 4 + T-cells, which affects the immune system. The symptoms in the primary stage of HIV are fever, night sweats, cough, weight loss, headache, diarrhea, rash, joint pain and muscle aches, swollen lymph glands, and sore throat. At this stage the load of viruses in the blood is high, and infection easily spreads in the body compared to other forthcoming stages. Furthermore, HIV infects via bodily fluids (blood, tears, urine, spit, etc.).
It is noted that CD 4 + T-cells play a role in the human body and are the target cells of HIV; moreover, these CD 4 + T-cells perform a key role in the adjustment of the immune order, so the declination of these cells can have widespread influence on causing damage to a working immune system. The level of concentration of these cells is used to determine the stage of HIV infection; therefore, it is valuable to determine its importance using a mathematical framework. Numerous mathematical models have been constructed to visualize the transmission phenomena of CD 4 + T-cells with HIV. In the literature, Perelson et al. [7] developed a simple mathematical system that describes the dynamics of HIV infection in CD 4 + T-cells. Furthermore, Perelson and Nelson [8] developed another model based on [7], which consists of four variables: the population of free HIV particles, latently infected, actively infected, and uninfected cells. The model demonstrates many characteristics observed clinically of AIDS, e.g., the reduction of CD 4 + T-cells, and the low levels and long latency time of the free virus in a human host body. After that, Raun and Calshaw [9] reduced the model developed by Perelson et al. [8] into one consisting of three independent variables: HIV-free virus particles, and infected and uninfected T-cells. The study by Bushnaq et al. [10] investigated the stability existence of steady states of HIV/AIDS mode with fractional derivatives and explored the role of memory in the dynamics of HIV infection. The authors in [11,12,13,14] implemented different techniques and methods to show the path-tracking behavior of the dynamics of HIV. The main objective of the current study is to formulate the dynamics of HIV with a variable source term instead of a constant term for the supply of new CD 4 + T-cells from the thymus. Furthermore, the mass action term k V T is introduced, which shows the infection caused by the virus V , interacting with uninfected CD 4 + T-cells, causing the loss of free virus at the rate k V T , where k is the rate of infection.
It is well known that fractional-order systems offer more reliable, deeper, more valuable and more accurate information about the dynamics of different diseases than classical integer-order systems [15,16,17,18]. Hereditary properties and a description of memory make them superior to integer-order models [19,20,21,22,23]; moreover, fractional-order models can easily explore and demonstrate the dynamics between two points. These new ideas have been effectively used in modeling real-world problems in physics, engineering, biology, economics, and several other areas [24,25,26,27,28,29,30,31]. In advance studies and in recent research, fractional calculus produced more accurate information about the dynamics of a system; in particular the dynamical behavior of infectious diseases can be more accurately highlighted through fractional calculus. In fractional calculus, Riemann–Liouville, Caputo, Hilfer and some other operators using the power-law kernel and encounter limitations in modeling natural phenomena. On the other hand, the Caputo–Fabrizio operator possess an exponential decay law kernel which produce more accurate results for natural phenomena. This new developed fractional operator contains stronger properties as compared to the other fractional operators. Therefore, we opt to investigate the dynamical behavior of HIV in the framework of the fractional derivative under the Caputo–Fabrizio operator to capture more accurate and more reliable information about the infection.
Motivated by the above accurate results of the fractional derivative, we opt to explore the fractional dynamics of HIV with a source term for the supply of new CD4+ T-cells depending on the viral load using the Caputo–Fabrizio derivative. We represent the basic results of fractional Caputo–Fabrizio (CF) derivative in the first section of the article. In the second section, a fractional model is formulated for HIV with a source term for the supply of new CD4+ T-cells depending on the viral load in the framework of Caputo–Fabrizio derivative. The proposed model of HIV infection is then investigated through mathematical skills. In the third section, a novel numerical scheme is presented to visualize the tracking path of HIV in the fractional framework. The dynamics of HIV infection is then investigated with the effect of fractional order numerically in section four. We observe that the input parameter ϑ has a significant effect on the dynamics of HIV infection and decreases the level of infection. In addition, we show the chaotic behavior of the system with variation of input parameters to conceptualize the chaotic pattern of HIV that justifies the suitability of the scheme. Our analysis predicts that fractional order can be used as a control parameter, and should be suggested to policymakers. Finally, concluding remarks and suggestions are presented in last section.

2. Structure of HIV Dynamics

To study HIV/AIDS, the most important requirement is to conceptualize the population dynamics of CD 4 + T-cells. These cells are produced in the bone marrow, and premature cells are shifted to the thymus, where they receive extra differentiation and mature into immunocompetent CD 4 + T-cells. In humans, the thymus reaches its maximum weight at the time of puberty and then slowly becomes complex. Thymus ejection from an adult has small effect, even though the thymus in adults is in a working state and its few lymphocytes perform as precursors of T-cells and immunocompetent T-cells. The focus of the presented model is on CD 4 + T-cells. The succession of HIV infection can be calculated by the number of CD 4 + T-cells, which tells us about initial symptoms. The interest of the present study is in observing the HIV infection model of CD 4 + T-cells, which performs path-tracking damped oscillations. The mathematical model presented here depicts the behavior of the infected rate of CD 4 + T-cells. The above assumptions of [13] lead us to the following system:
d T d t = s μ T T + r T ( 1 T + I T max ) k V T , d I d t = k V T μ I I , d V d t = N μ I I μ V V ,
where the three terms T ( t ) , I ( t ) , and V ( t ) illustrate the quantity of healthy CD 4 + T-cells in the blood, infected CD 4 + T-cells, and free HIV particles, respectively. The values of parameters, constants, and initial values of state variables and their explanation of the above-mentioned model are given in Table 1. In the upcoming section, the model will be extended in a fractional framework.

3. Fractional Dynamics of HIV Infection

In this section, a variable source term s ( V ) = s exp ( κ V ) is incorporated into the above system of HIV infection instead of constant the terms used in [13]. The new source term expresses the generated number of new healthy T-cells from the thymus, which depends on the concentration of the viral load. The production of healthy T-cells is observed to decrease due to the increase in the size of viral load; therefore, the source term is considered to be variable rather than constant.
The numerical solution of the new model using standard parameters with κ = 1 mm 3 , performs the key role in the depletion of healthy T-cells in such a way that their depletion is more gradual than in the previous model. Here, s indicates the concentration of new healthy T-cells that are generated from sources within the system, such as the thymus. Healthy T-cells can be increased rapidly through existing healthy T-cells, and κ is the constant. The infection model can be solved by a technique in which the system is extended by introducing a term k V T . Here, k behaves as a constant rate of infection. This kind of term has logical importance in such a way that the virus necessarily infects the healthy T-cells, so the chance of attacking healthy T-cells at a low rate can be assumed to be proportional to the product of V and T. For this reason, the infection takes place when the virus communicates with uninfected T-cells. First, the virus is induced in heathy T-cells, which makes these heathy T-cells irresistible at a rate of k V T and then produces infected T-cells at a rate of k V T . During HIV infection, the rate of the virus never exceeds the number of healthy T-cells. Infection can also take place via the communication of cells with each other, i.e., by meeting infected T-cells with healthy T-cells. In the presented model, V and T depicts the virus rate and measure of healthy T-cell count in the blood, as, obviously, there is no compulsion for the infected virus to just stay in the blood. In fact, a very large number of healthy T-cells is in lymphoid tissues. However, according to the given data, the calculated number of viruses and healthy T-cells in the blood adequately indicates their volume in the whole body [8,32,33], just as one would assume for a system to be in equilibrium. This extended model gives more accurate and more realistic forecasting of the healthy T-cell depletion for a given time, and can be written in the following way in the fractional framework:
0 C F D t ϑ T = s exp ( κ V ) μ T T + r T ( 1 T + I T max ) k V T , 0 C F D t ϑ I = k V T μ I I , 0 C F D t ϑ V = N μ I I μ V V k V T ,
where 0 C F D t ϑ indicates Caputo–Fabrizio fractional derivative of order ϑ . In the next part of the paper, we will discuss and represent the rudimentary knowledge of the Caputo–Fabrizio fractional derivative, which will be used for the analysis of our proposed HIV infection model.

Fundamental Knowledge and Concept

Here, we will represent the basic definitions and results important for the analysis of HIV with the source term for the supply of new CD4+ T-cells depending on the viral load in a fractional Caputo–Fabrizio(CF) framework. The fundamental knowledge of Caputo–Fabrizio fractional derivative is given as follows:
Definition 1.
Let us suppose h H 1 ( a , b ) , where b is greater than a, then the CF derivative [31] of order ϑ is given by
D t ϑ ( h ( t ) ) = U ( ϑ ) 1 ϑ a t h ( x ) exp ϑ t x 1 ϑ d x ,
where ϑ [ 0 , 1 ] and U ( τ ) denotes normality with U ( 0 ) = U ( 1 ) = 1 [31]. In this case, when h H 1 ( a , b ) , then the following fractional derivative is obtained
D t ϑ ( h ( t ) ) = ϑ U ( ϑ ) 1 ϑ a t ( h ( t ) h ( x ) ) exp ϑ t x 1 ϑ d x .
Remark 1.
Let us take α = 1 ϑ ϑ [ 0 , ) and ϑ = 1 1 + α [ 0 , 1 ] , then Equation (4) can be written in the following form
D t ϑ ( h ( t ) ) = M ( α ) α a t h ( x ) e [ t x α ] d x , M ( 0 ) = M ( ) = 1 .
In addition to this,
lim α 0 1 α exp t x α = δ ( x t ) .
In the next step, we will define the fractional integral, which was introduced in [34].
Definition 2.
Let h be a given function then the fractional integral is defined in the following manner
I t ϑ ( h ( t ) ) = 2 ( 1 ϑ ) ( 2 ϑ ) U ( ϑ ) h ( t ) + 2 ϑ ( 2 ϑ ) U ( ϑ ) 0 t h ( u ) d u , t 0 .
where 0 < ϑ < 1 and is the order of the above fractional integral.
Remark 2.
Further analysis of the above Definition 2 gives that
2 ( 1 ϑ ) ( 2 ϑ ) U ( ϑ ) + 2 ϑ ( 2 ϑ ) U ( ϑ ) = 1 ,
which gives U ( ϑ ) = 2 2 ϑ , 0 < ϑ < 1 . A new Caputo derivative of order ϑ was introduced by Nieto and Losada in [34] by using Equation (8) and is given by
D t ϑ ( h ( t ) ) = 1 1 ϑ 0 t h ( x ) exp ϑ t x 1 ϑ d x , 0 < ϑ < 1 .

4. Novel Numerical Scheme for Fractional Model

In this section, the main aim is to illustrate the dynamical behavior of the proposed fractional model (2) of HIV infection, numerically. Several numerical approaches have been introduced in the literature for Caputo–Fabrizio fractional systems to highlight their dynamics [35,36,37]. Here, we will use the scheme of [37] to represent the dynamics of our fractional system (2) of HIV infection, which is more reliable, easily implementable, and stable. The values of the parameter presented in Table 1 will be used for simulation purposes. First, we represent our proposed fractional model of HIV infection in Volterra type and then apply the idea of a fundamental theorem of integration for further simplification. To obtain the numerical scheme, the first equation of our proposed model of HIV infection implies:
z 1 ( t ) z 1 ( 0 ) = 1 ϑ U ( ϑ ) H 1 ( t , z 1 ) + ϑ U ( ϑ ) 0 t H 1 ( ω , z 1 ) d ω .
Next, we take the time t = t m + 1 , m = 0 , 1 , , and get the following
z 1 ( t m + 1 ) z 1 ( 0 ) = 1 ϑ U ( ϑ ) H 1 ( t m , z 1 ( t m ) ) + ϑ U ( ϑ ) 0 t m + 1 H 1 ( t , z 1 ) d t .
and
z 1 ( t m ) z 1 ( 0 ) = 1 ϑ U ( ϑ ) H 1 ( t m 1 , z 1 ( t m 1 ) ) + ϑ U ( ϑ ) 0 t m H 1 ( t , z 1 ) d t .
Here, the difference in the successive terms of the system are obtained as follows:
z 1 m + 1 z 1 m = 1 ϑ U ( ϑ ) H 1 ( t m , z 1 m ) H 1 ( t m 1 , z 1 m 1 ) + ϑ U ( ϑ ) m t m + 1 H 1 ( t , z 1 ) d t .
Furthermore, we approximate the above-mentioned function H 1 ( t , z 1 ) in the time interval [ t k , t k + 1 ] with the help of interpolation polynomial and get
P k ( t ) H 1 ( t k , z k ) h ( t t k 1 ) H 1 ( t k 1 , z k 1 ) h ( t t k ) ,
in which h is the time spent and h = t m t m 1 . The expression of P k ( t ) is used to determine the value of the following integral
m t m + 1 H 1 ( t , z 1 ) d t = m t m + 1 H 1 ( t m , z 1 m ) h ( t t m 1 ) H 1 ( t m 1 , z 1 m 1 ) h ( t t m ) d t , = 3 h 2 H 1 ( t m , z 1 m ) h 2 H 1 ( t m 1 , z 1 m 1 ) .
Here, substituting the value of (15) in Equation (13), we get the below
z 1 m + 1 = z 1 m + 1 ϑ U ( ϑ ) + 3 ϑ h 2 U ( ϑ ) H 1 ( t m , z 1 m ) 1 ϑ U ( ϑ ) + ϑ h 2 U ( ϑ ) H 1 ( t m 1 , z 1 m 1 ) ,
which is the required scheme for first equation of model (2) of HIV infection. Following the same procedure, we can determine the required scheme for the second and third equations of the proposed model (2) of HIV infection given by
z 2 m + 1 = z 2 m + 1 ϑ U ( ϑ ) + 3 ϑ h 2 U ( ϑ ) H 2 ( t m , z 2 m ) 1 ϑ U ( ϑ ) + ϑ h 2 U ( ϑ ) H 2 ( t m 1 , z 2 m 1 ) ,
and
z 3 m + 1 = z 3 m + 1 ϑ U ( ϑ ) + 3 ϑ h 2 U ( ϑ ) H 3 ( t m , z 3 m ) 1 ϑ U ( ϑ ) + ϑ h 2 U ( ϑ ) H 3 ( t m 1 , z 3 m 1 ) .
This method is a two-step Adams–Bashforth scheme for a Caputo–Fabrizio operator, which takes into account the nonlinearity of the kernel including exponential decay law for the Caputo–Fabrizio operator. In [37], the stability and convergence of the proposed scheme are presented in detail. We perform several simulations to understand the intricate dynamics of HIV, with the source term for the supply of new CD4+ T-cells depending on the viral load in the fractional framework. The values of input parameters and initial values of state variables given in Table 1 are used for simulation purposes. In Figure 1, Figure 2, Figure 3 and Figure 4, we demonstrate the dynamical behavior of healthy CD4+ T-cells, infected CD4+ T-cells, and HIV-free virus by varying the fractional-order ϑ , i.e., ϑ = 1.0 , 0.85 , 0.65 , 0.45 to observe the influence of fractional-order ϑ on the dynamics of HIV infection. We noticed that the index of memory or fractional order is an effective input parameter and can play an important role in visualizing a more accurate picture of the system; furthermore, it can be used as a preventive parameter for HIV infection. In Figure 5, Figure 6, Figure 7 and Figure 8, we demonstrate the chaotic behavior of our fractional model (2) of HIV infection with different values of fractional-order ϑ . The chaotic behavior of a system is very important in many scientific and engineering applications. It is well known that there is a strong tendency towards the formulation and visualization of chaotic behaviors of a system. The chaotic patterns justify the suitability and applicability of the proposed numerical scheme, and can be further applied to novel chaotic systems. From our simulations, we observe that the input parameter ϑ highly contributes to and can be used as an effective parameter for preventive measures.
In our analysis, we mainly focused on the solution pathway and chaotic behavior of the HIV infection. In addition to this, we showed the influence of fractional parameter ϑ on the system to observe whether the fractional order can be used as a preventive measure or not. It should be noted in simulations that the proposed numerical scheme is less time-consuming and a more easily implementable scheme. However, the investigation of effectiveness of this scheme such as stability, convergence, consistency, accuracy, and computational cost are left for future research.

5. Conclusions

In this article, we formulated the fractional dynamics of HIV with a source term for the supply of new CD4+ T-cells depending on the viral load via a Caputo–Fabrizio derivative. In our model, we incorporated a variable source term for the supply of new CD 4 + T-cells from the thymus depending on the concentration of the viral load. A rudimentary knowledge of fractional calculus in a Caputo–Fabrizio framework has been presented for the analysis of the model. We presented a novel numerical scheme for the solution of the Caputo–Fabrizio fractional model of HIV infection. The dynamics of HIV infection are demonstrated numerically with different values of fractional-order ϑ . The influence of fractional order on the chaotic behavior of our proposed model is illustrated. It is predicted that the index of memory ϑ has a positive effect on the system and can be used as a control parameter. Based on our analysis, we suggest that policymakers should not ignore the role of fractional order in the control of infections. In future research, we will compare the proposed numerical method with other existing numerical techniques and the dynamics of HIV infection will be presented with a time delay visualizing the effect of time delay on HIV infection.

Author Contributions

Conceptualization, Z.S. and R.J.; methodology, R.J.; software, Z.S. and M.S.; validation, P.K., W.D. and R.J.; formal analysis, R.J.; investigation, Z.S.; resources, P.K.; writing—original draft preparation, R.J. and W.D.; writing—review and editing, M.S. and R.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

The authors acknowledge the financial support provided by the Center of Excellence in Theoretical and Computational Science (TaCS-CoE), KMUTT. Moreover, this research project is supported by Thailand Science Research and Innovation (TSRI) Basic Research Fund: Fiscal year 2021 under project number 64A306000005.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, L.; Li, M.Y. Mathematical analysis of the global dynamics of a model for HIV infection of CD4+ T cells. Math. Biosci. 2006, 200, 44–57. [Google Scholar] [CrossRef]
  2. Antar, A.A.; Jenike, K.M.; Jang, S.; Rigau, D.N.; Reeves, D.B.; Hoh, R.; Krone, M.R.; Keruly, J.C.; Moore, R.D.; Schiffer, J.T.; et al. Longitudinal study reveals HIV-1–infected CD4+ T cell dynamics during long-term antiretroviral therapy. J. Clin. Investig. 2020, 130. [Google Scholar] [CrossRef] [Green Version]
  3. Jan, R.; Xiao, Y. Effect of partial immunity on transmission dynamics of dengue disease with optimal control. Math. Method. Appl. Sci. 2019, 42, 1967–1983. [Google Scholar] [CrossRef]
  4. Tamilalagan, P.; Karthiga, S.; Manivannan, P. Dynamics of fractional order HIV infection model with antibody and cytotoxic T-lymphocyte immune responses. J. Comput. Appl. Math. 2021, 382, 113064. [Google Scholar] [CrossRef]
  5. Jan, R.; Xiao, Y. Effect of pulse vaccination on dynamics of dengue with periodic transmission functions. Adv. Differ. Equ. 2019, 2019, 368. [Google Scholar] [CrossRef] [Green Version]
  6. Raza, A.; Ahmadian, A.; Rafiq, M.; Salahshour, S.; Naveed, M.; Ferrara, M.; Soori, A.H. Modeling the effect of delay strategy on transmission dynamics of HIV/AIDS disease. Adv. Differ. Equ. 2020, 2020, 663. [Google Scholar] [CrossRef] [PubMed]
  7. Perelson, A.S.; Kirschner, D.E.; De Boer, R. Dynamics of HIV infection of CD4+ T cells. Math. Biosci. 1993, 114, 81–125. [Google Scholar] [CrossRef] [Green Version]
  8. Perelson, A.S.; Nelson, P.W. Mathematical analysis of HIV-1 dynamics in vivo. Siam. Rev. 1999, 41, 3–44. [Google Scholar] [CrossRef]
  9. Culshaw, R.V.; Ruan, S. A delay-differential equation model of HIV infection of CD4+ T-cells. Math. Biosci. 2000, 165, 27–39. [Google Scholar] [CrossRef]
  10. Bushnaq, S.A.M.I.A.; Khan, S.A.; Shah, K.; Zaman, G. Existence theory of HIV-1 infection model by using arbitrary order derivative of without singular kernel type. J. Math. Anal. 2018, 9, 16–28. [Google Scholar]
  11. Ghoreishi, M.; Ismail, A.M.; Alomari, A.K. Application of the homotopy analysis method for solving a model for HIV infection of CD4+ T-cells. Math. Comput. Model. 2011, 54, 3007–3015. [Google Scholar] [CrossRef]
  12. Ongun, M. The laplace adomian decomposition method for solving a model for HIV infection of CD4+ 220 T-cells. Math. Comput. Model. 2011, 53, 597–603. [Google Scholar] [CrossRef]
  13. Vazquez-Leal, H.; Hernandez-Martinez, L.; Khan, Y.; Jimenez-Fernandez, V.M.; Filobello-Nino, U.; Diaz-Sanchez, A.; Herrera-May, A.L.; Castaneda-Sheissa, R.; Marin-Hernandez, A.; Rabago-Bernal, F.; et al. Multistage HPM applied to path tracking damped oscillations of a model for HIV infection of CD4+ T cells. J. Adv. Math. Comput. Sci. 2014, 1035–1047. [Google Scholar] [CrossRef] [PubMed]
  14. Abdel-Aty, A.H.; Khater, M.M.; Dutta, H.; Bouslimi, J.; Omri, M. Computational solutions of the HIV-1 infection of CD4+ T-cells fractional mathematical model that causes acquired immunodeficiency syndrome (AIDS) with the effect of antiviral drug therapy. Chaos. Soliton. Fract. 2020, 139, 110092. [Google Scholar] [CrossRef] [PubMed]
  15. Arshad, S.; Defterli, O.; Baleanu, D. A second order accurate approximation for fractional derivatives with singular and non-singular kernel applied to a HIV model. Appl. Math. Comput. 2020, 374, 125061. [Google Scholar] [CrossRef]
  16. Jan, A.; Jan, R.; Khan, H.; Zobaer, M.S.; Shah, R. Fractional-order dynamics of Rift Valley fever in ruminant host with vaccination. Commun. Math. Biol. Neurosci. 2020, 2020, 79. [Google Scholar]
  17. Baleanu, D.; Mohammadi, H.; Rezapour, S. Analysis of the model of HIV-1 infection of CD4+ CD4+ T-cell with a new approach of fractional derivative. Adv. Differ. Equ. 2020, 2020, 71. [Google Scholar] [CrossRef] [Green Version]
  18. Jan, R.; Khan, M.A.; Kumam, P.; Thounthong, P. Modeling the transmission of dengue infection through fractional derivatives. Chaos. Soliton. Fract. 2019, 127, 189–216. [Google Scholar] [CrossRef]
  19. Podlubny, I. Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications; Elsevier: Cambridge, MA, USA, 1998. [Google Scholar]
  20. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives; Yverdon-les-Bains, Gordon and Breach Science Publishers: Yverdon, Switzerland, 1993; Volume 1. [Google Scholar]
  21. Baleanu, D.; Aydogn, S.M.; Mohammadi, H.; Rezapour, S. On modelling of epidemic childhood diseases with the Caputo-Fabrizio derivative by using the Laplace Adomian decomposition method. Alex. Eng. J. 2020, 59, 3029–3039. [Google Scholar] [CrossRef]
  22. Sene, N. SIR epidemic model with Mittag-Leffler fractional derivative. Chaos. Soliton. Fract. 2020, 137, 109833. [Google Scholar] [CrossRef]
  23. Bushnaq, S.; Khan, S.A.; Shah, K.; Zaman, G. Mathematical analysis of HIV/AIDS infection model with Caputo-Fabrizio fractional derivative. Cogent. Math. Stat. 2018, 5, 1432521. [Google Scholar] [CrossRef]
  24. Brandibur, O.; Kaslik, E.; Mozyrska, D.; Wyrwas, M. Stability Results for Two-Dimensional Systems of Fractional-Order Difference Equations. Mathematics 2020, 8, 1751. [Google Scholar] [CrossRef]
  25. Srivastava, H.M.; Saad, K.M. A comparative study of the fractional-order clock chemical model. Mathematics 2020, 8, 1436. [Google Scholar] [CrossRef]
  26. Abdulhameed, M.; Vieru, D.; Roslan, R. Magnetohydrodynamic electroosmotic flow of Maxwell fluids with Caputo-Fabrizio derivatives through circular tubes. Comput. Math. Appl. 2017, 74, 2503–2519. [Google Scholar] [CrossRef]
  27. Fatmawati, F.; Jan, R.; Khan, M.A.; Khan, Y.; Ullah, S. A new model of dengue fever in terms of fractional derivative. Math. Biosci. Eng. 2020, 17, 5267–5287. [Google Scholar] [CrossRef]
  28. Kumar, S.; Ahmadian, A.; Kumar, R.; Kumar, D.; Singh, J.; Baleanu, D.; Salimi, M. An efficient numerical method for fractional SIR epidemic model of infectious disease by using Bernstein wavelets. Mathematics 2020, 8, 558. [Google Scholar] [CrossRef] [Green Version]
  29. Muresan, C.I.; Birs, I.R.; Dulf, E.H. Event-Based Implementation of Fractional Order IMC Controllers for Simple FOPDT Processes. Mathematics 2020, 8, 1378. [Google Scholar] [CrossRef]
  30. Goufo, E.D. A biomathematical view on the fractional dynamics of cellulose degradation. Fract. Calc. Appl. Anal. 2015, 18, 554–564. [Google Scholar]
  31. Caputo, M.; Fabrizio, M. A new definition of fractional derivative without singular kernel. Progr. Fract. Differ. Appl. 2015, 1, 1–13. [Google Scholar]
  32. Haase, A.T.; Henry, K.; Zupancic, M.; Sedgewick, G.; Faust, R.A.; Melroe, H.; Cavert, W.; Gebhard, K.; Staskus, K.; Zhang, Z.Q.; et al. Quantitative image analysis of HIV-1 infection in lymphoid tissue. Science 1996, 274, 985–989. [Google Scholar] [CrossRef] [Green Version]
  33. Perelson, A.S.; Essunger, P.; Cao, Y.; Vesanen, M.; Hurley, A.; Saksela, K.; Markowitz, M.; Ho, D.D. Decay characteristics of HIV-1-infected compartments during combination therapy. Nature 1997, 387, 188–191. [Google Scholar] [CrossRef] [PubMed]
  34. Losada, J.; Nieto, J.J. Properties of a new fractional derivative without singular kernel. Progr. Fract. Differ. Appl. 2015, 1, 87–92. [Google Scholar]
  35. Liu, Y.; Fan, E.; Yin, B.; Li, H. Fast algorithm based on the novel approximation formula for the Caputo-Fabrizio fractional derivative. Aims. Math. 2020, 5, 1729–1744. [Google Scholar] [CrossRef]
  36. Liu, Z.; Cheng, A.; Li, X. A second-order finite difference scheme for quasilinear time fractional parabolic equation based on new fractional derivative. Int. J. Comput. Math. 2018, 95, 396–411. [Google Scholar] [CrossRef]
  37. Atangana, A.; Owolabi, K.M. New numerical approach for fractional differential equations. Math. Model. Nat. Phenom. 2018, 13, 3. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Illustration of the solution pathway of (a) healthy CD4+ T-cells, (b) infected CD4+ T-cells, (c) free HIV virus particles and (d) chaotic behavior of the proposed fractional model (2) of HIV infection with fractional-order ϑ = 1.
Figure 1. Illustration of the solution pathway of (a) healthy CD4+ T-cells, (b) infected CD4+ T-cells, (c) free HIV virus particles and (d) chaotic behavior of the proposed fractional model (2) of HIV infection with fractional-order ϑ = 1.
Molecules 26 01806 g001
Figure 2. Illustration of the solution pathway of (a) healthy CD4+ T-cells, (b) infected CD4+ T-cells, (c) free HIV virus particles and (d) chaotic behavior of the proposed fractional model (2) of HIV infection with fractional-order ϑ = 0.85.
Figure 2. Illustration of the solution pathway of (a) healthy CD4+ T-cells, (b) infected CD4+ T-cells, (c) free HIV virus particles and (d) chaotic behavior of the proposed fractional model (2) of HIV infection with fractional-order ϑ = 0.85.
Molecules 26 01806 g002
Figure 3. Illustration of the solution pathway of (a) healthy CD4+ T-cells, (b) infected CD4+ T-cells, (c) free HIV virus particles and (d) chaotic behavior of the proposed fractional model (2) of HIV infection with fractional-order ϑ = 0.65.
Figure 3. Illustration of the solution pathway of (a) healthy CD4+ T-cells, (b) infected CD4+ T-cells, (c) free HIV virus particles and (d) chaotic behavior of the proposed fractional model (2) of HIV infection with fractional-order ϑ = 0.65.
Molecules 26 01806 g003
Figure 4. Illustration of the solution pathway of (a) healthy CD4+ T-cells, (b) infected CD4+ T-cells, (c) free HIV virus particles and (d) chaotic behavior of the proposed fractional model (2) of HIV infection with fractional-order ϑ = 0.45.
Figure 4. Illustration of the solution pathway of (a) healthy CD4+ T-cells, (b) infected CD4+ T-cells, (c) free HIV virus particles and (d) chaotic behavior of the proposed fractional model (2) of HIV infection with fractional-order ϑ = 0.45.
Molecules 26 01806 g004
Figure 5. Illustration of the (a) chaotic behavior, (b) phase portrait of healthy CD4+ T-cells verses infected CD4+ T-cells, (c) phase portrait of infected CD4+ T-cells verses free HIV virus particles and (d) phase portrait of healthy CD4+ T-cells verses free HIV virus particles of the fractional model (2) of HIV infection with fractional-order ϑ = 1.
Figure 5. Illustration of the (a) chaotic behavior, (b) phase portrait of healthy CD4+ T-cells verses infected CD4+ T-cells, (c) phase portrait of infected CD4+ T-cells verses free HIV virus particles and (d) phase portrait of healthy CD4+ T-cells verses free HIV virus particles of the fractional model (2) of HIV infection with fractional-order ϑ = 1.
Molecules 26 01806 g005
Figure 6. Illustration of the (a) chaotic behavior, (b) phase portrait of healthy CD4+ T-cells verses infected CD4+ T-cells, (c) phase portrait of infected CD4+ T-cells verses free HIV virus particles and (d) phase portrait of healthy CD4+ T-cells verses free HIV virus particles of the fractional model (2) of HIV infection with fractional-order ϑ = 0.8.
Figure 6. Illustration of the (a) chaotic behavior, (b) phase portrait of healthy CD4+ T-cells verses infected CD4+ T-cells, (c) phase portrait of infected CD4+ T-cells verses free HIV virus particles and (d) phase portrait of healthy CD4+ T-cells verses free HIV virus particles of the fractional model (2) of HIV infection with fractional-order ϑ = 0.8.
Molecules 26 01806 g006
Figure 7. Illustration of the (a) chaotic behavior, (b) phase portrait of healthy CD4+ T-cells verses infected CD4+ T-cells, (c) phase portrait of infected CD4+ T-cells verses free HIV virus particles and (d) phase portrait of healthy CD4+ T-cells verses free HIV virus particles of the fractional model (2) of HIV infection with fractional-order ϑ = 0.6.
Figure 7. Illustration of the (a) chaotic behavior, (b) phase portrait of healthy CD4+ T-cells verses infected CD4+ T-cells, (c) phase portrait of infected CD4+ T-cells verses free HIV virus particles and (d) phase portrait of healthy CD4+ T-cells verses free HIV virus particles of the fractional model (2) of HIV infection with fractional-order ϑ = 0.6.
Molecules 26 01806 g007
Figure 8. Illustration of the (a) chaotic behavior, (b) phase portrait of healthy CD4+ T-cells verses infected CD4+ T-cells, (c) phase portrait of infected CD4+ T-cells verses free HIV virus particles and (d) phase portrait of healthy CD4+ T-cells verses free HIV virus particles of the fractional model (2) of HIV infection with fractional-order ϑ = 0.4.
Figure 8. Illustration of the (a) chaotic behavior, (b) phase portrait of healthy CD4+ T-cells verses infected CD4+ T-cells, (c) phase portrait of infected CD4+ T-cells verses free HIV virus particles and (d) phase portrait of healthy CD4+ T-cells verses free HIV virus particles of the fractional model (2) of HIV infection with fractional-order ϑ = 0.4.
Molecules 26 01806 g008
Table 1. List of parameters and state variables with corresponding values used in the fractional model of HIV infection.
Table 1. List of parameters and state variables with corresponding values used in the fractional model of HIV infection.
SymbolsInterpretation of Parameters and State VariablesValues
kRate constant for healthy T-cells which become infected by free virus 2.4 × 10 5 days 1
T 0 Population of healthy T-cellsAssumed
I 0 Population of infected T-cellsAssumed
sThe supply rate of healthy T-cells from precursors 0.1 mm 3
rGrowth rate of the healthy T-cells population 3 day 1
μ T Death rate of healthy T-cells 0.02 day 1
μ I Death rate of latently infected T-cells 0.3 day 1
μ V Death rate of free virus of HIV infection 2.4 day 1
NNumber of virus produced by infected T-cellsAssumed
T max Maximum population level of healthy T-cells 1500 mm 3
V 0 Population of free HIV virusAssumed
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shah, Z.; Jan, R.; Kumam, P.; Deebani, W.; Shutaywi, M. Fractional Dynamics of HIV with Source Term for the Supply of New CD4+ T-Cells Depending on the Viral Load via Caputo–Fabrizio Derivative. Molecules 2021, 26, 1806. https://doi.org/10.3390/molecules26061806

AMA Style

Shah Z, Jan R, Kumam P, Deebani W, Shutaywi M. Fractional Dynamics of HIV with Source Term for the Supply of New CD4+ T-Cells Depending on the Viral Load via Caputo–Fabrizio Derivative. Molecules. 2021; 26(6):1806. https://doi.org/10.3390/molecules26061806

Chicago/Turabian Style

Shah, Zahir, Rashid Jan, Poom Kumam, Wejdan Deebani, and Meshal Shutaywi. 2021. "Fractional Dynamics of HIV with Source Term for the Supply of New CD4+ T-Cells Depending on the Viral Load via Caputo–Fabrizio Derivative" Molecules 26, no. 6: 1806. https://doi.org/10.3390/molecules26061806

APA Style

Shah, Z., Jan, R., Kumam, P., Deebani, W., & Shutaywi, M. (2021). Fractional Dynamics of HIV with Source Term for the Supply of New CD4+ T-Cells Depending on the Viral Load via Caputo–Fabrizio Derivative. Molecules, 26(6), 1806. https://doi.org/10.3390/molecules26061806

Article Metrics

Back to TopTop