Skip to main content
  • Research article
  • Open access
  • Published:

Understanding the HIV coreceptor switch from a dynamical perspective



The entry of HIV into its target cells is facilitated by the prior binding to the cell surface molecule CD4 and a secondary coreceptor, mostly the chemokine receptors CCR5 or CXCR4. In early infection CCR5-using viruses (R5 viruses) are mostly dominant while a receptor switch towards CXCR4 occurs in about 50% of the infected individuals (X4 viruses) which is associated with a progression of the disease. There are many hypotheses regarding the underlying dynamics without yet a conclusive understanding.


While it is difficult to isolate key factors in vivo we have developed a minimal in silico model based on the approaches of Nowak and May to investigate the conditions under which the receptor switch occurs. The model allows to investigate the evolution of viral strains within a probabilistic framework along the three stages of disease from primary and latent infection to the onset of AIDS with a a sudden increase in viral load which goes along with the impairment of the immune response. The model is specifically applied to investigate the evolution of the viral quasispecies in terms of R5 and X4 viruses which directly translates into the composition of viral load and consequently the question of the coreceptor switch.


The model can explain the coreceptor switch as a result of a dynamical change in the underlying environmental conditions in the host. The emergence of X4 strains does not necessarily result in the dominance of X4 viruses in viral load which is more likely to occur in the model after some time of chronic infection. A better understanding of the conditions leading to the coreceptor switch is especially of interest as CCR5 blockers have recently been licensed as drugs which suppress R5 viruses but do not seem to necessarily induce a coreceptor switch.


Although many studies aim to understand the evolution of the human immunodeficiency virus (HIV) in its host [13] there are still a lot of open questions related to the mechanisms driving the intra-host dynamics. One of these puzzles is associated with the HIV coreceptor switch [4]: To facilitate cell entry and subsequent replication HIV binds to the cell surface molecule CD4 as well as a chemokine coreceptor, most commonly the CCR5 or CXCR4 coreceptors. In most patients viruses that use the CCR5 coreceptor (R5 viruses) dominate in early stages of disease. This preference changes however in about 50% of patients during the course of disease towards viruses using the CXCR4 coreceptor (X4 viruses). This switch in coreceptor usage is associated with a worsened prognosis which makes a better understanding of the switch dynamics of direct clinical importance. At the same time there are several explanatory hypotheses but not yet a conclusive understanding of the underlying processes in terms of empirical evidence or model support. The transmission mutation hypothesis [4] assumes that R5 viruses are favoured in transmission and are in consequence more often found in early infection. X4 viruses with higher fitness can emerge from R5 viruses via intermediate mutants of lower fitness which will finally become dominant [4]. While there is some evidence for a fitness loss of intermediate mutants [57], this hypothesis strongly relies on the assumption that X4 viruses can hardly be transmitted in infectious doses as these would otherwise immediately become dominant due to their assumed selectional advantages in the host. There might be some strain selection during transmission but there is evidence that both R5 and X4 strains can be transmitted in infectious doses, however, with X4 viruses seemingly having a selectional disadvantage thereafter [8, 9]. The effective replication rate of X4 viruses exceeds R5 replication in some in vitro assays [10], there is, however, no conclusive picture of the in vivo situation with different environmental constraints [4, 11]. Taken this together [4, 1214], it is highly questionable whether a fragile setting relying solely on the transmission mutation hypothesis adequately describes the major processes in place or robustly represents the observed dynamics.

An alternative and less fragile hypothesis assumes that environmental conditions in the host change during the course of disease in favour of X4 viruses. Such changes may be the consequence of immune pressure or the availability of adequate target cells for replication in terms of (co-)receptor expression or replication efficiency by stage of cell development (naive/memory cells) [4, 15]. In a recent study Sguanci and coworkers [16] investigate the changing environment by assuming that the susceptibility of target cells (or infectivity of viruses) increases in favour of X4 cells as soon as these emerge in an auto-catalytic fashion, however, staying with the assumption that the initial infection is with R5 virus only. Earlier work by Wodarz et al. [17] and Callaway et al. [18] allows more realistically for a mixed infection and investigates how the shift in target cells and HIV specific T cells can dynamically induce a shift from R5 to X4 viruses. In a similar fashion, Ribeiro et al. investigate how the coreceptor switch can be facilitated due to preferential replication of X4 and R5 viruses in naive and memory T cells, respectively, with differences in efficiency of viral reproduction [15]. These approaches support the hypothesis of environmental change as a driving force to the coreceptor switch but do not consider the stochastic nature of virus evolution and the emergence of new strains and viral diversity which have been associated with the coreceptor switch. While the relative importance of the above hypotheses for the coreceptor switch can ultimately only be answered experimentally we will here focus on the more robust hypothesis of environmental change for further analysis. We want to identify key processes resulting in the observed phenomenon by a restriction to a minimal model representing the dynamics of the immune system and HIV [19]. Therefore, we develop an approach based on the well studied model of Nowak and May [20] that allows to study the dynamical processes associated with the HIV coreceptor switch in combination with viral evolution. While we do not go yet into the details of cellular subpopulations whose roles are empirically still under debate we rather take advantage of the relative simplicity of the model to study the effects of the stochastic nature of HIV evolution neglected in the models outlined above. The stochastic framework will allow to analyse the distribution of evolutionary pathways associated with the coreceptor switch as well as survival times in a more systematic way.

