Tutorials › Quantitative Systems Pharmacology (QSP) › Target-Cell Depletion Versus Immune Control of Viral Infection

Infectious Disease and Viral Dynamics · Tutorial 394 of 799

Target-Cell Depletion Versus Immune Control of Viral Infection

Learn to formulate and compare target-cell depletion and immune-mediated viral control models, calculate their local growth thresholds, and interpret which measured trajectories are consistent with each mechanism.

Advanced 12 min read

What You'll Learn

  • Formulate target-cell-limited and immune-mediated viral dynamics models with consistent units.
  • Derive the local viral growth condition and the target-cell threshold for control.
  • Calculate how effector-mediated infected-cell killing shifts that threshold.
  • Compare predicted viral-load and target-cell trajectories under the two mechanisms.
  • Identify measurements and diagnostic checks that help distinguish mechanisms.
  • Recognize why viral-load decline alone cannot establish its cause.

Two Mechanisms Can Produce Similar Viral Declines

A falling viral load is an outcome, not a mechanism. It can occur because infection has depleted the cells needed for viral production, because immune effectors are removing infected cells, or because both processes are operating. Distinguishing these explanations matters when interpreting treatment or disease trajectories: the mechanisms imply different target-cell preservation and different responses to perturbations.

The target-cell-limited model introduced in earlier material provides a baseline. Here, the comparison is between that model and a model in which an abstract immune effector population increases infected-cell loss. The effector is deliberately not assigned to a particular lineage or molecular pathway. The focus is on model structure and predicted trajectories, not on the detailed biology of a specific immune response.

Model conventions: Let \(T\) be uninfected target-cell concentration (cells/mL), \(I\) infected-cell concentration (cells/mL), and \(V\) viral concentration (virions/mL). Time is measured in days. Parameters and state variables are nonnegative. The models treat \(T\), \(I\), and \(V\) as well-mixed concentrations and do not represent spatial structure or cell heterogeneity.

Constructing the Competing Models

For target-cell-limited infection, use

$$ \frac{dT}{dt}=s_T-d_TT-\beta VT,\qquad \frac{dI}{dt}=\beta VT-\delta I,\qquad \frac{dV}{dt}=pI-cV. $$

Here \(s_T\) is target-cell supply (cells/mL/day), \(d_T\) is the uninfected-cell loss rate (day\(^{-1}\)), \(\beta\) is the infection rate constant (mL virion\(^{-1}\) day\(^{-1}\)), \(\delta\) is the infected-cell loss rate (day\(^{-1}\)), \(p\) is viral production per infected cell (virions cell\(^{-1}\) day\(^{-1}\)), and \(c\) is viral clearance (day\(^{-1}\)). In the absence of infection, the target-cell steady state is \(T_0=s_T/d_T\). During infection, the term \(\beta VT\) removes target cells and adds infected cells at the same rate.

An immune-mediated control model retains those processes and adds an effector concentration \(E\), measured in cells/mL:

$$ \frac{dT}{dt}=s_T-d_TT-\beta VT,\qquad \frac{dI}{dt}=\beta VT-(\delta+kE)I,\qquad \frac{dV}{dt}=pI-cV, $$
$$ \frac{dE}{dt}=s_E+aE\frac{I}{K_I+I}-d_EE. $$

The parameter \(k\) is an effector-mediated killing rate constant (mL cell\(^{-1}\) day\(^{-1}\)); thus \(kE\) is a rate in day\(^{-1}\). The effector source \(s_E\) has units cells/mL/day, \(a\) and \(d_E\) have units day\(^{-1}\), and \(K_I\) is an infected-cell concentration (cells/mL). This simple expansion term is one possible phenomenological representation of effector dynamics, not a claim that all immune responses follow this equation.

Both models assume that infection is proportional to \(V T\), viral production is proportional to \(I\), and viral clearance is first order in \(V\). The immune model additionally assumes that effector-mediated killing is proportional to \(EI\). These are useful starting assumptions, but they should be questioned if measured data show saturation, delays, or changing rates that the equations cannot reproduce.

Deriving the Local Control Threshold

To see how the mechanisms affect viral growth, hold \(T\) fixed over a short interval and, in the immune model, hold \(E\) fixed as well. The infected-cell and virus equations then form a linear system. Write \(q=\delta\) for the target-cell-limited model and \(q=\delta+kE\) for the immune model. A perturbation proportional to \(e^{\lambda t}\) satisfies

$$ \begin{pmatrix} \lambda+q & -\beta T\\ -p & \lambda+c \end{pmatrix} \begin{pmatrix}I\\V\end{pmatrix} = \begin{pmatrix}0\\0\end{pmatrix}. $$

A nonzero solution requires the determinant to equal zero. Therefore,

$$ (\lambda+q)(\lambda+c)-\beta pT=0. $$

At the boundary between local growth and decline, \(\lambda=0\), so the threshold target-cell concentration is

$$ T_{\mathrm{crit}}=\frac{cq}{\beta p}. $$

If \(T>T_{\mathrm{crit}}\), the dominant eigenvalue is positive and small infected-cell/virus perturbations grow under these fixed-state assumptions; if \(T<T_{\mathrm{crit}}\), they decline. Equivalently, the local reproduction ratio is \(\mathcal{R}_T=\beta pT/(cq)\), with growth when \(\mathcal{R}_T>1\). Immune killing raises \(q\), and therefore raises \(T_{\mathrm{crit}}\): control can occur while more target cells remain. This is a local result, not a guarantee about the full trajectory when \(T\) and \(E\) are changing.

The dimensions provide a useful check. \(\beta pT\) and \(cq\) both have units day\(^{-2}\), so \(\mathcal{R}_T\) is dimensionless. The ratio \(cq/(\beta p)\) has units cells/mL, as required for \(T_{\mathrm{crit}}\).

Worked Example: Calculating the Target-Cell Threshold

Suppose a hypothetical infection has \(\beta=2.0\times10^{-7}\) mL virion\(^{-1}\) day\(^{-1}\), \(p=100\) virions cell\(^{-1}\) day\(^{-1}\), \(c=10\) day\(^{-1}\), and \(\delta=0.50\) day\(^{-1}\). Let the uninfected target-cell concentration be \(T_0=1.0\times10^6\) cells/mL.

With no added immune killing, \(q=\delta=0.50\) day\(^{-1}\). The threshold is

$$ T_{\mathrm{crit}}=\frac{(10\ \mathrm{day}^{-1})(0.50\ \mathrm{day}^{-1})} {(2.0\times10^{-7}\ \mathrm{mL\ virion}^{-1}\mathrm{day}^{-1}) (100\ \mathrm{virions\ cell}^{-1}\mathrm{day}^{-1})} =\frac{5.0\ \mathrm{day}^{-2}}{2.0\times10^{-5}\ \mathrm{mL\ cell}^{-1}\mathrm{day}^{-2}} =2.50\times10^5\ \mathrm{cells/mL}. $$

At \(T_0\), \(\mathcal{R}_T=4.00\), so the frozen-state system predicts growth. At \(T=1.0\times10^5\) cells/mL, \(\mathcal{R}_T=0.400\), so it predicts decline. The result says where the local growth boundary lies; it does not establish that depletion caused a particular observed decline.

Predicted Patterns and Their Interpretation

In a target-cell-limited trajectory, infection initially expands while the target-cell pool is abundant. As infection consumes targets, \(T\) may approach or cross \(T_{\mathrm{crit}}\); infected cells and virus then stop growing and can decline. A key prediction is that viral control is associated with substantial target-cell reduction, unless target-cell supply or other processes maintain the pool.

In the immune-mediated model, a rising \(E\) increases the infected-cell loss rate \(\delta+kE\). Viral decline can therefore begin while \(T\) remains above the target-cell-only threshold. Because fewer cells become infected, target-cell loss may also be slower, so target cells can be relatively preserved. This does not mean every immune-controlled trajectory must preserve targets: depletion may have occurred before effector activity increased, or the effector response may be too weak to prevent it.

These are patterns of consistency, not unique fingerprints. The earlier material on biphasic viral decline and long-lived infected cells shows why a change in viral slope need not identify a single process. Here, the same caution applies: observing a viral decline alone cannot distinguish a shrinking target-cell pool from increased infected-cell loss.

Worked Example: Checking a Target-Depletion Explanation

Consider hypothetical state estimates from a target-cell-limited model using the parameters in the previous example. To make a quick local comparison, assume viral clearance is fast relative to changes in \(I\), so \(V\) is approximately \(pI/c\). Substitution gives the approximate infected-cell per-capita rate \(\beta pT/c-\delta\).

Target cells \(T\) (cells/mL)\(\mathcal{R}_T\)Approximate rate \(\beta pT/c-\delta\) (day\(^{-1}\))
\(1.0\times10^6\)4.00\(2.0-0.50=1.50\)
\(2.5\times10^5\)1.00\(0.50-0.50=0\)
\(1.0\times10^5\)0.400\(0.20-0.50=-0.30\)

For example, at \(1.0\times10^5\) cells/mL, \(\beta pT/c=(2.0\times10^{-7})(100)(1.0\times10^5)/10=0.20\) day\(^{-1}\), so the approximate net rate is \(-0.30\) day\(^{-1}\). A decline occurring after target cells fall through roughly \(2.5\times10^5\) cells/mL is consistent with target limitation in this model. The fast-clearance approximation is most appropriate when \(c\) is large compared with the changing infected-cell timescale; otherwise the full two-state system should be fitted or simulated. The exact threshold still follows from the eigenvalue condition above.

Worked Example: Quantifying Immune-Mediated Control

Use the same hypothetical parameters, but now let \(E=0.080\) cells/mL and \(k=25\) mL cell\(^{-1}\) day\(^{-1}\). The added killing rate is

$$ kE=(25)(0.080)=2.00\ \mathrm{day}^{-1},\qquad q=\delta+kE=0.50+2.00=2.50\ \mathrm{day}^{-1}. $$

The immune-adjusted threshold becomes

$$ T_{\mathrm{crit}}=\frac{(10)(2.50)}{(2.0\times10^{-7})(100)} =1.25\times10^6\ \mathrm{cells/mL}. $$

At \(T=1.0\times10^6\) cells/mL, \(\mathcal{R}_T=(2.0\times10^{-7})(100)(1.0\times10^6)/(10\times2.50)=0.800\). The local system therefore predicts decline even though \(T\) is above the target-cell-only threshold of \(2.50\times10^5\) cells/mL. For example, if \(I=100\) cells/mL, \(K_I=100\) cells/mL, \(s_E=0.004\) cells/mL/day, \(a=0.60\) day\(^{-1}\), and \(d_E=0.20\) day\(^{-1}\), then

$$ \frac{dE}{dt}=0.004+(0.60)(0.080)\frac{100}{100+100}-(0.20)(0.080) =0.012\ \mathrm{cells/mL/day}. $$

Thus the specified effector concentration is increasing at that instant in the example model. The calculation illustrates how immune activity can shift control without requiring target cells to have fallen below the target-only threshold; it does not establish that this abstract effector dynamic describes a particular infection.

Diagnostics, Alternatives, and Limitations

Mechanistic comparison requires more than matching a viral-load curve. A practical diagnostic is to examine whether one parameterization can account for viral load, target-cell abundance, and the timing of any proposed increase in immune activity. If the target-cell-only model predicts control only after a large target-cell decline that is not present in the data, its mechanism is weakened. Conversely, a decline aligned with target loss is compatible with target limitation, but does not rule out concurrent immune killing.

The main alternatives can be compared as follows:

FeatureTarget-cell-limited modelImmune-mediated model
Source of changing controlDeclining \(T\) reduces infection input.Increasing \(E\) raises infected-cell loss.
Threshold\(T_{\mathrm{crit}}=c\delta/(\beta p)\).\(T_{\mathrm{crit}}=c(\delta+kE)/(\beta p)\), for fixed \(E\).
Pattern consistent with mechanismViral control follows a substantial fall in target cells.Viral control occurs as effector activity rises, potentially with targets relatively preserved.
Useful additional measurementsTarget-cell abundance over time.Target-cell and effector measurements over time.

When both models fit viral load similarly, that fit alone cannot choose between them. More informative data may include serial target-cell counts, independent measures of effector abundance or activity, and observations before and during the transition to decline. Perturbations that selectively alter target supply or immune activity can be more discriminating than passive trajectories, provided their other effects are considered.

Several failure modes deserve attention. First, an apparent target-cell decrease may reflect sampling or changes in the measured cell population rather than depletion of the true susceptible pool. Second, the immune model may be weakly identifiable if \(k\) and \(E\) are observed only through their product \(kE\); viral data may constrain the total killing rate without determining its two components. Third, the assumed bilinear infection and killing terms may be inadequate if infection saturates or effector action depends on delays. Finally, parameter uncertainty in \(\beta,p,c,\delta\), or target-cell measurements shifts the inferred threshold. A threshold calculated from point estimates should not be treated as a precisely known biological boundary.

Common Mistake

  • Wrong reasoning: “Viral load fell, so the virus must have depleted its target cells.” A decline is compatible with that explanation but also with increased infected-cell loss or other omitted processes.
  • Correction: Compare the timing of viral decline with measured target-cell and immune trajectories. Use the model thresholds as conditional predictions, check units and assumptions, and describe the mechanism as consistent with the evidence rather than proven by viral load alone.
Key takeaway: Target depletion lowers infection input by reducing \(T\); immune-mediated control raises infected-cell loss through \(kE\). Both can lower viral load, but their predicted target-cell and effector patterns differ. Mechanistic interpretation depends on measuring those states, not on viral decline alone.

Check Your Understanding

Try each problem, then compare your reasoning with the solutions.

  1. With \(\beta=1.0\times10^{-7}\) mL virion\(^{-1}\) day\(^{-1}\), \(p=200\) virions cell\(^{-1}\) day\(^{-1}\), \(c=8\) day\(^{-1}\), and \(\delta=0.40\) day\(^{-1}\), calculate the target-cell-only threshold.
  2. For the parameters in Question 1, let \(T=4.0\times10^5\) cells/mL. Is local growth predicted in the target-cell-only model?
  3. In the immune model, what happens to \(T_{\mathrm{crit}}\) when \(kE\) increases, with other parameters fixed? Explain.
  4. A viral load declines while measured target cells remain nearly constant. Does this prove immune control? Give a model-based interpretation.

Solutions: 1. \(T_{\mathrm{crit}}=c\delta/(\beta p)=(8)(0.40)/[(1.0\times10^{-7})(200)]=3.2/(2.0\times10^{-5})=1.60\times10^5\) cells/mL. 2. \(\mathcal{R}_T=(\beta pT)/(c\delta)=(2.0\times10^{-5})(4.0\times10^5)/(3.2)=8.0/3.2=2.50\), so local growth is predicted under the frozen-state assumptions. 3. It increases because \(T_{\mathrm{crit}}=c(\delta+kE)/(\beta p)\); stronger killing permits decline at a larger target-cell concentration. 4. No. The pattern is consistent with immune-mediated control because decline occurs without target depletion, but other mechanisms or model misspecification may also explain it. Additional immune measurements and model checks are needed.

Clinical Biostats