- Research article
- Open access
- Published:

# Modelling imperfect adherence to HIV induction therapy

*BMC Infectious Diseases*
**volumeÂ 10**, ArticleÂ number:Â 6 (2010)

## Abstract

### Background

Induction-maintenance therapy is a treatment regime where patients are prescribed an intense course of treatment for a short period of time (the induction phase), followed by a simplified long-term regimen (maintenance). Since induction therapy has a significantly higher chance of pill fatigue than maintenance therapy, patients might take drug holidays during this period. Without guidance, patients who choose to stop therapy will each be making individual decisions, with no scientific basis.

### Methods

We use mathematical modelling to investigate the effect of imperfect adherence during the inductive phase. We address the following research questions: 1. Can we theoretically determine the maximal length of a possible drug holiday and the minimal number of doses that must subsequently be taken while still avoiding resistance? 2. How many drug holidays can be taken during the induction phase?

### Results

For a 180 day therapeutic program, a patient can take several drug holidays, but then has to follow each drug holiday with a strict, but fairly straightforward, drug-taking regimen. Since the results are dependent upon the drug regimen, we calculated the length and number of drug holidays for all fifteen protease-sparing triple-drug cocktails that have been approved by the US Food and Drug Administration.

### Conclusions

Induction therapy with partial adherence is tolerable, but the outcome depends on the drug cocktail. Our theoretical predictions are in line with recent results from pilot studies of short-cycle treatment interruption strategies and may be useful in guiding the design of future clinical trials.

## Background

Currently, 33 million people worldwide are infected with HIV/AIDS, of whom 2.7 million were infected in 2007 [1]. HIV is a disease that is accompanied by a profound depletion in the number of CD4^{+} T cells and can be transmitted by blood or other body fluids [2]. Most patients with HIV/AIDS are prescribed a triple-drug cocktail with either three nucleoside-analogue reverse transcriptase inhibitors (RTIs), or two RTIs and one protease inhibitor (PI) [3]. However, PI-sparing cocktails have been shown to have equivalent potency to PI-containing cocktails [4] and may reduce the risk of metabolic and potential cardiovascular consequences of PI-containing therapy, while providing similar or improved virologic control and durability of effect [5].

The importance of adherence to HIV drug regimens presents challenges that arise from the biology of HIV, the magnitude of the required therapeutic effort and the changing demography of HIV infection [6]. In order to determine regimens for partial adherence, a number of mathematical models have attempted to quantify how drug concentration levels in the body of an HIV patient affect viral replication [7â€“13].

Adherence to drug therapy is necessary in order to control HIV, but sometimes-overwhelming side effects, as well as the inconvenience of following a strict regimen, deter patients from taking their drugs [14]. Imperfect or partial adherence can facilitate the emergence of drug-resistant mutations [15].

Induction therapy is a HIV/AIDS treatment regime that hopes to benefit patients by decreasing drug resistance and reducing the overall number of drugs that must be taken. In order to minimise drug resistance, induction-maintenance (IM) therapy strategies begin with a period of intensified antiretroviral therapy (induction phase), followed by a simplified, long-term regimen (maintenance phase) [16â€“19].

Previous work with induction therapy failed due to uncalculated latently infected cells and imperfect adherence [17, 19]. Recently, however, Curlin *et al*. [20] have shown that a longer induction phase decreases the probability that viruses resistant to maintenance therapy will emerge. Their studies have shown that the probability of success (maintaining a suppressed, circulating, free-virus population for a period of at least 3 years after the end of induction therapy) varies with the length and time of the induction phase [20]. Using a stochastic model, it was shown that induction therapy would have to last at least 180 days for cocktails containing two RTI-like drugs and a PI-like drug [20, 21].

Imperfect adherence has led to failure in suppressing viral replication and often mutations develop before or during induction therapy [6]. Since induction therapy has a significantly higher chance of pill fatigue than maintenance therapy, it is likely the patients will take some holidays during this period. Scientific literature cautions patients against taking any holidays while on therapy [22], but many patients are underadherent or nonadherent [23, 24]. Without guidance, patients who choose to stop therapy will each be making individual decisions, with no scientific basis. Recently, the question of short-term holidays (such as weekends) have been examined. Patients were highly adherent to five days on/two days off (FOTO) therapy. When asked about their preference for this type of therapy versus continuous HAART (Highly Active Antiretroviral Therapy), on a 10-point scale, the mean response was 9.7 [25].

Here, we examine the effects of imperfect adherence during the induction phase using a mathematical model of impulsive differential equations. We use the model to address the following research questions: 1. Can we determine the maximal length of a drug holiday and the number of subsequent doses that must be taken to avoid resistance? 2. How many drug holidays can be taken during the induction phase?

## Methods

### Modelling drug therapy

When modelling drug therapy and trying to approximate the number of doses a patient can miss without gaining drug resistance, it is important to have a reliable threshold that will guarantee that viral replication will not exceed a safe limit and so that the mutant strain will not appear. The inhibition of viral replication, *s*, can be described by

where *R*(*t*) is the drug and *IC*
_{50} is the concentration of drug which inhibits viral replication by 50% [7].

Thus, when *s* â‰ˆ 0, the drug has no effect, while if *s* â‰ˆ 1, the drug completely inhibits viral replication. See Figure 1.