It has specifically become of major importance to understand the processes leading to the dominance of X4 viruses since new drugs have been licensed that block the CCR5 receptor [12, 2123] and in consequence can shift the selective advantage towards X4 viruses [24], although not necessarily inducing a coreceptor switch. Considering the association between X4 viruses and a worsened prognosis [4, 25] it is of vital importance to understand whether X4 can only flourish in the environment of a patient weakened by side effects of a chronic infection or whether X4 viruses themselves destabilise the patient and cause the worsened prognosis [13].

Results and Discussion

The Model

Following the ansatz of Nowak and May [20, 26] we study the time evolution of viral load v i of a set of n strains (i {1, ...,n}). These are controlled by a specific immune response with strength x i and as well as a cross-reactive immune response with strength z. Their time evolution is determined by the following set of equations:


with variables and parameters as summarised in table 1.

Table 1 Variables and parameters in equations (1)

The viral load of strain i grows at a rate r i and is diminished by the specific and cross-reactive immune response at a rate -p i x i - q i z. Vice versa, a specific immune response is stimulated by each strain as is the cross-reactive immunity by the total viral load. Immune responses decay at a rate b and further at a rate proportional to the viral load which is to account for the fact that HIV infected T helper cells are depleted (and in consequence impair the immune response mediated by other T and B cells). Among the 2npossible equilibrium solutions of a set of equations (1) for n strains at most one is stable [20], that is, for any number of strains the viral load either diverges or converges to a well defined equilibrium viral load. Strains are distinguished by their sensitivity towards the specific immune response, i.e. several related genomic sequences might represent the same epitope and consequently correspond to one strain in the model. The model does not yet address the question of different immune cell populations as for example CD4 vs. CD8 cells, naive vs. memory cells or resting vs. replicating cells. The distinction between x i and z takes only into account that there is a strain specific and an unspecific component of the immune response. The latter can be a composition of the response towards conserved parts of the virus, of cross-reactivity from earlier specific immune responses or innate immunity. Bare of any strain-specific responsiveness, z is a measure for the general activation of the immune system going along with an overall turnover of CD4 cells (irrespective of their specificity) seen in HIV infections [27, 28].

To focus on the evolution of R5 and X4 viruses in a patient we distinguish only between two parameter sets corresponding to either type of virus as listed in table 2. The values of initial viral growth rates are estimated from [29, 30]. However, viral growth rates or more generally viral fitness as the balance between viral replication and decay is not an intrinsic feature of the virus but is only well defined in the context of the viruses environmental conditions. Differences in viral fitness may therefore arise among viruses and over time with a change of environmental conditions such as the availability of target cells or the strength of an immune response. Here, we assume that X4 viruses are better recognised by the specific immune response than R5 viruses [31] and increase their replication rate proportional to the stimulation and activation of the immune system [15, 27, 28] represented by the cross-reactive immunity z. While the cumulative immune activation z will likely impact on both R5 and X4 viral replication [15, 32] we are only interested in the increase of X4 replication over R5 replication and therefore keep the R5 replication rate fixed. The cumulative immune activation z increases slightly during the course of disease and breaks only down with the collapse of the immune system at the onset of AIDS. The chosen parameter set corresponds to a situation in which the viral load can initially be controlled by the immune response, i.e. a situation in which limitations in target cell supplies or other saturation phenomena do not have to be taken into account. The model's dynamics is robust with respect to the choice of further parameters as long as they correspond to the model's regime of HIV-like dynamics (cf. equation (4)).

Table 2 Parameter settings

The model system is initialised with a strain of R5 virus and a strain of X4 virus and evolves according to equations (1). A new mutant arises at a hazard rate p m v, i.e. at each integration step of length dt with probability . In the current setting, X4 and R5 mutants are assumed to arise with equal probability irrespective of the composition of the viral load due to the close genetic neighbourhood of these viral subtypes [33, 34] (for further discussion cf. Additional file 1, Figure S1). After the emergence of each new mutant the system reaches a new equilibrium with nR5 R5 virus strains and nX4 X4 virus strains and viral load


The system is stable, i.e. there is a stable equilibrium solution only if

  1. 1.

    ru <kq, i.e. the viral load is fully controlled by the cross reactive immune response z alone or

  2. 2.

that is, only as long as the numbers of R5 and X4 virus strains do not exceed the critical numbers given by condition (3) the immune system can control the infection.

We are specifically interested in dynamical regimes of the model which remain only stable for finite numbers of R5 and X4 virus strains (as opposed to case 1.) but do not show immediate divergence (cpR5 - (ru - kq) > 0) which corresponds to


The first inequality in (4) implies that the cross-reactive immune response cannot fully control the viral infection as can be seen from equation (2). This means that each viral strain needs to be suppressed both by a specific and a cross-reactive immune response. The second inequality implies that the immune system can control only a limited number of strains, i.e. the limitations of the cross-reactive immune response can be compensated by the specific immune response if the critical number of strains given by equation (3) is not exceeded. As soon as this number of strains is exceeded this will lead to an uncontrolled increase in viral load in the model corresponding to the onset of AIDS in the patient [20]. The detailed dynamics in this regime are not longer described by equations (1) as saturation effects have to be considered due to limited resources.

