A hepatitis C virus infection model with time-varying drug effectiveness: solution and analysis

Abstract

Simple models of therapy for viral diseases such as hepatitis C virus (HCV) or human immunodeficiency virus assume that, once therapy is started, the drug has a constant effectiveness. More realistic models have assumed either that the drug effectiveness depends on the drug concentration or that the effectiveness varies over time. Here a previously introduced varying-effectiveness (VE) model is studied mathematically in the context of HCV infection. We show that while the model is linear, it has no closed-form solution due to the time-varying nature of the effectiveness. We then show that the model can be transformed into a Bessel equation and derive an analytic solution in terms of modified Bessel functions, which are defined as infinite series, with time-varying arguments. Fitting the solution to data from HCV infected patients under therapy has yielded values for the parameters in the model. We show that for biologically realistic parameters, the predicted viral decay on therapy is generally biphasic and resembles that predicted by constant-effectiveness (CE) models. We introduce a general method for determining the time at which the transition between decay phases occurs based on calculating the point of maximum curvature of the viral decay curve. For the parameter regimes of interest, we also find approximate solutions for the VE model and establish the asymptotic behavior of the system. We show that the rate of second phase decay is determined by the death rate of infected cells multiplied by the maximum effectiveness of therapy, whereas the rate of first phase decline depends on multiple parameters including the rate of increase of drug effectiveness with time.

Figure 1. Example of a biphasic decline of HCV, following a short delay, after initiation of interferon- therapy at .

Fit of Neumann et al. model (solid line) to data for Patient 1E (dots) from .

Figure 2. HCV viral load undergoes biphasic decay upon initiation of silibinin treatment at time .

The transition time between the first and second phases, , is calculated by maximizing the curvature in equation (14), and is marked by a vertical dashed line. VE model fit of Canini et al. (solid line) and HCV viral load data (dots) for (a) Patient 46, with transition time days, and for (b) Patient 48, with transition time days.

Figure 3. Approximate and analytic solution of VE model.

(a) Comparison of analytic solution (equation (11)) and the approximation (equation (16)) assuming sibilinin treatment (see Table 1 for parameters) and initial viral load of . (b) Relative error in of approximation.

Figure 4. Different exponential terms in approximate solution (16) compared with the exact solution and for silibinin treatment parameters, for which (see Table 1).

(a) Exponential terms from (16) plotted separately. (b) Exponential terms from (16) plotted in combined form.

Figure 5. Approximate and analytic solution of the VE model under danoprevir () or telaprevir () treatment with patient data.

(a,c) Approximate solution (16) compared to the analytic solution (11) for (a) danoprevir or (c) telaprevir treatment. (b,d) Different exponential terms in approximate solution compared with the exact solution, with decay phases indicated, for (b) danoprevir or (d) telaprevir treatment. Danoprevir treatment: data from patient 04-94XD (dosing 200 mg tid) in with associated parameter estimates for VE model , , , , , , and [unpublished]. Telaprevir treatment: data from patient 6 in with associated parameter estimates , , , , , , and .

Figure 6. Approximation to viral dynamics compared to exact dynamics under mericitabine treatment, 750 mg qd, .

(a) For patient 92102 from , characterized as “flat”. (b) For patient 92103 from , characterized as “non-flat”. (c) Different exponential terms in approximate solution (16) compared with the exact solution for patient 92103, characterized as “non-flat”. Parameter estimates from : For patient 92102, , , , , , , and ; for patient 92103, , , , , , , and .

Figure 7. Truncated series solutions for the VE model compared with the exact solution (11) under silibinin treatment (; see Table 1 for parameters).

Legend: (i) Series terms with exponents , , , and terms, included in the approximation (16), from the series solution (19); (ii) Series terms with exponents from (i) and also the and terms missing from the approximation.