Modelling Hepatotoxicity of Antiretroviral Therapy in the Liver during HIV Monoinfection

Liver related complications are currently the leading cause of morbidity and mortality among human immunodeficiency virus (HIV) infected individuals. In HIV monoinfected individuals on therapy, liver injury has been associated with the use of antiretroviral agents as most of them exhibit some degree of toxicity. In this study we proposed a mathematical model with the aim of investigating hepatotoxicity of combinational therapy of antiretroviral drugs. Therapy efficacy and toxicity were incorporated in the model as dose-response functions. With the parameter values used in the study, protease inhibitors-based regimens were found to be more toxic than nonnucleoside reverse transcriptase inhibitors-based regimens. In both regimens, the combination of stavudine and zidovudine was the most toxic baseline nucleoside reverse transcriptase inhibitors followed by didanosine with stavudine. However, the least toxic combinations were zidovudine and lamivudine followed by didanosine and lamivudine. The study proposed that, under the same second line regimens, the most toxic first line combination gives the highest viral load and vice versa.


Introduction
Among people infected with human immunodeficiency virus (HIV), liver related complications have become the leading causes of morbidity and mortality [1].There are a number of factors identified that are known to contribute to liver related mortality in HIV infected people.These are, for instance, coinfections with viral hepatitis B (HBV) and hepatitis C (HCV) as well as the use of antiretroviral therapy (ART) [2][3][4][5][6].Although liver disease can occur solely due to HIV infection [7][8][9][10], the use of ART is highly associated with end stage liver disease in HIV infected people [6,11].
A human body identifies every substance it absorbs, including drugs as a foreign substance.All substances that are identified as foreign are subjected to chemical processes, [12], to make them suitable for elimination.Although all tissues in the body have some ability to metabolise, the liver is the central metabolic clearing organ of all chemicals from a human body.This process is carried out by metabolising enzymes cytochrome P450 which are found in hepatocytes [12].The role played by hepatocytes in handling toxic drug substances makes them vulnerable to drug induced injury (hepatotoxicity) and, consequently, cell death [12][13][14].Druginduced liver injury is responsible for 5% of all hospital admissions and 50% of all acute liver failures [5].
According to [15], generally all antiretroviral drugs exhibit some level of toxicity to the liver.Hepatocytes that are responsible for metabolizing toxic antiretroviral drugs [12,13] also support all stages of HIV infection and replication [9,10,16,17].This double distress to hepatocytes would possibly partly explain the progression of liver disease in HIV infected people on ART.Using ART in combination makes it even more complicated to assess which drug leads to higher levels of hepatotoxicity than the other [18].
According to [8,10,19], HIV can directly infect human hepatocytes.Considering HIV infection and replication in hepatocytes [16,17,20] and CD4+ cells [21][22][23] as well as macrophages, [11,19], this study used a mathematical model 2 International Journal of Computational Mathematics  [35] to analyse therapeutic and toxic effects of HIV antiretroviral therapy.Defining hepatotoxicity generally as liver damage, there are a number of conditions that would fall under hepatotoxicity, and they include hepatitis, hepatic necrosis, and hepatic steatosis [24,25].This study considered only hepatic necrosis, which is the death of liver cells.Drugs studied include zidovudine (AZT), emtricitabine (FTC), didanosine (DDI), lamivudine (3TC), stavudine (d4T), efavirenz (EFV), atazanavir (ATV), and nelfinavir (NFV).Various drug combinations as recommended by World Health Organisation were considered.The study accounted for therapeutic and toxic effects of all classes of ART in form of dose-response functions [26,27].

Mathematical Model Development.
Recent study by [26] as well as [28] asserts that all classes of ART exhibit doseresponse curves.Finding a dose that gives 50% of maximal response is one method of determining how effective the drug is; therefore, [26] recommends that it would be of great contribution if efficacy of ART would be modelled as doseresponse.
Much as antiretroviral drugs are taken in doses at specific time intervals and the effectiveness and toxicity of the drug are largely dependent on the pharmacokinetics of the drug taken, our study preferred to use the Hill equation as recommended by [26], to model drug efficacy as opposed to Michaelis-Menten kinetics [29].Reactions involving a single substrate are often assumed to follow Michaelis-Menten kinetics, irrespective of the model's underlying assumptions.Since we are dealing with possible multiple drug reactions, the use of Michaelis-Menten equation would thus be inappropriate.The study has also adopted the use of the Hill equation because we are investigating therapeutic and toxic effect of therapy at steady state.
Therapeutic response function is defined as a Hill equation (1) to describe the effectiveness of the drug [26].Effectiveness of reverse transcriptase inhibitors (RTIs) and protease inhibitors (PIs) are represented, respectively, by drug efficacies  1 and  2 , where where  is the drug dose concentration, IC 50 is the drug concentration that leads to 50% of the maximal viral inhibition, and  is the gradient of the dose-response curve.Response is the ability of the drug to inhibit viral replication [28].
Gradients for dose-response curves of HIV drugs are given by [30] and are shown in Table 1.
In the process of metabolising toxic ART, some healthy as well as HIV infected hepatocytes get injured/die [12][13][14].The toxic effect of ART on hepatocytes is also assumed to be dose-dependent.Thus, infected and healthy hepatocytes are assumed to die during drug metabolism at a rate  for 0 ≤  ≤ 1 depending on the drug dose, [26], where  is defined as is the drug dose, TD 50 is the dose at which toxicity occurs in 50% of exposed cases, and  is the gradient of the doseresponse curve.We assume that both toxic and therapeutic effects of antiretroviral drugs exhibit dose-response curves of relative gradients since both are dose-dependent [31].Due to scarcity of literature regarding dose-response analysis of ART toxicity, we assume that all the gradients might not be the same but relative.This is based on the assumption that since nonnucleoside reverse transcriptase inhibitors (NNRTIs) and PIs are the most efficacious and have higher gradients than NRTIs, [28,30], this is consistent with toxicity and hence relative gradients.
Model parameters are as follows: CD4+ cells and hepatocytes are produced from within the body at rates  1 and  2 and die naturally at rates  1 and  3 , respectively.We assume that, among all cells in the liver, HIV has high affinity for CD4+ cells and hepatocytes [20].Thus, from the free viral population in the liver, if a virus is to infect a cell, there is a probability  that it will infect a hepatocyte at rate  2 and a probability 1− that it infects a CD4+ T cell at rate  1 .Infected CD4+cells die at rate  2 , where  2 >  1 , and are cleared by HIV-specific cytotoxic T-lymphocytes (CTLs) at a rate  1 .
When a hepatocyte is infected, there is a probability  that it becomes productively infected (viral replication will take place after successful reverse transcription) and probability (1 − ) that the cell will become latently infected, such that there is no viral production until cell activation.Decay rates for productive hepatocytes and latently infected hepatocytes are, respectively,  5 and  4 , where  5 >  4 [36].Productively infected hepatocytes are killed by HIVspecific CTLs at rate  2 and until activated, latently infected hepatocytes will not trigger the action of CTLs.This study assumes that latently infected hepatocytes will either get activated to become infectious or die.There is no possibility that they will become uninfected again [37].Latently infected hepatocytes are activated at rate , and it is assumed that this rate is reduced by the the efficacy of PIs.This is because our study assumed that reverse transcription has already taken place at this stage so reverse transcriptase inhibitors will have no effect.
With or without any pathogen in the body CTLs proliferate naturally at rate  and in the presence of HIV infection, they proliferate at rate  3 and are cleared at rate  6 .Let  1 and  2 represent the rates of HIV production per infected CD4+ cells and productively infected hepatocytes, respectively.In addition to CD4+ cells and hepatocytes, HIV productively infects other cells and macrophages like Kupffer cells in the liver [11,19,45].These cells produce virions at rate .The study assumed that the effect of medication is translated generally into minimal viral load.Thus, viral production from macrophages is also inhibited by both RTIs and PIs.As in [46], we assumed a synergy additivity of PIs and RTIs as Φ = (1 −  1 )(1 −  2 ) [39].Virions die at a per capita rate  7 .
This study assumes that some infectious hepatocytes that die due to drug metabolism are able to release viral particles, provided that at the point of hepatocyte's death all the stages of viral production have been attained. is the per capita rate of virions production by each infectious hepatocyte that dies due to drug metabolism.Viral production due to hepatotoxicity was inhibited by resultant efficacy of both RTIs and PIs.
Equation (10) defines the level of enzyme alanine aminotransferase (ALT) in the blood system.Among other enzymes, hepatocytes contain enzyme ALT, and when the cells die by any means, ALT leaks into the blood where it is clinically detected.According to [3], when there is no infection, the level of ALT in the blood system is generated from naturally dying hepatocytes at rate .As described in [42],  4 is the rate at which ALT is generated from hepatocytes that die due to HIV infection.( 5 +  +  2 ) and ( 4 + ) are the total death rate of productively infected hepatocytes and latently infected hepatocytes, respectively, attributed to HIV infection.The contribution to ALT by healthy hepatocytes is only due to drug metabolism.ALT is cleared from the blood naturally at rate  8 .
From the assumptions and description above we have the following system of ordinary differential equations: Given that  2 >  1 and  5 >  4 >  3 as in [36] and that all variables and parameters are positive (otherwise the model would not be biologically feasible) then Thus Implying hence where  is the total number of hepatocytes at the time of infection.Equation (15) indicates that there are two cases to consider: (1) consider the upper bound of  ℎ increases and approaches asymptotically the value  2 /( +  3 ); (2) consider the upper bound of  ℎ decreases and approaches asymptotically the value  2 /( +  3 ).
Similarly, it can also be shown that   () is bounded by If the total number of CD4+ cells (  ) and hepatocytes ( ℎ ) are bounded by  1 / 1 and  2 /( 3 + ), then the subclasses (  ) and ( ha and  hl ) respectively are bounded.Thus, letting Using ( 9) and considering   <  2 and   <  1 , then Using the same method as in (15) where  0 is initial viral load.
Using the same argument as in (15), we can deduce that  is bounded by  0 .
Using (8) and taking only infection-dependent proliferation, on assumption that antigen-independent proliferation is naturally bounded, it can be shown that where Equation ( 10) is defined by bounded functions; therefore,  is bounded by some  5 , where  5 = max{ 1 ,  3 ,  4 }.
Suppose the feasible solution of CD4+ cells in the liver lies in the region Θ  ; then, Similarly, assuming the feasible solution of hepatocytes in the liver to lie in the region Θ ℎ , then Letting  3 = / 6 , it can be deduced that the feasible solution for the model ( 3)-( 10) is