The emergence of viral strains in a patient that accumulate to the critical number at the onset of AIDS can also be understood in a probabilistic framework in order to study survival distributions. Therefore we derive the probability P(nR5, nX4, t) that a patient harbours nR5 strains of R5 virus and nX4 strains of X4 virus at time t. If the viral load equilibrates much faster after the emergence of a new mutant than the next mutant arises we can assume that new mutants will approximately arise at a rate proportional to the equilibrium viral load v(nR5, nX4). In the regime where and v(nR5, nX4) <v max the evolution of the number of viral strains is then described by a Markov process with the master equation [35]


For those states in which the critical number of strains given by equation (3) or a maximal viral load v max is exceeded, i.e. at the onset of AIDS, the last term in equation (5) is defined to vanish making these states absorbing states of the Markov process. v max limits the diverging viral load v(nR5, nX4) as the total viral load in the patient is limited. The distribution of times until absorption, i.e. the survival distributions, are described by phase type distributions [36].

The course of disease

First insight into the evolutionary dynamics of R5 and X4 viruses can be gained by investigating the time course of viral load attributed to each subtype in the model. Such a typical course of disease after a mixed infection with R5 and X4 viruses in the model is shown in Fig. 1.

Figure 1
figure 1

The course of disease in the model. Viral load of R5 viruses (blue) and X4 viruses (red) in a simulation of the model described by equations (1) and parameters as in table 2. The viral load is given in arbitrary units of viral load (VL).

The system shows a period of low viral load after the initial phase of disease with is interspersed with small outbreaks when a new viral mutant occurs - until eventually the viral load rises again. The system shows a quasi-stable behaviour as long as inequality (3) can be fulfilled for a non-zero but limited number of strains, i.e. if inequality (4) holds.

X4 viruses are underrepresented in viral load in the early stages of disease as the model assumes that they are more strongly suppressed by the specific immune response than R5 viruses (pX4 > pR5) [31] while having initially identical growth rates [11]. On the other hand the X4 replication rate in the model increases proportional to the overall activation of the immune system represented by the cross-reactive immune response. This shifts fitness advantages from R5 viruses to X4 viruses resulting in the phenomenology of a coreceptor switch.

The model captures the association between ongoing immune activation and the apparent shift in evolutionary advantages from R5 to X4 viruses during an HIV infection. There is however still controversy about the mechanisms underlying this association. Earlier assumptions about a shift in target cells favouring X4 over R5 viruses [4] are now questioned [37]. An alternative hypothesis links the efficiency of viral production to the division rate of target cells [15]. With the differential increase in target cell division rates during chronic infection [32] viral replication rates are differentially increased which eventually shifts the selectional advantage from R5 viruses to X4 viruses [15, 27, 28]. Another interpretation of the model equations is that X4 viruses experience a stronger specific immune response than R5 viruses but that they face a reduced cross-reactive immune response. In consequence X4 viruses will be more difficult to control as the number of strains grows. Any of these hypotheses allows for a dynamical change in environmental conditions in the host that become more favourable for X4 viruses after chronic infection while selection is in favour of R5 viruses in early stages of disease. Although, there is some evidence [4, 11, 15, 31, 32] in accordance with the model's assumptions on parameter settings it is difficult to get conclusive results from currently available empirical data. Many studies with one or only a few HIV or SIV cases point in different directions [5, 6, 7, 38]: Any comparative measurement of growth rates and the effect of an immune response among strains in vivo assesses quantities derived from a large variety of underlying processes that might vary among cases. Therefore, a more conclusive picture can only be expected with the integration of more data describing the interactions between HIV and the immune system in many patients. The observation that some viral strains may be present in a patient for a long time before they appear in detectable numbers favours the idea of a dynamical change in the environmental conditions in favour of X4 viruses [39, 40] at a higher level of abstraction. Reflecting these findings, the current model provides a coarse grained picture of possible scenarios for environmental shifts in the patient for further investigation. While details of these dynamics will have to be addressed within a refined model the current model has the advantage to be accessible to further analysis of the basic underlying stochastic processes as demonstrated in the following sections.

The routes of evolution

The probabilistic framework complementarily allows to study the probability to find a setting with nR5 R5 strains and nX4 X4 strains in a patient at any time during the course of disease - including the probability to have reached the final stage of disease, i.e. having exceeded the critical number of strains.

Fig. 2 shows that a patient is most likely to be found either with a few strains or in a situation in which the final stage of disease has already been reached. This is due to the fact that the viral load and in consequence the probability to generate new mutants increases with a growing burden of viral strains. As a result the mean residence time with a given number of strains will decrease as the number of strains grows, i.e. the course of disease speeds up with the accumulation of viral strains. The distribution of strain compositions shown in Fig. 2 is asymmetric with respect to the numbers of R5 and X4 strains in the presence of a moderate number of strains, i.e. in a stage of disease that is not yet progressed. A moderate number of strains will not result in strong stimulation of the immune system and X4 replication and viral load remain at a low level compared to the contribution of R5 strains. X4 strains will consequently not drive the course of disease as much ahead as R5 strains in early stages of disease. Therefore the model predicts that non-progressors may even have more X4 strains than R5 strains without experiencing a coreceptor switch. The intricate dependency of the coreceptor switch on the numbers and composition of strains is analysed in detail in the following section.

Figure 2
figure 2

The distribution of strains. The probability to find nR5 R5 strains and nX4 X4 strains in a patient at 1200 days since infection in the model with parameters as in table 2 and v max = 500 units of viral load, the final stage of disease (absorbing state) is marked by the accumulation of probability density at the diagonal border.

