Abstract

This paper aims to explore the dynamic characteristics of the post treatment human immunodeficiency virus (HIV) type-1 model by proposing the theoretical frameworks. Distinct from the previous works, this study explores the effect of effector cells, loss of functional effector cells, and two types of anti-retroviral therapies such as reverse transcriptase inhibitors (RTIs) and protease inhibitors (PIs) and also the effect of intracellular time delay. Based on the Routh—Hurwitz criterion and eigenvalue analysis, the stability of the proposed HIV-1 model is analyzed. To reveal the significance of time delay, the Hopf-type bifurcation analysis is performed. The optimal control algorithm is designed by choosing the antiviral therapies such as RTI and PI as control parameters. Numerical simulations are performed to validate the effectiveness of the proposed theoretical frameworks.

1. Introduction

According to the World Health Organization (WHO), 680,000 individuals died from HIV-1-related diseases worldwide in 2020, ranging from 480,000 to 1.0 million, while 1.5 million people were newly infected with HIV-1. Antiretroviral therapy (ART) has drastically reduced the number of people infected with HIV during the 1990s, with 27.5 million (approximately) people undergoing treatment in 2020. Because of antiretroviral medication, the HIV infection rate has decreased by 49% over the last two decades, from 2000 to 2020. Despite the researcher’s valiant efforts in terms of treatment options and drugs, a cure for HIV-1 remains a pipe dream, necessitating lifelong treatment. Mathematical models have been demonstrated to be a useful tool for comprehending the dynamics of disease progression, identifying key determinants, and evaluating the efficacy of antiretroviral therapy.

During the 1990s, HIV was thought to be a lethal disease, similar to other lentiviruses, because HIV remains within the host without causing symptoms and progresses to a chronic stage known as acquired immunodeficiency virus (AIDS), with a nearly ten-year delay between HIV and AIDS. In this case, mathematical modeling of HIV can help in estimating the lifespan of infected cells, evaluating therapeutic efficacy, and realizing that new virions require a host with deoxyribonucleic acid (DNA) to replicate. Currently, a variety of therapy options are available to help people with HIV get better. However, controlling the virus is the only option available and curing the disease still requires seamless efforts in the research domain.

Mathematical models have been created to investigate the dynamic properties of cell populations using parameters including latent reservoirs, immunological responses, total carrying capacity, and time-delayed fractional differential operators based on literature reviews (see [15]). Immature infected cells also known as latent stage are those that have been infected but are not yet infectious, which is considered in the present study. Immune responses such as CD4+ T-cells/CD8+ T-cells are activated when a foreign agent enters the body, causing the body’s alarm system to go into overdrive [6, 7].

When systems are in motion, there will always be a degree of lag time. Time delays are inevitable because it has an ability to cause a significant impact on cell populations. When modeling the kinetics of HIV-1 infection, two types of delays are taken into account. One situation is that the uninfected cells interact with infected/free virions and there is an intracellular temporal delay. Besides, after a foreign agent has been ingested, immune response cells must be activated, which results in a delay in the immunological response. This study aims to explore the effect of intracellular time delay with respect to infected cell population [8, 9]. Initially, ART is given to every primarily infected HIV person, but based on the stage of infection the level of drugs may be redefined. However, it is a challenging task with respect to the immune boosters which may vary in the individual. Besides, it can be seen that, if there is a change in the period of drugs provided, then it will reflect in the virion populations. Hence, this study explores the effect of antiretroviral therapies which are suggested as posttreatment for a long period of time. The proposed model’s stability analysis will provide some insights into disease progression in relation to the system parameters such as infection rate and time delays. Followed by, bifurcation analysis is used to determine the threshold value of the significant parameter that has the potential to cause fluctuations in the cell population. The optimal control algorithm is designed by selecting cell populations, which aids in better understanding and extraction of system parameters, resulting in cell populations with stable staining [1013]. Recently, the models were proposed on Zika virus, HIV, SARS-CoV-2, and other viral dynamics models such as maize streak virus in maize spread by leaf hopper mostly in Africa, canine distemper virus, and rabies epidemics in red fox with respect to significant factors such as vaccination parameter involved models and optimal control strategies (for more details, refer [1422]). Distinct from the existing models, the present study focuses on considering the effect of time delays and two kinds of antiretroviral therapies. The overall contribution is listed in the following:(1)This paper models the dynamics of HIV-1 infection by considering the factors such as healthy CD4+ cells, latent reservoirs, infected cells, free virions, and immune responses also considering the effect of time delay.(2)Positivity and boundedness of solutions of the differential model are proved. The reproduction number is determined through the next-generation matrix, which helps to identify the community spread.(3)Conditions for the existence of Hopf bifurcation are proved by choosing the intracellular time delay as a bifurcation parameter.(4)Optimal control algorithm is designed to ensure the stabilization of the proposed model.(5)Numerical simulations are performed to validate the proposed theoretical frameworks.

