# Rational Design and Adaptive Management of Combination Therapies for Hepatitis C Virus Infection

Ruian Ke, Claude Loverdo, Hangfei Qi, Ren Sun, James O. Lloyd-Smith

► To cite this version:

Ruian Ke, Claude Loverdo, Hangfei Qi, Ren Sun, James O. Lloyd-Smith. Rational Design and Adaptive Management of Combination Therapies for Hepatitis C Virus Infection. PLoS Computational Biology, 2015, 11 (6), pp.e1004040. 10.1371/journal.pcbi.1004040. hal-01224118

HAL Id: hal-01224118

https://hal.sorbonne-universite.fr/hal-01224118

Submitted on 4 Nov 2015

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

CrossMark
click for updates

RESEARCH ARTICLE

# Rational Design and Adaptive Management of Combination Therapies for Hepatitis C Virus Infection

Ruian Ke¹*, Claude Loverdo¹,², Hangfei Qi³, Ren Sun³,⁴,⁵,⁶, James O. Lloyd-Smith¹,⁷*

1 Department of Ecology and Evolutionary Biology, University of California, Los Angeles, Los Angeles, California, United States of America, 2 CNRS/UPMC Univ Paris 06, UMR 8237, Laboratoire Jean Perrin LJP, Paris, France, 3 Department of Molecular and Medical Pharmacology, University of California, Los Angeles, Los Angeles, California, United States of America, 4 The Molecular Biology Institute, University of California, Los Angeles, Los Angeles, California, United States of America, 5 Department of Infectious Diseases, Novartis Institutes for BioMedical Research, Emeryville, California, United States of America, 6 Zhejiang University, Hangzhou, China, 7 Fogarty International Center, National Institutes of Health, Bethesda, Maryland, United States of America

□ Current address: MS-K710, T-6 Theoretical Biology and Biophysics, Los Alamos National Laboratory, Los Alamos, New Mexico, United States of America

* rke.work@gmail.com (RK); jlloydsmith@ucla.edu (JOLS)

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Ke R, Loverdo C, Qi H, Sun R, Lloyd-Smith JO (2015) Rational Design and Adaptive Management of Combination Therapies for Hepatitis C Virus Infection. PLoS Comput Biol 11(6): e1004040. doi:10.1371/journal.pcbi.1004040

**Editor:** Roger Dimitri Kouyos, University of Zurich, SWITZERLAND

**Received:** August 2, 2014

**Accepted:** November 13, 2014

**Published:** June 30, 2015

**Copyright:** © 2015 Ke et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This work was supported by National Science Foundation (grant number EF-0928690 to JOLS). JOLS is grateful for support from the De Logi Chair in Biological Sciences, and from the RAPIDD program of the Science & Technology Directorate, Department of Homeland Security, and the Fogarty International Center, National Institutes of Health. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

## Abstract

Recent discoveries of direct acting antivirals against Hepatitis C virus (HCV) have raised hopes of effective treatment via combination therapies. Yet rapid evolution and high diversity of HCV populations, combined with the reality of suboptimal treatment adherence, make drug resistance a clinical and public health concern. We develop a general model incorporating viral dynamics and pharmacokinetics/ pharmacodynamics to assess how suboptimal adherence affects resistance development and clinical outcomes. We derive design principles and adaptive treatment strategies, identifying a high-risk period when missing doses is particularly risky for *de novo* resistance, and quantifying the number of additional doses needed to compensate when doses are missed. Using data from large-scale resistance assays, we demonstrate that the risk of resistance can be reduced substantially by applying these principles to a combination therapy of daclatasvir and asunaprevir. By providing a mechanistic framework to link patient characteristics to the risk of resistance, these findings show the potential of rational treatment design.

## Author Summary

Hepatitis C virus (HCV) affects approximately 170 million people world-wide and chronic infections can lead to cirrhosis and liver cancer. New combination therapies of direct acting antivirals have achieved remarkably high cure rates in clinical trials. However, high mutation rates and high diversity of HCV populations, combined with the reality of suboptimal treatment adherence, make drug resistance a clinical and public health concern. By constructing a mechanistic framework to assess the risk of drug resistance, we provide
**Competing Interests:** The authors have declared that no competing interests exist.

guidelines for rational design and adaptive management of these promising new therapies. In particular, we identify a high-risk period when missing doses is particularly risky, and quantify the number of extra doses needed to compensate when doses are missed. This framework is a step towards developing a tool for clinicians to design combination therapies and adaptively manage treatment regimens to achieve favorable clinical outcomes.

Introduction

Hepatitis C virus (HCV) affects approximately 170 million people and chronic infections can lead to cirrhosis and hepatocellular carcinoma [1,2]. Recently, development of direct acting antivirals (DAAs) against HCV infection has revolutionized the field of HCV treatment, because of their high potency, broad applicability and mild side effects [3,4]. Combination therapies of DAAs have achieved remarkably high rates of sustained virological response in clinical trials [5–10]. However, most DAAs have relatively low genetic barriers [11–13], with the exceptions of a few pan-genotypic, yet high-cost DAAs [6]. Because of the high intrinsic mutation rate of HCV [14,15] and the resulting high viral diversity [1,16,17], combined with the reality of suboptimal treatment adherence [18,19], viral resistance is still a clinical and public health concern [13,20]. This is especially true for high-risk groups such as patients with psychiatric disorders or depression [21], and in resource-limited settings where patients have limited access to clinical cares and cannot afford the expensive pan-genotypic DAAs with high genetic barriers [22,23]. If treatment is not properly managed, resistance could quickly develop to combination therapies and render these new DAAs useless, as observed for other antimicrobial treatments, squandering the potential health gains from these recent breakthroughs [24–26].

Suboptimal patient adherence to dosing regimens is a crucial risk factor for resistance development in both HIV and HCV treatments [18,19,27,28]. Although high rates of sustained virological response have been achieved in clinical trials, adherence levels may vary substantially among the vast population of infected patients, owing to long treatment periods, complicated regimens associated with DAA combination therapies and limited access to health care [18,19,29–31]. Rational design of combination therapy that achieves viral eradication in patients and maximizes the durability of available DAAs in the presence of suboptimal adherence is a research priority [18,32–34]. In addition, theories that guide individualized regimens based on the genetic composition of a patient's infection and real-time adjustments for missed doses are needed to avoid resistance. Mathematical models are well suited to address this problem. Previous modeling studies for HIV infections have illuminated potential mechanisms underlying treatment failure and explained puzzling clinical observations [35,36]. However, HCV is a curable disease and its infection, goal of treatment and mechanism of resistance differ from HIV in many respects [37], including no known latent reservoir and a finite treatment period to eradicate the virus. Here, by integrating pharmacokinetics/pharmacodynamics (PK/PD) and viral dynamics into mathematical models, we develop the first general theory to assess the impacts of suboptimal adherence on the outcome of DAA-based therapies for HCV infection. We derive design principles that can be generalized to therapies involving different classes and different numbers of drugs. Using large-scale data from in vitro resistance assays and human clinical trials, we apply this framework to a combination therapy of daclatasvir and asunaprevir [38], and derive evidence-based adaptive treatment strategies for treatment protocols over time according to resistance profiles and adherence patterns.
## Results

Resistance to antiviral treatments can develop through selection of preexisting mutants or *de novo* generation of new mutants. A core principle for designing effective combination therapy is that, if patients fully adhere to the treatment regimen, the treatment must suppress all preexisting mutants and *de novo* resistance should be unlikely [39]. Missing doses, however, can lead to suboptimal drug concentrations, allowing growth of some preexisting mutants with partially resistant phenotypes. Growth of these mutants allows the viral population to survive longer, possibly generating further mutations that contribute *de novo* resistance against the full combination therapy. For example, consider a combination therapy of two DAAs, A and B, as shown in Fig 1A. If missed doses and pharmacokinetics lead to a drop in the concentration of drug A, this allows growth of the preexisting mutant, AB', (which is already resistant to drug B), thus opening opportunities to generate the fully resistant mutant, A'B'. Therefore, the dynamics of the subset of preexisting mutants that have a high level of resistance against single DAAs determine resistance evolution and treatment outcomes for combination therapies. In the following, we denote these mutants as 'partially resistant' mutants.

## The effective viral fitness, $R_{eff}(t)$

The fitness of a particular strain in a treated patient is determined by the PK/PD of the drug, the level of resistance of the strain, and the availability of target cells, i.e. uninfected hepatocytes for HCV (Fig 1B). We can integrate all these factors (for any class of DAA therapy) into a single number, the effective reproductive number under treatment, $R_{eff}(t)$ (Fig 1C). $R_{eff}(t)$ is the average number of cells infected by viruses produced by a single infected cell. It acts as a measure of viral fitness, and can be calculated as:

$$R_{eff}(t) = (1 - \varepsilon(t)) \cdot R_0 \cdot h(t) \quad (1)$$

where $t$ is time since treatment starts, $\tau$ is the time since last dose, $\varepsilon(\tau)$ is the efficacy of the drug combination at time $\tau$ during the dosing cycle, $R_0$ is the reproductive number of the virus in the absence of treatment, and $h(t)$ is the normalized abundance of target cells (see S1 Text). Under effective treatment, the availability of target cells, $h(t)$, increases quickly to reach the infection-free level [40], and therefore, the overall viral fitness increases over time as $h(t)$ increases under effective treatment (Fig 1B and 1C). When adherence is optimal (i.e. no missing doses), the value of $R_{eff}$ for a partially resistant mutant is always less than 1 (i.e. viral suppression); however, if doses are missed, drug concentration declines exponentially and $R_{eff}$ can become greater than 1 (i.e. viral growth) (Fig 1C). Note that, although we consider the fitness of a single 'partially resistant' mutant here, the competition between different mutants is described implicitly in Eq 2 by the target cell availability, $h(t)$: if a 'partially resistant' mutant rises to a high abundance due to missing doses, then $h(t)$ will decrease to a low level again, leading to decreases in fitness for all viral mutants.

## The growth of partially resistant mutants and the need for extended treatment

We now consider how suboptimal adherence impacts the dynamics of partially resistant mutants. As an illustration, we contrast simulations assuming perfect adherence versus suboptimal adherence. Missing doses leads to rapid decreases in drug concentration, and thus increases in $R_{eff}$ of a partially resistant mutant (Fig 2A–2C). This means that extra doses are needed to compensate for the missed doses to suppress the mutant to extinction (Fig 2D), and also that the
A

<figure><img src="image_4.png" /></figure>

B

<figure><img src="image_5.png" /></figure>

C

<figure><img src="image_6.png" /></figure>

**Fig 1.** The impacts of suboptimal adherence on viral fitness. (A) A schematic illustrating how a non-preexisting mutant, A'B', fully resistant to a combination therapy involving two drugs, A and B, can be generated when adherence is suboptimal. Each black circle represents a mutant on the parameter space of resistance levels to A and B. AB, A'B and AB' are preexisting mutants that are non-resistant, resistant to A only and resistant to B only, respectively. Colored areas denote parameter regimes where mutants are fully resistant to the therapy (red), can grow if doses are missed (pink), and do not grow (blue). Note that the pink area can grow or shrink on the parameter space depending on the number of consecutively missed doses and drug PK/PD, and mutants lying in the pink area are 'partially resistant mutants'. (B) The dynamics of viral strains under treatment are determined by several factors: drug concentration, [D], which decreases with an increasing number of missed doses, *m* (upper panel); how viral replication is affected by drug (l-ε; middle panel); and the relative number of target cells, *h*(*t*) (lower panel). Upon effective treatment, *h*(*t*) increases to the infection-free level. (C) We integrate all these factors into a single fitness parameter, *R*<sub>eff</sub>(*t*). Viral fitness increases as drug concentration drops (indicated by shades of green) and as target cell abundance rises (the blue arrow). Values of *R*<sub>eff</sub>(*t*) can exceed 1, i.e. positive growth, if doses are missed after a period of effective treatment.

doi:10.1371/journal.pcbi.1004040.g001

number of newly infected cells rises substantially, which increases the opportunity for *de novo* resistance (Fig 2E).
<figure><img src="image_7.png" /><figcaption>**Fig 2.** Suboptimal adherence prolongs treatment time needed to eliminate partially resistant mutants and increases the risk of de novo evolution of fully resistant strains. Two simulations assuming perfect adherence (dashed black lines) and imperfect adherence (solid red lines) are shown. In the simulation assuming imperfect adherence, single doses are missed at day 5 and 15, and 3 consecutive doses are missed during days 22–24. (A) Drug concentration over time normalized by the maximum drug concentration $C_{max}$. (B) The abundance of target hepatocytes, which rebounds after the initiation of combination therapy. (C) Viral fitness of the partially resistant mutant under consideration. Missing doses increases the value of $R_{eff}$, especially when multiple consecutive doses are missed or when $h(t)$ has increased to high levels. (D) The dynamics of cells infected by the partially resistant mutant (on $log_{10}$ scale). The number of infected cells declines almost exponentially when doses are taken. Missed doses allow the number of infected cells to rebound. This means that an additional period of treatment is needed to suppress the mutant below the extinction threshold level, i.e. to achieve viral elimination. (E) The number of cells newly infected by the partially resistant mutants. Missed doses lead to substantial numbers of additional new infections, especially when 3 consecutive doses are missed.</figcaption></figure>

doi:10.1371/journal.pcbi.1004040.g002

We approximate the time-varying values of $R_{eff}(t)$ during periods when doses are missed, by calculating the average effective reproductive number, $R_{ave,m}$, as (see Materials and Methods):

$$R_{ave,m}(t) \approx (1 - \varepsilon_{ave,m}) \cdot R_0 \cdot h(t) \quad (2)$$

where $t$ is the time when the patient starts to miss doses, $m$ is the number of consecutive doses missed and $\varepsilon_{ave,m}$ is the average drug inhibition during the period when $m$ consecutive doses are missed. This allows us to generalize our theory to any DAA combinations for which $\varepsilon_{ave,m}$ can be either estimated from pharmacokinetics/pharmacodynamics data or calculated from mutant resistance profiles [36].
We then ask, if *m* consecutive doses are missed beginning at time *t*, how many extra doses, *N*<sub>*m*</sub>, are needed to compensate? This number, which we denote ‘compensatory doses’, can be approximated as (see *Materials and Methods*):

$$N_m(t) \approx m \cdot R_{ave,m}(t) \approx m \cdot (1 - \varepsilon_{ave,m}) \cdot R_0 \cdot h(t) \quad (3)$$

This allows us to estimate the total duration of treatment needed to clear infection for a given adherence pattern. Furthermore, since *h*(*t*) increases over time under effective treatment [40], Eq 3 shows that a higher number of extra doses are needed to eliminate the infection if doses are missed later in treatment.

## De novo generation of fully resistant mutants

To assess the risk that a partially resistant lineage will give rise to full resistance, we calculate the expected number of target cells, $\Phi_m$, that become infected by fully resistant mutant viruses due to *de novo* mutation during a period when *m* consecutive doses are missed. This quantity is the product of the cumulative number of cells newly infected by a partially resistant mutant and the effective mutation rate from that mutant to the fully resistant mutant, $\mu_{eff}$ (see *Materials and Methods*):

$$\Phi_m(t) \approx \mu_{eff} \cdot I(t) \cdot \underbrace{\frac{R_{ave,m}(t)}{R_{ave,m}(t) - 1} \cdot (e^{(R_{ave,m}(t)-1)\cdot\delta \cdot m \cdot T} - 1)}_{\text{cumulative number of cells newly infected by a partially resistant mutant}} \quad (4)$$

where *I*(*t*) is the number of cells infected by the partially resistant mutant at time *t* when the first dose is missed, and $\Theta(t)$ represents the potential to generate new infections. $\delta$ is the death rate of infected hepatocytes, and *T* is the scheduled interval between two doses. $\Phi_m$ quantifies the risk that a fully resistant mutant infects target cells, but whether it emerges and becomes established within the host depends on its fitness and the stochastic dynamics of invasion [41–43].

The strong dependence of $\Phi_m$ on $\mu_{eff}$ predicts that designing combination therapies to increase the genetic barrier to full resistance, e.g. using DAAs with higher genetic barrier or adding an extra drug into the combination, can reduce $\Phi_m$ by orders of magnitude or more, thus it would lead to drastic reductions in the probability of generating full resistance (compare trajectories a and b in Fig 3A).

Eq 4 also allows us to assess when during treatment it is most risky to miss doses, which can inform treatment guidelines. Changes in two quantities, $I(t)$ and $\Theta(t)$, determine changes in $\Phi_m$ over the course of a treatment regimen. For as long as adherence is perfect, $I(t)$ decreases exponentially, while $\Theta(t)$ increases over time since $R_{ave,m}(t)$ increases as the abundance of target cells rises over time (Fig 3B). Thus the value of $\Phi_m$ first increases (due to rapid recovery of target cells) and then decreases exponentially (due to decrease of infected cells). This leads to a high-risk window period, during which missing doses is especially risky for generating full resistance (Fig 3A). This qualitative finding is robust to changes in model parameters, though quantitative predictions of the risk of full resistance depend on the fitness of the mutant ($R_0$), the half-life of infected cells ($\delta$), and the rate at which the target hepatocytes become available upon treatment (S1 Fig).
<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /><figcaption>**Fig 3.** There is a high-risk window early in treatment when missing doses is more likely to cause *de novo* resistance. (A) The changes in the risk of *de novo* resistance, Φ<sub>m</sub>, generated by a partially resistant mutant over time. The two sets of trajectories, A and B, differ in that the value of μ<sub>eff</sub> for trajectory B is smaller by a factor of 10<sup>-5</sup> (representing one additional nucleotide mutation) than the value set for trajectory A. Each set of trajectories shows the risk when the number of doses missed (m) is 1, 2, or 3. (B) Dynamics of the two time-varying quantities in Eq 4, i.e. the number of cells infected by the partially resistant mutant relative to the initial number before treatment (I(t)/I(0); blue dashed line), and the value of Θ(t), green dotted lines, as shown in Eq 4. Under effective treatment, the number of infected cells I(t) decreases exponentially, while the number of target cells rebounds to the infection-free level quickly, causing an increase in R<sub>ave,m</sub> and thus Θ(t). Together these changes cause Φ<sub>m</sub> to increase initially and then to decrease exponentially at longer times (as seen in panel A).</figcaption></figure>

doi:10.1371/journal.pcbi.1004040.g003

Design principles and adaptive treatment strategy for DAA combination
therapy

These results suggest principles for designing combination therapies and rational optimization
of treatment outcomes. First, the genetic barrier to full resistance to a therapy is an important
determinant of the risk of resistance. Assessment of treatment readiness has been a low-cost
routine clinical practice for HIV treatment [44]. Similar strategies can be implemented for
HCV treatment. Based on the assessment, for patients who are predicted to maintain high ad-
herence, combinations of DAA that ensure the fully resistant mutants are not pre-existing
would be sufficient. For patients with risk factors for low adherence, therapies should be de-
signed by selecting drug combinations that impose a higher genetic barrier than required to
suppress all pre-existing mutants. Second, we have shown there exists a high-risk window peri-
od where the risk of *de novo* resistance is high. Intervention efforts to ensure a high level of ad-
herence during the high-risk window period (indicated by the value of Φ<sub>m</sub>) would reduce the
risk of resistance and treatment failure. Third, because of the exponential growth of ‘partially
resistant mutants' when doses are missed, missing a number of doses consecutively leads to a much higher risk of *de novo* resistance than missing the same number of doses separately [36]. Thus, missing a block of doses should be avoided.

*Adaptive treatment strategies* could be developed based on the theoretical findings shown above. If doses are missed during treatment, the patient should be treated with extra doses, computed as the maximum value of the $N_m$ values calculated for all partially resistant mutants. For the lowest risk of *de novo* resistance, the prescribed number of compensatory doses ($N_m$) should be taken, uninterrupted, immediately after doses are missed. Otherwise the infected cell population may rebound to a high level, which can make further missed doses very risky for resistance.

## Case study: Combination therapy of daclatasvir and asunaprevir

To demonstrate the practical applicability of our theory, we consider a recently developed interferon-free combination therapy based on an NS5A inhibitor, daclatasvir, and an NS3 protease inhibitor, asunaprevir [38]. In clinical trials, a large proportion of patients infected with HCV genotype-1b achieved sustained virological response (i.e. viral eradication) when treated with daclatasvir and asunaprevir for 24 weeks, although viral breakthrough and viral relapse occurred in a small fraction of patients [45,46].

We first consider patients with the wild-type virus at baseline, i.e. the wild-type virus is the dominant strain before treatment. Using the PK/PD data for each drug [47–49] and the resistance profiles data measured for genotype-1b HCV [50,51], we predicted which mutants are potentially fully-resistant to this combination therapy and calculated the values of $N_m$ and $\Phi_m$ for each of the partially resistant mutants (Fig 4A and 4B) (see Supplementary Materials for more detail). Choosing the highest values of $N_m$ and $\Phi_m$ among all the partially resistant mutants allows us to project the overall risk arising from missed doses over the course of treatment, and we found required numbers of compensatory doses were modest and the risk of *de novo* resistance is low (S2A Fig). To demonstrate that the theoretical approximations represent the full viral dynamics accurately, we simulated a multi-strain viral dynamics model (see Materials and Methods), assuming 1–3 day blocks of consecutive doses are missed randomly within a treatment regimen lasting 24 weeks. The model predicts that relapse of L31M+Y93H or L31W would be observed when overall adherence is less than 90% (Fig 4C and 4D). Indeed, the L31M+Y93H mutant has already been detected in one relapse patient in a clinical trial [46]. There is excellent agreement between simulation results and theoretical predictions (based on Eqs 3 and 4) for the number of cells infected by different mutants after 24 weeks of treatment and the cumulative number of cells infected by partially resistant mutants over the treatment period (Figs 4D and S3).

We then simulated outcomes when the doses are guided by the adaptive treatment strategy (*guided dosing*; see Methods for detailed simulation procedure). Because the risk of *de novo* resistance when doses are missed is low, there is no high-risk period for *de novo* resistance in this case (Fig 4B). If patient dosing is guided, i.e. all the required doses and the extra doses to compensate for the missed doses are taken, the infection can be cleared successfully (Fig 4E). Again, we find excellent agreement between simulation results and theoretical predictions (Fig 4F).

Many patients bear the Y93H mutation at baseline and this mutation reduces the genetic barrier to full resistance by one nucleotide[46]. Our theory suggests that reducing the genetic barrier to full resistance will drastically increase the risk of treatment failure. We repeated our analysis for patients with Y93H at baseline, to test how our adaptive treatment strategy works when the risk of resistance is high. As predicted, many more days of treatment are needed to compensate for missed doses, and the risks of generating full resistance *de novo* are high
<figure><img src="image_10.png" /><figcaption>**Fig 4.** Adaptive treatment strategy improves treatment outcome substantially—a case when the risk of *de novo* resistance is low (with wild-type genotype-1b HCV at baseline). (A) Theoretical prediction of the treatment duration needed to eliminate each partially resistant mutant under perfect adherence (gray bar), and the maximum number of additional doses needed to compensate for missing doses, *N*<sub>m,max</sub> (colored bars). Green, blue and red denote results when 1, 2 and 3 consecutive daily doses are missed, respectively. The longer the colored bars, the greater the impact of missing doses. The symbols next to the bars for *N*<sub>m,max</sub> show the type of mutant investigated in panels (B,D,F). (B) Theoretical prediction of the risk of *de novo* resistance, *Φ*<sub>m</sub>, over time (as shown in Fig 3D), for the three mutants with highest risks of generating fully resistant mutants. The dashed black line shows *Φ*<sub>m</sub> = 0.01. (C) Treatment outcomes 1–3 days of doses are missed randomly. Colored areas denote the fractions of simulations with outcomes of viral relapse without full resistance (light blue) and viral clearance (gray). (D) Comparison between theory predictions and simulations of the number of cells infected by different mutants after 24 weeks of treatment. (E) Treatment outcomes if adaptive treatment strategy is followed. The area above the black dashed line denotes the fraction of patients where virus is not cleared after 24 weeks' treatment. After 24 weeks, patients take the prescribed number of make-up doses without missing further doses. White areas denote adherence levels that are not allowed by the adaptive treatment strategy. (F) Same comparison as in panel (D) for the guided dosing simulation.</figcaption></figure>

doi:10.1371/journal.pcbi.1004040.g004

(>0.01) during the first 3 weeks of effective treatment if 2 consecutive doses are missed (or first 4 weeks if 3 doses are missed; Fig 5A and 5B and S2B Fig). *De novo* full resistance is likely if doses are missed randomly and adherence is less than 90% (dark red area in Fig 5C). The predicted number of infected cells agrees well with simulation, except when adherence is very low
<figure><img src="image_11.png" /><figcaption>**Fig 5.** Adaptive treatment strategy prevents *de novo* resistance and improves treatment outcome substantially—a case when the risk of *de novo* resistance is high. Theoretical prediction and simulation for patients with the Y93H mutant virus (genotype-1b) at baseline under combination therapy of daclatasvir and asunaprevir. Thus, the mutants considered here all have the Y93H mutation. The theoretical predictions and simulation results are plotted in the same way as in Fig 4. Dark red areas in panel (C,E) denote the fraction of patients with *de novo* full resistance to the combination therapy. Note that the fraction of patients with *de novo* resistance in the guided dosing scenario is very small (<0.1%). When doses are guided, so that mutant viral load does not rebound to the pre-treatment level, the theoretical prediction agrees well with simulation as shown in panel (F).</figcaption></figure>

doi:10.1371/journal.pcbi.1004040.g005

such that viral load rebounds back close to the pre-treatment level (Fig 5D and S4–S6 Figs). In stark contrast, when doses are guided, the risk of *de novo* resistance becomes much lower (compare Fig 5C with 5E). Again, for patients who do not clear infection after 24-week treatment, extended periods of treatment as predicted by our theory (using Eq 3) can clear infection with low risk of resistance. The efficacy of the adaptive treatment strategy is robust across different parameter values (S7–S12 Figs and S1 Text). Therefore, our treatment strategy can improve clinical outcomes substantially by adjusting on-going treatment based on patient adherence patterns.

## Discussion

In this study, we integrate PK/PD parameters and viral dynamics into a unified framework to assess the impacts of suboptimal treatment adherence on the risk of treatment failure. Using
simulations incorporating PK/PD and resistance profile data collected previously [48,50,51], we showed that treatment outcomes of combinations therapies of daclatasvir and asunaprevir can be improved by this adaptive treatment strategy, especially when the Y93H mutant is the dominant strain before treatment begins.

We have identified several factors that influence the risk of *de novo* resistance to a combination therapy. Among these factors, the genetic barrier to full resistance plays a dominant role. Thus, for patients with risk of low adherence, combinations that impose a higher genetic barrier are recommended. This is especially important in resource-limited settings where patients have limited access to health care and adherence is not closely monitored. The recently developed HCV entry inhibitors [52], which inhibit host factors that are required for viral entry (instead of viral factors), may offer a promising direction for HCV combination therapy, because of their high genetic barriers to resistance, and their synergistic interactions with other classes of DAAs. For situations where therapies with low genetic barriers to resistance are used, we have identified a high-risk window period during which *de novo* resistance is likely if doses are missed. Intervention efforts should focus on enhancing patients’ adherence during this period. Additional complementary strategies could further reduce the risk of treatment failure. First, if doses are missed during the high-risk window, the immediate addition of another drug with a different mechanism of action from existing drugs may eliminate any low level of fully resistant mutants that has arisen. Alternatively, a patient could be treated preemptively using additional drugs during the entire high-risk period and switched to fewer drugs afterwards. Another important factor is the number of consecutively missed doses as shown previously [36]. consecutively missed doses lead to exponential growth of ‘partially resistant’ mutants, and thus substantially increase the risk of *de novo* resistance.

Our theory also predicts the number of compensatory doses ($N_m$) needed to compensate for missed doses, in order to eliminate preexisting mutants. Interestingly, clinical trials have shown that adherence levels tend to decrease over time [19,31]; we show that more doses are needed to compensate for missed doses that occur later in treatment because of the rebound of target cells. While many previous studies have focused on average adherence [18,19,29–31,36], we emphasize that the timing of the missed doses is also a critical determinant of treatment outcome and the risk of resistance.

There exist substantial heterogeneities among patients owing to variation in HCV genotypes, patient viral loads, death rates of infected cells [40,53] and effectiveness of drug penetration [47]. Our analysis has identified several factors that influence the impact of suboptimal adherence, particularly the rebound rate of target cells under treatment, the half-life of infected cells and the overall viral fitness, $R_0$. We used the best available estimates of these parameters, but further empirical work is needed. If resistance profiles and viral parameters could be measured directly from a specific patient, then our framework linking these factors could be tailored to give patient-specific guidelines.

Certain model assumptions reflect uncertainties in our current knowledge of HCV infection. First, our prediction about time to viral extinction should be treated cautiously. We predict the time of extinction (as in other models [54–56]) by assuming that infected cells decline at a rate set by their death rate, and infection is cleared when the number of infected cells is below one. However, factors such as pressures from the immune system and infections in different tissue compartments may influence the extinction threshold. Furthermore, if DAA treatment causes intracellular viral RNA to decay with negligible replication [57], the decline of infected cells may result from a combination of cell recovery and death of infected cells. Indeed, sustained virological response has been observed in clinical trials of DAA combination therapies with shorter durations of treatment [5,6]. Our model can be adjusted easily once the decay dynamics of infected cells are understood better. Second, our model captures the main features
of pharmacodynamics and viral dynamics by assuming quasi-equilibrium for viral populations and drug penetration into liver cells. Further work that incorporates detailed intracellular interactions [57] and different body compartments may improve model accuracy, once pertinent parameters are measured. However, a more detailed model may become analytically intractable.

This quantitative framework is a step towards developing a tool (for example, see Ref. [58]) for clinicians to design combination therapies and adaptively manage treatment regimens to achieve favorable clinical outcomes. It highlights the importance of characterizing resistance profiles of HCV, assessing readiness for treatment, and monitoring adherence patterns during treatment, so that treatment can be designed and adjusted in an evidence-based manner. This framework can be adapted easily to combination therapies based on interferon, entry inhibitors [52] or other DAA candidates, or treatments of other curable diseases without a latent reservoir.

# Materials and Methods

## HCV model and viral fitness in the presence of drug, $R_{\text{eff}}(t)$

To analyze the dynamics of the virus, we constructed an ordinary differential equation (ODE) model to describe the long-term within-host dynamics of a single HCV strain under drug treatment, based on an established model developed by Neumann et al. [53] (see Supplementary Material). In the model, $\varepsilon$ represents the proportion by which the therapy reduces viral growth ($\varepsilon$ is in the range of 0 and 1). Then, the fitness of the virus, $R_{\text{eff}}(t)$, is the product of the complement of the therapy's efficacy (1- $\varepsilon(\tau)$), the reproductive number of the virus, $R_0$, and the availability of target cells, $h(t)$ (Eq 1).

## Average effective viral fitness when m doses are missed, $R_{\text{ave},m}$

To approximate the time-varying viral fitness, $R_{\text{eff}}(t)$, during the period when $m$ consecutive doses are missed, we assume that the abundance of target cells stays constant. This is a good approximation, because the length of the period when consecutive doses are missed tends to be short compared to the time scale of target cell rebound. Then the only time-varying quantity in Eq 1 is $\varepsilon(\tau)$. We can calculate the average level of drug inhibition during the period when $m$ doses are missed, $\varepsilon_{\text{ave},m}$, by incorporating parameters for pharmacokinetics and pharmacodynamics (for example, see Wahl and Nowak [36]). Then the time-average effective reproductive number, $R_{\text{ave},m}(t)$, for a mutant when $m$ consecutive doses are missed starting at time $t$ can be expressed as Eq 2. In practice, because the precise number of target cells at time $t$ is hard to estimate, we can approximate $R_{\text{ave},m}(t)$ by setting $h(t) = 1$, and then $R_{\text{ave},m}(t)$ becomes $R_{\text{ave},m}(t) \approx (1 - \varepsilon_{\text{ave},m}) \cdot R_0$. Because $h(t) \le 1$, this always overestimates the viral fitness and thus is a conservative estimate in terms of guiding treatment. Note that the assumption that $h(t) = 1$ is valid only when the viral load at time $t$ is much lower than it was before treatment, which is the case if adherence is not too low. Otherwise, $h(t)$ would decrease significantly due to large amount of infection.

## The number of compensatory doses needed ($N_m$)

To calculate $N_m$ for each mutant, we make the simplifying assumption that the dynamics of the viral populations are at quasi-equilibrium, because changes in the viral populations occur much faster than changes in infected hepatocytes. Then, the dynamics of the number of cells
infected by mutant viruses, $I(t)$, are described by:

$$ \frac{dI(t)}{dt} = (R_{\text{eff}}(t) - 1) \cdot \delta \cdot I(t) \qquad (5) $$

where $\delta$ is the death rate of infected hepatocytes. If we approximate $R_{\text{eff}}(t)$ using the constant $R_{\text{ave},m}$ for the period when doses are missed, Eq 5 can be solved analytically. Then, the number of infected cells after missing $m$ consecutive doses starting at time $t_0$ can be expressed as:

$$ I(t_0 + m + T) \approx I(t_0) \cdot \exp((R_{\text{ave},m}(t_0) - 1) \cdot \delta \cdot m \cdot T) \qquad (6) $$

We now consider the situation when $m$ consecutive doses are missed, and ask how many uninterrupted doses (compensatory doses) must be taken so that the number of cells infected by the mutant is suppressed to a same number as if the $m$ doses had not been missed. We first calculate the number of infected cells if the $m$ consecutive doses are taken, i.e. if no doses is missed:

$$ I_{\text{optimal}}(t_0 + m \cdot T) \approx I(t_0) \cdot \exp((R_{\text{ave},0}(t_0) - 1) \cdot \delta \cdot m \cdot T) \qquad (7) $$

where $I(t_0)$ is the number of cells infected by the mutant at time $t_0$, $R_{\text{ave},0}$ is the average effective reproductive number of the mutant when all doses are taken, and $T$ is the scheduled interval between doses.

We then analyze the situation where a patient skips $m$ consecutive doses, starting at time $t_0$, and then takes $N_m$ compensatory doses immediately afterwards. In this case, assuming the number of target cells does not change much during this period, we can approximate the number of cells infected by the mutant at the end of the $N_m$ doses as:

$$ I_{\text{suboptimal}}(t_0 + m \cdot T + N_m \cdot T) \approx I(t_0) \cdot \exp((R_{\text{ave},m}(t_0) - 1) \cdot \delta \cdot m \cdot T) \cdot I(t_0) \cdot \exp((R_{\text{ave},0}(t_0) - 1) \cdot \delta \cdot N_m \cdot T) \qquad (8) $$

By equating the right hand sides of Eqs 7 and 8 and solving the equation, we derive the expression for $N_m$:

$$ N_m(t_0) \approx \frac{R_{\text{ave},m}(t_0) - R_{\text{ave},0}(t_0)}{1 - R_{\text{ave},0}(t_0)} \cdot m \qquad (9) $$

For potent therapies, usually $R_{\text{ave},0}(t_0) \approx 0$. Then we get Eq 3.

In the derivation above, we have assumed that the target cell abundance stays constant during the period under consideration. This would be a good approximation if only a few days of doses are missed or if the target cell has already rebounded to the infection-free level. If the abundance of target cells changes considerably during the period under consideration, an alternative, conservative approach would be to assume $h(t) = 1$ and take $N_{m,max}(t_0) \approx m \cdot (1 - \varepsilon_{\text{ave},m}) \cdot R_0$ compensatory doses after missing $m$ consecutive doses of treatment.

## The number of doses to eradicate a mutant ($N_{\text{erad}}$) and the number of cells infected by a mutant ($I(t)$)

One important application of $N_m$ is to predict the number of remaining doses needed to eradicate a mutant, $N_{\text{erad}}$, in a patient during treatment. This number can be calculated as follows. If adherence is perfect, the number of infected cells declines exponentially at a rate set approximately by the death rate of infected cells, $\delta: (t) \approx I_0 \cdot \exp(-\delta \cdot t)$, where $I_0$ is the number of cells infected by a mutant of interest before treatment. If we assume that a mutant goes extinct if the expected number of infected cells in a patient goes below 1, the number of doses needed to
eradicate a mutant before treatment (assuming adherence is perfect), $N_{erad,0}$, is calculated as:
$$N_{erad,0} \approx \frac{\log(I_0)}{\delta \cdot T}$$

When doses are missed during treatment, it is clear from the calculation of $N_m$ above that $N_m-m$ extra doses of treatment are needed to eradicate the virus. Therefore, if a patient has taken a total of $x$ doses and has had $k$ instances of missing doses before time $t$, with $m_i$ days of doses missed in the $i^{th}$ instance ($i = 1,2,..,k$), then the number of remaining doses needed to eradicate the mutant is calculated as:

$$N_{erad} = N_{erad,0} - x + \sum_{i=1}^{k} (N_{m,i} - m_i) \quad (10)$$

We can use Eq 10 to predict the number of cells infected by a mutant as: $I(t) \approx \exp(\delta \cdot N_{erad}(t) \cdot T)$. In our model, and a patient is cleared of infection when all mutants are driven to extinction. The accuracy of this approximation is shown in Figs 4D and 4F and 5D and 5F.

## The risk of full resistance if doses are missed ($\Phi_m$)

To calculate the risk of full resistance during the period when $m$ doses are missed, we first calculate the number of cells newly infected by a partially resistant mutant when $m$ doses are missed, $\Omega_m(t)$. Again, we use $R_{ave,m}(t)$ to approximate $R_{eff}(t)$, the total number of cells infected by the mutant virus, starting at time $t$. $\Omega_m(t)$ can be expressed as an integration of new infections during the period of missing doses (according to Eq 5):

$$\Omega_m(t) \approx R_{ave,m}(t) \cdot \delta \cdot \int_{t}^{t+m \cdot T} I(x) dx = I(t) \cdot \frac{R_{ave,m}(t)}{R_{ave,m}(t) - 1} \cdot (e^{(R_{ave,m}(t)-1) \cdot \delta \cdot m \cdot T} - 1) \quad (11)$$

The expected number of target cells that become infected by fully resistant mutant viruses, $\Phi_m$, is a product of the effective mutation rate from the partially resistant mutant to the fully resistant mutant ($\mu_{eff}$) and the total number of cells infected by the partially resistant mutant ($\Omega_m$): $\Phi_m(t) = \mu_{eff} \cdot \Omega_m(t)$, as shown in Eq 4.

Note that we track the population of newly infected cells to assess the risk of *de novo* generation of full resistance. This assumes implicitly that the fully resistant mutant is selected only when it enters a cell. This is a good assumption for DAAs that act on intracellular stages of the viral life-cycle, such as viral genome replication or assembly. However, in situations where the drug blocks viral entry into the cell, the mutant virus may have a selective advantage for entering a cell. Then the viral population should be tracked instead, but the results presented here still can be applied to drugs that block cell entry by multiplying with a simple scaling factor [59].

## Stochastic-deterministic hybrid simulation of multiple strains of HCV

We constructed a simulation model considering the dynamics of the baseline virus and all the potentially partially resistant mutants (see Supplementary Material). This simulation model follows a hybrid approach used previously to simulate the evolutionary dynamics of HIV [60]. It considers the dynamics of multiple strains of HCV deterministically (using ODEs) while treating the extinction and mutation processes as stochastic events (see Supplementary Material for detail).

In the simulation, a patient is treated for a total period of 24 weeks. We generate two types of dosing patterns: *random dosing* and *guided dosing*. For the *random dosing* pattern, doses are missed in blocks of 1–3 days at times chosen randomly with equal probability during the treatment period. This probability is set as a constant in each run, but varied across runs such that different overall levels of adherence are generated. In each simulation, we assume that at least
one-day treatment is taken immediately after each dose-skipping event (i.e. 1, 2 or 3 consecutive missed doses), to ensure that two dose-skipping events do not occur consecutively (otherwise, longer blocks of doses would be missed than intended). For *guided dosing*, we ensure that doses are always taken during the high-risk window period predicted by our theory. After this high-risk window period, we set a constant probability of missing doses in blocks of 1–3 days. Immediately after a block of doses is missed, we ensure a sufficient number of uninterrupted doses (calculated as $N_m$) are always taken. If the virus is not eradicated after the 24-week treatment period, the patient is treated with an uninterrupted number of doses as predicted by our theory. The outcome of the simulation at the end of the procedure is reported.

**Supporting Information**

**S1 Text. Supplementary materials.**

(DOCX)

**S1 Fig.** Sensitivity analysis of the risk of *de novo* resistance to variations of key parameter values, $R_{0,mut}$ (panels A,B), $\delta$ (panel C) and the rate of recovery of target cells upon treatment (panel D). In each panel, the trajectories show how the risk of *de novo* resistance ($\text{Log}_10\Phi_m$) changes over time if adherence is perfect. Figures are plotted using the same parameter settings as trajectories ‘a’ in Fig 3 in the main text, except that $R_{0,mut} = 5$ in panel A, $R_{0,mut} = 15$ in panel B, $\delta = 0.5$ in panel C and $\alpha = 1.95 \times 10^4$, $d = 0.015$ in panel D. In the main results, i.e. Fig 3, the parameter values used are $R_{0,mut} = 10$, $\delta = 0.15$, $\lambda = 1.95 \times 10^5$, $d = 0.15$. (TIFF)

**S2 Fig.** The predicted number of additional days of doses needed to compensate for the first instance of missing 1 (green lines), 2 (blue lines) or 3 (red lines) consecutive days of doses (maximum $N_m$ for all partially resistant mutants), and the high-risk window period of *de novo* resistance (shaded area: $\Phi_m > 0.01$ for any of the partially resistant mutants as shown in Fig 4B). (A) Predictions for patients with the wild-type virus at baseline before treatment. The areas below the curves are white, indicating that the risk of *de novo* resistance is always low. (B) Predictions for patients with the Y93H mutant virus at baseline. The initial increases of the number of days of compensating doses are due to the increase of the number of target cells upon treatment, and the sudden drops during later periods of treatment are due to the elimination of particular partially resistant mutant lineages. Note that these curves are calculated under the assumption that adherence is perfect except for the 1–3 days of missed doses being considered, i.e. it is a prediction for the first instance of missed doses. For cases where multiple instances of missed doses have occurred, one needs to calculate the values of $N_m$ and $\Phi_m$ for each mutant based on the adherence pattern, and then integrate them together by choosing the highest values of $N_m$ and $\Phi_m$ for those mutants. (TIFF)

**S3 Fig.** Theory correctly predicts the number of cells infected by two mutants (L31V and L31M+Y93H) generated in the hybrid model simulation when doses are missed randomly (panels A,B) or guided by adaptive treatment theory (panels C,D) for patients with the wild-type virus at baseline. L31V and L31M+Y93H are the two most likely mutants that generate full resistance. The axes are the theory prediction (x-axis) and model simulation (y-axis) of the $\text{Log}_10$ of the number of mutants, which are calculated as the cumulative numbers of $\text{Log}_10 \Phi_m(t)/\mu_{mut}$ for all missed doses. (TIFF)
S4 Fig. The theory correctly predicts the number of cells infected by mutant viruses at the end of 24-weeks' treatment in the hybrid model simulation when adherence is greater than 70% (vertical dashed lines) and doses are missed randomly, for patients with the Y93H mutant virus at baseline. The y-axis shows the log₁₀ difference between the theory prediction and the model simulation at the end of the 24-weeks' treatment. Note that when adherence is lower than 70%, the population of infected cells grows to high levels close to the pre-treatment level, where further growth is curtailed by target cell limitation. As a result, the theoretical prediction overestimates the number of cells infected by the mutant virus significantly because we assume the number of target cells is not limited.

(TIFF)

S5 Fig. Comparison between theory prediction and simulation of the numbers of cells infected by three mutants (L31M+Y93H, L31V+Y93H, L31V+Q54H+Y93H) generated when doses are missed randomly, for patients with the Y93H mutant virus at baseline. (A,B,C)

The axes are the theory prediction (x-axis) and model simulation (y-axis) of the Log₁₀ of the number of cells infected by different mutants, which are calculated as the cumulative numbers of Log₁₀ Φₘ/μₘₐt for all missed doses (L31M+Y93H in panel A; L31V+Y93H in panel B; L31V+Q54H+Y93H in panel C). (D,E,F) Our theory prediction is accurate for adherence greater than 70%, but overestimates the number of cells infected by the mutant virus significantly when adherence is lower than 70%, for the same reason as explained in the legend of S4 Fig.

(TIFF)

S6 Fig. Theory correctly predicts the number of cells infected by three mutants (L31M +Y93H, L31V+Y93H, L31V+Q54H+Y93H) generated in the hybrid model simulation when doses are guided by adaptive treatment theory, for patients with the Y93H virus at baseline. (A,B,C) The axes are the theory prediction (x-axis) and model simulation (y-axis) of the Log₁₀ of the number of mutants (L31M+Y93H in panel A; L31V+Y93H in panel B; L31V+Q54H+Y93H in panel C), which are calculated as the cumulative numbers of Log₁₀ Φₘ(t)/μₘₐt for all missed doses. (D,E,F) the Log₁₀ differences between theory prediction and model simulation as shown in panels (A,B,C). Note that our theory agrees very well for mut-ants L31M+Y93H and L31V+Y93H. For mutant L31V+Q54H+Y93H, the stochastic extinc-tion and appearance of this mutant generates stochastic deviations of the simulation from theory prediction.

(TIFF)

S7 Fig. The impact of lower viral fitness ($R_0 = 5$) on treatment outcomes and adaptive treatment strategies of combination therapy with daclatasvir and asunaprevir, with the wild-type virus at baseline. Panels A-F show the same plots as Fig 4 in the main text, except that the fitness parameter $R_0$ for the wild-type virus is assumed to be 5. The treatment outcome im-proves for all scenarios for this lower viral fitness (compare with Fig 4 in the main text). Using the adaptive treatment strategy prevents viral relapse and *de novo* resistance if overall adher-ence is greater than 60% (panel E). Panel F is empty because all patients are cleared of infection after 24 weeks.

(TIFF)

S8 Fig. The impact of lower viral fitness ($R_0 = 5$) on treatment outcomes and adaptive treatment strategies of combination therapy with daclatasvir and asunaprevir, with the Y93H virus at baseline. Panels A-F show the same plots as Fig 5 in the main text, except that the fit-ness parameter $R_0$ for the wild-type virus is assumed to be 5. The treatment outcome improves for all scenarios for this lower viral fitness (compare with Fig 5 in the main text). Using adap-tive treatment strategy reduced the risk of *de novo* resistance (panels E). Our theory correctly
predicts the number of infected cells in a patient at the end of 24 weeks' treatment.
(TIFF)

**S9 Fig.** The impact of higher viral fitness ($R_0 = 15$) on treatment outcomes and adaptive treatment strategies of combination therapy with daclatasvir and asunaprevir, with the wild-type virus at baseline. Panels A-F show the same plots as Fig 4 in the main text, except that the fitness parameter $R_0$ for the wild-type virus is assumed to be 15. The treatment outcome improves for all scenarios for this lower viral fitness (compare with Fig 4 in the main text). Our theory correctly predicts the number of infected cells in a patient at the end of 24 weeks' treatment.
(TIFF)

**S10 Fig.** The impact of higher viral fitness ($R_0 = 15$) on the treatment outcomes and adaptive treatment strategies of combination therapy of daclatasvir and asunaprevir with the Y93H virus at baseline. Panels A-F show the same plots with Fig 5 except that in the analytical derivation and model simulation, the fitness parameter $R_0$ for the Y93H mutant virus is assumed to be 15. The risks of viral relapse and *de novo* resistance become higher when the viral fitness, $R_0$, is higher. Using adaptive treatment strategy can prevent *de novo* resistance and improve treatment outcomes (panels E and F). Our theory correctly predicts the number of infected cells in a patient at the end of 24 weeks' treatment when doses are guided (panels F). Our theory does not predict the number of infected cells at the end of treatment well, when doses are missed randomly and the adherence is low. This is because, when adherence is low, the viral load often rebounds back to the pre-treatment level, where it is limited by target cell availability. This phenomenon is not included in our theory, which overestimates the number of viruses as a result.
(TIFF)

**S11 Fig.** The impact of higher viral clearance rate ($\delta = 0.5$) on the treatment outcomes and adaptive treatment strategies of combination therapy of daclatasvir and asunaprevir with the wild-type genotype 1b virus at baseline. Panels A-F show the same plots with Fig 4 except that in the analytical derivation and model simulation, the viral clearance rate, $\delta$, is assumed to be 0.5 instead of 0.15 in Fig 4 (but note that $R_0$ for the viruses is kept the same). When the viral clearance rate increases, it takes less time to eradicate the virus from a patient. However, when doses are missed, the population of mutant viruses expands more quickly, because the half-life of the infected cells is shorter and thus it undergoes a higher number of replication generations during the period of missed doses. Using the adaptive treatment strategy can prevent viral relapse and *de novo* resistance and improve treatment outcome (panels E). Our theory correctly predicts the number of infected cells in a patient at the end of 24 weeks' treatment when doses are guided (panel D).
(TIFF)

**S12 Fig.** The impact of higher viral clearance rate ($\delta = 0.5$) on the treatment outcomes and adaptive treatment strategies of combination therapy of daclatasvir and asunaprevir with the Y93H virus at baseline. Panels A-F show the same plots with Fig 5 except that in the analytical derivation and model simulation, we assume the viral clearance rate, $\delta$, is 0.5 instead of 0.15 (but note that $R_0$ for the viruses is kept the same). As seen in S11 Fig for the scenario with the wild-type virus at baseline, it takes less time to eradicate the virus from a patient for this higher viral clearance rate. However, when doses are missed, the population of mutant viruses expands more quickly, increasing the risk of viral relapse and *de novo* resistance. Using adaptive treatment strategy can prevent viral relapse, *de novo* resistance and improve treatment outcome (panels E and F). Our theory does not predict the number of infected cells at the end of
treatment well, when doses are missed randomly and adherence is low. This is because, during the time period when doses are missed, the rebound of the viruses is quicker when $\delta$ is higher (because the viral generation time is shorter). When adherence is low, the viral load often re-bounds back to the pre-treatment level, where it is limited by target cell availability. This phenomenon is not included in our theory, which overestimates the number of viruses as a result. (TIFF)

## Acknowledgments

We thank Alan Perelson, Paul Fenimore and members in the Lloyd-Smith group and the Sun Group for helpful discussions.

## Author Contributions

Conceived and designed the experiments: RK RS JOLS. Analyzed the data: RK CL. Wrote the paper: RK CL HQ RS JOLS.

## References

1. Thomas DL (2013) Global control of hepatitis C: where challenge meets opportunity. Nat Med 19: 850–858. doi: 10.1038/nm.3184 PMID: 23836235

2. Lavanchy D (2009) The global burden of hepatitis C. Liver Int 29 Suppl 1: 74–81. doi: 10.1111/j.1478-3231.2008.01934.x PMID: 19207969

3. Scheel TK, Rice CM (2013) Understanding the hepatitis C virus life cycle paves the way for highly effective therapies. Nat Med 19: 837–849. doi: 10.1038/nm.3248 PMID: 23836234

4. Gane E (2011) Future hepatitis C virus treatment: interferon-sparing combinations. Liver Int 31 Suppl 1: 62–67. doi: 10.1111/j.1478-3231.2010.02383.x PMID: 21205140

5. Zeuzem S, Jacobson IM, Baykal T, Marinho RT, Poordad F, et al. (2014) Retreatment of HCV with ABT-450/r-ombitasvir and dasabuvir with ribavirin. N Engl J Med 370: 1604–1614. doi: 10.1056/NEJMoa1401561 PMID: 24720679

6. Kowdley KV, Gordon SC, Reddy KR, Rossaro L, Bernstein DE, et al. (2014) Ledipasvir and sofosbuvir for 8 or 12 weeks for chronic HCV without cirrhosis. N Engl J Med 370: 1879–1888. doi: 10.1056/NEJMoa1402355 PMID: 24720702

7. Feld JJ, Kowdley KV, Coakley E, Sigal S, Nelson DR, et al. (2014) Treatment of HCV with ABT-450/r-ombitasvir and dasabuvir with ribavirin. N Engl J Med 370: 1594–1603. doi: 10.1056/NEJMoa1315722 PMID: 24720703

8. Afdhal N, Zeuzem S, Kwo P, Chojkier M, Gitlin N, et al. (2014) Ledipasvir and sofosbuvir for untreated HCV genotype 1 infection. N Engl J Med 370: 1889–1898. doi: 10.1056/NEJMoa1402454 PMID: 24725239

9. Pol S, Ghalib RH, Rustgi VK, Martorell C, Everson GT, et al. (2012) Daclatasvir for previously untreated chronic hepatitis C genotype-1 infection: a randomised, parallel-group, double-blind, placebo-controlled, dose-finding, phase 2a trial. Lancet Infect Dis 12: 671–677. doi: 10.1016/S1473-3099(12)70138-X PMID: 22714001

10. Lok AS, Gardiner DF, Lawitz E, Martorell C, Everson GT, et al. (2012) Preliminary study of two antiviral agents for hepatitis C genotype 1. N Engl J Med 366: 216–224. doi: 10.1056/NEJMoa1104430 PMID: 22256805

11. Qi H, Olson CA, Wu NC, Ke R, Loverdo C, et al. (2014) A quantitative high-resolution genetic profile rapidly identifies sequence determinants of hepatitis C viral fitness and drug sensitivity. PLoS Pathog 10: e1004064. doi: 10.1371/journal.ppat.1004064 PMID: 24722365

12. Robinson M, Tian Y, Delaney WEt, Greenstein AE (2011) Preexisting drug-resistance mutations reveal unique barriers to resistance for distinct antivirals. Proc Natl Acad Sci U S A 108: 10290–10295. doi: 10.1073/pnas.1101515108 PMID: 21646519

13. Pawlotsky JM (2011) Treatment failure and resistance with direct-acting antiviral drugs against hepatitis C virus. Hepatology 53: 1742–1751. doi: 10.1002/hep.24262 PMID: 21374691

14. Ribeiro RM, Li H, Wang S, Stoddard MB, Learn GH, et al. (2012) Quantifying the diversification of hepatitis C virus (HCV) during primary infection: estimates of the in vivo mutation rate. PLoS Pathog 8: e1002881. doi: 10.1371/journal.ppat.1002881 PMID: 22927817
15. Powdrill MH, Tchesnokov EP, Kozak RA, Russell RS, Martin R, et al. (2011) Contribution of a mutational bias in hepatitis C virus replication to the genetic barrier in the development of drug resistance. Proc Natl Acad Sci U S A 108: 20509–20513. doi: 10.1073/pnas.1105797108 PMID: 22135458

16. Simmonds P (2004) Genetic diversity and evolution of hepatitis C virus—15 years on. J Gen Virol 85: 3173–3188. PMID: 15483230

17. Pybus OG, Charleston MA, Gupta S, Rambaut A, Holmes EC, et al. (2001) The epidemic behavior of the hepatitis C virus. Science 292: 2323–2325. PMID: 11423661

18. Lieveld FI, van Vlerken LG, Siersema PD, van Erpecum KJ (2013) Patient adherence to antiviral treatment for chronic hepatitis B and C: a systematic review. Ann Hepatol 12: 380–391. PMID: 23619254

19. Lo Re V 3rd, Teal V, Localio AR, Amorosa VK, Kaplan DE, et al. (2011) Relationship between adherence to hepatitis C virus therapy and virologic outcomes: a cohort study. Ann Intern Med 155: 353–360. doi: 10.7326/0003-4819-155-6-201109200-00003 PMID: 21930852

20. Sarrazin C, Zeuzem S (2010) Resistance to direct antiviral agents in patients with hepatitis C virus infection. Gastroenterology 138: 447–462. doi: 10.1053/j.gastro.2009.11.055 PMID: 20006612

21. Mathes T, Antoine SL, Pieper D (2014) Factors influencing adherence in Hepatitis-C infected patients: a systematic review. Bmc Infectious Diseases 14.

22. Chung RT, Baumert TF (2014) Curing chronic hepatitis C—the arc of a medical triumph. N Engl J Med 370: 1576–1578. doi: 10.1056/NEJMp1400986 PMID: 24720678

23. Ford N, Singh K, Cooke GS, Mills EJ, von Schoen-Angerer T, et al. (2012) Expanding access to treatment for hepatitis C in resource-limited settings: lessons from HIV/AIDS. Clin Infect Dis 54: 1465–1472. doi: 10.1093/cid/cis227 PMID: 22431808

24. Smith RJ, Okano JT, Kahn JS, Bodine EN, Blower S (2010) Evolutionary dynamics of complex networks of HIV drug-resistant strains: the case of San Francisco. Science 327: 697–701. doi: 10.1126/science.1180556 PMID: 20075214

25. Roberts RR, Hota B, Ahmad I, Scott RD 2nd, Foster SD, et al. (2009) Hospital and societal costs of antimicrobial-resistant infections in a Chicago teaching hospital: implications for antibiotic stewardship. Clin Infect Dis 49: 1175–1184. doi: 10.1086/605630 PMID: 19739972

26. DiMasi JA, Hansen RW, Grabowski HG (2003) The price of innovation: new estimates of drug development costs. J Health Econ 22: 151–185. PMID: 12606142

27. Bangsberg DR, Perry S, Charlebois ED, Clark RA, Roberston M, et al. (2001) Non-adherence to highly active antiretroviral therapy predicts progression to AIDS. AIDS 15: 1181–1183. PMID: 11416722

28. Paterson DL, Swindells S, Mohr J, Brester M, Vergis EN, et al. (2000) Adherence to protease inhibitor therapy and outcomes in patients with HIV infection. Ann Intern Med 133: 21–30. PMID: 10877736

29. Gordon SC, Yoshida EM, Lawitz EJ, Bacon BR, Sulkowski MS, et al. (2013) Adherence to assigned dosing regimen and sustained virological response among chronic hepatitis C genotype 1 patients treated with boceprevir plus peginterferon alfa-2b/ribavirin. Aliment Pharmacol Ther 38: 16–27. doi: 10.1111/apt.12342 PMID: 23710734

30. Evon DM, Esserman DA, Bonner JE, Rao T, Fried MW, et al. (2013) Adherence to PEG/ribavirin treatment for chronic hepatitis C: prevalence, patterns, and predictors of missed doses and nonpersistence. J Viral Hepat 20: 536–549. doi: 10.1111/jvh.12079 PMID: 23808992

31. Weiss JJ, Brau N, Stivala A, Swan T, Fishbein D (2009) Review article: adherence to medication for chronic hepatitis C—building on the model of human immunodeficiency virus antiretroviral adherence research. Aliment Pharmacol Ther 30: 14–27. doi: 10.1111/j.1365-2036.2009.04004.x PMID: 19416131

32. zur Wiesch PA, Kouyos R, Engelstadter J, Regoes RR, Bonhoeffer S (2011) Population biological principles of drug-resistance evolution in infectious diseases. Lancet Infect Dis 11: 236–247. doi: 10.1016/S1473-3099(10)70264-4 PMID: 21371657

33. Gelman MA, Glenn JS (2010) Mixing the right hepatitis C inhibitor cocktail. Trends Mol Med.

34. Spellberg B, Guidos R, Gilbert D, Bradley J, Boucher HW, et al. (2008) The epidemic of antibiotic-resistant infections: a call to action for the medical community from the Infectious Diseases Society of America. Clin Infect Dis 46: 155–164. doi: 10.1086/524891 PMID: 18171244

35. Rosenbloom DI, Hill AL, Rabi SA, Siliciano RF, Nowak MA (2012) Antiretroviral dynamics determines HIV evolution and predicts therapy outcome. Nat Med 18: 1378–1385. PMID: 22941277

36. Wahl LM, Nowak MA (2000) Adherence and drug resistance: predictions for therapy outcome. Proc Biol Sci 267: 835–843. PMID: 10819155

37. Soriano V, Perelson AS, Zoulim F (2008) Why are there different dynamics in the selection of drug resistance in HIV and hepatitis B and C viruses? J Antimicrob Chemother 62: 1–4. doi: 10.1093/jac/dkn175 PMID: 18441341

38. Suzuki Y, Ikeda K, Suzuki F, Toyota J, Karino Y, et al. (2012) Dual Oral Therapy with Daclatasvir and Asunaprevir for Patients with HCV Genotype 1b Infection and Limited Treatment Options. J Hepatol.
39. Ribeiro RM, Bonhoeffer S (2000) Production of resistant HIV mutants during antiretroviral therapy. Proc Natl Acad Sci U S A 97: 7681-7686. PMID: 10884399

40. Rong L, Dahari H, Ribeiro RM, Perelson AS (2010) Rapid emergence of protease inhibitor resistance in hepatitis C virus. Sci Transl Med 2: 30ra32. doi: 10.1126/scitranslmed.3000544 PMID: 20445200

41. Loverdo C, Lloyd-Smith JO (2013) Intergenerational phenotypic mixing in viral evolution. Evolution 67: 1815-1822. doi: 10.1111/evo.12048 PMID: 23730772

42. Loverdo C, Park M, Schreiber SJ, Lloyd-Smith JO (2012) Influence of viral replication mechanisms on within-host evolutionary dynamics. Evolution 66: 3462-3471. doi: 10.1111/j.1558-5646.2012.01687.x PMID: 23106710

43. Alexander HK, Bonhoeffer S (2012) Pre-existence and emergence of drug resistance in a generalized model of intra-host viral dynamics. Epidemics 4: 187-202. doi: 10.1016/j.epidem.2012.10.001 PMID: 23351371

44. Nordqvist O, Sodergard B, Tully MP, Sonnerborg A, Lindblad AK (2006) Assessing and achieving readiness to initiate HIV medication. Patient Educ Couns 62: 21-30. PMID: 16766244

45. Kosaka K, Imamura M, Hayes CN, Abe H, Hiraga N, et al. (2014) Emergence of resistant variants detected by ultra-deep sequencing after asunaprevir and daclatasvir combination therapy in patients infected with hepatitis C virus genotype 1. J Viral Hepat.

46. Karino Y, Toyota J, Ikeda K, Suzuki F, Chayama K, et al. (2013) Characterization of virologic escape in hepatitis C virus genotype 1b patients treated with the direct-acting antivirals daclatasvir and asunaprevir. J Hepatol 58: 646-654. doi: 10.1016/j.jhep.2012.11.012 PMID: 23178977

47. Ke R, Loverdo C, Qi H, Olson CA, Wu NC, et al. (2014) Modelling clinical data shows active tissue concentration of daclatasvir is 10-fold lower than its plasma concentration. J Antimicrob Chemother 69: 724-727. doi: 10.1093/jac/dkt423 PMID: 24169581

48. Nettles RE, Gao M, Bifano M, Chung E, Persson A, et al. (2011) Multiple ascending dose study of BMS-790052, a nonstructural protein 5A replication complex inhibitor, in patients infected with hepatitis C virus genotype 1. Hepatology 54: 1956-1965. doi: 10.1002/hep.24609 PMID: 21837752

49. Eley T, Pasquinelli C, Wendelburg P, Villegas C, He B, et al. (2010) Safety and Pharmacokinetics of BMS-650032 Following Multiple Ascending Doses for 14 Days in Healthy Subjects. ICAAC. Boston.

50. Fridell RA, Qiu D, Wang C, Valera L, Gao M (2010) Resistance analysis of the hepatitis C virus NS5A inhibitor BMS-790052 in an in vitro replicon system. Antimicrob Agents Chemother 54: 3641-3650. doi: 10.1128/AAC.00556-10 PMID: 20585111

51. Fridell RA, Wang C, Sun JH, O'Boyle DR 2nd, Nower P, et al. (2011) Genotypic and phenotypic analysis of variants resistant to hepatitis C virus nonstructural protein 5A replication complex inhibitor BMS-790052 in humans: in vitro and in vivo correlations. Hepatology 54: 1924-1935. doi: 10.1002/hep.24594 PMID: 21809362

52. Xiao F, Fofana I, Thumann C, Mailly L, Alles R, et al. (2014) Synergy of entry inhibitors with direct-acting antivirals uncovers novel combinations for prevention and treatment of hepatitis C. Gut.

53. Neumann AU, Lam NP, Dahari H, Gretch DR, Wiley TE, et al. (1998) Hepatitis C viral dynamics in vivo and the antiviral efficacy of interferon-alpha therapy. Science 282: 103-107. PMID: 9756471

54. Dahari H, Shteingart S, Gafanovich I, Cotler SJ, D'Amato M, et al. (2014) Sustained virological response with intravenous silibinin: individualized IFN-free therapy via real-time modelling of HCV kinetics. Liver Int.

55. Guedj J, Perelson AS (2011) Second-phase hepatitis C virus RNA decline during telaprevir-based therapy increases with drug effectiveness: implications for treatment duration. Hepatology 53: 1801-1808. doi: 10.1002/hep.24272 PMID: 21384401

56. Snoeck E, Chanu P, Lavielle M, Jacqmin P, Jonsson EN, et al. (2010) A comprehensive hepatitis C viral kinetic model explaining cure. Clin Pharmacol Ther 87: 706-713. doi: 10.1038/clpt.2010.35 PMID: 20463660

57. Guedj J, Dahari H, Rong L, Sansone ND, Nettles RE, et al. (2013) Modeling shows that the NS5A inhibitor daclatasvir has two modes of action and yields a shorter estimate of the hepatitis C virus half-life. Proc Natl Acad Sci U S A 110: 3991-3996. doi: 10.1073/pnas.1203110110 PMID: 23431163

58. Garg AX, Adhikari NK, McDonald H, Rosas-Arellano MP, Devereaux PJ, et al. (2005) Effects of computerized clinical decision support systems on practitioner performance and patient outcomes: a systematic review. JAMA 293: 1223-1238. PMID: 15755945

59. Perelson AS, Nelson PW (1999) Mathematical analysis of HIV-1 dynamics in vivo. SIAM review 41: 42.

60. Ke R, Lloyd-Smith JO (2012) Evolutionary analysis of human immunodeficiency virus type 1 therapies based on conditionally replicating vectors. PLoS Comput Biol 8: e1002744. doi: 10.1371/journal.pcbi.1002744 PMID: 23133349