The Basic Reproductive Number.
When there is no HIV infection in the liver, the system of ( 3)-( 10) settles to a disease-free equilibrium   defined by Taking the system of ( 3)-( 10) without medication, that is,  1 =  2 =  = 0, the basic reproductive number calculated using the next generation method as in [47] is given by where 0hl is the number of secondary infections in latently infected hepatocytes compartment. 0ha is the number of secondary infections in productively infected hepatocytes. 0 is the number of secondary infections produced by one virus producing CD4+ cell. 0 is the total number of secondary infections in the liver.The total number of secondary infections is directly proportional to the clearance rate of CTLs and inversely proportional to the clearance rate of virions.Generally,  0 is dependent on antigen-independent CTLs proliferation rate () and independent of antigen-dependent proliferation rate ( 3 ).This indicates that if the CTLs are boosted prior to infection, then the body can handle infection better than when they (CTLs) proliferate in the presence of infection.High CTL clearance rate implies higher numbers of secondary infection.Probably, if many CTLs are cleared from the body by any means, then there are fewer cells left to fight the infection hence faster progression of the infection.Using Theorem 2 of [47] we establish the following result that the disease-free equilibrium   is locally asymptotically stable when  0 < 1 and unstable when  0 > 1.

Effective Reproductive Number.
Analysing the system of ( 3)- (10) with the medication ( 1 =  2 =  ̸ = 0), we calculated the effective reproductive number (  ) using the next generation method as in [47].  is the actual average number of secondary cases per primary case that reflects the impact of therapy on infection.Technically, the number of secondary infections during therapy (  ) should be less than those without therapy ( 0 ).  is defined as where ,  ha , and  hl are defined the same way as  0 ,  0ha , and  0hl , respectively, in (25) except that the latter are functions of therapy efficacy and toxicity.Secondary infections in either type of cells largely depend on the drug efficacy.It can be seen that if the drug is 100% effective ( 1 =  2 = 1), then there would be no secondary infections in either type of cell.However, if the drug is totally ineffective ( 1 =  2 = 0), it would be expected that the number of secondary infections would be equivalent to those produced when no therapy is administered.This is however not the case, because whether the drug is effective or not, its toxicity will affect the cells involved in drug metabolism.At present there is no HIV antiretroviral drug that is 100% effective, [48]; likewise, there is no HIV therapy that is free of toxicity [18].When  1 =  2 = 0, it is not straightforward to deduce whether the number of secondary infections in hepatocytes during therapy is less or greater than the number of secondary infections when therapy is not administered.We therefore assess the toxic function that would result in the same number of secondary infections as those when no therapy is administered.We suppose that this is the critical toxic effect  crit of the drug that should not be exceeded if the therapy should reduce the number of secondary infections in the liver.Since toxicity is a function of drug dose, the dose that corresponds to the critical toxic effect can consequently be computed.Using  0hl +  0ha =  hl +  ha , the critical toxic effect which is a function of drug efficacy is then calculated and satisfies the equation where It can be shown that  0 > 0 and  2 < 0 and in this case we have two possibilities for a positive  crit .
Case 1.When  1 > 0, that is, then Case 2. When  1 < 0, that is, then Comparing the two scenarios, we take max{ crit 1 ,  crit 2 } and hence deduce the critical drug toxicity that should not be exceeded by HIV medication if it should help in reducing the number of secondary infections, given the appropriate choice of parameters, as 2.2.4.Sensitivity Analysis.We carried out sensitivity analysis of the toxic function in the effective reproductive number of the hepatocytes.The study ignored CD4+ cells component of the basic reproductive number in toxic analysis, on assumption that toxic effect of the medication does not have direct effect on CD4+ cells.Following [49], we determined the effect of  on   first by calculating the difference between  0ℎ and  ℎ as Δ.If Δ > 0 then the toxic effect will slow down the progress of the infection; otherwise it will speed it up: where Since 0 ≤  ≤ 1, the right hand side of ( 36) is greater than the left hand side provided (Φ +  2 (1 −  2 )) < 1.This is not likely because of the values of  and  2 (100 and 1000, resp.).Therefore Δ < 0; this would imply that  0ℎ <  ℎ , hence indicating that the toxic effect of medication speeds up the progression of infection.
Differentiating  ℎ with respect to   ℎ  = (( + ) ( +  6 ) ( − ) where From (37), if ( − ) > 0 and ( −  6 ) > 0, that is, then According to [49], in order to slow down the infection rate of the virus when treatment is implemented, the conditions Δ > 0 and  ℎ / < 0 should be satisfied.Based on the condition derived from (36), the effective reproductive number of hepatocytes is an increasing function of the toxic effect of therapy .This implies that increase in toxicity of medication leads to increase in secondary infections in the liver, provided the relationship between drug efficacies as stated in conditions (39) is fulfilled.
Carrying out sensitivity analysis on the drug efficacies  1 and  2 , we considered the effective reproductive number   of both CD4+ cells and hepatocytes.Starting with reverse transcriptase inhibitors (RTIs), where Since , ,  > 0, then   / 1 < 0 indicating that   is a decreasing function of RTIs.Thus, increase in efficacy of RTIs implies a reduction in secondary infection.
With protease inhibitors (PIs), we have where It can be shown that   / 2 < 0. We can therefore conclude that   is a decreasing function of  1 and  2 as shown in Figure 1.Increasing the drug efficacy will result in a decrease in secondary infections.Referring to (1), increasing efficacy would necessitate increasing drug dose concentration so that it is much greater than its IC 50 .

Endemic Equilibrium.
Due to the complexity of the model, mathematical analysis is rather cumbersome.We therefore assume a free virus spread of infection and no  2.
cell-to-cell transfer of HIV, [50], ruling out the interaction between CD4+ cells and hepatocytes.We therefore split the model into two subpopulations of either cell.We further consider virions produced from either type of cells as the viral production in that subpopulation.We assume a fraction  of the total viral population from macrophages to contribute to the viral population in hepatocytes subpopulation while the remainder would contribute to the viral population in CD4+ cells subpopulation.We subdivide the system of ( 3)-( 10) into the following.CD4+ cells subpopulation infection dynamics: Hepatocytes subpopulation infection dynamics: We derive the endemic equilibrium state as where and  *  and  * ha are given by the following equations: International Journal of Computational Mathematics in which Differentiating (49) with respect to   to obtain the steady states values of   , it can be shown that given then there would be two equilibrium values of   .Otherwise CD4+ cells would have one positive endemic equilibrium point given by provided the relationship between the two drug efficacies was given by With hepatocytes subpopulation, however, there was no finite number of equilibrium states but rather a rage of values.Therefore, using Descartes rule of signs [51], ( 50) had at most one positive endemic equilibrium everywhere except in the region  satisfied by  = { 1 < 0,  2 > 0,  3 < 0}.

Reduced Model.
We reduce the model for further analysis with assumption that RTIs and PIs have an overall drug efficacy given by Φ, where Φ = 1 − (1 −  1 )(1 −  2 ) [52].The combined efficacy Φ reduces the rate at which the virus infects CD4+ cells ( 1 ) and hepatocytes ( 2 ).We further assumed that much as macrophages act as reserve source for HIV production throughout HIV infection [53], the level of replication is low as compared to other cells.
We therefore revised the model to exclude viral production from macrophages and subdivided it into CD4+ cells and hepatocytes subpopulation.The objective was to investigate the critical combined efficacy of antiretroviral therapy above which the infection would be managed as well as the optimal toxicity below which alanine aminotransferase level in the blood would be negligible.CD4+ cells subpopulation: International Journal of Computational Mathematics 9 Hepatocytes subpopulation: The corresponding endemic equilibrium values for the variables in CD4+ cells subpopulation are in which where It is clear that  2 > 0 and  1 > 0. Comparing with ( 28), (considering ) and assuming the effective reproductive number of CD4+ cells to be above unity in HIV infection during therapy, it can be shown that  0 < 0. Therefore, (58) has one positive root given by To find the critical combined efficacy of antiretroviral therapy, we assume this happens when  + = 0. Since all parameters are nonnegative, then  + can only be zero if  +  = 0. Thus, the critical efficacy of ART in CD4+ cells above which the infection would be managed would be given by If Φ > Φ critical then the medication will be able to keep the virus in a controlled state even if the infection remains endemic; that is, the disease-free state and the endemic steady state coexist [54].
In the same way we investigated the optimal toxicity of ART that would not elevate liver enzymes in the blood system even if HIV infection remained endemic.This was assumed to happen when  + = 0; we thus established the endemic equilibrium point in the hepatocyte subpopulation as in which where We previously analysed the critical toxicity that would result in the same level of ALT in the blood as when no therapy is administered and found that this was a function of therapy efficacy.Our aim now is to identify the optimal toxic value  opt of medication that would not lead to elevated alanine aminotransferase in the blood even if HIV remains endemic in the liver.We therefore compute  opt as Like our previous findings, system (65) also shows that the toxicity is a function of efficacy.Analytical solution of  opt from (65) could not be easily obtained since it is also a function of  + ha .We therefore investigated for  + = 0 when  + ha = 0.When  + ha = 0 then  6  opt ( 4 + )( 4 +  +  opt ) = 0. Since ( 4 +  +  opt ) ̸ = 0 and ( 4 + ) ̸ = 0 then the only possibility of having no elevated alanine aminotransferase in the blood when HIV is endemic in the liver and the infectious hepatocytes have been reduced to zero is when the medication is not toxic at all.This concurs with various researches that report high levels of liver enzyme in the blood due to toxic nature of all classes of antiretroviral drugs [11,[13][14][15].

Numerical Simulation.
In this section we present the numerical simulations of the model while both the therapeutic and toxic effects of the drugs are incorporated.Medications used in this study are listed in Table 1, while parameter values used in simulations are in Table 2. Parameter values used in calculation of therapeutic and toxic functions in ( 1) and ( 2) are shown in Table 1.
Simulating the dynamics of HIV in the liver, three cases were considered: (1) no treatment, (2) with treatment but without toxic effect of the drugs, and (3) with treatment plus toxic effect.Figure 2 depicts the results.Clearly, therapy reduced viral load resulting in fewer number of latent and infectious cells.In case of CD4+ cells, the medication consequently resulted into increased number of uninfected cells.Healthy CD4+ cells were many when the toxic therapy was administered as opposed to when it was not used.However, the dynamics were different in hepatocytes.When the toxic medication was administered, the number of uninfected cells was much less than when medication was not applied at all.This suggested that much as HIV infects hepatocytes, the rate of progression is slower in these types of cells as compared to CD4+ cells [20].However, hepatocytes seemed to be more affected by toxicity than HIV infection.This was consequently seen in the level of enzymes in the blood that was a lot higher in the case of drug toxicity included than when it was not considered.However, these findings are subject to parameter values and our results are based on parameters in Tables 1 and 2.
The World Health Organisation recommends that antiretroviral therapy be used in combination of NRTIs together with NNRTIs or PIs [55].This study considered three drugs from NRTIs, one drug from NNRTIs and two drugs from PIs.All combinations from the sample drugs were studied.It is further recommended that atazanavir (ATV) should be used with another PI in a PI-based regimen [55], so it was combined with nelfinavir (NFV).Therapeutic effect of the drugs was modeled using a Hill function as shown in (1).In NNRTI-based regimen, Figure 3 shows that 3TC + DDI + EFV gave the minimal viral load while AZT + d4T + EFV gave the maximal viral load.In PI-based combinations, Figure 4 shows that DDI + 3TC + ATV + NFV gave the least viral load and AZT + d4T + ATV + NFV gave the most.This would suggest that, within the parameter values in Table 2, among   2.
the baseline NRTIs studied, DDI + 3TC is the most efficacious and AZT + d4T is the least.We compared Figures 2, 3, and 4 and established that the level of ALT in the blood when no medication was used was higher than the level of ALT when medication was used.This was contradicting with a number of researches [2][3][4][5][6], indicating that the use of ART increases the level of ALT in the blood stream.
All antiretrovirals are associated with some level of toxicity [5,56].During drug metabolism, the toxic nature of ART causes liver cells stress and consequently cell death.This study assumes that the lower the ALT level in the blood the less toxic the therapy.With toxic effect (2) incorporated in the model, simulation results are shown in Figure 5 for NNRTI-based regimen and Figure 6 for PI-based regimen.2.
Comparing Figures 5 and 6 and their corresponding dynamics in toxic-free cases in Figures 3 and 4, respectively, it is noted that the level of ALT is highest in the former.It can also be noted that PI-based regimens are more toxic than NNRTI-based regimens among the studied drugs.In NNRTbased regimen, the least toxic combination was 3TC + DDI + EFV while AZT + d4T + EFV was the most toxic.In PIbased regimen, 3TC + DDI + ATV + NFV and AZT + d4T + ATV + NFV were the least and most toxic combinations, respectively.Combinations that contained 3TC were least toxic while combining d4T with either AZT or DDI made the combination toxic.This is consistent with a number of researches that assert that 3TC is the least toxic NRTI while d4T is the most toxic [18,25,48,55,57].
well as toxicity incorporated, to study various types of HIV drugs.
The study investigated for toxicity threshold between 0 and 1 inclusive, below which HIV infection in the liver can be controlled.Mathematical analysis showed that the critical toxicity threshold was a function of drug efficacy.The dependency of toxicity on efficacy would be due to the fact that both are dose-dependent [26,28].
With parameter values as given in Table 2, without including the toxic effect of medication, in NNRTI-based regimen, the combination of DDI + 3TC that contained EFV gave the minimal viral load while the combination of AZT + d4T gave the maximum.In PI-based regimen, the combination of ATV + NFV + DDI + 3TC gave the least viral load while ATV + NFV + AZT + d4T gave the highest.
With toxicity incorporated and parameter values as in Table 2, in NNRTIs-based combinations, DDI + EFV + 3TC produced the minimal level of ALT, while including d4T + EFV + AZT gave maximal ALT level in the blood system.In PI-based regimen the most toxic combination was d4T + ATV + NFV + AZT while the least toxic was DDI + ATV + NFV + 3TC.Our findings are consistent with [24,58] 2.
who recommend that d4T should not be combined with AZT because the combination is highly toxic.PI-based regimens were found to be the most toxic and at the same time the best in reducing viral load in the liver as compared to NNRTIbased regimens [24,25,58].

Conclusion
The findings regarding toxic and therapeutic effects of NRTI drugs were consistent with the already existing literature.
We therefore deduce that, with the drugs studied and parameter values used, the most toxic combination gave the highest viral load in the liver and vice versa.It is important to note that there was no variability in NNRTIs and PIs as a single drug from NNRTIs and same combination from PIs was used.The same method could possibly be used with various drugs in NNRI-based and PI-based regimens to explore the most therapeutic and toxic combinations in HIV therapy.2.

Figure 1 :
Figure 1: Increase in effective reproductive number with increase in drug efficacy with  = 0.5.Horizontal axes represent time in days and parameter values are as shown in Table2.

Figure 2 :
Figure 2: HIV infection dynamics with and without medication and with and without drug toxic effect.Horizontal axes represent time in days and parameter values are as shown in Table2.

Figure 3 :
Figure 3: HIV infection dynamics when NNRTI-based regimen is administered without considering the toxic effect of the therapy.Horizontal axes represent time in days and parameter values are as shown in Table2.

Figure 4 :
Figure 4: HIV infection dynamics when PI-based regimen is administered without considering the toxic function of the therapy.Horizontal axes represent time in days and parameter values are as shown in Table2.

Figure 5 :
Figure 5: HIV infection dynamics when NNRTI-based regimen is administered with the toxic function of the therapy considered.Horizontal axes represent time in days and parameter values are as shown in Table2.

Figure 6 :
Figure 6: HIV infection dynamics when PI-based regimen is administered with the toxic function of the therapy considered.Horizontal axes represent time in days and parameter values are as shown in Table2.

Table 1 :
Medication and corresponding parameter values. and IC 50 are from Supplementary Table 1 of [30].

Table 2 :
Parameters used in simulations. ,