Thus, the antiretroviral drug effect can be split into three regions: in Region 1, drug levels are insufficient to control either the wild-type or the mutant strain. In Region 2, drug levels are sufficient to control the wild-type strain but not a 10-fold mutant strain of the virus (ie a mutant strain that requires ten times the amount of drug to be controlled). In Region 3, drug levels are sufficient to control replication of both virus strains. These findings provide a threshold above which resistant viruses will be eradicated. We let *R*
_{2} be the threshold between Regions 2 and 3.

### The mathematical model

We adapt the mathematical model used from Smith & Wahl [26], to include latently infected cells [27]:

for *t* â‰ *t*
_{
k
}, where

In these equations, *V*
_{
I
}and *V*
_{
Y
}denote the wild-type and mutant virus respectively, *V*
_{
NI
}denotes the non-infectious virus, *T*
_{
S
}denotes the susceptible CD4^{+} T cells, *T*
_{
I
}denotes CD4^{+} T cells infected by the wild-type virus, denotes *T*
_{
LI
}denotes CD4^{+} T cells latently infected by the wild-type virus, *T*
_{
Y
}denotes CD4T^{+} T cells infected by the mutant virus, *T*
_{
LY
}denotes CD4^{+} T cells latently infected by the mutant virus, *T*
_{
RY
}denotes the noninfected CD4^{+} T cells which have absorbed enough drug so the wild-type strain is inhibited, but not enough to prevent infection from the mutant strain, *T*
_{
RY
}denotes the noninfected CD4^{+} T cells which have absorbed enough drug to prevent infection from both virus strains, *t* is the time in days, *n*
_{
I
}is the number of virions produced per infected cell per day, *Ï‰* is the fraction of virions produced per day by an infected CD4^{+} T cell, *d*
_{
V
}is the clearance rate of free virus, *r*
_{
I
}is the rate at which a susceptible cell becomes infected by the wild-type strain, *r*
_{
Y
}is the rate at which a susceptible cell becomes infected by the mutant strain, *d*
_{
S
}is the death rate of noninfected CD4^{+} T cells, *d*
_{
I
}is the death rate of infected CD4^{+}T cells, *Ïˆ* is the proportion of cells which become latently infected, *p*
_{
L
}is the rate at which latently infected cells become productive, *r*
_{
P
}is the rate at which the drug inhibits the wild-type T cells when drug concentrations are in Region 2, and *r*
_{
R
}and *r*
_{
Q
}are the rates at which the drug inhibits the wild- type and drug-resistant T cells, respectively, when drug concentrations are in Region 3. The constant *Î»* is the birth rate of CD4^{+} T cells, while *m*
_{
RI
}and *m*
_{
RY
}are the rates at which the drug is cleared from the intracellular compartment for intermediate and high drug concentrations, respectively. For parameter values and references, see [26].

The dynamics of a drug can be modelled using impulsive differential equations. The exponential decay can be written as a differential equation, where *R*(*t*) is the drug concentration during induction therapy. The dynamics of the drug are

with impulsive conditions, at times *t* = *t*
_{
k
},

The rate at which the drug is cleared is *d*
_{
r
}and *R*
^{i}is the dosage. Assuming a drug is taken at time *t*
_{
k
}, by the definition of an impulsive effect, we have

### Determining the Region 2 threshold

To find *R*
_{2}, the Region 2 threshold, we determined the time taken for resistance levels to reach a minimum. The drug levels at this time were evaluated from the antiretroviral effect curves and used as the *R*
_{2} threshold. This ensures that, when a drug holiday occurs, resistance levels are guaranteed to be low. Missing several doses increases resistance, but, by using the local minimum values, we ensure that resistance cannot emerge when patients are not taking a drug holiday.

To determine the threshold, note that

since *T*
_{
Y
}(0) = 0 at the beginning of infection. It follows that the viral load is initially decreasing. If the viral load reaches a minimum at time , then define *R*
_{2} = *R*(). This ensures that < 0 for 0 <*t* <. If the viral load decreases indefinitely, then we could define *R*
_{2} to be any value of *R* less than the trough value of the periodic orbit of the drug dynamics. However, this case is not realistic, since the virus does not clear on its own.

We define *R*
_{1} to be the value of *R* such

Thus, *R*
_{1} = 0.1*R*
_{2}. See Figure 1.

### Impulsive differential equations

The dynamics of both the wild-type and the resistant strains can be modelled using impulsive differential equations. Impulsive differential equations consist of a system of ordinary differential equations (ODEs), together with difference equations. Between "impulses", *t*
_{
k
}, the system is continuous, behaving as a system of ODEs. At the impulse points, there is an instantaneous change in state in some or all of the variables. This instantaneous change can occur when certain spatial, temporal or spatio-temporal conditions are met. We refer the interested reader to Bainov & Simeonov [28â€“30] and Lakshmikantham et al. [31] for more details on the theory of impulsive differential equations.

The change in drug concentration depends on whether a drug is taken or not. There is an instantaneous increase in the drug concentration immediately after a dose is taken and then an exponential decay while the drug is being absorbed in the body. The case of perfect adherence is illustrated in Figure 2A. However, as long as the drug concentration level does not drop below *R*
_{2}, there is a sufficient amount of drug to control both viral strains. We can thus determine the number of doses that can be missed and the number of doses subsequently taken in order to stay above the *R*
_{2} threshold. See Figure 2B.

The differential equations describing the virus and T cells depend on the dynamic behaviour of the drugs. Thus, for example, the rate of change of susceptible T cells decreases in Regions 2 or 3 (at different rates), but not in Region 1. The T cell and virus dynamics are continuous, but their derivatives are not, since those derivatives depend on the drugs, which are discontinuous. Since the drug equations decouple from the remaining equations, we develop theoretical results using the drug equations and apply those results numerically to the entire model.