The coreceptor switch

The coreceptor switch is a phenomenological observation for which a corresponding process has to be characterised in our model. As the simulations start with a dominance of R5 viruses due to the bias in immune recognition (i.e. pX4 > pR5) we define that a coreceptor switch has occurred if the viral load is dominated by X4 viruses at the onset of AIDS which is determined in the model by either exceeding a viral load v max or the critical numbers of strains nR5 and nX4 (cf. Fig. 2). From


it can be seen that it depends on the number and composition of viral strains in the system which type of virus dominates, i.e. when the trade-off between replication rate and sensitivity to the immune response favours X4 viruses. X4 viral load dominates in a patient who has not yet reached the terminal stage of disease if


which corresponds to the area shaded in red in Fig. 3. In early infection X4 strains do not impose a high viral load and the number of X4 strains required to dominate the viral population in the presence of R5 viruses is usually not found in a patient. This picture changes in late infection where a smaller fraction of X4 strains among all strains can be sufficient for dominance in in viral load. This can be seen in Fig. 3 which shows graphically what strain composition leads to a dominance in X4 viral load (red shaded area) and when the system is destabilised indicating the onset of AIDS (grey shaded area). The saturation behaviour of equation (7) implies that dominance of X4 viruses in viral load can be attained with a minority of X4 strains as soon as

Figure 3
figure 3

X4 dominance and the onset of AIDS. The figure shows how the dominance of X4 viral load and the onset of AIDS depends on the number of X4 and R5 strains present in the system according to equations (3) and (7), for parameters cf. table 2.

Note that this situation is only observed in the latent phase defined by equation (3) if X4 viruses gain sufficiently fast replicative fitness, i.e. .

Within the model the occurrence of a coreceptor switch depends on the probabilistically chosen evolutionary path in the patient (cf. Fig. 2) - i.e. whether a path is chosen that leads directly to destabilisation and AIDS or whether the final stage is reached after a switch in coreceptor usage. The probabilities for either choice are shown in Fig. 4. It shows that for any newly infected patient the probabilities of eventually reaching the final stage of disease with dominance of R5 viruses or X4 viruses are of the same order of magnitude with the parameter set chosen in in table 2. Dominance of X4 viruses is only possible with a certain amount of strains and immune activation being established and is therefore associated with an increased number of strains and a progressed stage of disease.

Figure 4
figure 4

Survival distribution and the coreceptor switch. The probability to have reached the final stage of disease (black) either with dominance in viral load of R5 viruses (blue) or X4 viruses (red) according to the master equation (5), the survival distribution sampled from 200 simulations of equation (1) (grey) agree well with the master equation's predictions i.e. giving support to the assumption of quasi-stationary viral populations, parameters according to table 2.

The results from the probabilistic approach considering equilibrium viral load according to equation (5) are compared with survival curves sampled from simulations of equations (1) based on the actual viral load. The good agreement shows that the flow of mutants is strongly determined by the equilibrium viral load (cf. Additional file 2 and 3, Figure S2 and S3).

As the average time to the onset of AIDS as well as the mean residence times in a given stage decrease with the accumulation of mutant strains (cf. equation (16)) patients with dominance in X4 viral load have a worsened prognosis in this respect. The higher burden in viral load associated with the higher number of strains leads to an accelerated progression of disease (cf. equation (16) and Figs. 2 and 3).

Note that suppression of R5 viruses as facilitated by coreceptor blockers [12, 2123] does not only reduce R5 viral load in the model but also diminishes the X4 viral load due to reduced immune activation (cf. Additional file 4, Figure S4).


A mechanistic assessment of the processes leading to the coreceptor switch in HIV infected patients is difficult - and maybe even inadequate - due to the large variability and stochasticity of the underlying processes in vivo. Quantities such as viral growth or clearance rates are often only known in very specific settings without detailed knowledge of possible hidden dependencies from the underlying dynamical processes - and can in consequence point in contradictory directions. The observation that certain HIV strains may be present in a patient for a long time before they grow to a detectable viral load [39, 40] hints however towards a picture of a dynamical change in the environmental conditions in the host. This holds specifically for X4 viruses [40].

The current model provides a coarse grained, probabilistic picture of possible dynamical scenarios of the coreceptor switch without going yet into the details of virus replication and immune response. This allows to study the basic dynamical features and conditions for a coreceptor switch. In the model, the viral quasispecies [41] is generated by a probabilistic process under the pressure of the immune system. Its composition depends on the randomly chosen evolutionary path that may show some biases due to environmental conditions. These include the viral replication efficiencies based on adequate target cells as well as constraints of immune recognition or therapy. The answer to the question whether a coreceptor switch has occurred is determined by the composition of the viral quasispecies found at the time when the host's immune system loses control and the viral load diverges, i.e. at the onset of AIDS in the model system. One prediction is that the suppression of R5 viruses in the host will not necessarily result in an expansion of X4 viruses. The control of R5 viruses can even delay or prevent the emergence of considerable X4 viral load - which strongly depends on the constitution of the host, i.e. in biases in the choice of the possible evolutionary paths. This is well in accordance with the studies done for a newly licensed drug which blocks the CCR5 receptor [22, 23, 40]. The model provides a stochastic framework to investigate R5 and X4 virus evolution in a dynamic environment, however staying with a qualitative description of the cellular processes. This is a strong simplification of the complex interaction network that interlinks HIV and the hosts immune system. The growing knowledge about their interactions [15, 27, 32, 42] will be implemented within a refined model. This should combine the advantages of this stochastic approach with more details on the involved immune cells populations - both with regard to their role in the immune response and their contribution to persistent immune activation and viral replication - to allow for more quantitative predictions.


Virus evolution

The differential equation for the evolution of viral load and immune response (1) have been solved numerically with routine rk4 from the Math::RungeKutta module for PERL. The system is initialised with one R5 viral strain and one X4 viral strain with vR5(0) = vX4(0) = 1V L, i.e. one unit of viral load for each subtype. During each integration step of duration dt = 0.02d a new mutant arises with probability dt × p m × v = 0.02d × 0.01d-1VL-1 × v proportional to the viral load v. R5 and X4 mutants occur with equal probability.

Equation (2) can be derived from equations (1) following [20] (chapt. 13) by assuming the stationary case, i.e. , i.e.


for {v i |v i ≠ 0}. With the parameters from table 2 this results in


with being the total viral load and and the viral load of one R5 or X4 strain respectively. The total viral load attributed to R5 and X4 viruses can be derived to


leading to equations (6) and (7) (v(nR5, nX4) = vR5(nR5, nX4) + vX4(nR5, nX4)). Note that vR5(nR5, nX4) < 0 corresponds to the case in which no R5 viruses can co-exist with X4 viruses, i.e. a total coreceptor switch [20].

In case of a reduction of the growth rate in R5 viruses from r to rCCR5 <r the above equations change to


leading to a condition for the onset of AIDS


that allows for more virus strains and the condition for dominance of X4 viruses becomes


Note that less R5 viruses can coexist with X4 viruses in case of rCCR5 <r, i.e. neither R5 infection nor X4 infection (due to lack of immune activation) can progress as fast as without therapy.

Probabilistic approach

Assuming fast equilibration of viral load relative to the time scale of emergence of new mutants allows to model the evolution in the number of viral strains as a Markov process [35] with constant transmission rates between states which are proportional to the equilibrium viral load. With P(nR5, nX4, t) being the probability to find nR5 R5 strains and nX4 X4 strains at time t the master equation [35] for this process can be written as


with p1 being the probability per time and viral load to mutate within the same subtype (R5 → R5, X4 → X4) and p2 being the probability per time and viral load to mutate between subtypes (R5 → X4, X4 → R5). For the data shown it holds p1 = p2 = = 0.005d-1VL-1 (for p1p2 cf. Additional file 1, Figure S1). The last term in equation (15) vanishes if the critical number of stains according to equation (3) or the maximal viral load v max = 500V L is exceeded leading to an absorbing boundary. The time evolution of P(nR5, nX4, t) was solved numerically with routine rk4 from the Math::RungeKutta module for PERL with the initial condition P(1, 1, 0) = 1 (zero otherwise) corresponding to the initialisation with one R5 and one X4 strain.

The distribution of the waiting times to the onset of AIDS on the basis of an initial distribution in strain numbers nR5, nX4can be described by a phase type distribution [36]. With p1 = p2 = , a patient's average residence time in a configuration with nR5 R5 virus strains and nX4 X4 virus strains is i.e. decreasing with a growing number of viral strains and in consequence viral load v(nR5, nX4). This allows to derive the mean time to the onset of AIDS in a patient with nR5 R5 virus strains and nX4 X4 virus strains to


The mean waiting time until the onset of AIDS is derived in equation (16) by averaging the waiting times arising along all possible evolutionary paths starting form a configuration with nR5 R5 virus strains and nX4 X4 virus strains (cf. Additional file 5, Figure S5).


  1. Shankarappa R, Margolick JB, Gange SJ, Rodrigo AG, Upchurch D, Farzadegan H, Gupta P, Rinaldo CR, Learn GH, He X, Huang XL, Mullins JI: Consistent viral evolutionary changes associated with the progression of human immunodeficiency virus type 1 infection. J Virol. 1999, 73 (12): 10489-10502.

    PubMed Central  CAS  PubMed  Google Scholar 

  2. Grenfell BT, Pybus OG, Gog JR, Wood JL, Daly JM, Mumford JA, Holmes EC: Unifying the epidemiological and evolutionary dynamics of pathogens. Science. 2004, 303 (5656): 327-332. 10.1126/science.1090727.

    Article  CAS  PubMed  Google Scholar 

  3. Ball CL, Gilchrist MA, Coombs D: Modeling within-host evolution of HIV: mutation, competition and strain replacement. Bull Math Biol. 2007, 69 (7): 2361-2385. 10.1007/s11538-007-9223-z.

    Article  PubMed  Google Scholar 

  4. Regoes RR, Bonhoeffer S: The HIV coreceptor switch: a population dynamical perspective. Trends Microbiol. 2005, 13 (6): 269-277. 10.1016/j.tim.2005.04.005.

    Article  CAS  PubMed  Google Scholar 

  5. Ho SH, Tasca S, Shek L, Li A, Gettie A, Blanchard J, Boden D, Cheng-Mayer C: Coreceptor switch in R5-tropic simian/human immunodeficiency virus-infected macaques. J Virol. 2007, 81 (16): 8621-8633. 10.1128/JVI.00759-07.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Coetzer M, Nedellec R, Salkowitz J, McLaughlin S, Liu Y, Heath L, Mullins JI, Mosier DE: Evolution of CCR5 Use Before and During Coreceptor Switching. J Virol. 2008, 82 (23): 11758-11766. 10.1128/JVI.01141-08.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Tasca S, Ho SH, Cheng-Mayer C: R5X4 viruses are evolutionary, functional, and antigenic intermediates in the pathway of a simian-human immunodeficiency virus coreceptor switch. J Virol. 2008, 82 (14): 7089-7099. 10.1128/JVI.00570-08.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Koning FA, Schols D, Schuitemaker H: No selection for CCR5 coreceptor usage during parenteral transmission of macrophagetropic syncytium-inducing human immunodeficiency virus type 1. J Virol. 2001, 75 (18): 8848-8853. 10.1128/JVI.75.18.8848-8853.2001.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. van't Wout AB, Kootstra NA, Mulder-Kampinga GA, Albrecht-van Lent N, Scherpbier HJ, Veenstra J, Boer K, Coutinho RA, Miedema F, Schuitemaker H: Macrophage-tropic variants initiate human immunodeficiency virus type 1 infection after sexual, parenteral, and vertical transmission. J Clin Invest. 1994, 94 (5): 2060-2067. 10.1172/JCI117560.

    Article  PubMed Central  PubMed  Google Scholar 

  10. Bjorndal A, Deng H, Jansson M, Fiore JR, Colognesi C, Karlsson A, Albert J, Scarlatti G, Littman DR, Fenyo EM: Coreceptor usage of primary human immunodeficiency virus type 1 isolates varies according to biological phenotype. J Virol. 1997, 71 (10): 7478-7487.

    PubMed Central  CAS  PubMed  Google Scholar 

  11. Arien KK, Gali Y, El Abdellati A, Heyndrickx L, Janssens W, Vanham G: Replicative fitness of CCR5-using and CXCR4-using human immunodeficiency virus type 1 biological clones. Virology. 2006, 347: 65-74. 10.1016/j.virol.2005.11.045.

    Article  CAS  PubMed  Google Scholar 

  12. Lewis M, Simpson P, Fransen S, Huang W, Whitcomb J, Mosley M, Dobertson DL, Mansfield R, Ciaramella G, Westby M: CXCR4-using Virus Detected in Patients Receiving Maraviroc in the Phase 3 Studies MOTIVATE 1 and 2 Originates form Pre-existing Minority of CXCR4-using Virus. Conference Reports for NATAP. 2007, []

    Google Scholar 

  13. Hoffmann C: The epidemiology of HIV coreceptor tropism. Eur J Med Res. 2007, 12 (9): 385-390.

    PubMed  Google Scholar 

  14. Keele BF, Giorgi EE, Salazar-Gonzalez JF, Decker JM, Pham KT, Salazar MG, Sun C, Grayson T, Wang S, Li H, Wei X, Jiang C, Kirchherr JL, Gao F, Anderson JA, Ping LH, Swanstrom R, Tomaras GD, Blattner WA, Goepfert PA, Kilby JM, Saag MS, Delwart EL, Busch MP, Cohen MS, Montefiori DC, Haynes BF, Gaschen B, Athreya GS, Lee HY, Wood N, Seoighe C, Perelson AS, Bhattacharya T, Korber BT, Hahn BH, Shaw GM: Identification and characterization of transmitted and early founder virus envelopes in primary HIV-1 infection. Proc Natl Acad Sci USA. 2008, 105 (21): 7552-7557. 10.1073/pnas.0802203105.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Ribeiro RM, Hazenberg MD, Perelson AS, Davenport MP: Naive and memory cell turnover as drivers of CCR5-to-CXCR4 tropism switch in human immunodeficiency virus type 1: implications for therapy. J Virol. 2006, 80 (2): 802-809. 10.1128/JVI.80.2.802-809.2006.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Sguanci L, Bagnoli F, Lio P: Modeling HIV quasispecies evolutionary dynamics. BMC Evol Biol. 2007, 7 (Suppl 2): S5-10.1186/1471-2148-7-S2-S5.

    Article  PubMed Central  PubMed  Google Scholar 

  17. Wodarz D, Lloyd AL, Jansen VAA, Nowak MA: Dynamics of macrophage and T cell infection by HIV. Journal of Theoretical Biology. 1999, 196: 101-113. 10.1006/jtbi.1998.0816.

    Article  CAS  PubMed  Google Scholar 

  18. Callaway DS, Ribeiro RM, Nowak MA: Virus phenotype switching and disease progression in HIV-1 infection. Proceedings of the Royal Society of London Series B-Biological Sciences. 1999, 266 (1437): 2523-2530. 10.1098/rspb.1999.0955.

    Article  CAS  Google Scholar 

  19. Alizon S, van Baalen M: Acute or chronic?: Within-host models with immune dynamics, infection outcome, and parasite evolution. Am Nat. 2008, 172 (6): E244-E256. 10.1086/592404.

    Article  PubMed  Google Scholar 

  20. Nowak MA, May RM: Virus Dynamics, Mathematical Principles of Immunology and Virology. 2000, New York: Oxford University Press

    Google Scholar 

  21. Heinzl S: CCR5-Hemmer: Neuer Angriffspunkt bei HIV Infektion. CTJ. 2007, 16 (6): 195-198.

    CAS  Google Scholar 

  22. Gulick RM, Lalezari J, Goodrich J, Clumeck N, DeJesus E, Horban A, Nadler J, Clotet B, Karlsson A, Wohlfeiler M, Montana JB, McHale M, Sullivan J, Ridgway C, Felstead S, Dunne MW, Ryst van der E, Mayer H: Maraviroc for previously treated patients with R5 HIV-1 infection. New England Journal of Medicine. 2008, 359 (14): 1429-1U27. 10.1056/NEJMoa0803152.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Fätkenheuer G, Nelson M, Lazzarin A, Konourina I, Hoepelman AI, Lampiris H, Hirschel B, Tebas P, Raffi F, Trottier B, Bellos N, Saag M, Cooper DA, Westby M, Tawadrous M, Sullivan JF, Ridgway C, Dunne MW, Felstead S, Mayer H, van der RE: Subgroup analyses of Maraviroc in previously treated R5 HIV-1 infection. N Engl J Med. 2008, 359 (14): 1442-1455. 10.1056/NEJMoa0803154.

    Article  PubMed  Google Scholar 

  24. Regoes R, Bonhoeffer S: HIV coreceptor usage and drug treatment. J Theor Biol. 2002, 217 (4): 443-457. 10.1006/jtbi.2002.3049.

    Article  CAS  PubMed  Google Scholar 

  25. Mosier DE: Changes in HIV-1 tropism: Clinical and prognostic consequences. European Journal of Medical Research. 2007, 12 (9): 371-374.

    PubMed  Google Scholar 

  26. Nowak MA, May RM, Anderson RM: The evolutionary dynamics of HIV-1 quasispecies and the development of immunodeficiency disease. AIDS. 1990, 4 (11): 1095-1103. 10.1097/00002030-199011000-00007.

    Article  CAS  PubMed  Google Scholar 

  27. Lawn SD, Butera ST, Folks TM: Contribution of immune activation to the pathogenesis and transmission of human immunodeficiency virus type 1 infection. Clin Microbiol Rev. 2001, 14 (4): 753-77. 10.1128/CMR.14.4.753-777.2001.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Davenport MP, Zaunders JJ, Hazenberg MD, Schuitemaker H, van Rij RP: Cell turnover and cell tropism in HIV-1 infection. Trends Microbiol. 2002, 10 (6): 275-278. 10.1016/S0966-842X(02)02370-3.

    Article  CAS  PubMed  Google Scholar 

  29. Perelson AS, Neumann AU, Markowitz M, Leonard JM, Ho DD: HIV-1 dynamics in vivo: virion clearance rate, infected cell life-span, and viral generation time. Science. 1996, 271 (5255): 1582-1586. 10.1126/science.271.5255.1582.

    Article  CAS  PubMed  Google Scholar 

  30. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Bunnik EM, Quakkelaar ED, van Nuenen AC, Boeser-Nunnink B, Schuitemaker H: Increased neutralization sensitivity of recently emerged CXCR4-using human immunodeficiency virus type 1 strains compared to coexisting CCR5-using variants from the same patient. J Virol. 2007, 81 (2): 525-531. 10.1128/JVI.01983-06.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Hazenberg MD, Stuart JW, Otto SA, Borleffs JC, Boucher CA, de Boer RJ, Miedema F, Hamann D: T-cell division in human immunodeficiency virus (HIV)-1 infection is mainly due to immune activation: a longitudinal analysis in patients before and during highly active antiretroviral therapy (HAART). Blood. 2000, 95: 249-255.

    CAS  PubMed  Google Scholar 

  33. Bozek K, Thielen A, Sierra S, Kaiser R, Lengauer T: V3 loop sequence space analysis suggests different evolutionary patterns of CCR5- and CXCR4-tropic HIV. PLoS ONE. 2009, 4 (10): e7387-10.1371/journal.pone.0007387.

    Article  PubMed Central  PubMed  Google Scholar 

  34. Chesebro B, Wehrly K, Nishio J, Perryman S: Macrophage-tropic human immunodeficiency virus isolates from different patients exhibit unusual V3 envelope sequence homogeneity in comparison with T-cell-tropic isolates: definition of critical amino acids involved in cell tropism. J Virol. 1992, 66 (11): 6547-6554.

    PubMed Central  CAS  PubMed  Google Scholar 

  35. Gardiner CW: Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences. 2008, Springer, 3:

    Google Scholar 

  36. Latouche G, Ramaswami V: Introduction to Matrix Analytic Methods in Stochastic Modeling. ASA-SIAM. 1999

    Book  Google Scholar 

  37. van Rij RP, Hazenberg MD, van Benthem BH, Otto SA, Prins M, Miedema F, Schuitemaker H: Early viral load and CD4+ T cell count, but not percentage of CCR5+ or CXCR4+ CD4+ T cells, are associated with R5-to-X4 HIV type 1 virus evolution. AIDS Res Hum Retroviruses. 2003, 19 (5): 389-398. 10.1089/088922203765551737.

    Article  CAS  PubMed  Google Scholar 

  38. Salemi M, Burkhardt BR, Gray RR, Ghaffari G, Sleasman JW, Goodenow MM: Phylodynamics of HIV-1 in lymphoid and non-lymphoid tissues reveals a central role for the thymus in emergence of CXCR4-using quasispecies. PLoS ONE. 2007, 2 (9): e950-10.1371/journal.pone.0000950.

    Article  PubMed Central  PubMed  Google Scholar 

  39. Madsen TV, Gerstoft J, Nielsen C, Jorgensen LB: Reappearance of an 11-year-old sequence in an HIV-1 infected patient during treatment interruption. Scand J Infect Dis. 2008, 40 (2): 174-182. 10.1080/00365540701558706.

    Article  PubMed  Google Scholar 

  40. Westby M, Lewis M, Whitcomb J, Youle M, Pozniak AL, James IT, Jenkins TM, Perros M, van der RE: Emergence of CXCR4-using human immunodeficiency virus type 1(HIV-1) variants in a minority of HIV-1-infected patients following treatment with the CCR5 antagonist Maraviroc is from a pretreatment CXCR4-using virus reservoir. J Virol. 2006, 80 (10): 4909-4920. 10.1128/JVI.80.10.4909-4920.2006.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  41. Kamp C: A quasispecies approach to viral evolution in the context of an adaptive immune system. Microbes Infect. 2003, 5 (15): 1397-1405. 10.1016/j.micinf.2003.10.001.

    Article  CAS  PubMed  Google Scholar 

  42. Catalfamo M, Di Mascio M, Hu Z, Srinivasula S, Thaker V, Adelsberger J, Rupert A, Baseler M, Tagaya Y, Roby G, Rehm C, Follmann D, Lane HC: HIV infection-associated immune activation occurs by two distinct pathways that differentially affect CD4 and CD8 T cells. Proc Natl Acad Sci USA. 2008, 105 (50): 19851-19856. 10.1073/pnas.0810032105.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references


I would like to thank B. Schnierle and B. Krause for many helpful discussions on HIV evolution during the preparation of the manuscript. Their feedback as well as that of colleagues from the section of Biostatistics at the Paul-Ehrlich-Institute, specifically of K.M. Hanschmann, have contributed much to the improvement of the manuscript. Finally I would like to thank the referees for their detailed remarks allowing me to elaborate the approach more transparently.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Christel Kamp.

Additional information

Authors' contributions

CK carried out the study and wrote the manuscript.

Electronic supplementary material


Additional file 1: Mutation between R5 and X4 viruses. The Figure shows a comparison between the situation of homogeneous and heterogeneous mutation rates among R5 and X4 viruses. While the former case (left panel) corresponds to the situation discussed in the main paper, i.e. mutation rates being identical among and between subtypes p1 = p2 = , the latter case (right panel) has the same total mutation rate p m = 0.01d-1 V L-1 but a threefold higher probability for intra-subtype mutation than inter-subtype mutation, i.e. p1 = , p2 = . The top row shows how this assumption shifts evolutionary paths towards routes with a higher fraction of R5 viruses. In consequence, X4 dominance is less often attained before the onset of AIDS leading to a lower fraction of coreceptor switches. This might however be shifted to the observed levels by a stronger coupling of X4 growth to immune activation. The survival curves decay slightly steeper in the case of heterogeneous mutation patterns because more viable R5 mutants (than not yet adapted X4 mutants) are generated in the earlier stages of disease. The exact mutation rates among R5 and X4 viruses are hard to estimate, but their sequence similarity suggests them to be of a similar order of magnitude (p1p2) [33]. (PDF 438 KB)


Additional file 2: Equilibrium viral load. Viral load of R5 viruses (blue) and X4 viruses (red) in a simulation of the model described by equations (1) as shown in Fig. 1 (parameters as in table 2). In addition the equilibrium viral load according to equations (11) is shown with good agreement with the simulation data. Deviations occur only at the breakdown of the system where no equilibration can be expected any more. (PDF 56 KB)


Additional file 3: Variation of mutation rate p m . Survival distribution for p m = 0.01 (cf. Fig. 4) and p m = 0.1 sampled from equations (1) (200 runs each) in comparison with the probabilistic approach (5), other parameters as in table 2. The ten-fold increase of mutation rate shrinks the time axis of the survival distribution by a factor of 10. The survival curves for p m = 0.1 are still close but the sampled curve decays faster than predicted by the master equation approach as mutants are already likely to be established in the initial viral peak, i.e. do not allow for equilibration. (PDF 48 KB)


Additional file 4: Suppression of R5 viruses and viral load. The density plots show the equilibrium viral load determined by equations (11,12) (as visualised in Additional file 2, Figure S2) depending on the numbers of R5 and X4 virus strains being present. The top row shows the equilibrium viral load for R5 viruses vR5(nR5, nX5) and X5 viruses vX4(nR5, nX4) for parameter values as depicted in Table 2. The bottom row shows the same situation, however, with the growth rate of R5 viruses halved to 1d -1 leading to a situation in which hardly any R5 viral load can be established. This however, results at the same time in a lower viral load from X4 viruses (cf. right panel). In consequence, the model predicts an indirect positive effect from the suppression by R5 viruses as induced by CCR5 blockers. (PDF 330 KB)


Additional file 5: Mean waiting time to the onset of AIDS. The Figure shows a graphical representation of equation (16) in the nR5-nX4-plane showing the probability to find nR5 R5 strains and nX4 X4 strains in a patient at 1200 days since infection (cf. Fig. 2). The mean waiting time to the onset of AIDS is determined by averaging the waiting times along all possible evolutionary paths. Therefore, decreases with a growing number of strains (here: T4,4 > T6,5). (PDF 11 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Kamp, C. Understanding the HIV coreceptor switch from a dynamical perspective. BMC Evol Biol 9, 274 (2009).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: