Abstract

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.

1. 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) [26]. Although liver disease can occur solely due to HIV infection [710], 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 [1214]. Drug-induced 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 [2123] as well as macrophages, [11, 19], this study used a mathematical model 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].

2. Material and Methods

2.1. Mathematical Model Development

Recent study by [26] as well as [28] asserts that all classes of ART exhibit dose-response 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 dose-response.

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   and , where , . If the drug is 100% effective, then , where is the drug dose concentration, is the drug concentration that leads to 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 [1214]. 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 depending on the drug dose, [26], where is defined as is the drug dose, is the dose at which toxicity occurs in 50% of exposed cases, and is the gradient of the dose-response 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.

In model formulation, we define eight variables as follows: uninfected CD4+ cells (), infectious CD4+ cells (), uninfected hepatocytes (), latently infected hepatocytes (), [9], productively infected hepatocytes (), HIV-specific cytotoxic T lymphocytes (), viral load (), and the level of enzyme alanine aminotransferase in blood ().

Model parameters are as follows: CD4+ cells and hepatocytes are produced from within the body at rates and and die naturally at rates and , 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 and a probability that it infects a CD4+ T cell at rate . Infected CD4+cells die at rate , where , and are cleared by HIV-specific cytotoxic T-lymphocytes (CTLs) at a rate .

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 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, and , where [36]. Productively infected hepatocytes are killed by HIV-specific CTLs at rate 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 and are cleared at rate . Let and 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 [39]. Virions die at a per capita rate .

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], is the rate at which ALT is generated from hepatocytes that die due to HIV infection. and 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 .

From the assumptions and description above we have the following system of ordinary differential equations:

2.2. Model Analysis
2.2.1. Positivity and Boundedness of Solutions

With initial conditions , , , , , , , and , the solutions for , , , , , , , and , respectively, remain positive and bounded provided .

We define and as the total number of CD4+ cells and hepatocytes, respectively, where Given that and 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)consider the upper bound of decreases and approaches asymptotically the value .In either case remains bounded.

Similarly, it can also be shown that is bounded by .

If the total number of CD4+ cells and hepatocytes are bounded by and , then the subclasses and ( and ) respectively are bounded. Thus, letting and then , , , , and .

Using (9) and considering and , then Letting , where since , . Using the same method as in (15) where is initial viral load.

Using the same argument as in (15), we can deduce that is bounded by .

Using (8) and taking only infection-dependent proliferation, on assumption that antigen-independent proliferation is naturally bounded, it can be shown that where . Letting , then .

Equation (10) is defined by bounded functions; therefore, is bounded by some , where .

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 , it can be deduced that the feasible solution for the model (3)–(10) is

2.2.2. 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, , the basic reproductive number calculated using the next generation method as in [47] is given by where is the number of secondary infections in latently infected hepatocytes compartment. is the number of secondary infections in productively infected hepatocytes. is the number of secondary infections produced by one virus producing CD4+ cell. 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, is dependent on antigen-independent CTLs proliferation rate () and independent of antigen-dependent proliferation rate (). 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 and unstable when .

2.2.3. Effective Reproductive Number

Analysing the system of (3)–(10) with the medication (), 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 (). is defined as where , , and are defined the same way as , , and , 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 (), then there would be no secondary infections in either type of cell. However, if the drug is totally ineffective (), 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 , 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 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 , the critical toxic effect which is a function of drug efficacy is then calculated and satisfies the equation where It can be shown that and and in this case we have two possibilities for a positive .

Case 1. When , that is, then

Case 2. When , that is, then Comparing the two scenarios, we take 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 and as . If then the toxic effect will slow down the progress of the infection; otherwise it will speed it up: where .

Since , the right hand side of (36) is greater than the left hand side provided . This is not likely because of the values of and (100 and 1000, resp.). Therefore ; this would imply that , hence indicating that the toxic effect of medication speeds up the progression of infection.

Differentiating with respect to where From (37), if and , that is, then According to [49], in order to slow down the infection rate of the virus when treatment is implemented, the conditions   and 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 and , we considered the effective reproductive number of both CD4+ cells and hepatocytes. Starting with reverse transcriptase inhibitors (RTIs), where Since , then 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 . We can therefore conclude that is a decreasing function of and 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 .

2.2.5. 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 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 are given by the following equations: 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 , .

2.2.6. Reduced Model

We reduce the model for further analysis with assumption that RTIs and PIs have an overall drug efficacy given by , where [52]. The combined efficacy reduces the rate at which the virus infects CD4+ cells () and hepatocytes (). 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: Hepatocytes subpopulation: The corresponding endemic equilibrium values for the variables in CD4+ cells subpopulation are in which where It is clear that and . 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 . Therefore, (58) has one positive root given by To find the critical combined efficacy of antiretroviral therapy, we assume this happens when . Since all parameters are nonnegative, then can only be zero if . Thus, the critical efficacy of ART in CD4+ cells above which the infection would be managed would be given by If 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 ; 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 of medication that would not lead to elevated alanine aminotransferase in the blood even if HIV remains endemic in the liver. We therefore compute as Like our previous findings, system (65) also shows that the toxicity is a function of efficacy. Analytical solution of from (65) could not be easily obtained since it is also a function of . We therefore investigated for when . When then . Since and 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, 1315].

3. Results

3.1. 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 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 [26], 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.

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 NNRT-based regimen, the least toxic combination was 3TC + DDI + EFV while AZT + d4T + EFV was the most toxic. In PI-based 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].

4. Discussion

Much as liver injury can occur solely due to HIV infection [710], liver related mortality and morbidity are being associated with the use of antiretroviral therapy, [26], because all antiretroviral therapy has some degree of toxicity, [15]. This study used a mathematical model with therapy efficacy as 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], 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 NNRTI-based regimens [24, 25, 58].

5. 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.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to acknowledge the Sida/SAREC bilateral research cooperation programme of Makerere University for funding this research. They also would like to thank the Center for International Mobility (CIMO), Finland, for funding research visits to Lappeenranta University of Technology.