HIV-1 CCR5 gene therapy will fail unless it is combined with a suicide gene.

Highly active antiretroviral therapy (ART) has successfully turned Human immunodeficiency virus type 1 (HIV-1) from a deadly pathogen into a manageable chronic infection. ART is a lifelong therapy which is both expensive and toxic, and HIV can become resistant to it. An alternative to lifelong ART is gene therapy that targets the CCR5 co-receptor and creates a population of genetically modified host cells that are less susceptible to viral infection. With generic mathematical models we show that gene therapy that only targets the CCR5 co-receptor fails to suppress HIV-1 (which is in agreement with current data). We predict that the same gene therapy can be markedly improved if it is combined with a suicide gene that is only expressed upon HIV-1 infection.

HIV infects human T cells by binding to the CD4 receptor and the CCR5 or the CXCR4 co-receptor. HIV-1 primarily targets effector/memory (T) CD4+ T cells because they express CCR512. Individuals carrying a CCR5 deletion mutation (CCR5Δ-32) are resistant to HIV-1 because HIV-1 does not bind efficiently to the target cells13. Recently, an HIV-1 infected individual, the “Berlin patient”, was cured after receiving a hematopoietic stem cell (HSC) transplant from a donor homozygous for the CCR5Δ-32 mutation4. Mathematical, experimental, and clinical studies have demonstrated that targeting the CCR5 co-receptor can reduce the HIV-1 viral load145.

Current ART therapies suppress HIV-1 and have drastically improved the survival of HIV-1 infected individuals, but must be administered lifelong. Gene therapy therefore offers an attractive alternative to ART. Current gene therapies aim to create a population of CD4+ T cells that carry at least partially dysfunctional CCR5 co-receptors. These genetically modified cells are less susceptible to HIV-1 infection6. Recently Tebas et al.1 used zinc finger based gene editing targeting CCR5 gene and modified a fraction (11–28%) of host CD4+ T cells1. They found that genetically modified T cells have a survival advantage over unmodified T cells upon ART interruption and that gene therapy can boost total CD4+ T cell count. However, their data suggest that gene therapy does not suppress the viral load (see Supplementary Table S1) possibly because the genetically modified T cells were still susceptible to HIV-1 infection. Other researchers are targeting the virus by knocking down CCR5 using si/shRNA constructs or by inhibiting viral fusion (using C46 peptides)3789. In summary, most of the current gene therapy trials try to decrease the infection rate by targeting viral entry via the CCR5 co-receptor1011. Because targeting viral entry reduces the susceptibility of a subset of the target CD4+ T cells, CD4+ T cell counts should increase and the viral load should decrease.

Since the data remain inconclusive (Supplementary Table S1), we ask the question: is CCR5 gene therapy expected to reduce the viral load over the long term, and to replace lifelong ART? Using several generic mathematical models we show that the viral load will not necessarily decrease, and propose to include a suicide gene in the genetically modified T cells to develop a much more efficient therapy.

Results and Discussion

Typical HIV infection

To model HIV-1 infection, we considered that the target uninfected T cells are largely self-renewing and model this with a logistic growth for uninfected dividing T cells, with a small constant influx from naive T cells becoming effector/memory cells21213. Uninfected T cells die at rate δ and become infected at rate β. Infected T cells (I) die at rate δ producing virus particles (V) at rate p per infected cell; viruses are cleared at rate c (Fig. 1a and Eqs 1, 2, 3). In short, we wrote a typical model for the main features of HIV-1 infection131415.

Figure 1

Models of the HIV infection.

(a) Schematic of the typical HIV model. The uninfected CD4+ T cells can self-renew (logistic growth), die at rate δ day−1, and get infected at rate β particle−1 day−1. To model the influx from the naive compartment, uninfected T cells additionally get replenished at rate λ cells day−1. Infected cells produce p virus particles cell−1 day−1 and die at rate δ day−1. Viruses are cleared at rate c day−1 (see section Models). (b) The infection rate in the typical HIV model (β) was reduced, and is expressed as a percentage of the initial infection rate (β0). Only once the infection rate is reduced below a certain threshold (β < β*; vertical gray dashed line) will the viral load start to decrease. Several HIV-1 models let uninfected T cells be replenished at a constant rate from thymus, bone marrow, naive, and memory T cells131423. A quantitavily similar model with a constant rate of replenishment (i.e. λ = 14.3 and r = 0 in Eq. 1) also has the property that β has to be decreased markedly (here more than 50%) to have a significant effect on the viral load14 (black dashed line). Modeling one μl of blood, the following parameter values were used: r = 0.06 day−1, K = 1500 cells, δ = 0.02 day−1, β0 = 3.6 × 10−6 particle−1 day−1, δ = 1 day−1, p = 2.14 × 104 virus particles cell−1 day−1, c = 23 day−1 (see section Models). The values of λ, r, and K were chosen to have 1000 T cells per μl blood in the virus-free steady state. (c) For the GT model, we model two populations of uninfected T cells (normal and genetically modified) sharing logistic growth. Normal uninfected T cells additionally get replenished at rate λ cells day−1. Normal and genetically modified T cells get infected at rate β and β particle−1 day−1 (where β ≤ β), respectively. Infected cells die at rate δ and δ day−1.

To derive analytical solution for the steady state, we make the reasonable assumption that the uninfected T cell population is largely self-renewing (i.e. λ = 0 in Eq. 1)212 to obtain the following steady state:

Our typical HIV model shows that decreasing the infection rate (β) increases the uninfected T cell count (Eq. 4a; Supplementary Fig. S1). However, decreasing the infection rate (β) can either increase or decrease the viral load (decreasing β increases the first term cr/pβ and decreases the term in the parentheses for Eqs 4b and 4c). Thus, if the infection rate does not decrease substantially gene therapy can actually have either no effect or even marginally increase the viral load (note the slight increase in the viral load for β ≥ β*; Fig. 1b). Obviously decreasing the infection rate markedly (β < β*; i.e. a nearly complete suppression of CCR5) should always decrease the viral load (Fig. 1b) as the infection can no longer sustain itself14. Note that we assumed λ = 0 in Eq. 1 to present a simple analytical solution, and that similar results are obtained when λ > 0 (Eq. 5 and Supplementary Fig. S4a). For all subsequent analysis, we use the full model with λ > 0.

Inefficient HIV-1 gene therapy

Current gene therapies generate modified T cells with either homozygous or heterozygous disruption of the CCR5 gene to reduce the expression of the CCR5 co-receptor on host T cells16. However, genetically modified T cells remain susceptible to infection by HIV-1 unless both alleles of CCR5 gene are completely disrupted3. In the example of Fig. 1b we find that the infection rate has to be reduced by more than 58% to reduce the viral load. Given that gene therapy cannot ensure complete disruption of the CCR5 gene, it is unclear if this substantial decrease in the infection rate is achieved by current gene therapies.

To study the efficacy of gene therapy, we make a gene therapy (GT) model from the typical model by allowing for two T cell populations: normal host T cells (T) and genetically modified T cells (T) (described in Fig. 1c and section Models). The GT model suggests that CCR5 gene therapy should always increase the total T cell count (T + T: Fig. 2a) but will always decrease the normal T cell count (T: Fig. 2b). The competition between the normal T cells and the genetically modified T cells can thus reduce the diversity of the T cell repertoire, as there will be fewer normal T cells.

Figure 2

The effects of gene therapy (i.e. decreasing the infection rate β) and the introduction of suicide gene (i.e. increasing the death rate, δ, of the genetically modified T cells) on the steady state T cell count and the steady state viral load.

We kept the infection rate of unmodified target cells constant (i.e. β = 3.6 × 10−6) and varied the infection rate of the genetically modified cells (β). The value of β is given as a percentage of the β value, where 100% means that β = β and 20% means that β = 0.2β. The effect of a change in the death rate of genetically modified cells (δ) is shown by different lines. Black lines represent the current CCR5 gene therapies (with δ = δ). The value for δ was changed as a fold increase of δ. (a) Decreasing the infection rate (β) of genetically modified cells increases the total T cell count (black line). (b) Decreasing the infection rate (β) of genetically modified cells decreases the normal uninfected T cell count (black line). Decreasing the infection rate (β) below a threshold rescues the normal uninfected T cell count. (c) Decreasing the infection rate (β) of genetically modified cells can increase the viral load slightly for a mildly effective gene therapy. A suicide gene (to increase δ) expressed in infected genetically modified cells increases the total T cell count (a), rescues the normal uninfected T cell count (b), and decreases the viral load (c) indicated by different colors (see legend in Fig. 2a). To initialize the GT model, we replaced 10% of the T cells by T cells in the infected steady state and we ran the model until the steady state is approached. Parameters: thymic influx for T (λ) = 1 cells μl−1 day−1, logistic growth parameter for T and T (r) = 0.057 day−1, and death rate of I (δ); the other parameters remain the same (Fig. 1). The values for r and λ were chosen to have T = 1000 cells per μl blood as virus-free steady state using parameter values given in222324. (d) The effect of variation of β along with δ in GT model (see section Models). The colors indicate the steady state viral load.

The GT model suggests that the viral load is only guaranteed to decrease if the infection rate of the genetically modified T cells falls below a certain threshold i.e., less than 40% of that of the normal T cells (β < 0.4β; see black line in Fig. 2c). Similar to the results shown in Fig. 1b, gene therapy does not decrease the viral load (see Fig. 2c,d, 0.4β ≤ β) if the infection rate of the genetically modified T cells does not decrease substantially.

Gene therapy cannot guarantee homozygous CCR5 gene disruption since the majority of the modified T cells will be heterozygous for the CCR5 gene3. The infection rate of CD4+ T cells from individuals genetically heterozygous for the CCR5Δ-32 mutation (where all cells have the CCR5 disruption) is reduced by a maximum of 65%, making them slow progressors3. Unlike the individuals carrying the CCR5Δ-32 mutation, gene therapy can modify only a fraction of T cells (11%–28%) most of which will still be susceptible to HIV-1 infection (because they will be heterozygous for the CCR5 gene). Since only a fraction of T cells is modified, current gene therapies are not expected to decrease the viral load. For current CCR5 gene therapies to be successful and to control viral transmission, the infection rate has to be strongly reduced, which can only be achieved if almost all cells that are genetically modified become homozygous for the CCR5 deletion. This perfectly explains why introducing genetically modified T cells1 failed to decrease the viral load in the absence of ART, except in one patient who was heterozygous for the CCR5Δ-32 mutation (see Supplementary Table S1).

Modifications to current gene therapy

Novel gene therapies targeting two or more steps in the HIV-1 life cycle have been proposed, and some are currently being tried10. Since we find that targeting CCR5 alone probably cannot replace ART, we studied whether or not combining CCR5 gene disruption with other factors could result in a more effective gene therapy. These other factors are: reducing the burst size (p) or increasing the virus clearance rate (c). Since these have the same effect as reducing the infection rate (β), they also result in an ineffective therapy (Supplementary Fig. S2).

HIV-1 suicide gene therapy

Since the uninfected T cell count increases linearly with the death rate, δ, of the infected cells (Eq. 4a) but the viral load decreases with δ (Eqs 4b, 4c, 5b, and 5c), we quantified the effect of increasing the death rate of the infected cells. In practice, the death rate of the infected genetically modified T cells can be increased by introducing a suicide gene. Introducing a suicide gene and thus increasing the death rate of the infected genetically modified T cells (δ) without CCR5 gene therapy (β = β) does not affect the viral load as the genetically modified T cells without CCR5 disruption go extinct (Fig. 2a).

Remarkably, introducing a suicide gene with CCR5 gene therapy strongly suppresses the viral load (Fig. 2c). For example, a 50% decrease in the infection rate of the genetically modified T cells (β) along with a 100% increase in the death rate of infected genetically modified T cells reduces the viral load by 50%. A 50% decrease in the infection rate (β) without changing the death rate (i.e. without a suicide gene) does not decrease the viral load (it even slightly increased the viral load in Fig. 2c). Moreover, introducing a suicide gene in CCR5-disrupted genetically modified T cells increases the total T cell count (Fig. 2a) and rescues the normal T cell count (Fig. 2b). Thus the expression of a suicide gene in genetically modified T cells upon infection can markedly increase the effectiveness of CCR5 gene therapies (Fig. 2c). Moreover, HIV-1 has high mutation rates16 and can adapt alternate mechanisms (like receptor switching) to infect the genetically modified T cells. Introduction of suicide gene in CCR5 modified cells will provide protection against HIV-1 variants that evolve usage of alternate co-receptors because HIV-1 variants which infect using alternate receptors will activate the suicide gene in modified T cells.

The GT model was designed as a simple model of HIV-1 dynamics to help us gain insight into various features of CCR5 gene therapy17. To test its robustness, we studied several extensions allowing: for latent cells (Supplementary Fig. S3a); for cytotoxic T cell (CTL) responses (Supplementary Fig. S3b); for genetic modification of hematopoietic stem cells (Supplementary Fig. S3c and d); and for different stages of infected cells (i.e. an eclipse phase; Supplementary Fig. S3e). Irrespective of the model used, we predict that incomplete gene therapy targeting only the CCR5 receptor fails to reduce the viral load. However, a concomitant increase in the death rate of infected cells always resulted in an effective therapy.

In the GT Model, we considered that growth of target uninfected T cells is largely via self-renewal of T cells with a small influx from the source (i.e. naive T cells becoming effector/memory cells). To test how this affects our results, we performed a parameter sensitivity analysis, and found that when contribution from the source is less than 10%, genetically modified (T) cells can survive and form a major fraction of the T cell population. Hence the introduction of a suicide gene along with CCR5 gene therapy can be beneficial only when the contribution from the source is less than 10% (Supplementary Fig. S4a and b). When the contribution from the source is high, T cells either go extinct or are present in such low numbers that any form of CCR5 based gene therapy will not be effective. Since the CCR5 receptor is found on T cells that are formed largely via self-renewal the contribution from the source should be low. Moreover, Tebas et al.1 found that the genetically modified cells can survive for long time, indicating that the source contribution is low. A similar parameter sensitivity analysis for the rate of infection (β) of normal T cells confirmed that the introduction of suicide gene always reduces the viral load and increases the T cell count compared to the CCR5 gene therapy without suicide gene (Supplementary Fig. S4c and d).

Conclusions

A great deal of effort is being made to improve HIV-1 gene therapy after seeing the “curative” effects of the CCR5Δ-32 mutation in the “Berlin patient”. Mitsuyasu et al.8 and Tebas et al.1 achieved an increase in CD4+ T cell counts within patients as predicted by Fig. 2a and Supplementary Fig. S1; however, neither gene therapy significantly reduced the viral load. Different suicide genes (e.g., HSV1-TK and iCasp9) are known to cause cell death when expressed at high concentration18. HIV-1 infected cells can be killed in a specific manner if a suicide gene (HSV1-TK or iCasp9) is successfully incorporated with the HIV-1 promoter and transactivation response element sequences181920. Our models show that the therapeutic effect of targeting viral entry is markedly increased when the infected modified T cells undergo rapid cell death when they are infected. It therefore seems feasible and promising to combine a suicide gene with zinc finger or si/shRNA based constructs to knock down CCR5.

Methods

Mathematical models for HIV-1 gene therapy

We constructed a suite of models, considering various extensions like latency, cytotoxic T cell (CTL) responses, different stages of infected cells, to test the predictions of our simple model. We found that in all models CCR5 gene therapy is predicted to fail unless it is combined with a suicide gene. Mathematical models were simulated in MATLAB R2014b (http://www.mathworks.com/products/matlab/).

For the typical HIV model, we considered that T cells get replenished via logistic growth (dividing T cells) and have a small contribution from another compartment (e.g. naive T cells). We consider that upon infection both uninfected (T) and infected T cells (I) compete in the logistic growth term (Eqs 1, 2, 3 and Fig. 1a). To obtain the simple steady state expression presented in the text (Eq. 4), we assumed that and . If we relax these assumptions the results do not change qualitatively as shown by the steady state expression obtained for the full typical HIV model:

The GT model considers two populations of uninfected and infected T cells: normal uninfected (T) and infected (I) T cells, and genetically modified uninfected (T) and infected (I) T cells. T cells get replenished via logistic growth and have a small contribution from another compartment (e.g. naive T cells). We allowed all actively dividing cells (T, T, I, and I) to compete in the logistic growth of T and T cells. The virus particles (V) are produced from I and I (Eqs 6, 7, 8, 9, 10 and Fig. 1b). Parameters are described in the legend of Fig. 2:

In the Latent model, we changed the dynamics for the infected cells by adding compartments for latently infected T cells (normal: L and genetically modified: L). The equations for T cells, T cells, and V (Eqs 6, 7, and 10) remain the same, and the equations for infected cells are replaced by Eqs 11, 12, 13, 14. We added the following parameters: f = 0.995, δ = 0.0001, and α = 0.1, the other parameters remained the same:

To model the effect of CTL killing, we consider that a fraction of infected T cells are killed by CTLs (Eqs 15 and 16, where δ = 0.9, δ = 0.1, δ = d × δ) keeping equations 6, 7, and 10, and the remaining parameters the same (Fig. 2):

If CTLs are largely responsible for killing the infected T cells (δ = 0.9 and δ = 0.1), a suicide gene in I cells should kill the infected cells much more rapidly to curb HIV transmission (Supplementary Fig. S3b). This is not unreasonable because Marin et al.18 showed that suicide genes can indeed cause rapid cell death. Adding a CTL population with a competitive saturating function for CTL growth, and a mass-action term for CTL-mediated killing (Eqs 17, 18, 19) enables that the viral load depends on many parameters of the model21. Considering such a CTL model does not change the results qualitatively (data not shown):

To study hematopoietic stem cell (HSC) transplantation as an alternative to conventional gene therapy with mature T cells, we allow HSC cells to contribute to the replenishment of genetically modified T cells. We replaced the T growth dynamics by equation 20 keeping Eqs 6, 8, 9, 10 and other parameters the same:

Althaus et al.22 recently described a within-host HIV model considering different sub-classes of HIV-1 infected cells, namely, latently infected cells (L1 and L2), persistently infected cells (M1 and M2), defectively infected cells (D), newly infected cells (I1; without integrated virus), infected cells with integrated proviral DNA (I2), infected cells with increasing HIV transcriptional activity (I3, I4, and I5) and virus producing cells (I6). The I1, I2, I3, I4, and I5 sub-classes of HIV-1 infected cells describe the intracellular eclipse phase of HIV-1 dynamics. We extended this model to allow for a population of normal T cells (i = 1) and genetically modified T cells (i = 2) (Eqs 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33), and extend the model to have logistic growth for T1 and T2 cells:

The following parameters were kept the same22: δ = 0.0165, β1 = 1.35 × 10−6, , δ = 8.3 × 10−5, δ = δ = δ, f = 0.141, f = 0.0035, f = 0.25, γ1 = 3, γ2 = 6, γ3 = γ4 = γ5 = 12, σ1 = 0.0789, σ2 = 1, κ1 = 0.103, κ2 = 1, α = 2.67 × 10−3, N = 2.14 × 104, and c = 23. We added the following parameters to make the virus-free steady state equivalent to that of the Althaus22 model: λ1 = 1, λ2 = 0, r = 0.0856, and K = 5000. To model the expression of suicide genes, we increased the death rate of genetically modified infected cells expressing viral RNA to , and .

Additional Information

How to cite this article: Pandit, A. and de Boer, R. J. HIV-1 CCR5 gene therapy will fail unless it is combined with a suicide gene. Sci. Rep. 5, 18088; doi: 10.1038/srep18088 (2015).