Mathematical modeling of escape of HIV from cytotoxic T lymphocyte responses

Vitaly V. Ganusov,Richard A. Neher,Alan S. Perelson
DOI: https://doi.org/10.1088/1742-5468/2013/01/P01010
2012-07-24
Abstract:Human immunodeficiency virus (HIV-1 or simply HIV) induces a persistent infection, which in the absence of treatment leads to AIDS and death in almost all infected individuals. HIV infection elicits a vigorous immune response starting about 2-3 weeks post infection that can lower the amount of virus in the body, but which cannot eradicate the virus. How HIV establishes a chronic infection in the face of a strong immune response remains poorly understood. It has been shown that HIV is able to rapidly change its proteins via mutation to evade recognition by virus-specific cytotoxic T lymphocytes (CTLs). Typically, an HIV-infected patient will generate 4-12 CTL responses specific for parts of viral proteins called epitopes. Such CTL responses lead to strong selective pressure to change the viral sequences encoding these epitopes so as to avoid CTL recognition. Here we review experimental data on HIV evolution in response to CTL pressure, mathematical models developed to explain this evolution, and highlight problems associated with the data and previous modeling efforts. We show that estimates of the strength of the epitope-specific CTL response depend on the method used to fit models to experimental data and on the assumptions made regarding how mutants are generated during infection. We illustrate that allowing CTL responses to decay over time may improve the fit to experimental data and provides higher estimates of the killing efficacy of HIV-specific CTLs. We also propose a novel method for simultaneously estimating the killing efficacy of multiple CTL populations specific for different epitopes of HIV using stochastic simulations. Lastly, we show that current estimates of the efficacy at which HIV-specific CTLs clear virus-infected cells can be improved by more frequent sampling of viral sequences and by combining data on sequence evolution with experimentally measured CTL dynamics.
Populations and Evolution
What problem does this paper attempt to address?
### Problems the paper attempts to solve The problems that the paper attempts to solve are: **How does the human immunodeficiency virus (HIV) establish a chronic infection under a strong cytotoxic T - lymphocyte (CTL) immune response?** Specifically, the paper focuses on the following aspects: 1. **HIV escape mechanisms**: - HIV can evade CTL recognition by mutating its proteins. This escape mechanism enables the virus to persist under the pressure of the immune system. - Typically, an HIV - infected patient will generate 4 - 12 CTL responses against specific regions of viral proteins (called epitopes). These CTL responses exert strong selective pressure on the parts of the viral sequence encoding these epitopes, causing the virus to mutate to avoid CTL recognition. 2. **Experimental data and mathematical models**: - The paper reviews the experimental data on the evolution of HIV under CTL pressure and discusses the mathematical models used to explain this evolution. - The researchers point out that the methods and assumptions for estimating the strength of CTL responses to specific epitopes will affect the effectiveness of the model in fitting the experimental data. 3. **Attenuation of CTL responses**: - Allowing CTL responses to decay over time can improve the fit of the model to the experimental data and provide a higher estimate of the HIV - specific CTL killing efficiency. 4. **Simultaneous estimation of multiple CTL populations**: - A new method is proposed to simultaneously estimate the killing efficiencies of multiple CTL populations against different HIV epitopes using stochastic simulation. 5. **The influence of data sampling frequency and depth**: - Current estimates of the efficiency of HIV - specific CTLs in clearing infected cells can be improved by sampling viral sequences more frequently and combining sequence evolution with experimentally measured CTL dynamics data. ### Key formulas 1. **Viral kinetics model**: \[ \frac{dT}{dt}=\lambda - dT-\beta T(V_w + V_m) \] \[ \frac{dI_w}{dt}=(1 - \mu)\beta T V_w-(\delta + k)I_w \] \[ \frac{dI_m}{dt}=\mu\beta T V_w+\beta T V_m-\delta I_m \] \[ \frac{dV_w}{dt}=p_w I_w - c_V V_w \] \[ \frac{dV_m}{dt}=p_m I_m - c_V V_m \] 2. **Dynamics of escape variant frequencies**: \[ \frac{dz(t)}{dt}=\mu r+z(t)(k - r(c - \mu)) \] \[ z(t)=z_0 e^{\epsilon t}+\frac{\mu r}{\epsilon}(e^{\epsilon t}-1) \] 3. **Estimation of escape rates**: \[ f(t)=\frac{f_0}{f_0+(1 - f_0)e^{-\epsilon t}} \] \[ f(t)=\frac{z(t)}{1 + z(t)}=\frac{f_0}{f_0+(1 - f_0)e^{-\epsilon t}\times(1 - e^{-\epsilon t})} \] 4. **Likelihood function**: \[ L=\prod_i\binom{n_i}{k_i}f(t_i)^{k_i}(1 - f(t_i))^{n_i - k_i} \] \[ \log L=\sum_i[k_i\ln(f(t_i))+(n_i - k_i)\ln(1 - f(t_i))]