2. Model of Posttreatment HIV-1 Viral Dynamics

The schematic representation of our model is given below.

In model (1), the first equation represents the rate of change in the susceptible cell populations, where denotes the rate of production, is decay rate and denotes the efficacy of antiretroviral therapies. Consider that , where , in ; if , where , then cent percent the treatment is effective which makes the infection zero, and if , , then, no progress in the therapy. represents the rate of infection between the infected and uninfected cell population. stands for intracellular time delay. The second equation describes the state of latent infection; that is, target cells are infected but not yet infectious. Suppose if the infected cells are matured enough to infect the susceptible cells. , , and are scalars. The third equation explains the rate of change in the infectious cell population with the death rate , migrated from latent infection , removing the infection by immune responses at the rate of . The fourth equation explores the rate of change in the free virions; the rate of proliferation from the infection is given by and the decay rate is . Finally, the production of effector cells can have the maximum proliferation in an infected cell with a maximum rate and is given by the term . The loss of functional effector cells is defined by , with an assumption .

2.1. Positivity and Boundedness

In order to ensure the convergence of the solutions of the model, it becomes necessary to ensure that solutions of the state variables in model (1) are positive and ultimately bounded.

Theorem 1. Assume that , , , , ) be a solution of the proposed model (1) with the initial conditions , , , , . It is positive and ultimately bounded for .

Proof. Consider the proof regarding for all . In this regard, assume that there exists , which implies , , such that . From (1), it is clear that , , which is a contradiction and leads to the proof that .
Similarly, the proof can be extended for all the remaining state equations.Now, to prove {, , , , . And consider that and defineIf , for and , , for  = 0, then we have , which is a contradiction.
If , for and {, , for with . However, from (1), one can have .
Similarly, for and is also a contradiction. Thus, , . To proceed with ensuring the boundedness of the solutions, we extend the results of the positivity of the solution for model (1).Taking the limits will lead toLet Define that  =  andTo prove the ultimate boundedness of free virions, the same approach is followed:Since cannot be negative, if , then . Similarly, for , we getApplying integration on both sideswhereHence, it is proved that all the state variables solutions are ultimately bounded.
Considerwhere , , and are given in the above equation. From the given theorem, it is clear that, within the region, is a positive invariant.

3. Equilibria

This section describes the derivation of equilibria of PTC HIV model based on three cases such as disease-free, immune-free equilibrium, and endemic equilibrium.(i)Model (1) without infection exhibits disease-free equilibrium , where .(ii)If , model (1) has an immune-free equilibrium with the coefficients(iii)The endemic equilibrium of the system , , , , is given bywith and .

3.1. Basic Reproduction Number

In epidemiology, a basic reproduction number is the average number of persons an affected person can transmit the secondary infection. The basic reproduction number is an indicator that helps to determine the community spread of infection, which can be calculated using the next generation matrix. For the proposed model, the basic reproduction number, say, , is determined for various situations and described in the following sections.

3.1.1. For Disease-Free Equilibrium

The basic reproduction number of the model without a viral latent reservoir is

Since , then

3.1.2. For Immune-Free Equilibrium

The basic reproduction number of the model with immune-free equilibrium is

Then, and help us to find the next generation matrix as calculated:

Now, the characteristic equation for the above matrix iswhere , then the basic reproduction number is given by

3.1.3. For Endemic Equilibrium

The basic reproduction number of the model with endemic equilibrium is given by

Then, and help us to find the next generation matrix as calculated below:

Now, the characteristic equation for the above matrix is

The basic reproduction number is calculated from the next generation matrix is given bywhere , , and are described above. The spectral radius or the largest eigenvalue of the next generation matrix is called the basic reproduction number.

4. Stability Analysis

Let be any arbitrary positive equilibrium of the system (1). Then, the Jacobian matrix was evaluated at positive equilibrium’s in the view of biological aspects. This leads us to the following characteristic polynomial. The plus signs are ignored for the convenience of calculations.where and

4.1. Stability Analysis of the Model without Time Delay

The characteristic polynomial without delay is given by

Routh–Hurwitz criterion: define root-based Routh–Hurwitz matrix as follows:where , if . When , the matrix is simplified into

Suppose all the roots of the characteristic polynomial with negative real part, then the determinant of Routh–Hurwitz matrices are positive and vice versa. That is, . This can be employed to verify the proof of Theorem 1. The necessary and sufficient condition to exist for the negative real part for equation (6) is , and .

The expansion of the coefficients is given in Appendix.

4.2. Stability Analysis of the Model with Time Delay

The characteristic polynomial of the Jacobian is

The coefficients , , and , , are given in Appendix. We rewrite the above equation as

Suppose some of the eigenvalues are purely imaginary, that is, , then the characteristic equation becomes

Equating real and imaginary parts in the above equation:

Squaring and adding equations (32) and (33), we get

Simplifying the above equations leads to the following:

Since this is the differential equation of order 10, we get almost “10” roots such as .

Eliminating from equations (32) and (33), we getwhere , and .

Hence, for , is asymptotically stable by Routh–Hurwitz criterion [1], remains stable for ; we choosewhich completes the proof.

5. Hopf Bifurcation Analysis

In general, any physical system will reflect the changes in the qualitative behavior subject to states and parameters changes; for instance, changes in the rate of infection will reflect in the cell populations. Hence, determining the significant parameters that have an ability to affect the stability of the systems is considered bifurcation parameters. By choosing bifurcation parameters, various kinds of solution nature can be realized; among that Hopf-type bifurcation explores the point where solution trajectories cross the origin, say, from negative to positive, in which the system has purely imaginary eigenvalues. The present model possesses the Hopf-type bifurcation while the bifurcation parameter exceeds the threshold value. The process of deriving the stability conditions and proving the existence of Hopf bifurcation are given below.

Theorem 2. Consider that the characteristic polynomial is of the formwhere and are continuously differentiable with respect to . One of the roots is , where is continuously differentiable with respect to , and satisfies and for a positive real number . Denotewhich results in

Proof. Consider equation (38), and by taking the derivative of with respect to , will lead toThen,Since , we haveNow, we turn to the left side of (38), calculating the derivative of both sides of with respect to , we obtainThus,Since , we haveTherefore,The proof is completed.

6. Optimal Control Design

Considering (1), based on two variables such as RTI and PI as controls, namely, and . Here, represents the drug reverse transcriptase, and it safeguards the healthy CD4+ from infection, so that the healthy immune cells are maintained in the right proportion. Also, is represented the drug protease inhibitors, which maintains the release of the free virions to burst which are active and fully infected. In general, the treatment is initiated with an antiretroviral drug or the combination of two or more drugs, which creates some side effects while using a regular basis. Since the decision of choosing the drug combination is complex and, in this regard, an optimal strategy is a useful tool to understand the situation and make decisions. Now, optimal control provides different options with respect to estimating the costs of the drugs used for the therapies and observing the drug’s effectiveness for the disease. The therapy should last long as the optimal control helps us to find the suitable drug combination for the disease. Formulating the optimization problem based on optimal pair, existence is discussed in the following sections.

6.1. The Optimization Problem

In order to state the optimization problem, we first consider and vary with time.

The optimization problem is designed in terms of maximizing the following objective function constructed from the model parameters:

Here, the upper bound denotes the period of the treatment and assumptions and , respectively, stand for benefit and treatment costs. The scalars and are bounded and Lebesgue integrable. The objective of the control is to increase the uninfected cell populations through immune cells and decrease the cell count of free virions and infected cells. Hence, is the control pair that needs to be investigated. We assume that the control pair is nonempty, convex, and closed and it is integrable in the objective functional.where U is the control set defined by measurable, .

6.2. Optimal System

In order to investigate the properties of the optimal system, the Pontryagin’s minimum principle given in [13] is utilized and it provides necessary stability conditions for the designed optimal control problem. The advantage of the principle is that it handles (47)–(49) in terms of maximizing a Hamiltonian through and with

Based on the above discussion, the following theorem can be derived.

Theorem 3. For any optimal control , , and solutions of the corresponding state system (1), there exists adjoint variables , , , , and satisfying the equationswith the transversality conditionsMoreover, the optimal control is given by

Proof. The adjoint equations and transversality conditions can be obtained by using Pontryagin’s minimum principle with delay, such that suppose if consider the equation in the vector for .The optimal control and can be solved from the optimality conditionsThat is,Consider the optimal controls and defined in (54) are bounded. If we substitute and in systems (47) and (52), we obtain the following optimality system:The significance of stability conditions, validating the existence of Hopf bifurcation properties, and effect of the optimal controller are explored in the following section.

7. Numerical Simulation

This section describes the numerical evaluations of the proposed model (1) by choosing the experimental range of parameter values provided in Table 1. The simulations are performed through Runge–Kutta fourth-order numerical approximation scheme. The outcomes of the cell simulations corresponding to cell populations are demonstrated through phase-space diagrams. Figures 1 and 2 illustrate the solution behavior of cell populations within the threshold and exceeding the threshold of time delay, respectively. Also, Figure 3 explores the effect of optimal control design, in which the objective is achieved by increasing the number of uninfected cell populations and diminishing the number of virions. In addition, Figures 47 provide insights by comparing the controlled and uncontrolled system states such as uninfected, infected, latent, and free virions, respectively. Figures 2 shows the interpretation of variables of the equilibrium when . The top left corner in the panel of Figure 2 evidence the rate of change in healthy CD4+ and infected cells and the bottom left depicts the uninfected cells and the free virions. Similarly, the phase portrait of all the cell populations such as uninfected and infected cells with free virions and immune cells is accordingly picturized. The drug parameter and the transmission rate have a high influence on the stability of the equilibrium. That is, for , the equilibrium is stable. For , the equilibrium remains stable. When the transmission rate gets higher or lesser, it loses the stability of equilibrium, where , , , and also affect the stability when it increases rapidly. The rest of the parameters are having less significance in terms of stability when compared with the remaining parameters.

Figure 1 represents the state trajectories where the equilibrium point loses its stability due to time delay that makes the transmission rate slow. If the drug is already administered into the body, the drug starts to function among the cells; there is a period for the drug to make progress and with respect to time delays, it becomes necessary to have the effect of time response, and if the drug efficacy is not up to expectation, then the endemic equilibrium will loose the stability.

Figure 8 shows that the uninfected cells start to bifurcate as crosses the value 0.93.

Figures 3 represents the uninfected cell population after the optimal control is established. The top figures show the increasing cell population of the uninfected cells and the reducing rate of the infected cells; it is visible after a certain stage the immune cells maintained in proportion, with no increasing number of cells. Similarly, the free virions are depleted in the process.

Figures 47 provide the comparisons between controlled and uncontrolled susceptible cell populations, infected cell populations, latent cell populations, and free virions. The graphs exhibit the controlled cell population after the implementation of optimal control.

8. Conclusion

The aim of the paper is to explore the dynamical analysis of posttreatment HIV-1 infection with respect to various significant parameters such as the effect of time delays, two different kinds of antiretroviral therapies, and loss of functional effector cells. By employing the Routh–Hurwitz criterion, the stability properties of the model with respect to discrete-type constant time-delay which is chosen as a bifurcation parameter have been presented. To ensure the effect of time delays, the existence of Hopf-type bifurcation in the behavior of solutions has been proved through proving the corresponding transversality conditions. To reveal the effect of combination of drug therapies, along with those parameters, the situation is modeled as an optimal control problem in terms of the objective function. Through objective function, the results explore the maximization in the number of uninfected cell population and minimization in the number of infected cells and corresponding results are picturized in the numerical section. In the future direction, the model can be extended in terms of fractional order differential operator, considering continuous-type time-delays, stochastic disturbances, and impulse in the antiviral therapy, which also plays a significant role in the system dynamics.

Appendix

The coefficient of system variables is provided in the following:

Data Availability

All the data are available within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest.