- Research
- Open Access
- Published:

# Improving the estimation of the death rate of infected cells from time course data during the acute phase of virus infections: application to acute HIV-1 infection in a humanized mouse model

*Theoretical Biology and Medical Modelling*
**volume 11**, Article number: 22 (2014)

## Abstract

### Background

Mathematical modeling of virus dynamics has provided quantitative insights into viral infections such as influenza, the simian immunodeficiency virus/human immunodeficiency virus, hepatitis B, and hepatitis C. Through modeling, we can estimate the half-life of infected cells, the exponential growth rate, and the basic reproduction number (*R*_{0}). To calculate *R*_{0} from virus load data, the death rate of productively infected cells is required. This can be readily estimated from treatment data collected during the chronic phase, but is difficult to determine from acute infection data. Here, we propose two new models that can reliably estimate the average life span of infected cells from acute-phase data, and apply both methods to experimental data from humanized mice infected with HIV-1.

### Methods

Both new models, called as the reduced quasi-steady state (RQS) model and the piece-wise regression (PWR) model, are derived by simplification of a standard model for the acute-phase dynamics of target cells, viruses and infected cells. By having only a limited number of parameters, both models allow us to reliably estimate the death rate of productively infected cells. Simulated datasets with plausible parameter values are generated with the standard model to compare the performance of the new models with that of the major previous model (i.e., the simple exponential model). Finally, we fit models to time course data from HIV-1 infected humanized mice to estimate the several important parameters characterizing their acute infection.

### Results and conclusions

The new models provided much better estimates than the previous model because they more precisely capture the *de novo* infection process. Both models describe the acute phase of HIV-1 infected humanized mice reasonably well, and we estimated an average death rate of infected cells of 0.61 and 0.61, an average exponential growth rate of 0.69 and 0.76, and an average basic reproduction number of 2.30 and 2.38 in the RQS model and the PWR model, respectively. These estimates are fairly close to those obtained in humans.

## Background

In most viral infections, the initial exponential growth phase is followed by a second exponential phase known as contraction. In hosts exposed to viruses such as influenza and coronaviruses (the causative agents of severe acute respiratory syndrome), the viral load continuously declines during the contraction phase [1, 2]. In contrast, in chronic viral infections, such as human immunodeficiency virus (HIV) and hepatitis C virus (HCV) infections, contraction slows down such that the viral load approaches a steady state, called the virological set point [3, 4]. In both infection types, the expansion and contraction of the viral load have been modeled as single exponential functions, with parameters determined by linear regression of the log transformed data [3, 5–7]. This simple approach is reasonable as long as the conditions, e.g., the availability of target cells or the immune response, hardly change within each phase. Using this approach, the initial growth rate, death rate of the infected cells, and the basic reproduction number (i.e., *R*_{0}) have been estimated, which has improved our understanding of particular virus infections, and has guided medical treatment [7–9]. For example, once the basic reproduction number is estimated, the critical inhibition, 1–1/*R*_{0}, induced by vaccines, or by antiviral drugs, to prevent primary virus infection can be calculated [6].

Knowledge of the death rate of infected cells is crucial for properly understanding viral dynamics because the average life time is required for calculating the basic reproduction number. In chronic viral infections such as HIV and HCV, the death rate is estimated from large perturbations of the set point viral load data instigated by potent anti-viral therapy [4, 10–13]. Shortly after effective treatment, the decay rate of viral load approaches the death rate of productively infected cells. However, estimating the death rate of infected cells during the acute phase remains a challenging task, and, as a consequence, it is difficult to accurately estimate the basic reproduction number from viral load data during the early stages of viral infection. In addition to calculation of the basic reproduction number, the death rate *per se* is useful for evaluating the efficacy of vaccine-induced cellular immune responses during the acute phase of virus infection [14–16]. Therefore, an improved method for estimating the death rate of infected cells during this phase is urgently required.

In this study, we first generated simulated datasets with biologically plausible parameter values, using a population dynamics model of virus population. The time evolution of target cell densities and viral load were modeled during the acute phase. The datasets describing acute infection were subsequently analyzed by two novel mathematical models to evaluate whether the new models could accurately estimate the known parameters. Our proposed models properly described the artificial datasets and delivered better estimates of the parameters and well calculated indices than conventional models (i.e., simple exponential models). Our methods proved especially effective for calculating the death rate of infected cells. We then applied our models to time course data from a human hepatopoietic stem cell-transplanted humanized mouse model infected with HIV type-1 (HIV-1) [17–20], to quantify the infection dynamics during the acute phase. To our knowledge, this is the first report quantifying the dynamics of acute HIV-1 infection in humanized mice. Finally, we discuss how our approach may be combined with animal experiments. Like previous simple exponential models [3, 5–7], our approach is quite general and can be used in several infection models.

## Methods

### Mathematical models describing the acute phase of virus infection

#### (I) Reduced quasi-steady state (RQS) model

The standard model for viral infection consists of three differential equations for target cells, *T*(*t*), infected cells, *I*(*t*), and viral particles, *V*(*t*) [7–9]. Since during acute infection the normal production and loss of target cells is much smaller than the loss due to viral infection and/or its side effects [14–16], the standard model can be reduced to

where the parameters *β*, *δ*, *p* and *c* represent the conventional rate constants for viral infection of target cells, the death rate of infected cells, the virus production rate in an infected cell, and the clearance rate of virus particles, respectively. The initial expansion of viral load in this model is well approximated by *V*(*t*) ≈ *V*(0) exp (*g*_{0}*t*) [3, 5–7], with an exponential growth rate, *g*_{0} (the Malthusian parameter), given by the positive root of the characteristic equation *g*_{0}^{2} + (*c* + *δ*)*g*_{0} + *cδ* - *pβT*(0) = 0, i.e., {\mathit{g}}_{0}=\left\{-\left(\mathit{c}+\mathit{\delta}\right)+\sqrt{{\left(\mathit{c}-\mathit{\delta}\right)}^{2}+4\mathit{p}\mathit{\beta}\mathit{T}\left(0\right)}\right\}/2.

This model can be simplified further by a quasi-steady state (QSS) approximation for the viral particles [7, 21–23]. Typical estimated half-lives of viruses (1/*c*) such as HIV, HCV and hepatitis B virus are of the order of minutes (or hours), whereas those of infected cells (1/δ) *in vivo* are of the order of days [4, 5, 7–13]. Since, the clearance rate of viral particles, *c*, is typically much larger than the death rate, *δ*, of the infected cells, we can make a QSS assumption, *V*′(*t*) = 0, and replace Eq. (3) by *V*(*t*) = *pI*(*t*)/*c*. Because we fit viral loads, *V*(*t*), rather than the number of infected cells, *I*(*t*), we also substitute *I*(*t*) = *cV*(*t*)/*p* into Eq. (2) to obtain

where *r* = *pβ*/*c* is the viral replication rate per target cell, and *δ* is the death rate of infected cells. Eqs. (1) and (4) together form our first model, that we here call the “reduced quasi-steady state” (RQS) model. The RQS model lumps the 8 parameters of the reduced standard model of Eqs. (1–3) into five parameters, i.e., *β, r, δ, T*(*0*), and *V*(*0*). Because there is no production of target cells the infection will ultimately be cleared. The five parameters together define several “observables”. First, the basic reproduction number is *R*_{0} = *rT*(0)/*δ*. Second, the initial exponential growth rate is *g*_{0} ≈ *rT*(0) - *δ* ( when *δ* and *g*_{0} ≪ *c*, one can directly calculate the same *g*_{0} from the characteristic equation). Third, the final level of target cells is given by the epidemiological “final size equation” [24, 25], as the solution of *f* = exp [ - *R*_{0}(1 - *f*)], where *f* = *T*_{
min
}/*T*(0) is the fraction of surviving target cells.

#### (II) Piece wise targets (PWT) model

In the standard model, and its simplifications outlined above, the dynamics of the target cells are coupled to the density of viral particles because target cells disappear by infection. Since target cell densities in the peripheral blood (PB) may also depend on other factors, like inflammation, activation and redistribution, we next write a model where the target cell dynamics are decoupled from the viral dynamics. Since, during the acute phase of several virus infections, such as HIV, SIV and SHIV, the decrease in number of target cells in PB is preceded by an initial flat phase [14–20, 26–28], we propose a phenomenological model for the target cells consisting of an initial flat phase, and a second phase of exponential loss (Figure 1). This basically implies that we assume that target cells only become a limiting factor when their density starts to decline. Thus, the dynamics of target cells is described as follows:

where the parameter Δ represents the daily rate of target cells loss following the initial flat phase, and *t*^{*} corresponds to the time at which the target cell densities begin to decrease. Eqs. (4–6) define our second model that we here call the “piece wise targets” (PWT) model. The PWT model has six parameters, i.e., *β, r, δ, T*(*0*), *V*(*0*), Δ, and *t*^{*}, and because it shares Eq. (4) with the RQS model, it has the same definitions for the replication rate, *r* = *pβ*/*c*, the Malthusian parameter, *g*_{0} = *rT*(0) - *δ*, and the reproduction number *R*_{0} = *rT*(0)/*δ*. In contrast to the partial depletion, *f*, in the RQS model, the target cells will ultimately be completely depleted, i.e., *T*(∞) → 0, in the PWT model. Thanks to the decoupled dynamics of the target cells the PWT model can be solved analytically:

where the replication phase (Eq. (7)) is identical to the initial phase of the standard model (see the remark above). Thus, at the price of one additional parameter, we can generalize the depletion of target cells to mechanisms (e.g., inflammation, activation and redistribution of target cells) other than infection only, and have a model with very similar parameters characterizing the acute viral infection.

#### (III) Piece wise regression (PWR) model

As a control method, we additionally consider a classical method that has been widely adopted in earlier literatures [3, 5–7]. Previously acute infection data have been quantified using piece wise linear regression of the log transformed viral loads before and after the peak in the viral load [3, 5–7]. The ascending and descending slopes of the log viral load roughly correspond to the exponential growth rate, *g*_{0}, and the death rate, *δ*, of infected cells, respectively [3, 5–7]. Since *g*_{0} = *rT*(0) - *δ* and *R*_{0} = *rT*(0)/*δ*, these two slopes suffice to estimate the basic reproduction number, *R*_{0} ≈ 1 + *g*_{0}/*δ*, and knowing the initial target cell density, *T*(0), the viral replication rate per target cell can be estimated from the *R*_{0}, i.e., *r* = *δR*_{0}/*T*(0) [3, 5–7]. Here we call this classical model as the “piece wise regression” model (PWR). It has been realized before that the down-slope will only reflect the death rate of infected cells, *δ*, if there is hardly any residual infection of target cells during the contraction phase, i.e., if target cells are markedly depleted [3, 5–7].

### Artificial data generated from a population dynamics model of virus infection

To estimate the accuracy of the parameters estimated by our two novel models, we created simulated time course data of target cell densities and viral load during the acute phase of viral infection (lasting approximately 21 days [14–20, 26–28]) assuming biologically plausible parameter values. The artificial datasets were generated with the reduced standard population dynamics model of viral infection, i.e., Eqs. (1–3), in which the target cell dynamics are coupled to the viral dynamics by the infection term.

We added stochastic variation to the “data” generated by this model by adding “observational” noise and/or by varying the parameter values. The log transformed data were perturbed by adding a normally distributed noise variable with zero mean and standard deviation *σ* (see Results). The datasets describing acute infection were subsequently analyzed by the two novel RQS and PWT models, and by the previous PWR model.

### HIV-1 infection in humanized mice

The dynamics of HIV-1 infection during acute infection were quantified in a human hepatopoietic stem cell-transplanted humanized mouse model (NOG-hCD34 mice) [17–20]. Five humanized mice were infected with the CCR5-tropic HIV-1 (strain AD8) [29], and 100 μl of peripheral blood (PB) was routinely collected under anesthesia through the retro-orbital venous plexus at 0, 3, 7, 14, and 21 days post-infection, as previously described [17–20]. The amount of viral RNA in 50 μl of plasma was quantified by RT-PCR (Bio Medical Laboratories, Inc). To estimate target cell densities, the number of memory CD4^{+} T cells was measured by hematometry and flow cytometry, as previously described [17–20]. Briefly, the number of human leukocytes in 10 μl of peripheral blood (PB) was measured in a Celltac α MEK-6450 hematology analyzer (Nihon Kohden, Co.), and the percentage of memory CD4^{+} T cells in human CD45^{+} leukocytes (i.e., CD45^{+} CD3^{+} CD4^{+} CD45RA^{-} cells) was quantified in a FACSCanto II (BD Biosciences) flow cytometer. In the flow cytometry analyses, APC-conjugated anti-CD4 antibody (RPA-4; Biolegend), APC-Cy7-conjugated anti-CD3 antibody (HIT3a; Biolegend), and PE-conjugated anti-CD45 antibody (HI30; Biolegend) were used.

All protocols involving human subjects were reviewed and approved by the Kyoto University institutional review board. Informed written consent from the human subjects was obtained in this study.

## Results

### Coverage probability of the mathematical models

In the Methods section we formulate two novel mathematical models describing the target cell densities and the viral load during acute infection. We created artificial data with target cell densities and virus loads during acute infection using the reduced standard model for viral infection (i.e., Eqs. (1–3)). The data was generated for one ml of PB with “typical” values of the parameters for HIV-1, i.e., an infection rate *β* = 10^{-8}per cell per day, a virus production rate *p* = 4000 particles per day, a death rate of infected cells *δ* = 0.7 per day, and a clearance rate of *c* = 23 per day [7–11]. We study whether our simplified models can describe the *in silico* data, and whether their (lumped) parameters are identifiable. The major biological observables of this model are the initial viral growth rate, {\mathit{g}}_{0}=\left\{-\left(\mathit{c}+\mathit{\delta}\right)+\sqrt{{\left(\mathit{c}-\mathit{\delta}\right)}^{2}+4\mathit{p}\mathit{\beta}\mathit{T}\left(0\right)}\right\}/2, the viral replication rate per target cell, *r* = *pβ*/*c*, the death rate of infected cells, *δ*, and the basic reproduction number, *R*_{0} = *rT*(0)/*δ*.

First, we created *in silico* data with “observational” noise by adding proportional random variation to each data point. Specifically, we drew random values from a Gaussian distribution with a mean of one and a standard deviation *σ* = 0.2, and added these values to the log transformed data. This seemed natural as we are also fitting the log transformed data, and on a log scale this corresponds to the measurement error of about 60%. The generated datasets were fitted to the numerical solutions of the RQS model (Eqs. (1) and (4)), the analytical solution of the PWT model (Eqs. (5–8)), and to the previous PWR model. The sum of squared residuals was minimized using the FindMinimum package of *Mathematica 9.0*, fitting the target cell and viral load data simultaneously. The typical behavior of the models using these best-fit parameter estimates is depicted in Figure 2, together with the simulated data. Other standard deviations of the parameters yielded similar results (results not shown). The two novel models reasonably describe the acute phase of viral infection. Note that target cell densities are partially depleted in the RQS model, and will ultimately approach “0” in the PWT model.

Second, we created different “patients” by randomly drawing the parameter values from normal distributions centered at their typical values. Thus, the infection rate, death rate of infected cells and virus production rate were assumed to be normally distributed as \begin{array}{cc}\hfill \mathit{\beta}~\mathit{N}\left({\mathit{\mu}}_{\mathit{\beta}},{\mathit{\sigma}}_{\mathit{\beta}}^{2}\right),\hfill & \hfill \mathit{\delta}~\mathit{N}\left({\mathit{\mu}}_{\mathit{\delta}},{\mathit{\sigma}}_{\mathit{\delta}}^{2}\right),\hfill \end{array} and \mathit{p}~\mathit{N}\left({\mathit{\mu}}_{\mathit{p}},{\mathit{\sigma}}_{\mathit{p}}^{2}\right), with \begin{array}{ccc}\hfill {\mathit{\mu}}_{\mathit{\beta}}={10}^{-8},\hfill & \hfill {\mathit{\mu}}_{\mathit{\delta}}=0.7,\hfill & \hfill {\mathit{\mu}}_{\mathit{p}}=4000,\hfill \end{array} respectively. The standard deviations were set as \begin{array}{ccc}\hfill {\mathit{\sigma}}_{\mathit{\beta}}={10}^{-9},\hfill & \hfill {\mathit{\sigma}}_{\mathit{\delta}}=0.3,\hfill & \hfill {\mathit{\sigma}}_{\mathit{p}}=400\hfill \end{array}. In this way, we obtained a distribution of the basic reproduction number centered around the true value (*R*_{0} = 2.48). We then randomly sampled one parameter set of *β*, *δ*, *p* from the distributions, and produced 200 different artificial datasets as explained above. Analyzing each dataset with the same three models (RQS, PWT and PWR), we calculated 95% confidence intervals (CI) for *g*_{0}, *r*, *δ* and *R*_{0}, and investigated whether the 95% CI successfully contained the true values of *g*_{0}, *r*, *δ* and *R*_{0} used to create the data set. This procedure was repeated for all 200 datasets, and Table 1 provides the frequency of datasets for which the 95% CI successfully contained the true values of *g*_{0}, *r*, *δ* and *R*_{0}, i.e., the coverage probability. Although the initial growth rate was well estimated by the previous PWR model, the novel RQS and PWT models estimated the viral replication rate, death rate of infected cells, and basic reproduction number much more accurately than the PWR model. Thus, the novel models can more accurately extract information from acute-phase viral infection data.

### Sensitivity of parameter estimation to target cell depletion

At the beginning of the infection all models are identical as they all predict exponential growth of the virus load. The models differ around the peak because the previous PWR model assumes an exponential contraction after the peak, whereas the new models allow the peak be formed by the loss of target cells (i.e., by the *βT*(*t*)*V*(*t*) term in the RQS model, and by the exponential loss of target cells in the PWT model). In both models this loss of target cells continues during the contraction phase. However, if target cells were depleted rapidly such that there would be hardly any infection of target cells during the contraction phase, this difference among the models would vanish, and the contraction phase of the new models would also be dominated by the death rate of infected cells.

The level to which the target cells become depleted in the reduced standard model of viral infection (Eqs. (1–3)) can be computed with the epidemiological “final size equation” [24, 25] (see the Methods section). Using the same equation, Petravic *et al.*[16] show that this final size of the target cell level provides a good description of the nadir, *T*_{
min
}, of the target cell density during an acute infection. Defining *f* = *T*_{
min
}/*T*(0) as the fraction of surviving target cells, we use the final size equation to compute different values of the infection rate *β* to vary the nadir of the target cells over the interval *f* ∈ [0.001, 0.200]. Doing so we again created different cases, each with a different level of target cell depletion. For simplicity this was done in the absence of noise (which does not affect these results). Using the same approach as explained above we fit the “data” generated by these cases with the 3 models (Figure 3).

As discussed previously [3, 5–7], the PWR model (orange symbols) fails to correctly estimate the death rate of infected cells in the presence of continued *de novo* infections, i.e., when many target cells survive; see Figure 3A (in which the dashed line denotes the correct death rate at *δ* = 0.7).

Conversely, the two novel models (RQS: blue symbols and PWT: green symbols) accurately estimate the death rate even if many target cells survive. Note that it is not surprising that the RQS model provides better parameter estimates than the PWT model, because the artificial data was generated by the reduced standard model (Eqs. (1–3)), from which the RQS model was derived by a reasonable quasi-steady state assumption. The parameters are accurately estimated by all methods when target cells are severely depleted (i.e., when *f* approaches 0), reconfirming that the previous PWR model can accurately quantify the infected cell death when the target cells are severely depleted. This, for example, occurs in CXCR4-tropic simian–human immunodeficiency virus (SHIV) infections (which deplete naïve and memory CD4 T cells during the acute phase) [14, 16, 26–28]. In cases where target cells are not as strongly depleted, the new models do much better than the PWR model in estimating the death rate (Figure 3A), the viral replication rate (Figure 3B) and the basic reproduction number (Figure 3C). All models perform similarly on estimating the exponential growth rate (Figure 3D), because they all assume (implicitly or explicitly) that the number of target cells remains constant during the earliest phases of viral infections.

### Application to acute HIV-1 infection in a humanized mouse model

Having established that the novel models outperform the old one, they were fitted to the 21-day time courses of viral loads and target cells observed in five virus-infected humanized mice. From the parameters estimated with the individual mouse datasets (see Table 2 for the RQS model and Table 3 for the PWT model), we obtained the similar *g*_{0} ranging from 0.43 to 1.07 and from 0.48 to 1.06 per day in the RQS and PWT models, respectively. These estimates are not in disagreement with the replication rate of HIV-1 in human patients, which has been estimated to be 1.01 ± 0.37 per day [6]. We estimate a death rate of HIV-1 infected cells in humanized mice as *δ* ranging from 0.30 to 0.76 and from 0.38 to 0.76 per day with the RQS model and the PWT model, respectively (see Tables 2 and 3). Again, this result is in concordance with estimates of viral death rate in treated HIV-1 infected patients [10, 11]. We then determined the basic reproduction number *R*_{0} of HIV-1 in humanized mice from the individual estimates of *T*(0), *r* and *δ*, obtaining *R*_{0} ranging from 1.58 to 3.88 and from 1.77 to 3.26 in the RQS and PWT models, respectively (see Tables 2 and 3). Since the mean of *R*_{0} corresponds to a predicted target cell nadir of *f* = 0.12 in the RQS model, the simple exponential model is not expected to perform equally well on this data (see Figure 3).

The estimated parameter values of each individual mouse are given in Tables 2 and 3. Using the best-fit parameter estimates the behavior of RQS and PWT models is depicted together with the individual data in Figures 4A and B, respectively, which confirms that both models reasonably describe the acute phase of HIV-1 infection in humanized mice.

## Discussion

We here propose two novel models to quantify the most important parameters characterizing acute viral infections. Both models are major improvements over the previous simple exponential model [3, 5–7] because that model has difficulties estimating the death rate of infected cells when target cells are not depleted after the viral load peak (see Figure 3). The novel models use the observed target cell densities when estimating the parameter values, and by using simulated data we have demonstrated that the new models typically outperform the previous model. Applying the new models to data obtained in humanized mice estimates that the rates at which the virus expands, and at which infected cells die, resemble those measured in humans.

The efficacy of vaccines eliciting cytotoxic immune responses [14–16] could be quantified by our new approach by comparing the estimated death rate of infected cells between normal and vaccinated individuals during acute infection. Additionally, by comparing the death rates estimated during acute infection with those estimated by perturbations with antiviral drugs during the chronic phase [4, 10–13], we might be able to compare the impact of cellular immunity during both phases of an infection. Both examples illustrate that the new models can markedly improve our understanding of viral infection dynamics *in vivo*.

Our novel approach overcomes the difficulty the previous PWR model had with estimating the death rate of infected cells in situations where target are not severely depleted. Additionally, the choice between the two different models that we propose here can be made on the final densities of the target cells. If the number of target cells reaches a nadir during acute infection, the RQS model seems most appropriate. If the target cells continue to decrease the PWT model should be better.

Thus, whenever one has sufficient time course data from an acute infection, the new models should allow one to estimate the death rate of infected cells, and hence the *R*_{0}, with reasonable accuracy. Indeed, in animal experiments using rhesus macaques [14–16, 26–28], ferrets [30], and mice [17–20, 31–34], both target cell densities and viral loads have been measured during the acute phase of viral infection. For example, the number of uninfected (and infected) cells in the lungs or respiratory tracts of macaques, ferrets and mice that were experimentally infected with influenza could be measured [30, 31, 34]. Using the CXCR4-tropic SHIV/macaque model, both target cell densities (naïve and memory CD4 T cells) and viral loads from PB have been measured [14, 16, 26–28]. The target cells of simian immunodeficiency virus (SIV), or CCR5-tropic SHIV infection (memory CD4 T cells expressing CCR5), have been measured from gastrointestinal mucosa samples [15]. Thus there are several infection models that can be analyzed with our new models.

## Conclusion

In this paper, we developed novel mathematical approaches to estimating parameters from acute viral infection data. We demonstrated that the new models outperform the previous model using simulated data. We quantified the dynamics of acute-phase HIV-1 infections by measuring their time course data in a humanized mouse model. Interestingly we find that the rates at which the virus expands, and at which infected cells die, are similar to those in humans.

## References

Carrat F, Vergu E, Ferguson NM, Lemaitre M, Cauchemez S, Leach S, Valleron AJ: Time lines of infection and disease in human influenza: a review of volunteer challenge studies. Am J Epidemio. 2008, 167 (7): 775-785. 10.1093/aje/kwm375.

Peiris JS, Chu CM, Cheng VC, Chan KS, Hung IF, Poon LL, Law KI, Tang BS, Hon TY, Chan CS, Chan KH, Ng JS, Zheng BJ, Ng WL, Lai RW, Guan Y, Yuen KY, HKU/UCH SARS Study Group: Clinical progression and viral load in a community outbreak of coronavirus-associated SARS pneumonia: a prospective study. Lancet. 2003, 361 (9371): 1767-1772. 10.1016/S0140-6736(03)13412-5.

Little SJ, McLean AR, Spina CA, Richman DD, Havlir DV: Viral dynamics of acute HIV-1 infection. J Exp Med. 1999, 190 (6): 841-850. 10.1084/jem.190.6.841.

Neumann AU, Lam NP, Dahari H, Gretch DR, Wiley TE, Layden TJ, Perelson AS: Hepatitis C viral dynamics in vivo and the antiviral efficacy of interferon-alpha therapy. Science. 1998, 282 (5386): 103-107.

Nowak MA, Lloyd AL, Vasquez GM, Wiltrout TA, Wahl LM, Bischofberger N, Williams J, Kinter A, Fauci AS, Hirsch VM, Lifson JD: Viral dynamics of primary viremia and antiretroviral therapy in simian immunodeficiency virus infection. J Virol. 1997, 71 (10): 7518-7525.

Ribeiro RM, Qin L, Chavez LL, Li D, Self SG, Perelson AS: Estimation of the initial viral growth rate and basic reproductive number during acute HIV-1 infection. J Virol. 2010, 84 (12): 6096-6102. 10.1128/JVI.00127-10.

Nowak MA, May RM: Virus Dynamics. 2000, New York, USA: Oxford University Press

Perelson AS: Modelling viral and immune system dynamics. Nat Rev Immunol. 2002, 2 (1): 28-36. 10.1038/nri700.

Simon V, Ho DD: HIV-1 dynamics in vivo: implications for therapy. Nat Rev Microbiol. 2003, 1 (3): 181-190. 10.1038/nrmicro772.

Perelson AS, Essunger P, Cao Y, Vesanen M, Hurley A, Saksela K, Markowitz M, Ho DD: Decay characteristics of HIV-1-infected compartments during combination therapy. Nature. 1997, 387 (6629): 188-191. 10.1038/387188a0.

Markowitz M, Louie M, Hurley A, Sun E, Di Mascio M, Perelson AS, Ho DD: A novel antiviral intervention results in more accurate assessment of human immunodeficiency virus type 1 replication dynamics and T-cell decay in vivo. J Virol. 2003, 77 (8): 5037-5038. 10.1128/JVI.77.8.5037-5038.2003.

Dahari H, Guedj J, Perelson AS, Layden TJ: Hepatitis C viral kinetics in the era of direct acting antiviral agents and IL28B. Curr Hepa Rep. 2011, 10 (3): 214-227. 10.1007/s11901-011-0101-7.

Iwami S, Koizumi Y, Ikeda H, Kakizoe Y: Quantification of viral infection dynamics in animal experiments. Front Microbiol. 2013, 4: 264-

Davenport MP, Zhang L, Shiver JW, Casmiro DR, Ribeiro RM, Perelson AS: Influence of peak viral load on the extent of CD4+ T-cell depletion in simian HIV infection. J Acquir Immune Defic Syndr. 2006, 41 (3): 259-265. 10.1097/01.qai.0000199232.31340.d3.

Wilson DP, Mattapallil JJ, Lay MD, Zhang L, Roederer M, Davenport MP: Estimating the infectivity of CCR5-tropic simian immunodeficiency virus SIV(mac251) in the gut. J Virol. 2007, 81 (15): 8025-8029. 10.1128/JVI.01771-06.

Petravic J, Ribeiro RM, Casimiro DR, Mattapallil JJ, Roederer M, Shiver JW, Davenport MP: Estimating the impact of vaccination on acute simian-human immunodeficiency virus/simian immunodeficiency virus infections. J Virol. 2008, 82 (23): 11589-11598. 10.1128/JVI.01596-08.

Nie C, Sato K, Misawa N, Kitayama H, Fujino H, Hiramatsu H, Heike T, Nakahata T, Tanaka Y, Ito M, Koyanagi Y: Selective infection of CD4+ effector memory T lymphocytes leads to preferential depletion of memory T lymphocytes in R5 HIV-1-infected humanized NOD/SCID/IL-2Rgammanull mice. Virology. 2009, 394 (1): 64-72. 10.1016/j.virol.2009.08.011.

Sato K, Izumi T, Misawa N, Kobayashi T, Yamashita Y, Ohmichi M, Ito M, Takaori-Kondo A, Koyanagi Y: Remarkable lethal G-to-A mutations in vif-proficient HIV-1 provirus by individual APOBEC3 proteins in humanized mice. J Virol. 2010, 84 (18): 9546-9556. 10.1128/JVI.00823-10.

Sato K, Misawa N, Fukuhara M, Iwami S, An DS, Ito M, Koyanagi Y: Vpu augments the initial burst phase of HIV-1 propagation and downregulates BST2 and CD4 in humanized mice. J Virol. 2012, 86 (9): 5000-5013. 10.1128/JVI.07062-11.

Sato K, Misawa N, Iwami S, Satou Y, Matsuoka M, Ishizaka Y, Ito M, Aihara K, An DS, Koyanagi Y: HIV-1 Vpr accelerates viral replication during acute infection by exploitation of proliferating CD4+ T cells in vivo. PLoS Pathog. 2013, 9 (12): e1003812-10.1371/journal.ppat.1003812.

Iwasa Y, Michor F, Nowak M: Some basic properties of immune selection. J Theor Biol. 2004, 229 (2): 179-188. 10.1016/j.jtbi.2004.03.013.

Thieme HR: Mathematics in Population Biology. 2003, Princeton, USA: Princeton University Press

Komarova NL, Levy DN, Wodarz D: Effect of synaptic transmission on viral fitness in HIV infection. PLoS One. 2012, 7 (11): e48361-10.1371/journal.pone.0048361.

Anderson RM: The Kermack-McKendrick epidemic threshold theorem. Bull Math Biol. 1991, 53 (1–2): 3-32.

Diekmann O, Heesterbeek JAP: Mathematical Epidemiology of Infectious Diseases: Model Building, Analysis and Interpretation. 2000, New York, USA: John Wiley & Sons

Igarashi T, Brown CR, Endo Y, Buckler-White A, Plishka R, Bischofberger N, Hirsch V, Martin MA: Macrophage are the principal reservoir and sustain high virus loads in rhesus macaques after the depletion of CD4+ T cells by a highly pathogenic simian immunodeficiency virus/HIV type 1 chimera (SHIV): Implications for HIV-1 infections of humans. Proc Natl Acad Sci U S A. 2001, 98 (2): 658-663. 10.1073/pnas.98.2.658.

Igarashi T, Endo Y, Nishimura Y, Buckler C, Sadjadpour R, Donau OK, Dumaurier MJ, Plishka RJ, Buckler-White A, Martin MA: Early control of highly pathogenic simian immunodeficiency virus/human immunodeficiency virus chimeric virus infections in rhesus monkeys usually results in long-lasting asymptomatic clinical outcomes. J Virol. 2003, 77 (20): 10829-10840. 10.1128/JVI.77.20.10829-10840.2003.

Kozyrev IL, Ibuki K, Shimada T, Kuwata T, Takemura T, Hayami M, Miura T: Characterization of less pathogenic infectious molecular clones derived from acute-pathogenic SHIV-89.6p stock virus. Virology. 2001, 282 (1): 6-13. 10.1006/viro.2000.0839.

Theodore TS, Englund G, Buckler-White A, Buckler CE, Martin MA, Peden KW: Construction and characterization of a stable full-length macrophage-tropic HIV type 1 molecular clone that directs the production of high titers of progeny virions. AIDS Res Hum Retroviruses. 1996, 12 (3): 191-194. 10.1089/aid.1996.12.191.

Imai M, Watanabe T, Hatta M, Das SC, Ozawa M, Shinya K, Zhong G, Hanson A, Katsura H, Watanabe S, Li C, Kawakami E, Yamada S, Kiso M, Suzuki Y, Maher EA, Neumann G, Kawaoka Y: Experimental adaptation of an influenza H5 HA confers respiratory droplet transmission to a reassortant H5 HA/H1N1 virus in ferrets. Nature. 2012, 486 (7403): 420-428.

Itoh Y, Shinya K, Kiso M, Watanabe T, Sakoda Y, Hatta M, Muramoto Y, Tamura D, Sakai-Tagawa Y, Noda T, Sakabe S, Imai M, Hatta Y, Watanabe S, Li C, Yamada S, Fujii K, Murakami S, Imai H, Kakugawa S, Ito M, Takano R, Iwatsuki-Horimoto K, Shimojima M, Horimoto T, Goto H, Takahashi K, Makino A, Ishigaki H, Nakayama M: In vitro and in vivo characterization of new swine-origin H1N1 influenza viruses. Nature. 2009, 460 (7258): 1021-1025.

Mercer DF, Schiller DE, Elliott JF, Douglas DN, Hao C, Rinfret A, Addison WR, Fischer KP, Churchill TA, Lakey JR, Tyrrell DL, Kneteman NM: Hepatitis C virus replication in mice with chimeric human livers. Nat Med. 2001, 7 (8): 927-933. 10.1038/90968.

Miao H, Hollenbaugh JA, Zand MS, Holden-Wiltse J, Mosmann TR, Perelson AS, Wu H, Topham DJ: Quantifying the early immune response and adaptive immune response kinetics in mice infected with influenza A virus. J Virol. 2010, 84 (13): 6687-6698. 10.1128/JVI.00266-10.

Morita M, Kuba K, Ichikawa A, Nakayama M, Katahira J, Iwamoto R, Watanebe T, Sakabe S, Daidoji T, Nakamura S, Kadowaki A, Ohto T, Nakanishi H, Taguchi R, Nakaya T, Murakami M, Yoneda Y, Arai H, Kawaoka Y, Penninger JM, Arita M, Imai Y: The lipid mediator protectin D1 inhibits influenza virus replication and improves severe influenza. Cell. 2013, 153 (1): 112-125. 10.1016/j.cell.2013.02.027.

## Acknowledgments

This research is partly supported by the Kyushu University Fund (to H.I.); Grants-in-Aid for Young Scientists B23790500 (to K.S.) and B25800092 (to S.I.) from the Japan Society for the Promotion of Science (JSPS); the Uehara Memorial Foundation (to K.S.); the Shimizu Foundation for Immunological Research Grant (to K.S.); Takeda Science Foundation (to K.S.); Sumitomo Foundation Research Grant (to K.S.); Senshin Medical Research Foundation (to K.S.); JST PRESTO program (to S.I.); the Kyushu University Interdisciplinary Programs in Education and Projects in Research Development (to S.I.); the Aihara Innovative Mathematical Modeling Project, JSPS through the “Funding Program for World-Leading Innovative R & D on Science and Technology (FIRST program)”, initiated by Council for Science and Technology Policy (to K.S., K.A., and S.I.).

## Author information

### Authors and Affiliations

### Corresponding author

## Additional information

### Competing interests

The authors declare that they have no competing interests.

### Authors’ contributions

Conceived and designed the study: HI, RDB, KA and SI. Analyzed the data: HI, SM and SI. Wrote the paper: HI, RDB, KA and SI. Contributed reagents/materials/analysis tools: KS, NM and YK. All 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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

## About this article

### Cite this article

Ikeda, H., de Boer, R.J., Sato, K. *et al.* Improving the estimation of the death rate of infected cells from time course data during the acute phase of virus infections: application to acute HIV-1 infection in a humanized mouse model.
*Theor Biol Med Model* **11**, 22 (2014). https://doi.org/10.1186/1742-4682-11-22

Received:

Accepted:

Published:

DOI: https://doi.org/10.1186/1742-4682-11-22

### Keywords

- Population dynamics model
- Parameter estimation
- Virus infection dynamics
- Death rate of infected cells
- HIV-1
- Humanized mouse model