## Results

### Theoretical results

We used our model to examine the effects of imperfect adherence on the induction phase of IM therapy. First, it is necessary to model perfect adherence to locate the impulsive periodic orbit in the drug levels. This provides a region where the drug concentration level must reach in order to sustain a low viral load. As can be seen in Figure 2, drug levels start at zero during induction therapy (since induction therapy starts at the beginning of drug therapy). Each time a drug is taken, the dose decays at a rate of , where is the value at which the drug starts to decay instantaneously after the drug is ingested. Since we assume perfect adherence, we get

where *Ï„* = *t*
_{
k+1 }- *t*
_{
k
}is the (fixed) time between doses for perfect adherence. We thus have

Furthermore,

as *n* â†’ âˆž.

Therefore, assuming perfect adherence, the impulsive orbit has endpoints

Knowing the values of the endpoints for the impulsive orbit after *n* = *n*
_{1} doses, we are able to incorporate imperfect adherence and see its effects. In order to avoid Region 2 after missing many doses and to maintain an average drug concentration level within Region 3, we impose conditions to ensure proper therapy. To guarantee successful induction therapy, after the first *n*
_{1} doses are taken, we will force the lower endpoint of the drug concentration to be within a tolerance *Îµ*
_{1} of the impulsive orbit. Thus, we require

Once the drug concentration level has reached the impulsive orbit, a patient may take a drug holiday. If *h*
_{1} doses are subsequently missed (see Figure 2B), then

In order to avoid Region 2 after *h*
_{1} doses are missed, we impose the condition . This will allow us to find the maximum number of doses a patient can miss after being *Îµ*
_{1} away from the impulsive orbit. This results in

After a patient has missed *h*
_{1} doses, in order to keep the viral replication low, they must take enough doses, *n*
_{2}, to return to the impulsive orbit. In the worst-case scenario, the exponential decay has reached Region 2; thus, starting at *R*
_{2}, we get

After *n*
_{2} doses are taken, we must impose a new condition that forces the drug concentration level to be *Îµ*
_{2} away from the impulsive orbit. We need

In order to determine the number of times a patient can miss a fixed amount of doses, we must verify if missing *h*
_{2} doses is the same as missing *h*
_{1} doses. After missing *h*
_{2} doses, we have

Patients are able to miss *h*
_{2} doses as long as their drug concentration levels do not drop below Region 2. Thus we repeat the same condition on *h*
_{2}:

At the end of induction therapy, *k* doses must be taken to ensure that, before the start of maintenance therapy, there is sufficient drug to control viral replication. After missing *h*
_{2} doses and assuming we are at Region 2, *k* subsequent doses are taken and the drug level becomes

As can be seen, because we started at the threshold after missing *h*
_{2} doses, as long as *n*
_{2} = *k*. Finally, after *k* doses, a patient needs to return to the periodic orbit. Thus, we impose

which is the same as the constraint for *n*
_{2} as long as *Îµ*
_{2} = *Îµ*
_{3}. If these conditions are satisfied, we are able to guarantee that the drug concentration levels do not enter Region 2 and significant drug resistance will not emerge.

### Imperfect adherence

The number of missable and subsequent doses that must be taken to avoid significant drug resistance for all FDA-approved drugs that are part of a PI-sparing cocktail is shown in Table 1. These are defined by (1) and (2), respectively. However, we stress that these results are theoretical and have not been tested clinically. In particular, it should be noted that pharmacokinetic parameters can vary from patient to patient.

There are fifteen FDA-approved PI-sparing triple-drug cocktails, for which we calculated (a) the initial number of doses that must be taken to be within a prescribed tolerance of perfect adherence, (b) the number of doses that could be missed without significant drug resistance emerging and (c) the number of doses that must be taken subsequently.

To determine the value of the prescribed tolerance, we examined two possibilities: a tolerance of 0.1 *Î¼M* and a tolerance of 0.01 *Î¼M*. That is, the number of doses is considered sufficient if the trough value of the periodic orbit of the drug dynamics is within 0.01 *Î¼M* of the trough value of therapy without drug holidays. We imposed a further condition: that the mean drug concentration be larger than the trough value of drugs when no drug holidays are taken. This is illustrated in Figure 3. This ensures that, over the length of the entire induction phase, drugs are maintained at sufficiently high levels (see [26] for more discussion). In Figure 3A, using a tolerance of 0.1 *Î¼M*, the overall mean drug concentration is below the trough value during therapy. Using a tolerance of 0.01 *Î¼M*, as shown in Figure 3B, shows that the overall mean drug concentration is above the trough value during therapy.

For the fifteen FDA-approved PI-sparing triple-drug cocktails, we identified the "weakest" drugs in each cocktail; ie, those for which the least number of doses can be missed. These drugs are Abacavir (ABC), Lamivudine (3TC), Stavudine (d4T), Emtricitabine (FTC), Zidovudine (ZDF), Didanosine (ddI) and Nevirapine (NVP). Thus, for each cocktail, the maximal number of missable doses is the same as that of its "weakest" drug. By combining the steps in (b) and (c) above, it was possible to theoretically calculate the number of drug holidays that could be taken during the inductive phase, based on the regimen for the "weakest" drug. See Table 2.

Since the minimum number of doses required to be taken and the maximum number of doses allowed to be missed follow a reliable pattern, we can extend this to fit into a baseline induction phase of 180 days [20]. This means, for example, that a patient taking the triple-drug cocktail FTC/TDF/EFV can theoretically have a 6 day holiday, as long as each holiday is followed by 17 days of perfect adherence; patients can take seven such holidays during the induction phase, and are thus able to miss a total of 42 days out of 180. A patient taking ABC/3TC/NVP can theoretically have sixteen drug holidays of 3 days each in a 180 day period, as long as each holiday is immediately followed by a 7 day period of strict adherence.

### Numerical simulations

In order to determine the long-term effects of taking the prescribed drug holidays, we simulated the worst-case scenario: monotherapy to the "weakest" drug in each combination from Table 2. This has the effect of overestimating the development of resistance: if no resistance is predicted to emerge during monotherapy, then it is unlikely to emerge during combination therapy. Conversely, if resistance does emerge during monotherapy, then there is no guarantee that it would emerge during combination therapy, due to the presence of the other two drugs.

We considered an extinction threshold of 2 Ã— 10^{-4} virions/mL. This corresponds to the concentration at which the virus falls below 1 per body. Thus, missing the maximum number of doses would theoretically lead to extinction of both strains (at least up to the level of detection), whereas missing more doses does not. However, it should be noted that we did not curtail the viral dynamics at this threshold.

We used the model in Section (describing the dynamic interaction between virus, T cells and drugs) and the calculations in Section (summarised in Tables 1 and 2) to illustrate our theoretical results. In order to demonstrate the effects of taking the prescribed drug holidays, we first ran simulations where patients missed the maximum number of doses and then took the required number of subsequent doses; this cycle was repeated for 180 days. Next, we ran the same simulations, with the same parameters, except that one additional dose of the drug was skipped during each drug holiday.

We performed these simulations for each of the "weakest" drugs identified in Table 2: ABC (Figure 4), 3TC (Figure 5), d4T (Figure 6), FTC (Figure 7), ZDV (Figure 8), ddI (Figure 9) and NVP (Figure 10). The first figure in each case illustrates the case of missing the maximal drug holiday and taking the minimum number of subsequent doses. The second figure in each case illustrates the same case, except that one additional dose was missed during each drug holiday. The exception is NVP, in which resistance did not emerge until three extra doses were missed (Figure 10B, inset).

For the first case, the wild-type virus oscillated at low levels during each drug holiday, but significant levels of resistance did not appear. Thus, taking the required number of doses successfully keeps the mutant strain at low levels. Conversely, missing one extra dose per holiday (three in the case of NVP) resulted in a significant buildup of resistance by the end of the induction phase.

In this case, there is a tremendous increase in the mutant strain by the end of the inductive phase, indicating that therapy has failed. Resistance to Abacavir increases from 10^{-3} to 10^{4}; resistance to Lamivudine increased from less than 10^{-3} to 10^{4}; resistance to Stavudine increased from 10^{-3} to 10^{3}; resistance to Emtricitabine increased from 10^{-3} to 10^{4}; resistance to Zidovudine increased from 10^{-3} to 10^{3}; resistance to Didanosine increased from 10^{-3} to 10^{2}; and resistance to Nevirapine increased from 10^{-3} to 10^{4}.

### Comparison with clinical results

A number of studies have attempted to characterise the safety of regular (and irregular) treatment interruptions, generally referred to as structured treatment interruptions (STIs). Pai *et al*. [22] summarised the to-date evidence of STIs in patients with chronic unsuppressed HIV infection due to drug-resistant HIV. They concluded that there were no significant virologic or immunologic benefit to STIs and that there is evidence that STIs have a prolonged negative impact on CD4 response and other disease events.

Subsequently, the SMART trial [32] examined CD4^{+} guided interruptions, of an average duration of 16 months. The DART trial [33] examined fixed 12 week interruptions. Both trials showed no benefit to these treatment interruptions. Indeed, the SMART trial was halted prematurely, due to significant morbidity and mortality among participants. Holkmann *et al*. [34] reported a two-fold risk of AIDS or death for patients who underwent treatment interruptions that lasted three months or longer.

It should be noted that all these trials involved lengthy periods of treatment interruption, of the order of weeks. Our results here recommend signficantly shorter periods of treatment interruption, of the order of days. Furthermore, our results predict significant increase in resistance if these periods are exceeded, consistent with the results from the majority of trials.

Shorter treatment interruptions have also been investigated. A study comparing interruptions of less than 7 days compared to longer interruptions showed that only 5% of men who discontinued HAART for short periods increased their HIV RNA. Conversely, men with longer interruptions had significantly higher rates (35.7 of HIV RNA increase [35]. Another study investigating cycles of 2-6 week fixed interruptions observed no clinically significant benefit with regard to viral suppression when off HAART, but also observed no evidence for an increase of viral resistance among patients undergoing repeated interruptions [36].

Recently, a pilot study examining five days on, two days off (FOTO) followed patients for 48 weeks [25]. Virologic suppression was maintained in 89.6% of patients. Combinations included 3TC/TDF/EFV, ABC/TDF/EFV, ddI/3TC/EFV and ABC/ddI/TDF/EFV; 100% of subjects on EFV- based regimens on the FOTO treatment schedule maintained virologic suppression at weeks 24 and 48. Combinations also included nevirapine- based regimens where one subject, on NVP/ABC/3TC/ZDV, had viral rebound at week 12 that was confirmed at week 16 on the FOTO schedule. It was also noted that 30% of the subjects on nevirapine-based regimens had blips of viral increase during therapy. Other combinations included TDF/3TC/NVP, ZDV/3TC/NVP, d4T/3TC/TDF/NVP. They also observed excellent adherence to the FOTO treatment schedule and a strong preference for this schedule compared to HAART. None of the observed rebounds in viral load were associated with the reported adherence of more than 2 days off therapy.

These preliminary results are in line with our theoretical recommendations. For regimens that include EFV-based regimens, all therapies included NRTIs and NNRTIs that we recommend a maximum of more than 2 days per drug holiday, followed by at least 5 days of subsequent therapy (Table 2). The NVP-based regimen with viral rebound included ZDV; our results predict that drug holidays on such a regimen should be no longer than 1.33 days (Figure 8). The three other NVP-based regimens with viral blips included ZDV and d4T; our results predict that neither would allow drug holidays as long as two days (Table 1).

### Sensitivity to variations

Since individual patients may respond differently to drugs, we explore the sensitivity of the number of missable doses to variations in parameters. The number of missable doses depends on the dosing interval, the drug decay rate, the drug concentration, the Region 2 threshold and the number of initial doses, which itself depends on the prescribed tolerance. Since we have already explored variations in the dosing interval and the prescribed decay rate, we now examine the variation with respect to the other parameters.

Figure 11 demonstrates the effect of variations in the drug decay rate, the Region 2 threshold and the drug concentration. Since the slope of the curves is low for the second and third figures, we conclude that the results are not highly sensitive to variations in the Region 2 threshold or the drug concentration, although small fluctuations may decrease the number of missable days (Figure 11B and Figure 11C). The outcome is more sensitive to variations in the drug decay rates, but is still not highly sensitive (Figure 11A).

## Discussion

It is vital to provide HIV patients with an effective drug regimen. Not only is it important that the drugs have a high effcacy, but it is also important that patients follow a regimen that will benefit both their mental and physical states. Since there are such a large number of patients who are unable to take their drugs regularly, it is important to understand the impact of drug holidays upon a patient's ability to control the virus. Induction therapy provides patients with the chance to submit to a very strict, but short, period of intense drug taking, followed by a long period of less-restrictive and more-relaxed therapy (maintenance therapy). We have demonstrated the effects of taking drug holidays during induction therapy. Instead of taking drugs two to three times a day for the entire length of the induction period, we were able to show that a patient can have drug holidays with sometimes as much as six days off each time. This form of treatment allows patients to take drug holidays with very little negative effect.

However, missing more doses than stated can highly affect the amount of resistant virus created. We have demonstrated that there is a large increase of mutant virus by simply missing one extra dose during each drug holiday (three for the TDF- FTC-NVP combination). Induction therapy with partial adherence works as long as a patient does not exceed the maximum length of the drug holiday; if they follow the prescribed regime, they can control the effects of drug resistance. It should be noted that Curlin *et al*. [20] showed that an induction phase on the order of 180-days was ideal for a triple-drug therapy including two RTI-like drugs and one PI-like drug. Since the results for a triple-drug therapy including three RTIs do not show a dramatic increase in resistant virus while taking the patterns suggested, we used an 180 day induction phase as a baseline.

These results apply to the fifteen FDA-approved, PI-sparing triple-drug cocktails, but simulations were only performed for the drugs with the least number of missable doses: Abacavir, Lamivudine, Stavudine, Emtricitabine, Zidovudine, Didanosine and Nevirapine. Missing one extra doses at the end of each drug holiday (three for Nevirapine) drastically increases the amount of resistant virus. However, it should be noted that the simulations were for monotherapy only and thus, in a triple-drug cocktail, the remaining two drugs inure against resistance.

Efavirenz and nevirapine only require a single mutation to confer resistance, and cross resistance affecting these three NNRTIs is common [37]. Both Lamivudine and Emtricitabine select for the M184V resistance mutation, which confers high-level resistance to both drugs, a modest decrease in susceptibility to Didanosine and Abacavir, and improved susceptibility to Zidovudine, Stavudine and Tenofovir [38]. It should be noted that our model assumes that the mutant is always present. By simulating the results for monotherapy, we illustrated the worst-case scenario; this is illustrated by Figure 10B, which shows that missing one extra dose per holiday is not disastrous; in this example, the mutant only takes hold when three extra doses are missed. Thus, our results are more conservative than is strictly necessary.

Double mutation happens less frequently; emergence of the M184V mutation is less frequent with Tenofovir/Emtricitabine than with Zidovudine/Lamivudine, while selection of the Lamivudine-associated M184V mutation to the Zidovudine/Lamivudine combination has been associated with increased susceptibility to Zidovudine [37]. It follows that, when the combinations are taken synchronously, the selection of mutants will be signficantly less likely than under monotherapy.

Discontinuous dosing is, of course, not realistic. There is a delay, the time-to- peak, between taking a drug and it reaching peak values in cells. Consequently, estimates based on maximal concentrations and terminal plasma half-lives could overestimate drug exposure. However, such delays can be approximated by an instantaneous change if the time-to-peak is sufficiently short, compared to the time between doses. This approximation has been shown to be robust, even for quite large delays [39].

Other limitations to our model are the assumption that the CD4^{+} pool of lymphocytes is the most significant source of HIV infection and that maintaining drug concentrations at clinical levels results in maximal control of virus replication. However, not all HIV-susceptible tissues are equally susceptible to antiretroviral drugs. For example, lymphoid cells in the gut are not completely suppressed [40]. These and other reservoirs will contribute to the long-term generation of virus particles, both during therapy and while undergoing a drug holiday. The relative rates of mutation or selection of resistant viruses for the various drugs are modelled via the choice of infection rate, *r*
_{
Y
}, compared to the infection rate, *r*
_{
I
}, for the wild-type strain. For numerical simulations, we used the intracellular half-life of each drug, if known; in the case of nucleosides, it is the cellular concentration of active nucleotide that is responsible for inhibition of viral reverse transcription. Furthermore, we assume that all tissues harbouring HIV are exposed to the same concentration of drug.

Previously [13], we showed how many doses can be missed for each PI- sparing drug, for only a single drug holiday during any given therapy. Here, we extend this to the case of more than one drug holiday.

Furthermore, all previous mathematical models of adherence considered therapy without an endpoint. Since induction therapy only occurs for a finite time, we have to consider the viral load when induction therapy ends. In particular, if a drug holiday coincided with the end of induction therapy, then the induction phase would functionally have ended at an earlier time and may thus be significantly less effective. Some of the key differences between our earlier work and the results provided here occur due to the fact that here we use 10-fold resistance, rather than 50- fold resistance; multiple holidays occur during a finite time interval; and the tolerance used was 0.01 *Î¼M* instead of 1% of the minimum value of periodic orbit; the tolerance we used here is more conservative.

Future work will investigate the effects of imperfect adherence to triple-drug cocktails involving protease inhibitors. We will also investigate the compounding effects of combination therapy in slowing the emergence of resistance and the effect of inter-individual variances in pharmacokinetics. Our modelling process could also be extended to additional treatment scenarios in which patients might be tempted to take drug holidays due to a high pill burden, such as booster therapies or the initial year of HAART.

## Conclusions

Using readily available pharmacokinetic data, we can theoretically determine the maximal length of drug holidays and the number of subsequent doses that must be taken. Since the induction phase lasts for a finite time, we can thus determine how many drug holidays can be taken within a 180-day induction period. Our theoretical results are in line with recent results concerning five-days-on/two-days-off (FOTO) for most cocktails, suggesting that drug holidays may be limited to very short breaks, rather than the longer holidays previously examined.

We thus conclude that induction therapy with partial adherence is tolerable, but the outcome depends on the drug cocktail. We have also demonstrated a robust method by which to determine therapy guidelines for patients who are unable or unwilling to adhere completely. Treatment interruptions, if they occur, must be short and followed by a strict period of dose taking. Thus, while continuous therapy is preferable, FOTO therapy is acceptable for all RTI cocktails except those containing ZDV, d4T or DLV, which can only tolerate extremely short drug holidays.

## References

Steinbrook R: The AIDS Epidemic - A Progress Report from Mexico City. N Engl J Med. 2008, 359: 885-887. 10.1056/NEJMp0805761.

McCune J: The dynamics of CD4

^{+}T-cell depletion in HIV disease. Nature. 2001, 410: 974-979. 10.1038/35073648.Yeni PG, Hammer SM, Hirsch M, Saag M, Schechter M, Carpenter CCJ, Fischi MA, Gatell JM, Gazzard BG, Jacobsen D, Katzenstein DA, Montaner JSG, Richman DD, Schooley RT, Thompson MA, Vella S, Volberding PA: Treatment for adult HIV infection: 2004 recommendations of the International AIDS Society - USA Panel. JAMA. 2004, 292: 251-265. 10.1001/jama.292.2.251.

Staszewski S, Morales-Ramirez J, Tashima KT, Rachlis A, Skiest D, Stanford J, Stryker R, Johnson P, Labriola DF, Farina D, Manion DJ, Ruiz NM: Efavirenz plus zidovudine and lamivudine, efavirenz plus indinavir and indinavir plus zidovudine and lamivudine in the treatment of HIV-1 infection in adults. N Engl J Med. 1999, 341: 1865-1873. 10.1056/NEJM199912163412501.

Moyle G: Protease inhibitor-sparing regimens: new evidence strengthens position. J Acq Immun Def Synd. 2003, 33 (Suppl 1): 17-25.

Atlice FL, Friedland GH: The era of adherence to HIV therapy. Ann Intern Med. 1998, 129: 503-505.

Wahl LM, Nowak MA: Adherence and drug resistance: predictions for therapy outcome. Proc R Soc B. 2000, 267: 835-843. 10.1098/rspb.2000.1079.

Philips AN, Youle M, Johnson M, Loveday C: Use of stochastic model to develop understanding of the impact of different patterns of antiretroviral drug use on resistance development. AIDS. 2001, 15: 2211-2220. 10.1097/00002030-200111230-00001.

Tchetgen E, Kaplan EH, Friedland GH: Public health consequences of screening patients for adherence to highly active antiretroviral therapy. J AIDS. 2001, 26: 118-129.

Huang Y, Rosenkranz S, Wu H: Modeling HIV dynamics and antiviral response with consideration of time-varying drug exposures, adherence and phenotypic sensitivity. Math Biosci. 2003, 184: 165-186. 10.1016/S0025-5564(03)00058-0.

Huang Y, Liu D, H W: Hierarchical Bayesian methods for estimation of parameters in longitudinal HIV dynamic system. Biometrics. 2004, 62: 413-423. 10.1111/j.1541-0420.2005.00447.x.

Ferguson NM, Donnelly CA, Hooper J, Ghani AC, Fraser C, Bartley L: Adherence to antiretroviral therapy and its impact on clinical outcome in HIV-infected patients. J R Soc Interface. 2005, 2: 349-363. 10.1098/rsif.2005.0037.

Smith RJ: Adherence to antiretroviral HIV drugs: how many doses can you miss before resistance emerges?. Proc R Soc B. 2006, 273: 617-624. 10.1098/rspb.2005.3352.

Krakovska O, Wahl LM: Optimal drug treatment regimens for HIV depend on adherence. J Theor Biol. 2007, 246: 499-509. 10.1016/j.jtbi.2006.12.038.

Friedland GH, Williams A: Attaining higher goals in HIV treatment: the central importance of adherence. AIDS. 1999, 13 (Suppl 1): 61-72.

Reijers MH, Weverling GJ, Jurriaans S, Wit FW, Weigel HM, Ten Kate RW, Mulder JW, Frissen PH, van Leeuwen R, Reiss P, Schuitemaker H, de Wolf F, Lange JM: Maintenance therapy after quadruple induction therapy in HIV-1 infected individuals: Amsterdam Duration of Antiretroviral Medication study. Lancet. 1998, 352: 185-190. 10.1016/S0140-6736(98)06193-5.

Havlir DV, Marschner IC, Hirsch MS, Collier AC, Tebas P, Bassett RL, Ioannidis JP, Holohan MK, Leavitt R, Coone G, Richman DD: Maintenance antiretroviral therapies in HIV infected patients with undetectable plasma HIV RNA after triple-drug therapy. AIDS Clinical Trials Group Study 343 Team. N Engl J Med. 1998, 339: 1261-1268. 10.1056/NEJM199810293391801.

Pialoux G, Raffi F, Brun-Vezinet F, MeiffrÃ©dy V, Flandre P, Gastaut J, Dellamonica P, Yeni P, Delfraissy J, Aboulker J: A randomized trial of three maintenance regimens given after three months of induction therapy with zidovudine, lamivudine, and indinavir in previously untreated HIV-1-infected patients. Trilege (Agence Nationale de Recherches sur le SIDA 072) Study Team. N Engl J Med. 1998, 339: 1269-1276. 10.1056/NEJM199810293391802.

Descamps D, Flandre P, Calvez V, Peytavin G, MeiffrÃ©dy V, Collin G, Delaugerre C, Robert-Delmas S, Bazin B, Aboulker JP, Pialoux G, Raffi F, Brun-VÃ©zinet F: Mechanisms of virologic failure in previously untreated HIV-infected patients from a trial of induction-maintenance therapy. JAMA. 2000, 283: 205-211. 10.1001/jama.283.2.205.

Curlin M, Iyer S, Mittler J: Optimal timing and duration of induction therapy for HIV-1 infection. PLoS Comput Biol. 2007, 3 (7): e133-10.1371/journal.pcbi.0030133.

Havlir D, Hellmann NS, Petropoulos CJ, Whitcomb JM, Collier AC, Hirsch MS, Tebas P, Sommadossi JP, Richman DD: Drug Susceptibility in HIV infection after Viral Rebound in Patients Recieving Indinavir-Containing Regimens. JAMA. 2000, 283: 299-234. 10.1001/jama.283.2.229.

Pai N, Lawrence J, Reingold AL, Tulsky J: Structured treatment interruptions (STI) in chronic unsuppressed HIV infection in adults. Cochrane Database of Systematic Reviews. 2006, 3: CD006148-

Chesney MA, Morin M, Sherr L: Adherence to HIV combination therapy. Social Science & Medicine. 2000, 50: 1599-1605. 10.1016/S0277-9536(99)00468-2.

Curlin ME, Wilkin T, Mittler J: Induction-maintenance therapy for HIV-1 infection. Future HIV Therapy. 2008, 2: 175-185. 10.2217/17469600.2.2.175.

Cohen CJ, Colson AE, Sheble-Hall AG, McLaughlin KA, Morse GD: Pilot Study of a Novel Short-Cycle Antiretroviral Treatment Interruption Strategy: 48-Week Results of the Five-Days-On, Two-Days-Off (FOTO) Study. HIV Clin Trials. 2007, 8: 19-23. 10.1310/hct0801-19.

Smith RJ, Wahl LM: Drug resistance in an immunological model of HIV-1 infection with impulsive drug effects. Bull Math Biol. 2005, 67: 783-813. 10.1016/j.bulm.2004.10.004.

Smith RJ, Aggarwala BD: Can the viral reservoir of latently infected CD4+ T cells be eradicated with antiretroviral drugs?. J Math Biol. 2009, 59: 697-715. 10.1007/s00285-008-0245-4.

Bainov DD, Simeonov PS: Systems with Impulsive Effect. 1989, Chichester: Ellis Horwood Ltd

Bainov DD, Simeonov PS: Impulsive differential equations: periodic solutions and applications. 1993, Burnt Mill: Longman Scientific and Technical

Bainov DD, Simeonov PS: Impulsive Differential Equations: Asymptotic Properties of the Solutions. 1995, Singapore: World Scientific

Lakshmikantham V, Bainov DD, Simeonov PS: Systems with Impulsive Effect. 1989, Singapore: World Scientific

Strategies for Management of Antiretroviral Therapy (SMART) Study Group, El-Sadr WM, Lundgren JD: CD4+ count-guided interruption of antiretroviral treatment. N Engl J Med. 2006, 355 (22): 2283-2296. 10.1056/NEJMoa062360.

DART Trial Team: Fixed duration interruptions are inferior to continuous treatment in African adults starting therapy with CD4 cell counts < 200 cells/microl. AIDS. 2008, 22 (2): 237-247. 10.1097/QAD.0b013e3282f2d760.

Holkmann Olsen C, Mocroft A, Kirk O, Vella S, Blaxhult A, Clumeck N, Fisher M, Katlama C, Phillips A, Lundgren J: Interruption of combination antiretroviral therapy and risk of clinical disease progression to AIDS or death. HIV Med. 2007, 8: 96-104. 10.1111/j.1468-1293.2007.00436.x.

Li X, Margolick J, Conover C, Badri S, Riddler SA, Witt MD, Jacobson LP: Interruption and Discontinuation of Highly Active Antiretroviral Therapy in the Multicenter AIDS Cohort Study. JAIDS. 2005, 38: 320-328.

Papasavvas E, Kostman JR, Mounzer K, Grant RM, Gross R, Gallo C, Azzoni L, Foulkes A, Thiel B, Pistilli M, Mackiewicz A, Shull J, Montaner LJ: Randomized, Controlled Trial of Therapy Interruption in Chronic HIV-1 Infection. PLoS Medicine. 2004, 1 (3): e64-10.1371/journal.pmed.0010064.

Panel on Antiretroviral Guidelines for Adults and Adolescents. Guidelines for the use of antiretroviral agents in HIV-1-infected adults and adolescents. Department of Health and Human Services. 2008, 1-139. [http://www.aidsinfo.nih.gov/ContentFiles/AdultandAdolescentGL.pdf]

Ait-Khaled M, Stone C, Amphlett G, Clotet B, Staszewski S, Katlama C, Tisdale M, CNA3002 International Study Team: M184V is associated with a low incidence of thymidine analogue mutations and low phenotypic resistance to zidovudine and stavudine. AIDS. 2002, 16: 1686-1689. 10.1097/00002030-200208160-00017.

Smith RJ, Schwartz EJ: Predicting the potential impact of a cytotoxic T-lymphocyte HIV vaccine: how often should you vaccinate and how strong should the vaccine be?. Math Biosci. 2008, 212: 180-187. 10.1016/j.mbs.2008.02.001.

Chun TW, Nickle DC, Justement JS, Meyers JH, Roby G, Hallahan CW, Kottilil S, Moir S, Mican JM, Mullins JI, Ward DJ, Kovacs JA, Mannon PJ, Fauci AS: Persistence of HIV in gut-associated lymphoid tissue despite long-term antiretroviral therapy. J Infect Dis. 2008, 197: 714-720. 10.1086/527324.

Scott L, Perry C: Delavirdine - A review of its use in HIV infection. Drugs. 2000, 60: 1411-1444. 10.2165/00003495-200060060-00013.

Perry C, Noble S: Didanosine - An update review of its use in HIV infection. Drugs. 1999, 58: 1099-1135. 10.2165/00003495-199958060-00009.

Bhana N, Ormrod D, Perry C, Figgitt D: Zidovudine - A review of its use in management of vertically-aquired pediatric HIV infection. Drugs. 2002, 4: 515-553.

Hervey P, Perry C: Abacavir - A review of its clinical potential in patients with HIV infection. Drugs. 2000, 60: 447-479. 10.2165/00003495-200060020-00015.

Frampton J, Perry C: Emtricitabine - A review of its use in management of HIV infection. Drugs. 2006, 65: 1427-1448. 10.2165/00003495-200565100-00008.

Moore KHP, Barrett JE, Shaw S, Pakes GE, Churchus R, Kapoor A, Barry MG, D B: The pharmacokinetics of lamivudine phosphorylation in peripheral blood mononuclear cells from patients infected with HIV-1. AIDS. 1999, 13: 2239-2250. 10.1097/00002030-199911120-00006.

Fletcher C, Brundage R, Remmel R, Page L, Weller D, Calles N, Simon C, Kline M: Pharmacologic characteristics of indinivir, didanosine, stavudine, in human immunodeficiency virus-infected children receiving combination therapy. Antimicrobial Agents and Chemotherapy. 2000, 44: 1029-1034. 10.1128/AAC.44.4.1029-1034.2000.

Hawkins T, Veikley W, St Claire R, Guyer B, Clark N, Kearney B: Intracellular pharmacokinetics of tenofovir diphosphate, carbovir triphosphate, and lamivudine triphosphate in patients receiving triple-nucleoside regimens. J Acquir Defic Syndr. 2001, 39: 406-411. 10.1097/01.qai.0000167155.44980.e8.

### Pre-publication history

The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/10/6/prepub

## Acknowledgements

The authors are grateful to Frithjof Lutscher, Lindi Wahl and Jing Li for technical discussions, and to Robert Stengel, Selwyn Hurwitz and John Mittler for their careful reading of the manuscript and constructive comments that greatly improved the final version. RJS is supported by an NSERC Discovery grant, an Early Researcher Award and funding from MITACS.

## Author information

### Authors and Affiliations

### Corresponding author

## Additional information

### Competing interests

The authors declare that they have no competing interests.

### Authors' contributions

Both authors wrote the manuscript and performed numerical simulations. RJS designed the study, while REM performed the mathematical analysis. Both authors read and approved the final manuscript.

## Authorsâ€™ original submitted files for images

Below are the links to the authorsâ€™ original submitted files for images.

## Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## About this article

### Cite this article

Miron, R.E., Smith, R.J. Modelling imperfect adherence to HIV induction therapy.
*BMC Infect Dis* **10**, 6 (2010). https://doi.org/10.1186/1471-2334-10-6

Received:

Accepted:

Published:

DOI: https://doi.org/10.1186/1471-2334-10-6